An Adaptive Interpolation Scheme for Wideband Frequency Sweep in Electromagnetic Simulations

by   Kai Zhu, et al.
Beihang University

An adaptive interpolation scheme is proposed to accurately calculate the wideband responses in electromagnetic simulations. In the proposed scheme, the sampling points are first carefully divided into several groups based on their responses to avoid the Runge phenomenon and the error fluctuations, and then different interpolation strategies are used to calculate the responses in the whole frequency band. If the relative error does not satisfy the predefined threshold in a specific frequency band, it will be refined until the error criteria is met. The detailed error analysis is also presented to verify the accuracy of the interpolation scheme. At last, two numerical examples including the antenna radiation and the filter simulation are carried out to validate its accuracy and efficiency.



There are no comments yet.


page 1

page 2

page 3

page 4


Fast Green Function Evaluation for Method of Moment

In this letter, an approach to accelerate the matrix filling in method o...

Assessing Prediction Error at Interpolation and Extrapolation Points

Common model selection criteria, such as AIC and its variants, are based...

Error constant estimation under the maximum norm for linear Lagrange interpolation

For the Lagrange interpolation over a triangular domain, we propose an e...

A shape preserving quasi-interpolation operator based on a new transcendental RBF

It is well-known that the univariate Multiquadric quasi-interpolation op...

Spatially Adaptive Inference with Stochastic Feature Sampling and Interpolation

In the feature maps of CNNs, there commonly exists considerable spatial ...

A Novel Computational Method for Band Structures of Dispersive Photonic Crystals

We propose a new method to compute band structures of dispersive photoni...

Rational-based model order reduction of Helmholtz frequency response problems with adaptive finite element snapshots

We introduce several spatially adaptive model order reduction approaches...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

I Introduction

With the fast development of high-speed electronic systems in the past decades, the demand for accurate and efficient electromagnetic simulations keeps increasing. The finite element method (FEM) [1]

and the method of moment (MOM)


in the frequency domain have been widely used in electrical engineering applications, such as antenna design

[3], passive component design [4]. Unlike the finite-difference time-domain (FDTD) [5] methods, which can obtain a wideband response through a single simulation, the frequency sweep is required if the wideband responses are desired, such as microwave imaging [6], parameter extraction in integrated circuit design [7].

However, the wideband frequency responses usually show large and sharp variations and fluctuations, as shown in Fig. 1(a), which is the S parameter of a waveguide from 1.7 GHz to 1.9 GHz. To obtain the wideband responses, an accurate interpolation scheme is preferred since only a small number of discrete frequency points are required and significant efficiency improvement can be obtained. However, the traditional interpolation approaches usually suffer from numerical instability since many discontinuities may exist in the interested frequency band, and the increasing degree of the interpolation polynomial leads to the so-called Runge phenomenon [8]. It can be suppressed if the Chebyshev sampling scheme is applied [9]. However, the Chebyshev points with high non-uniformity may not be desirable in many applications.

The model order reduction (MOR) technique can reduce the computational complexity in electromagnetic simulations. One typical approach is the asymptotic waveform evaluation (AWE) [10] based on the Taylor series and the Pade approximation. In [11]

, the multipoint MOR was proposed to reduce the number of basis functions. However, the MOR methods suffer from accuracy issue and ill-conditioned matrix calculation


In this letter, an adaptive interpolation scheme is proposed to accurately and efficiently predict the wideband responses in electromagnetic simulations. A small number of frequency points in the frequency band is first selected, and their responses are calculated through field solvers. Then, we adaptively cast the sampling points into several groups to avoid numerical instability due to discontinuities. The Lagrange interpolation polynomial in each group is then constructed. If the relative error is larger than the predefined threshold in a specific frequency band, more sampling points are adaptively selected in the corresponding band until convergent results are obtained. The main merit of the proposed scheme is independent of the complexity of the matrix reduction and applicable for extremely complicated wideband responses as shown in the numerical results.

This paper is organized as follows. In Section II, the proposed adaptive interpolation scheme is detailed presented. In Section III, its rigorous error analysis is discussed. Two numerical examples are carried out to demonstrate its effectiveness in Section IV . Finally, some conclusions are drawn in Section V.

Ii Methodology

Ii-a Numerical Issues in the Wideband Interpolation

To obtain accurate responses in a wide frequency band, the sampling points can be divided into several groups to obtain an effective interpolation scheme as shown in Fig. 1(a). Since p5 is an extreme point, an ideal division is that p1-p5 is casted into a group, and p5-p9 as another group. Then, the polynomial interpolation can be used to calculate the responses in the whole frequency band. However, it is not always applicable for complex scenarios. There are three issues in the practical applications: (1) If there are too many points in a group and the polynomial interpolation is used, the Runge phenomenon may happen. (2) Several extreme points are adjacent to each other when insufficient sampling occurs. (3) The extreme points may not be sampled and other points may be erroneously regarded as the extreme points. To overcome these issues, an adaptive sampling scheme is proposed in this section.

To make our explanation clear, the th group denotes , which has sampling points, and the distance between the th and th groups is . is the desired frequency point and its response is . The non-single group is defined as a group including more than one points, and the single group contains only one isolated or extreme point. and are the number of non-single and single groups, respectively. The single group closest to is , and denote the latter and former groups, respectively. is the previous point of . The superscript indicates the th point in the group.

(a) (b)
Fig. 1: Illustration of selecting the sampling points to calculate the S parameter in the frequency band for a waveguide, (a) the ideal group subdivision scheme and (b) the proposed group subdivision scheme.

Ii-B The Proposed Adaptive Interpolation Scheme

The proposed interpolation scheme can be divided into four phases. First, the responses of the uniform sampling points are calculated at the beginning, which are then divided into several groups according to the extreme points. Second, the Lagrange interpolation is applied in each group, and the frequency response in the whole frequency range are analyzed. Third, the frequency band is divided into several parts, the relative error in each part is evaluated. Finally, more sampling points are carefully selected in the specific frequency bands, where the relative errors are larger than the predefined values, and then they are adaptively casted into several groups according to their responses. This procedure will repeat until the convergent results are obtained.

Ii-B1 Adaptive Sampling and Classification

The interested frequency band is first divided into several parts for evaluating the errors. Then, we select a small number of sampling points and find the extreme points according to their responses. They are then divided into several groups. In order to reduce the error in the whole frequency band, the maximum number of points in each group is limited to no more than 3. Fig. 1(b) presents an example with four groups including two non-single and single groups.

Ii-B2 Interpolation

In our implementation, the Lagrange linear and quadratic interpolation are used. presents the procedure of calculating the approximate response at a desired frequency point. The group type should be checked before the interpolation. Since non-single groups do not include extreme points, the interpolation results are usually more accurate and stable. In order to ensure the convergence, the points with the distance not greater than at the endpoints are considered to be in the group .

3:for  do
4:     if  then
5:          isfind = true;
6:           = LagrangeInterpolation();
7:     end if
8:end for
9:if isfind = false then
10:     for  do
11:          if  then
12:               mindistance = ;
13:               ;
14:          end if
15:     end for
16:      = AdaptiveInterpolation();
17:end if
Algorithm 1 Interpolation

For these points, which are not in the non-single groups, is used to construct the adaptive interpolation scheme, which are demonstrated in . presents the interpolation strategy for the scenario, where contains only an extreme point and is in its positive direction. For the points in its negative direction, a similar strategy can be used with only replacing by and changing with . is used to record the chosen strategy for subsequent optimization. When does not include an extreme point, the procedure is introduced in . In this scenario, the former or latter group includes an extreme point.

3:if  then
4:     if  then
5:          if  then
6:                = LagrangeInterpolation(
7:        );
8:          else
9:                = LagrangeInterpolation(
10:        );
11:                = 2;
12:          end if
13:     else
14:          if  then
15:                = LinearInterpolation();
16:                = 3;
17:          else
18:                = LinearInterpolation();
19:          end if
20:     end if
21:end if
Algorithm 2 Scenario #1: Adaptive Interplation
3:if  then
4:     if  then
5:          if  then
6:                = LagrangeInterpolation(
7:        );
8:          else
9:                = LagrangeInterpolation(
10:        );
11:          end if
12:     else
13:          if  then
14:               if  then
15:                     = LagrangeInterpolation(
16:        );
17:               else
18:                     = LagrangeInterpolation(
19:        );
20:               end if
21:          else
22:                = LinearInterpolation(,,);
23:          end if
24:     end if
25:end if
Algorithm 3 Scenario #2: Adaptive Interpolation

Another essential step is to optimize the previously calculated results since larger errors occur near the extreme points. The responses in the adjacent non-single group are used to optimize the results near the extreme points to improve the accuracy. For instance, if is in a non-single group and , we can recalculate the Lagrange polynomials with to get the corrected value at . As for the linear interpolation, take for example, if , is calculated through the linear extrapolation of and . Otherwise, the linear extrapolation of and .

Ii-B3 Error Estimate

If we obtain the responses of frequency points, to quantitatively measure the performance, the relative error is defined as , where is the previously approximated results, is the refined results. The error is used to control the iteration procedure in the proposed scheme.

Ii-B4 Error Control

After completing one interpolation process, the relative errors in each frequency range are calculated. In order to increase the convergence and reduce the computational overhead, only the frequency bands, which have larger relative error than the predefined threshold, are refined until the convergent results are obtained.

Iii Error Analysis

If we have the sampling points together with frequency responses . The Lagrange’s interpolation formulation can be expressed as


where is the Lagrange polynomial associated to , which can be expressed as .

Using the Taylor series analysis, the truncation error can be expressed by


where is the th derivative of at , and . If , the truncation error reduces to zero. We denote . If the sampling points are even and the interval is , is bounded by [13].

By using the differentiation theorem for the Fourier transforms

[14] and applying the Schwartz inequality, we obtain


where , and denotes the cut-off frequency. The first term is bounded by


The second term is known as the finite energy of , which can be represented by . By substituting (4) into (III), we obtain


By combining (5) and (2), we obtain


Therefore, if sufficient frequency points are selected, where the interval meets , and the energy of is finite, the convergent results can be obtained. In addition, as the number of sampling points increases, the error becomes smaller.

Iv Numerical Results and Discussion

Iv-a Antenna Radiation

The first numerical example is to calculate the responses of the Potter Horn Antenna, which is a dual-mode feedhorn and provides excellent radiation pattern with suppressed sidelobes and symmetrical beamwidths [15]. The frequency band is 27 GHz-33 GHz. If the frequency band is divided into 70 parts and the predefined relative error in each part is set to 0.03, 130 frequency points are selected to calculate and the responses in other frequency points can be evaluated through the proposed interpolation scheme. The average error in the whole frequency region is converged to 0.0019. Fig. 2 shows the results obtained from the HFSS, the proposed scheme in this case. It is easy to find that the results obtained from the proposed interpolation scheme show excellent agreement with those obtained from the HFSS both in the discrete and interpolation sweep modes. The frequency interval is 0.01 GHz and the number of points is 601 in the discrete mode. However only 130 points need to be calculated directly in the proposed scheme, which accounts for 21.6% of the number of frequency points in the discrete mode.

Fig. 2: The S Parameter obtained from the HFSS in the discrete mode, the HFSS in the interpolation mode, and the proposed interpolation scheme.

Fig. 3 (a) shows the interpolation error and the number of points versus the number of parts, the error in each part is set to 0.03. As the frequency band divided into different parts, different number of points are needed, with more points calculated, the smaller errors are usually obtained. In this numerical example, if the frequency range is divided into 70 parts, it can be found that the smallest error is achieved, with 130 points are used. Fig. 3 (b) shows the interpolation error and the number of points versus the error limit. It can be found that as the predefined error increases, the number of points required to be calculated decreases.

(a) (b)
Fig. 3: (a) The error and the number of points with the error limit in each part set to 0.03, (b) the error and the number of points with the frequency band divided into 70 parts.

Iv-B The S Parameters of A Filter

In order to further validate the accuracy of the proposed scheme with large variations, we consider an interdigital bandpass filter with a 1 GHz bandwidth. The S parameter from 0.6 GHz to 2.4 GHz is considered, and 601 points with an interval of 0.003 GHz is used in our simulations. When the frequency band is divided into 70 parts and the relative error in each part is limited to 0.05, only 129 discrete points are required to be calculated in the proposed interpolation scheme. It accounts for only 21.5% of the desired frequency points. The global error in the whole frequency band converges to 0.012. Fig. 4 shows the results obtained from the HFSS and the proposed interpolation scheme. Similar to the previous numerical example, as a result of the adaptive sampling strategy, the extreme points can be selected and the proposed scheme shows excellent agreement with those from the HFSS in the discrete and interpolation modes.

Fig. 4: The S Parameter obtained from the HFSS in the discrete mode, the HFSS in the interpolation mode, and the proposed interpolation scheme.

Fig. 5 (a) shows the interpolation error and the number of points with the frequency band divided into different parts when the relative error in each part is limited to 0.03. In Fig. 5 (b), the interpolation results are listed with the frequency band divided into 70 parts and the relative error in each part limited to be different values. As we can see, the observations are similar to the previous example. As the relative error increases, less points are required to be calculated.

(a) (b)
Fig. 5: (a) The error and the number of points with the error limit in each part set to 0.03, (b) the error and the number of points with the frequency band divided into 70 parts.

V Conclusion

An adaptive interpolation scheme with a controllable error for fast frequency sweep in a wide frequency band is proposed. In the proposed scheme, we only need to calculate a small number of frequency points, and then adaptively use different interpolation strategies for other required frequency points. Numerical examples show that the proposed scheme can produce very accurate results and reveals its fast frequency sweep capability. Moreover, the proposed scheme can be used to calculate other complex responses in a wide variable range without being limited to the specific problems.


  • [1] J. Jin, The finite element method in electromagnetics. New York, NY, USA:Wiley, 2015.
  • [2] W. Gibson, The Method of Moments in Electromagnetics. CRC press, 2015.
  • [3] A. Barka, A. Jouadé, and D. Jacquinot, “Simulation of Active Reflection Coefficient Phenomena of Large Antenna Array Using FEM Domain Decomposition Methods,” IEEE Antennas Wireless Propag. Lett., vol. 19, no. 10, pp. 1789-1792, Oct. 2020.
  • [4] T. Lu and J. Jin, “Electrical-Thermal Co-Simulation for Analysis of HighPower RF/Microwave Components,” IEEE Trans. Electromagn. Compat., vol. 59, no. 1, pp. 93-102, Sep. 2016.
  • [5] A. Taflove and S. Hagness, Computational Electrodynamics: the Finite- Difference Time-Domain Method. Artech house, 2005.
  • [6] R. Chandra and H. Zhou, “On the Opportunities and Challenges in Microwave Medical Sensing and Imaging,” IEEE Trans. Biomed. Eng., vol. 62, no. 7, pp. 1667-1682, Jul. 2015.
  • [7] Z. Gao and K. Kang, “Analysis and Equivalent-Circuit Model for CMOS On-Chip Multiple Coupled Inductors in the Millimeter-Wave Region,” IEEE Trans. Electron Devices, vol. 62, no. 12, pp. 3957-3964, Dec. 2015.
  • [8] M. Grasselli and D. Pelinovsky, Numerical Mathematics, CA:Jones & Bartlett Publishers, Inc., 2008.
  • [9] J. P. Boyd and J. R. Ong, “Exponentially-convergent strategies for defeating the Runge phenomenon for the approximation of non-periodic functions, part I: single-interval schemes,” Commun. Comput. Phys., vol. 5, no. 2-4, pp. 484-497, 2009.
  • [10] Y. Jeong, I. Hong, and K. Lee, “Fast frequency sweep using asymptotic waveform evaluation technique and thin dielectric sheet approximation,” IEEE Trans. Antennas Propag., vol. 64, no. 5, pp. 1800-1806, May 2016.
  • [11] S. V. Polstyanko, R. Dyczij-Edlinger, and J. F. Lee, “Fast frequency sweep technique for the efficient analysis of dielectric waveguides,” IEEE Trans. Microw. Theory Tech., vol. 45, no. 7, pp. 1118–1126, Jul. 1997.
  • [12] Z. Song, D. Su, and F. Duval, “Model order reduction for PEEC modeling based on moment matching,” Progress In Electromagnetics Research, vol. 114, pp. 285-299, 2011.
  • [13] R. Radzyner and P. Bason, “An error bound for Lagrange interpolation of low-pass functions,” IEEE Trans. Inf. Theory, vol. 18, no. 9, pp. 669–671, Sep. 1972.
  • [14] J. S. Walker, Fast fourier transforms. CRC press, 2017.
  • [15] H. Wang, M. Yi, and G. Liu, “Potter horn and compact orthomode transducer at 150 GHz,” IEEE Trans. Antennas Propag., vol. 62, no. 10, pp. 5403–5405, Oct. 2014.