Multi-Fidelity Local Surrogate Model for Computationally Efficient Microwave Component Design Optimization

In order to minimize the number of evaluations of high-fidelity (“fine”) model in the optimization process, to increase the optimization speed, and to improve optimal solution accuracy, a robust and computational-efficient multi-fidelity local surrogate-model optimization method is proposed. Based on the principle of response surface approximation, the proposed method exploits the multi-fidelity coarse models and polynomial interpolation to construct a series of local surrogate models. In the optimization process, local region modeling and optimization are performed iteratively. A judgment factor is introduced to provide information for local region size update. The last local surrogate model is refined by space mapping techniques to obtain the optimal design with high accuracy. The operation and efficiency of the approach are demonstrated through design of a bandpass filter and a compact ultra-wide-band (UWB) multiple-in multiple-out (MIMO) antenna. The response of the optimized design of the fine model meet the design specification. The proposed method not only has better convergence compared to an existing local surrogate method, but also reduces the computational cost substantially.


Introduction
Internet of Things (IoT) technology allows a large number of electronic devices interconnected wirelessly [1,2]. It demands higher data transmission rate and higher operating frequency than the current communication technology (e.g., 5G) can offer. The microwave components for IoT applications require performances such as wide frequency bands, high power, and small size. Therefore, the structure of these components have become increasingly complex [3][4][5][6][7]. Reliable and computationally efficient design for high performance microwave components is one of the major challenges in IoT technology.
By encoding given design specifications into an objective function (typically related to S-parameters evaluated over certain frequency bands), the microwave component design task can be formulated as a nonlinear optimization problem, and solving this problem yields a better design. In the complex microwave component optimization problem, high-fidelity electromagnetic (EM) simulation allows accurate response evaluation but its computation is expensive. It is impracticable to use high-fidelity EM simulation models directly for numerical optimization in many cases, because conventional

Multi-Fidelity Local Surrogate Model Optimization
In this section, we present the multi-fidelity surrogate model based on local response surface approximations. In the optimization process, local region modeling and optimization are performed iteratively. The judgment factor is introduced to adaptively change local region size, so that the algorithm converges quickly and accurately to the optimal solution. The last surrogate model is refined by space mapping to obtain the optimal design with higher accuracy.

Microwave Component Design Problem
A microwave component design task can be formulated as a nonlinear minimization problem where x ∈ n is a vector of n design variables The vector x * is the optimal design to be found. R f ∈ m represents a fine model response vector obtained by EM simulation at high computational costs (e.g., the transmission coefficient |S 21 | of a component evaluated at m different frequencies). U(·) is a scalar merit function related to a design specification. In many cases, the objective function can be defined in the norm sense as U R f (x) = R f (x) − R * T , where R * T is a target response (or design specification). It is impractical to directly solve the above problem using conventional optimization techniques (e.g., quasi-Newton algorithm and genetic algorithm, etc.). These techniques usually require a large number of the fine model evaluations, resulting in prohibitive computational cost.
The SBO approach aims at addressing the aforementioned difficulties [13]. It speeds up the design process by transferring the optimization burden to inexpensive yet reasonably accurate surrogate models. Thus, the original design problem (1) is replaced by an iterative procedure The procedure produces a sequence of designs x (i) , i = 0, 1, . . .. Each one is the optimal design for the ith surrogate model R (i) s respectively. Assuming that the accuracy of surrogate models are increasing in optimization process, the design sequence can eventually converge to the optimal design of the fine model.

Multi-Fidelity Local Surrogate Models
The design optimization technique proposed here utilizes a series of coarse models (coarse-discretization EM models) {R c.j }, j= 1, . . . , K, evaluated by the same EM solver as fine model but faster. It is recommended that the number of models K is less than 4. The fidelity of these models can be adjusted by controlling the mesh density through model parameters (specific to a given EM solver, e.g., cells per wavelength or edge refinement, etc.). Therefore, the total number of simulated meshes of a component and its computational cost can be reduced. The first coarse model R c.1 , is the fastest and the least accurate. The last one R c.K , is closest to the fine model accuracy. There is an important trade-off between the computational complexity and accuracy of the coarse models, which greatly affects the performance of the surrogate-based optimization process. In particular, it is essential that the coarse model capture all important features present in the fine model (e.g., the general "shape" of responses etc.).
The surrogate models are constructed from simulation data of coarse models at various accuracy using an approximation-based response surface method. This method can be applied to the cases where the computational cost of coarse models cannot be ignored, such as the coarse-discretization EM model of complex microwave devices or antenna structures.
T be a certain design obtained in the ith surrogate model optimization built with the jth coarse model R c.j . x (0,1) is denoted as the initial design of the 1st surrogate model built using the 1st coarse model R c.1 . Instead of optimization R c.j directly, the ith surrogate model is built using the 2nd-order response surface model in the region containing x (i,j) , . Initial region size determined by d (0,1) usually equal 2% to 10% of x (0,1) , and its region center z (0,1) coincides with x (0,1) .
We construct the surrogate model using 2n+1 sample points of evaluated R c.j at z (i,j) and at boundary of the region [z (i,j) − d (i,j) , z (i,j) + d (i,j) ]. The coordinates of these samples are z where λ l.r , l = 1, . . . , m, r = 0, 1, . . . , 2n, is the rth coefficient of the 2nd-order response surface model. The model coefficients are obtained by solving the liner regression problem where R is an m-dimensional coarse model frequency response function vector over a specific frequency band (e.g., the return loss |S 11 |). Therefore, in building a surrogate, the coefficients are solved by (4) separately at each frequency point.

Surrogate Model Optimization with Adaptive Region Updating
This algorithm is an iterative process of local modeling and optimization. Therefore, the size and moving direction of each local region are two important factors, which affect the accuracy of the optimal design and the convergence speed of the algorithm. First, a local surrogate model s (0,1) is established using R c.1 and optimized in the region containing the initial design similar to [12], until the new local optimal design lies within the interior of the current region. At this time, we regard as the optimal design of the surrogate model at the lowest accuracy is found.
Then, this optimal design is used as the regional center of the next higher-accuracy local surrogate model (related to R c.2 ) and the new surrogate model is established and optimized using the above process. The process of establishing and optimizing the local surrogate model at the same accuracy is defined as one iteration of our algorithm.
The fidelity of the local surrogate increases with the number of iterations. In addition, the current optimal design serves as a starting point for the next iterative process. The size of local region at beginning of each iteration is adaptively updated according to a judgment factor. The judgment factor of the (i+1)th surrogate can be expressed as δ where β l can adjust the weight of the judgment factor at different frequencies if needed (here, β l = 1, l= 1, . . . , m). r (i+1) l representing the similarity between the surrogate model and coarse model at the lth frequency is defined for, where R l c.j is the lth component of the coarse model response vector R c.j . The judgment factor δ (i+1) r represents the overall degree of similarity (for all frequency points) between the current local surrogate model and its corresponding coarse model. If δ (i+1) r ≥ δ 1 (here, δ 1 =0.7), it means the current surrogate model is accurate, and the next local region size can be expanded , it indicates the current surrogate may not be accurate enough to reflect the true response of the current coarse model. Therefore, the next local region size is reduced (d (i+1,j) = d (i,j) /2). When δ 2 ≤ δ (i+1) r ≤ δ 1 , the next region size is kept the same as the current one and the starting point of the (i+1)th iteration is taken form the local optimal design found by the ith iteration.
The multi-fidelity surrogate models are built on various fidelity of the coarse models. Each iteration in the algorithm corresponds to a coarse model R c. j at specific fidelity. The algorithm converges in K iterations and K coarse models of different fidelity are used. Figure 1 shows the conceptual illustration of the process. The process in the figure exploits four coarse models:R c.1 , R c.2 , R c.3 , and R c.4 . For each iteration a particular color is used to mark local regions, their centers, and local optimal designs. The dashed line is used to mark the local regions in which the surrogate model is established. The solid dot z (i,j) is the center of the local region and the hollow circles is the current local optimal design x (i,j) . The red dot represents the optimal design of the surrogate model obtain after K = 4 iterations. accurate enough to reflect the true response of the current coarse model. Therefore, the next local ≤ , the next region size is kept the same as the current one and the starting point of the (i+1)th iteration is taken form the local optimal design found by the ith iteration. The multi-fidelity surrogate models are built on various fidelity of the coarse models. Each iteration in the algorithm corresponds to a coarse model is the center of the local region and the hollow circles is the current local optimal design ( , ) i j x . The red dot represents the optimal design of the surrogate model obtain after K=4 iterations. 1) x (1,1) z (2,2) =x (2,1) z (3,3) =x (3,2) x (4,3) z (5,4) =x (5,3) x (7,4) z (4,3) x (6,4) z (1,1) z (6,4) ). It is optimized to yield a new design (1,1) x , where the next surrogate model (1,1) s is established and centered at (1,1) z . After (2,2) is found, the algorithm enters the second iteration and the next surrogate model is built on

Refinement of Surrogate with Fine Model
The optimized design * s x of the last local surrogate model based on .
c K R (closest to the fine model accuracy) is a good estimation of the optimal design of the fine (high-fidelity) model. To improve the accuracy of the optimal design for the fine model, the last surrogate model is further refined. The refinement process is performed in the local region centered at * s x . The surrogate model (3) and (4).
In order to improve the surrogate model * ( ) s  , an input space mapping method [14] is applied.
A vector * c is obtained to improve the matching between * s and fine model R in the parameter . It is optimized to yield a new design x (1,1) , where the next surrogate model s (1,1) is established and centered at z (1,1) . After z (2,2) = x (2,1) is found, the algorithm enters the second iteration and the next surrogate model is built on R c.2 with an enlarged local region size assuming a judgement factor δ r > 0.7. The process continues and terminates after 4 iterations at x * s = x (7,4) .

Refinement of Surrogate with Fine Model
The optimized design x * s of the last local surrogate model based on R c.K (closest to the fine model accuracy) is a good estimation of the optimal design of the fine (high-fidelity) model. To improve the accuracy of the optimal design for the fine model, the last surrogate model is further refined. The refinement process is performed in the local region centered at x * s . The surrogate model s * (·) is established using coarse model R c.K within this region following (3) and (4).
In order to improve the surrogate model s * (·), an input space mapping method [14] is applied. A vector c * is obtained to improve the matching between s * and fine model R f in the parameter extraction process Then an output space mapping [15] correction term is added to the surrogate model to reduce the zeroth-order error [16]. The refined optimal design is then found as where d is the same as the last local region in Section 2.3. The refinement for local surrogate requires only one evaluation of R f .

Summary of Multi-Fidelity Local Surrogate Model Optimization
Based on the above description, Figure 2 summarizes the complete flowchart of the proposed algorithm.
Sensors 2019, 19, x FOR PEER REVIEW 6 of 13 Then an output space mapping [15] correction term is added to the surrogate model to reduce the zeroth-order error [16]. The refined optimal design is then found as ( ) where d is the same as the last local region in Section 2.3. The refinement for local surrogate requires only one evaluation of f R .

Summary of Multi-Fidelity Local Surrogate Model Optimization
Based on the above description, Figure 2 summarizes the complete flowchart of the proposed algorithm.

Verification Examples and Comparison
In this section, we present the numerical verifications of the multi-fidelity local surrogate model optimization algorithm in Section 2. Our proposed algorithm is applied to design EM-based microwave component models. A microstrip filter, an UWB monopole antenna, and a compact UWB MIMO antenna are considered. We compare the performance of the proposed method with that of the original algorithm [12] and direct fine model optimization (genetic algorithm). To make a fair comparison, we used a moderate mesh size for c R when testing the original algorithm. For example, the second coarse model

Verification Examples and Comparison
In this section, we present the numerical verifications of the multi-fidelity local surrogate model optimization algorithm in Section 2. Our proposed algorithm is applied to design EM-based microwave component models. A microstrip filter, an UWB monopole antenna, and a compact UWB MIMO antenna are considered. We compare the performance of the proposed method with that of the original algorithm [12] and direct fine model optimization (genetic algorithm). To make a fair comparison, we used a moderate mesh size for R c when testing the original algorithm. For example, the second coarse model R c.2 of the proposed algorithm is used as R c for the original algorithm in our comparison. To determine the mesh size of the coarse model R c.j in the proposed algorithm, we performed a preliminary test. It is found that when the mesh size of the coarse model is set so that the number of cells per wavelength are less than 10, the simulation response of the coarse model cannot represent the response characteristics of a fine model, therefore cannot be used. Generally, three to five iterations (or number of coarse models) are needed in the proposed algorithm. The mesh sizes of coarse models for R c.1 , R c.2 , and R c.3 are set to 10, 15, and 20 cells per wavelength respectively in our examples.

A Microstrip Filter Optimization
The first example is a second-order ring resonator microstrip bandpass filter, which is proposed in [17]. The geometry of the microstrip filter is shown in Figure 3. There are five designable parameters:

A Microstrip Filter Optimization
The first example is a second-order ring resonator microstrip bandpass filter, which is proposed in [17]. The geometry of the microstrip filter is shown in Figure 3. There are five designable parameters: x =[L1, L2, S1, W1, W2] T mm. The initial design is where w1, w2, and w3 are weighting factors for different frequency bands (here, 1 2 =1, =10, w w and is the fine model frequency response value at ω GHz. The filter structure is simulated using Computer Simulation Technology (CST) Microwave Studio transient solver with Intel(R) Core(TM) i7-4790 processor and 16 GB of RAM. The fine model f R has 253,561 meshes and its evaluation time is 5 min. To optimize the filter using multi-fidelity local surrogate algorithm of Section 2, we use two low-fidelity models: where w 1 , w 2 , and w 3 are weighting factors for different frequency bands (here, w 1 = 1, w 2 = 10, and w 3 = 5). And R f (x,ω) is the fine model frequency response value at ω GHz.
The filter structure is simulated using Computer Simulation Technology (CST) Microwave Studio transient solver with Intel(R) Core(TM) i7-4790 processor and 16 GB of RAM. The fine model R f has 253,561 meshes and its evaluation time is 5 min. To optimize the filter using multi-fidelity local surrogate algorithm of Section 2, we use two low-fidelity models: R c.1 (40,480 meshes at x (0,1) , evaluation time 0.6 min) and R c.2 (96,330 meshes, 1.3 min). As a comparison, the filter is also optimized using the original algorithm, in which the low-fidelity model R c is the same as R c.2 with 96,330 meshes and evaluation time of 1.3 min. Figure 4 shows the responses of fine model R f at x (0,1) and at its optimal designs x * = [19.81, 23.48, 0.32, 0.88, 1.48] T mm. In addition, the figure also compares the response of the optimal designs of the original algorithm (optimal design x 1 * = [18. 16    The optimization costs of the three algorithms are compared in Table 1. The original algorithm terminated after 5 local surrogates' optimization without finding the design satisfying the specifications. The direct optimization (genetic algorithm) obtains a design almost satisfying the specification, but it consumes a lot of computational cost (397 fine model evaluations). Our proposed algorithm obtains a design satisfying the specification and the total optimization cost is equivalent to only about 16 fine model f R evaluations. The computational saving of the proposed method is 95% or even more as compared to the direct optimization of the fine model.

Compact UWB MIMO Antenna Design
The second example is to design a compact UWB MIMO antenna using the proposed algorithm. In the design process, two optimization steps. In the first step, a compact UWB monopole antenna is The optimization costs of the three algorithms are compared in Table 1. The original algorithm terminated after 5 local surrogates' optimization without finding the design satisfying the specifications. The direct optimization (genetic algorithm) obtains a design almost satisfying the specification, but it consumes a lot of computational cost (397 fine model evaluations). Our proposed algorithm obtains a design satisfying the specification and the total optimization cost is equivalent to only about 16 fine model R f evaluations. The computational saving of the proposed method is 95% or even more as compared to the direct optimization of the fine model.

Compact UWB MIMO Antenna Design
The second example is to design a compact UWB MIMO antenna using the proposed algorithm. In the design process, two optimization steps. In the first step, a compact UWB monopole antenna is optimized using our proposed method. In the second step, a MIMO antenna (constructed using two optimized UWB monopole antennas) is optimized using a simplify version of our proposed method.

UWB Monopole Antenna Optimization
First, the monopole antenna structure shown in Figure 5 is optimized. The structure consists of a rectangular radiating plate fed through a microstrip line and a ground plane with a small rectangular slot and a long ground stub. According to Reference [18], the small rectangular slot with a size of s 1 × s 2 is for better impedance matching at high frequencies and the long stub for isolation enhancement is bent to reduce the overall antenna area. It is designed on a Rogers substrate, RO4350B (ε r = 3.55, tan δ = 0.0037, and h = 0.8 mm).
Sensors 2019, 19, x FOR PEER REVIEW 9 of 13 optimized using our proposed method. In the second step, a MIMO antenna (constructed using two optimized UWB monopole antennas) is optimized using a simplify version of our proposed method.

UWB Monopole Antenna Optimization
First, the monopole antenna structure shown in Figure 5 is optimized. The structure consists of a rectangular radiating plate fed through a microstrip line and a ground plane with a small rectangular slot and a long ground stub. According to Reference [18], the small rectangular slot with a size of 12 ss  is for better impedance matching at high frequencies and the long stub for isolation enhancement is bent to reduce the overall antenna area. It is designed on a Rogers substrate, RO4350B ( The fine model f R is evaluated using the CST (1168614 mesh cells, evaluation time is 6.5 min). We optimized this example using three low-fidelity models:   The optimization goal is to meet the following specifications: reflection coefficients 11 12dB S  for 3.1 GHz to 10.6 GHz. In addition, the antenna footprint is also considered while maintaining the above electrical performance. The size reduction is selected as the main optimization problem, and good electrical performance is achieved in the form of a penalty function. Thus, the local region optimization according to the formula (2) is defined as where  is a penalty factor (here, =1000  ) and ( , ) () ij s is the ith multi-fidelity surrogate model built using After the optimal solution of the final local surrogate model is found, the surrogate model is refined. The optimal design * [11.2, 0.7, 8.85, 8.77, 13.8, 5.06, 1.01, 2.3, 1.9, 6.96, 1.71] T  x obtained by In this case, 11 design variables x = [l g , g, a 1 , a 2 , l 1 , l 2 , w 1 , s 1 , s 2 , o 1 , o 3 ] T mm are considered. The feedline width is fixed at w 0 = 1.8 mm. The initial design is x (0,1) = [7, 1, 10, 10, 11, 3, 1, 3, 1, 6, 1] T mm. The size of the initial local region is defined as d (0,1) = [1, 1, 0.1, 0.1, 1, 1, 0.1, 1, 0.1, 1, 0 .1] T mm. The fine model R f is evaluated using the CST (1168614 mesh cells, evaluation time is 6.5 min). We optimized this example using three low-fidelity models: R c.1 (134560 meshes, 1 min), R c.2 (187110 meshes, 1.5 min), and R c.3 (338680 meshes, 2.8 min). As a comparison, this example is also optimized using the original algorithm, in which the low-fidelity model R c is the same as R c.2 with 187110 meshes and evaluation time of 1.5 min.
The optimization goal is to meet the following specifications: reflection coefficients |S 11 | ≤ −12 dB for 3.1 GHz to 10.6 GHz. In addition, the antenna footprint S(x) = (o 1 + 0.5 · a 2 + o 3 ) × (l g + l 1 + w 1 ) is also considered while maintaining the above electrical performance. The size reduction is selected as the main optimization problem, and good electrical performance is achieved in the form of a penalty function. Thus, the local region optimization according to the formula (2) is defined as where β is a penalty factor (here, β = 1000) and s (i,j) (·) is the ith multi-fidelity surrogate model built using R c. j , (j = 1, 2, 3). And c is a penalty function defined as c(s (i,j) (x)) = max max |S 11 | 3.1GHz to 10.6GHz + 12, 0 .
After the optimal solution of the final local surrogate model is found, the surrogate model is refined. The optimal design x * = [11.2, 0.7, 8.85, 8.77, 13.8, 5.06, 1.01, 2.3, 1.9, 6.96, 1.71] T obtained by the refined local surrogate model (8) is very close to that of the fine model and they both satisfy the specification, as shown in Figure 6. the refined local surrogate model (8) is very close to that of the fine model and they both satisfy the specification, as shown in Figure 6.   ， is also shown in Figure 6 (three runs of genetic algorithm are performed and none of the results satisfies the design specification). The optimization costs of UWB monopole antenna using the three algorithms are compared in Table 2. The original algorithm terminated after 4.2 h without finding the design satisfying the  Figure 7 shows the response of fine model R f at initial design x (0,1) and the optimal design obtained x * by the proposed algorithm. The corresponding antenna size is 13 mm × 26 mm (338 mm 2 ). In addition, the figure also shows results obtained by the original algorithm and the direct optimization algorithm (genetic algorithm). It can be seen that the optimal design x *  Figure 6 (three runs of genetic algorithm are performed and none of the results satisfies the design specification). the refined local surrogate model (8) is very close to that of the fine model and they both satisfy the specification, as shown in Figure 6.   ， is also shown in Figure 6 (three runs of genetic algorithm are performed and none of the results satisfies the design specification). The optimization costs of UWB monopole antenna using the three algorithms are compared in Table 2. The original algorithm terminated after 4.2 h without finding the design satisfying the The optimization costs of UWB monopole antenna using the three algorithms are compared in Table 2. The original algorithm terminated after 4.2 h without finding the design satisfying the specifications. The direct optimization (genetic algorithm) is not able to obtain a design satisfying the specification after an average of 53.8 h (CST genetic algorithm terminates after 497 fine model evaluations by default). Our proposed algorithm obtains a design satisfying the specification and the total optimization cost is 5.6 h.

Compact UWB MIMO Antenna
We construct a UWB MIMO antenna using two radiating patches (UWB monopole antenna element) placed perpendicularly. Figure 8 shows the geometry of the compact UWB MIMO antenna. A thin ground strip is used to electrically connect the two ground planes together. The ground strip width and length affect lower cutoff frequency and the isolation. evaluations by default). Our proposed algorithm obtains a design satisfying the specification and the total optimization cost is 5.6 h. We construct a UWB MIMO antenna using two radiating patches (UWB monopole antenna element) placed perpendicularly. Figure 8 shows the geometry of the compact UWB MIMO antenna. A thin ground strip is used to electrically connect the two ground planes together. The ground strip width and length affect lower cutoff frequency and the isolation.
The optimal design x * in the Section 3.  The MIMO antenna is complex and only the high-fidelity EM simulation is reliable. The number of design variables is only two. For these reasons, a simplified version of our multi-fidelity local surrogate model optimization algorithm is used. Instead of multiple coarse models, the algorithm directly use the response of the fine model to build the local surrogate models, The initial design is  Figure 9 shows the response of f R at the optimal design * x . As can be seen, the return loss 11 S (or 22 S ) of the MIMO antenna is less than −12 dB The optimal design x * in the Section 3.2.1 is used as the geometric parameter of the two UWB monopole antenna elements in the MIMO antenna. At this time, the design variable x = [l c w c ] T of the MIMO antenna are optimized.
The MIMO antenna is complex and only the high-fidelity EM simulation is reliable. The number of design variables is only two. For these reasons, a simplified version of our multi-fidelity local surrogate model optimization algorithm is used. Instead of multiple coarse models, the algorithm directly use the response of the fine model to build the local surrogate models, R c.j = R f , j = 1, . . . , K. The initial design is x (0,1) = [3, 3] T mm. The specifications are mutual coupling |S 21 | ≤ −15 dB for 3.1 GHz to 10.6 GHz. The optimal design x * = [4.68 2.21] T mm is obtained by establishing and optimizing three local surrogate models iteratively. Figure 9 shows the response of R f at the optimal design x * . As can be seen, the return loss |S 11 | (or |S 22 |) of the MIMO antenna is less than −12 dB for 3.1 GHz to 10.6 GHz. Since the geometry parameter values of the UWB monopole antenna elements are the same, the return loss response curves at the two ports are similar. The mutual coupling of less than −15 dB between the input ports is achieve (the compact UWB MIMO antenna obtained high isolation). for 3.1 GHz to 10.6 GHz. Since the geometry parameter values of the UWB monopole antenna elements are the same, the return loss response curves at the two ports are similar. The mutual coupling of less than −15 dB between the input ports is achieve (the compact UWB MIMO antenna obtained high isolation).

Conclusions
A novel multi-fidelity local surrogate model optimization approach for microwave component design is presented. In this optimization process, multiple local surrogate models are established and optimized iteratively. It exploits a series of coarse models with variable fidelity to construct the local surrogate models with sequentially increased accuracy. The proposed algorithm introduces the judgment factor to implement adaptive local surrogate model optimization. The algorithm converges quickly and efficiently to the optimal design of fine model. This is demonstrated through a microstrip filter and a compact UWB MIMO antenna design examples Compared with the original local surrogate model optimization algorithm, the accuracy of the optimal design is improved when applying the proposed algorithm. Compared with the direct optimization algorithm, the proposed algorithm finds the optimal design efficiently.

Conclusions
A novel multi-fidelity local surrogate model optimization approach for microwave component design is presented. In this optimization process, multiple local surrogate models are established and optimized iteratively. It exploits a series of coarse models with variable fidelity to construct the local surrogate models with sequentially increased accuracy. The proposed algorithm introduces the judgment factor to implement adaptive local surrogate model optimization. The algorithm converges quickly and efficiently to the optimal design of fine model. This is demonstrated through a microstrip filter and a compact UWB MIMO antenna design examples Compared with the original local surrogate model optimization algorithm, the accuracy of the optimal design is improved when applying the proposed algorithm. Compared with the direct optimization algorithm, the proposed algorithm finds the optimal design efficiently.