EM-Driven Multi-Objective Optimization of a Generic Monopole Antenna by Means of a Nested Trust-Region Algorithm

Antenna structures for modern applications are characterized by complex and unintuitive topologies that are difficult to develop when conventional, experience-driven techniques are of use. In this work, a method for the automatic generation of antenna geometries in a multi-objective setup has been proposed. The approach involves optimization of a generic spline-based radiator with an adjustable number of parameters using a nested, trust region-based algorithm. The latter iteratively increases the dimensionality of the radiator in order to gradually improve its performance. The method has been used to generate a set of nine antenna designs, representing a trade-off between minimization of reflection within 3.1 GHz to 10.6 GHz and a reduction of size. The properties of the optimized designs vary along the Pareto set from −10 dB to −20 dB and from 230 mm2 to 757 mm2 for the first and second objectives, respectively. The presented design approach has been validated against a genuine, population-based optimization routine. Furthermore, the smallest Pareto-optimal design has been compared to the antennas from the literature.


Introduction
Ultra-wideband (UWB) antennas belong to the key components of modern location and tracking systems capable of providing centimeter-level accuracy [1][2][3][4][5][6]. For the mentioned services, the radiators may find application not only to data exchange, but also as energy harvesting devices that provide the necessary power for sensors [7][8][9][10]. A variety of application scenarios for tracking and location devices impose a range of potential requirements on antenna performance (in terms of field and electrical properties), as well as its dimensions (small footprint or specific shape) [10][11][12][13][14][15]. Consequently, a coherent system often comprises a range of antenna solutions.
Other approaches to the generation of an antenna structure involve representation of its geometry using a set of characteristic points interconnected by lines or appropriately defined curves [26][27][28][29][30]. Alternatively, the topology can be defined using a binary matrix (often with overlapping components with variable dimensions) [31][32][33], [34] where zeros and ones correspond to the areas with etched or continuous metallization [32]. Contrary to conventional radiators, antenna topologies described using generic models are predominantly determined using optimization methods [35][36][37][38][39]. Notwithstanding, the applicability of numerical algorithms to the generation of unconventional geometries is often demonstrated based on relatively simple, single-objective problems [28,29]. Furthermore, the optimization of generic models is normally undertaken using population-based routines (e.g., particle swarm optimization, genetic algorithms, or ant colony methods), which require thousands of evaluations to obtain the final solution [31,35,38,40]. To put that into perspective, in [41], the particle swarm algorithm required a total of 248 EM simulations to solve a set of 12 problems represented using only one input parameter. In [42], the same method required over 1400 model evaluations to find an optimal design of frequency selective surface represented using three adjustable variables. Due to the high evaluation cost (being in the order of minutes to hours), the application of conventional, populationbased routines for the EM-driven optimization of antennas defined using generic models is impractical from a computational standpoint.
The design of modern antennas is a multi-objective (MO) problem, where a set of performance or geometry requirements is to be simultaneously fulfilled [38,43,44]. Conventional approaches to MO tasks involve conversion of the performance requirements into a single-objective function using aggregation techniques, penalty components, or design constraints [35,39,43]. The design task, formulated in this manner, can be solved using conventional single-objective routines [35,43]. On the other hand, aggregation techniques provide only one solution to the problem per algorithm run, whereas the quality of the result depends on the weights assigned to individual figures of interest [35]. Alternatively, MO optimization can be performed using the genuine MO algorithms [39,43]. The latter ones are especially useful when the designer's priorities are not clearly defined, as they provide a set of solutions representing trade-offs between the requirements. Such cognitive studies are useful from the perspective of gaining a detailed insight into the capabilities of the antenna, with respect to the desired application scenarios or environments [45]. Although multi-objective problems could be solved using population-based metaheuristic methods, their practical applications are limited to rather simple structures due to the prohibitive computational cost [35,37]. Less numerically demanding routines involve the execution of evolutionary algorithms on the approximation models. Due to the curse of dimensionality [46], the usefulness of approximation models and metaheuristics for reliable MO optimization is limited either to low-dimensional structures or narrow search spaces [43,47]. Notwithstanding, the mentioned methods proved to be useful for solving problems characterized by over a dozen design parameters [43]. Alternatively, MO design can be handled using dedicated heuristic approaches oriented toward the sequential generation of Pareto-optimal solutions [43,47]. Nonetheless, neither of the mentioned methods are suitable for handling antenna models characterized by a flexible number of adjustable parameters.
The above considerations indicate that contemporary UWB antennas can be categorized into a few classes [16][17][18][19][20][21][22][23][24]. The topologies from each group are similar except for features such as slots or stubs that are normally introduced as a result of manual, greedy search-based studies to alter the functionality of the structure. In this regard, representation of the antenna using a generic model might be useful for the automated development of its topology regarding the imposed performance and geometry requirements. Although "universal" antenna models have been reported in the literature, they lack flexibility in terms of the number of degrees of freedom used for their representation. Consequently, they are either too simple or too complex (parameter-wise) to obtain compromise designs. The application of generic topology (or its derivatives) for a range of application scenarios should reduce the cost of system development, as well as its time to market. Unfortunately, the available MO routines are inadequate for handling models with flexible parameterization due to the lack of mechanisms for the adjustment of problem dimensionality in the course of the optimization process. The motivation of this work is to address the mentioned challenges.
In this paper, a generic EM model of a planar antenna with flexible parameterization is presented. The structure is optimized in a bi-objective setup concerning minimization of in-band reflection within a 3.1 GHz to 10.6 GHz band and size. A set of designs, representing the trade-off between the requirements, is obtained using a novel framework that sequentially generates candidate solutions. The latter are refined in a nested trust region-based optimization loop that enables one to gradually increase the number of parameters. The proposed model and optimization framework are demonstrated through the generation of nine Pareto-optimal designs. Their geometries are represented using 30 to 46 input parameters. The variability of the antenna's in-band performance and size along the Pareto set is from −10.3 dB to −20.5 dB (over 10 dB) and from 230 mm 2 to 757 mm 2 (almost 70 percent), respectively. The presented algorithm has been compared against other population-based metaheuristic routines in terms of optimization cost and the quality of obtained solutions. Numerical validation of the generated designs and comparisons of the smallest obtained UWB structure with the state-of-the-art compact antennas are also provided.

Generic UWB Monopole Antenna
Consider a planar monopole antenna, shown in Figure 1. The structure consists of a driven element in the form of a spline-parameterized radiator fed through a microstrip line. The ground plane edge oriented toward the radiating component is also implemented in the form of a spline curve. The presented topology generates closely located resonant modes, which contribute to the broadband behavior [48,49]. Furthermore, the bandwidthrelated performance is enhanced using an L-shaped ground plane stub, which improves the impedance matching of miniaturization-oriented designs [24]. The structure is implemented on a Rogers RO4003 dielectric substrate with a permittivity of 3.58, thickness of 0.813 mm, and loss tangent of 0.0027. The vector of the antenna input parameters is y = [y a y g y r ] T , where y a = [X l f l 1 l 2r w 1 o r ] T , y g = [y g.1 . . . y g.L ] T and y r = [y r.1 . . . y r.L ] T (L = 50 represents the number of points used to represent the ground plane and radiator in the antenna model). The relative variables are l 2 = (X − w 1 )l 2r , Y = l 1 + w 1 , l fr = min{0.5X, 0.5(Y − l f )}, o = 0.5X + o r . Dimension w f = 1.8 remains fixed to ensure 50 Ohm input impedance. The unit for all variables except dimensionless l 2r is mm. It should be noted that vectors y g = Y·z g , and y r = S·z r , where S = min{0.5(X − o r ), 0.5(Y − l f )} represent the characteristic points (so-called knots) of the spline curves used to define the ground plane and radiator, respectively (see Figure 1). To ensure the structure consistency, the gap between each knot y g and y r is even. Hence, the vectors denote the relative distance of the curve from the reference point (edge of the substrate for y g and center of the radiator for y r , respectively; Figure 1). It is worth emphasizing that y r represents axial dimensions, as the parameters of the radiator are defined in a cylindrical coordinate system. The EM model of the antenna structure R is implemented in CST Microwave Studio and evaluated using its time domain solver. The presented antenna supports flexible parameterization. In other words, the number of adjustable variables used for representation of the structure geometry can be changed in the course of the design process. The rationale behind this feature is that-when the generation of compromise design solutions is considered-the initial search space ranges are set to be large to ensure that the optimal designs are captured within them. The application of conventional (population-based) methods for exploration of the search space with wide parameter ranges is numerically impractical due to the tremendous number of iterations required for convergence [38,43]. Having a model with adjustable dimensionality, one can narrow down the search space to the region of interest using a relatively small number of parameters. The promising region can be further exploited through optimizations of the model characterized by gradually increasing dimensionality. Flexible parameterization has been implemented through implicit handling of the input variables. Let . . , L represent the d-knot spline curve used to represent the antenna ground plane and radiator, respectively. Note that when d < 4, x g.d = x g.1 and x r.d = x r.1 , because the implementation of the spline used in this work requires at least four knots. The change of the antenna dimensionality is realized through the interpolation of x g.d and x r.d to z g = z g.L and z r = z r.L vectors, respectively. The process can be defined as follows: where sp is the function that realizes L-point spline-based interpolation of the d-dimensional input vector. The conversion process is transparent to the design optimization algorithm. The lower and upper bounds for the structure optimization are l b = [6 4 10 0.05 0.

Design Methodology
In this section, the proposed design methodology is explained. To make the paper self-contained, a short discussion of the multi-objective design problem, followed by a brief description of the utilized optimization algorithm, are provided here. Next, the introduced multi-objective optimization method that permits a gradual increase of the number of antenna parameters is described. The section is concluded with the summary of the presented design framework. Numerical results are provided in Section 4.

Multi-Objective Design Problem
The multi-objective design problem can be defined as follows [39]: where x = x d (see Section 2) denotes the vector of the design parameters, X is the feasible region of the search space, and x * represents the optimal design. The objective functions are given as is a vector of responses obtained from EM evaluation of the antenna model at the design x. The responses may represent the electrical performance (e.g., reflection over a frequency range or isolation between radiators), field (e.g., single-direction gain over the frequency or radiation pattern at the given frequency), or geometry (e.g., footprint or volume) of the structure. If Q > 1, the conventional notion of performance improvement known from single-objective optimization (i.e., the one based on direct comparison of the objective function responses obtained for candidate solutions) is not applicable [39]. Instead, designs can be compared using a domination relation [39]. One can say that the design x 1 dominates over for all q and F q (x 1 ) < F q (x 2 ) for at least one objective. A design x q * is optimal in the Pareto sense if there is no other design x r * (q = r), such that x r * dominates over x q * . The solution to the multi-objective problem is a set of Pareto-optimal designs. For a more detailed description of the MO optimization concept, see [35,39,43].

Nested Optimization with a Variable Number of Design Parameters
Although conventional global-optimization algorithms proved to be useful for solving complex problems, they suffer from an unacceptable numerical cost, associated with the tremendous number of model evaluations required to find the optimal designs. On the other hand, deterministic algorithms (such as gradient or pattern search methods [36,50]) are often characterized by a lower computational cost yet limited capability in terms of escaping from local minima. In this work, the problem is mitigated by performing a series of optimizations in a nested loop. The latter gradually increases the number of design variables, first to identify and then to exploit the region of the search space of interest. The analogy to this approach, from the perspective of the design objective space, is observation of the landscape from a high altitude where the amount of visible details is limited. With each iteration of the nested loop (that is associated with an increasing the number of parameters), the "altitude" decreases, which allows the algorithm to "identify" more details of the landscape and improve the design response. Let . . , L, n = 1, . . . , N represent the set defining the number of spline parameters (positive integer values) used to represent the ground plane and radiator, respectively. In each step, the final d n dimensional design x dn * , obtained in the nth nested iteration, is used as a starting point for the n + 1 step. Assuming that the vector d is defined, n is set, and the initial design x dn * is available, the nested optimization can be summarized as follows: Find design x dn+1 * using x dn+1 (0) as a starting point for optimization; 3.
If n = N, END; otherwise, set n = n + 1 and go to 1.
It should be noted that the size N of vector d substantially affects the algorithm cost, as the optimization process is reset at each iteration n = 1, . . . , N. On the other hand, the number of iterations and step size between d n and d n+1 can be adjusted so as to mitigate the risk of getting stuck in a poor local optimum while managing the optimization cost at the acceptable level. Conceptual visualization of the objective function (a so-called feature) and design spaces for the structure, represented using the increasing number of parameters, is shown in Figure 2.
The red line represents the Pareto front [39]. Note that the solutions that are attainable for the structure represented using, for example, d 1 variables, represent a subset of the responses that can be found when the dimensionality is increased.

Generation of Trade-Off Designs
The presented bi-objective optimization algorithm involves two stages: (i) generation of the extreme (or corner) designs and (ii) sequential determination of the Pareto set using constrained optimizations. The first corner solution is obtained by solving the following minimization problem: Note that Equation (4) involves optimization of the structure in regard to a single objective. The second design is obtained through constrained optimization, oriented toward minimization of the second performance criterion: The goal of Equation (5) is to minimize F 2 while ensuring that the response of the optimized design x 2.dn * does not violate the defined threshold F 1.max . The reasoning behind the introduction of a design constraint is to permit minimization of the second objective while maintaining the response regarding F 1 (x dn ) at an acceptable level. Generation of the extreme points is conceptually illustrated in Figure 3.
In the second stage, the algorithm sequentially identifies the Pareto designs. In each step, the new points are obtained between the already existing optimal solutions. Let x p.dn * (p = 1, . . . , l k ; l k = 2 k−1 + 1) represent the Pareto designs obtained in the k − 1 iteration, where k = 1, . . . , K. It should be noted that l 0 = 2 refers to two corner designs-x 1.dn * , x 2.dn * -found in the first design stage. The candidate solutions are determined as follows: where c = l k + b and b = 1, . . . , l k − 1. It should be noted that the reference designs in Equation (6) are approximated using d m variables. Hence, the cth candidate design can be generated even when using x * p.dn and x * p+1.do , d n = d o (o = 1, . . . , N). This is ensured by interpolating the reference designs' dimensionality as ← sp(x p+1.do * , d m ) (see Section 2). The obtained candidate solutions are utilized as starting points for constrained optimization: The goal of Equation (7) is to minimize the objective function F 1 while maintaining the structure response with respect to the second objective below the specified threshold. The latter is set to F 2.c.max = 0.5(F 2 (x p.dn * ) + F 2 (x p+1.dn * )). In order to facilitate generation of the Pareto set, the constraint on F 2 is implemented in the form of a penalty coefficient included with the objective function. The process in Equation (7) is performed in a nested optimization loop as shown in Section 3.2. As already mentioned, such an approach allows one to gradually increase the quality of the Pareto designs by appropriate adjustment of the antenna dimensionality.  Figure 4 for a block diagram) can be summarized as follows: 1.
Set the initial designs x 1 . Generate x c.dn * by solving Equation (7) in a nested optimization loop; 7.
If p = l k , go to 7; otherwise set p = p + 1 and go to 5; 8.
If k = K, END; otherwise, set k = k + 1 and go to 4. In the first stage of the presented framework, the initial designs x 1.d1 (0) and x 2.d1 (0) are obtained using a greedy search-based topology adjustment. The goal of the procedure is to find reasonable starting points for the optimization process. Upon identification of the initial designs, the generation of the Pareto set is handled by the above-outlined routine.
It is worth mentioning that the concept of generating intermediate designs is similar to the one presented in [47]. Here, however, the bisections are used only for the generation of candidate designs, which are further refined by means of constrained optimization. Furthermore, the algorithm of [47] does not support increasing the dimensionality of the problem in the course of structure optimization. Owing to the identification of new designs between the already existing ones, the algorithm allows for obtaining a Pareto set characterized by relatively even separation between the consecutive solutions (in the objective function space). Generation of the Pareto designs using the discussed algorithm is conceptually illustrated in Figure 5.

Optimization Engine
The optimization engine used in this work for solving Equations (4), (5), and (7) is a local gradient method [36]. The algorithm is embedded in a trust region (TR) framework [37]. The latter generates a series (i = 1, 2, 3, . . . ) of local models around what is currently the best design x (i) in order to gradually approximate the optimal solution to be found. The TR-based optimization is defined as [37] x (i+1) = arg min where G (i) is a first-order Taylor expansion model, given by [51] where J(x (i) ) represents the Jacobian generated around x (i) using a large-step finite differentiation [51]. The TR radius r (i) (i.e., the region around x (i) for which the model is considered valid) is adjusted based on the calculated gain ratio [51]: The candidate design x (i+1) is accepted only if ρ > 0. Otherwise, the iteration is reset. The radius is updated according to the standard rules [37]. The initial value r (0) is controlled by the nested optimization loop (i.e., for n = 1, the radius is set to r (0) = 1, whereas for n > 1, r (0) = 1/d n ) (see Section 3.2). The algorithm is terminated either when the Euclidean distance between two consecutive designs or the radius are smaller than a specified threshold (i.e., ||x (i+1) − x (i) || ≤ ε, or r (i+1) ≤ ε (here, ε = 10 −2 )). The algorithm can be summarized as follows: 1.
Set i = 0 and x (0) , then set r (0) according to the iteration of the nested loop in Section 3.2; 2.
Generate the Jacobian J(x (i) ) and construct the local approximation model G (i) ; 3.
If the termination condition is met, END; otherwise, set i = i + 1 and go to 2.
Owing to utilization of the Taylor model, the computational cost of the algorithm is low, as it requires d n + 1 simulations of the R model per successful iteration. Additional evaluations are required for the steps that do not improve the objective function value. The justification for decreasing the initial radius when n > 1 is that with each iteration of the nested loop, the goal of the algorithm is to exploit the region identified in the previous step. For more comprehensive discussion on TR optimization methods, see [37,51].

Numerical Results
In this section, the presented antenna structure and design optimization algorithm are numerically validated. The design problem involves bi-objective optimization of the radiator for size reduction and minimization of reflection within the frequency band of interest. The obtained EM simulation results were verified based on evaluations performed using a frequency domain solver. Furthermore, the presented routine has been compared against a conventional multi-objective optimization algorithm in terms of computational cost and quality of the obtained Pareto designs. Finally, the smallest obtained antenna solution has been benchmarked against other structures from the literature. All computations have been performed using an Intel Xeon E5-2650 processor with 32 GB of RAM.

Bi-Objective Antenna Optimization
The design objectives for the considered antenna optimization problem included the following: the minimization of reflection F 1 , such that|S 11  It is worth noting that, due to being represented using only one variable, the radiators and ground planes of the initial designs were in the forms of circles and rectangles, respectively. Both antenna geometries are shown in Figure 6. The objective function for optimization with respect to the first performance requirement was defined as where Here, |S 11 | fv represents the reflection at M frequency points f v ∈ f, the value of which was higher than the S max = −15 dB threshold. Note that the definition of Equation (12) in the least squares sense provided smooth changes of the objective function values in the course of the optimization process [51]. On the other hand, E 1 (x) = 0 when the antenna response was below S max for all f v points. In this case, the control over optimization was taken over by the min-max component in Equation (13). The scaling coefficient was set to α 1 = 10 4 in order to ensure a smooth transition between the E 1 (x) and E 2 (x) functions. The second design objective was given as where α 2 = 0.01 in order to ensure that the contribution of E 1 (x) to F 2 (x) was noticeable for a relatively small violation of S max = F 1.max = −11 dB. It should be noted that for the majority of the antenna structures, a typical design requirement is to ensure that the reflection is below the −10 dB threshold. However, in practice, minimization of Equation (14) often leads to violation of the condition [43]. Setting a lower-than-required threshold S max allowed us to ensure that |S 11  Note that the dimensionality of the extreme points was different (12 relative variables per feature-30 total-for x 1 * vs. 20-46 total for x 2 * ), which was due to the lack of performance improvement from increasing the number of x 1 variables above 12. The obtained solutions were used as starting points for generation of the Pareto set. In the second design stage, the Pareto designs were obtained using an objective of the following form: The additive component of Equation (15) is a penalty coefficient activated when the threshold A c = 0.5(A(x m.d * ) + A(x m+1.d * )) is violated. In F 1 (x), the parameter S max = F 1.c.max (see Section 3.3). The number of algorithm iterations was set to K = 3. The candidate solutions for optimization were generated using Equation (6). The Pareto set, composed of two corner and seven intermediate designs {x 1.20 * x 2.20 * x 3.20 * x 4.20 * x 5.20 * x 6.20 * x 7.20 * x 8.12 * x 9.12 * }, and the frequency responses of the selected optimal solutions are shown in Figure 8. The obtained results indicate that the footprints of the optimized antennas ranged from 230 mm 2 (with an in-band reflection of −10.3 dB) to 757 mm 2 (reflection of −20.5 dB). Hence, the change in antenna size and frequency responses along the Pareto set were 527 mm 2 and 10.2 dB, respectively. It is worth noting that the miniaturization rate of the smallest antenna with reference to the largest one was almost 70%. The gain pattern responses were obtained for the selected antenna designs (i.e., x 1.20 * , x 9.12 * (corner designs), and  Figure 9 for a visualization of the topologies). The characteristics shown in Figure 10 indicate that the x 1.20 * design was characterized by slightly more omnidirectional radiation patterns in the x-y plane and lower gain compared with the x 5.20 * and x 9.12 * geometries.

Comparisons and Discussion
To confirm the validity of the numerical results, the x 1.20 * , x 5.20 * , and x 9.12 * designs were evaluated using a frequency domain solver (finite elements method with adaptive meshing) implemented in CST Studio [52]. The response characteristics (shown in Figure 11) obtained using the time and frequency domain methods were similar, which indicates that the EM models used in the optimization process were valid. Furthermore, the plots of Figure 11 show that, at the lower corner frequency of 3.1 GHz, the x 1.20 * , x 5.20 * , and x 9.12 * designs featured radiation efficiencies of 76%, 89%, and 99%, respectively. Note that for the x 1.20 * design, an efficiency of over 90% was achieved above the 3.4 GHz frequency, which is considered acceptable, keeping in mind its small footprint. Overall, the considered antennas were characterized by average in-band efficiencies of 93%, 97%, and 98%, respectively. Clearly, the high values were due to the utilization of a low-loss substrate for implementation of the structures.  Figure 12a shows a comparison of the optimal solutions obtained at different iterations of the nested loop (represented using a varying number of parameters). As can be seen, for the majority of the designs, increasing the dimensionality resulted in improvement in the structure responses. However, for certain solutions, a larger number of variables led to performance degradation. This was due to insufficient precision of the spline approximation when the number of knots increased, which led to undesirable frequency shifts at the edges of the operational bandwidth. Notwithstanding this, for most of the designs, the imperfections were corrected by the optimization algorithm.  Table 1. It should be noted that the average numerical expenses associated with the identification of a single Pareto-optimal design corresponded to about 488 R model evaluations, which is low when having in mind the dimensionality of the problem (up to 46 independent parameters). It is also worth noting that optimization of the designs represented using more than 30 variables amounted to 57% of the overall cost. For the sake of comparison, the antenna at hand was also optimized using a standard population-based multi-objective evolutionary algorithm (MOEA) with fitness sharing, mating restrictions, and Pareto dominance tournament selection [53]. The setup of the benchmark routine was as follows. The population size and the number of generations were set to 200 and 50 (a total of 10,000 EM simulations), respectively, whereas the probability of mutation and crossover were fixed at 0.2 and 0.5, respectively. Since the discussed algorithm did not support adjusting the antenna dimensionality, the number of parameters was set to 30 (i.e., d 4 = 12). Furthermore, the corner designs obtained using the method in Section 3 were included into the initial population. The latter was generated between the extreme points in order to narrow down the search space to the region of interest. Figure 12b shows a comparison of the designs obtained using both routines (the plots present the responses obtained for the antenna model, represented using 30 parameters). As can be seen, the algorithm of Section 3 outperformed the benchmark method in terms of performance of the resulting designs. More importantly, the number of R simulations required to obtain the final designs using the proposed algorithm amounted to 1889, which was 81% less compared to the cost of MOEA operation. Even when the optimization was performed in a nested loop, the cost of the presented algorithm (4389 EM simulations) was over two times lower compared with the benchmark routine (while providing noticeably better responses characterized by a more even spread between the extreme designs and better in-band reflection). One should emphasize that the dimensions of the smallest Pareto-optimal antenna design were only 12.5 mm × 18.4 mm (footprint = 230.4 mm 2 ). The structure was compared against other compact UWB monopoles from the literature in terms of size and bandwidth. For the sake of fair comparison, the dimensions of the considered structures were expressed in terms of the guided wavelength λ g calculated for the lower corner frequency and the electrical performance of the substrates on which the considered radiators were implemented. The results gathered in Table 2 indicate that the proposed structure was characterized by a competitive size while providing acceptable performance. Small dimensions make the design's potential use not only for a range of mobile UWB sensors and IoT-based localization services, but also for body of area network applications [54].

Conclusions
In this work, a spline-based UWB monopole antenna with an adjustable number of parameters has been presented. The structure has been optimized in a bi-objective setup using an algorithm that enables sequential generation of the Pareto designs between the extreme solutions. Each step of the design process involves a trust region-based optimization of the structure in a nested loop. The latter gradually increases the dimensionality of the antenna. The trade-off solutions are generated through the optimization of candidate designs, determined from bisections between the already existing solutions. The proposed method was used to obtain a total of nine antenna designs, representing compromises between a minimization of reflection within 3.1 GHz to 10.6 GHz band and size reduction. The geometries of the optimized radiators were represented using 30 to 46 independent parameters. The performance and sizes of the final designs along the Pareto front varied from −10.3 dB to −20.5 dB and from 230 mm 2 to 757 mm 2 , respectively. The method has been favorably validated against a genuine multi-objective optimization algorithm in terms of computational cost and the quality of the obtained solutions. Furthermore, the smallest Pareto-optimal antenna was compared against the state-of-the-art monopoles from the literature. The optimized radiator is characterized not only by competitive performance, but also a small footprint, which makes it of potential use for UWB sensors, mobile localization services, or body area networks. It is worth noting that the availability of the Pareto set provides the designer a detailed insight into the capabilities of the antenna for a range of application scenarios, which might be of high importance when designing heterogeneous sensors or IoT networks dedicated to operating in diverse environments is considered.
Future work will focus on the application of spline-parameterized antenna models to problems defined using more complex sets of requirements. The development of routines oriented toward reduction of the optimization cost will also be considered.