Estimation of resonant frequencies and quality factors from time domain computations

A new frequency domain subspace algorithm for the estimation of resonant frequencies and quality factors given the output from time domain computations is developed and tested. It applies the discrete Fourier transform to the time domain signal and the estimates are computed from the frequency data within the frequency band of interest. Thus, out-of-band disturbances are suppressed through a band-pass filtering operation in frequency domain, which is free from the transient effects that are inevitable for standard time domain filtering techniques. We also propose an efficient procedure to generate short time domain excitations with an a priori specified frequency spectrum. By combining such excitations with the frequency domain subspace algorithm, we estimated to high accuracy the lowest 167 resonant frequencies of a three-dimensional bow-tie microwave resonator given two separate time domain computations with, in total, 23,000 time steps. A free-space computation for a cavity with an aperture shows that the estimation algorithm can be used to efficiently and accurately extrapolate time domain signals. The tests show that it is possible to reduce the overall computation time by several orders of magnitude for systems with undamped or weakly damped modes. It operates as a post-processing step and, thus, decouples the estimation of frequencies and amplitudes from the time domain solver which offers a considerable flexibility.

[1]  Mats Viberg,et al.  Subspace-based methods for the identification of linear time-invariant systems , 1995, Autom..

[2]  Steven Kay,et al.  Modern Spectral Estimation: Theory and Application , 1988 .

[3]  R. O. Schmidt,et al.  Multiple emitter location and signal Parameter estimation , 1986 .

[4]  Boualem Boashash,et al.  Estimating and interpreting the instantaneous frequency of a signal. I. Fundamentals , 1992, Proc. IEEE.

[5]  Jian-Ming Jin,et al.  The Finite Element Method in Electromagnetics , 1993 .

[6]  Jin-Fa Lee,et al.  Time-domain finite-element methods , 1997 .

[7]  John Litva,et al.  Fast finite-difference time-domain analysis of resonators using digital filtering and spectrum estimation techniques , 1992 .

[8]  J. Nédélec.  Mixed finite elements in ℝ3 , 1980 .

[9]  Josef A. Nossek,et al.  Unitary ESPRIT: how to obtain increased estimation accuracy with a reduced computational burden , 1995, IEEE Trans. Signal Process..

[10]  K. Yee.  Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media , 1966 .

[11]  Allen Taflove,et al.  Computational Electrodynamics the Finite-Difference Time-Domain Method , 1995 .

[12]  Raj Mittra,et al.  A combination of FD-TD and Prony's methods for analyzing microwave integrated circuits , 1991 .

[13]  Andreas C. Cangellaris,et al.  A Reflectionless Sponge Layer Absorbing Boundary Condition for the Solution of Maxwell's Equations with High-Order Staggered Finite Difference Schemes , 1998 .

[14]  Sun-Yuan Kung,et al.  A new identification and model reduction algorithm via singular value decomposition , 1978 .

[15]  Petre Stoica,et al.  Introduction to spectral analysis , 1997 .

[16]  William H. Press,et al.  The Art of Scientific Computing Second Edition , 1998 .

[17]  Jian-Ming Jin,et al.  Fast and Efficient Algorithms in Computational Electromagnetics , 2001 .

[18]  A. Vegas,et al.  Computation of resonant frequencies and quality factors of open dielectric resonators by a combination of the finite-difference time-domain (FDTD) and Prony's methods , 1992, IEEE Microwave and Guided Wave Letters.

[19]  T. McKelvey,et al.  A robust frequency domain subspace algorithm for multi-component harmonic retrieval , 2001, Conference Record of Thirty-Fifth Asilomar Conference on Signals, Systems and Computers (Cat.No.01CH37256).

[20]  Thomas Kailath,et al.  ESPRIT-estimation of signal parameters via rotational invariance techniques , 1989, IEEE Trans. Acoust. Speech Signal Process..

[21]  Patrick Dewilde,et al.  Subspace model identification Part 1. The output-error state-space model identification class of algorithms , 1992 .

[22]  R. Mittra,et al.  Efficient computation of resonant frequencies and quality factors of cavities via a combination of the finite-difference time-domain technique and the Pade approximation , 1998 .

[23]  Krishna Naishadham,et al.  ARMA-based time-signature estimator for analyzing resonant structures by the FDTD method , 2001 .

[24]  Feng Liu,et al.  Full-wave analysis and modeling of multiconductor transmission lines via 2-D-FDTD and signal-processing techniques , 2002 .

[25]  Thomas Rylander,et al.  Stable FEM-FDTD hybrid method for Maxwell's equations , 2000 .

[26]  Chao Yang,et al.  ARPACK users' guide - solution of large-scale eigenvalue problems with implicitly restarted Arnoldi methods , 1998, Software, environments, tools.

[27]  L. Ljung,et al.  Subspace-based multivariable system identification from frequency response data , 1996, IEEE Trans. Autom. Control..

[28]  Gene H. Golub,et al.  Matrix computations , 1983 .