Fast solution of electromagnetic integral equations using adaptive wavelet packet transform

The adaptive wavelet packet transform is applied to sparsify the moment matrices for the fast solution of electromagnetic integral equations. In the algorithm, a cost function is employed to adaptively select the optimal wavelet packet expansion/testing functions to achieve the maximum sparsity possible in the resulting transformed system. The search for the best wavelet packet basis and the moment matrix transformation are implemented by repeated two-channel filtering of the original moment matrix with a pair of quadrature filters. It is found that the sparsified matrix has above-threshold elements that grow only as O(N/sup 1.4/) for typical scatterers. Consequently the operations to solve the transformed moment equation using the conjugate gradient method scales as O(N/sup 1.4/). The additional computational cost for carrying out the adaptive wavelet packet transform is evaluated and discussed.

[1]  I. Daubechies.  Orthonormal bases of compactly supported wavelets , 1988 .

[2]  R. Coifman,et al.  Fast wavelet transforms and numerical algorithms I , 1991 .

[3]  Jorge Nocedal,et al.  Global Convergence Properties of Conjugate Gradient Methods for Optimization , 1992, SIAM J. Optim..

[4]  Ronald R. Coifman,et al.  Entropy-based algorithms for best basis selection , 1992, IEEE Trans. Inf. Theory.

[5]  Charles K. Chui,et al.  An Introduction to Wavelets , 1992 .

[6]  Ronald R. Coifman,et al.  Wavelet-Like Bases for the Fast Solution of Second-Kind Integral Equations , 1993, SIAM J. Sci. Comput..

[7]  Hao Ling,et al.  On the application of fast wavelet transform to the integral‐equation solution of electromagnetic scattering problems , 1993 .

[8]  Y. Leviatan,et al.  On the use of wavelet expansions in the method of moments (EM scattering) , 1993 .

[9]  Linda P. B. Katehi,et al.  Analysis of integrated millimeter-wave and submillimeter-wave waveguides using orthonormal wavelet expansions , 1994 .

[10]  Weng Cho Chew,et al.  A ray‐propagation fast multipole algorithm , 1994 .

[11]  M. Victor Wickerhauser,et al.  Adapted wavelet analysis from theory to software , 1994 .

[12]  L. B. Montefusco.  Parallel Numerical Algorithms with Orthonormal Wavelet Packet Bases , 1994 .

[13]  Andrew K. Chan,et al.  Linear frequency modulated signal detection using wavelet packet, ambiguity function and Radon transform , 1995, IEEE Antennas and Propagation Society International Symposium. 1995 Digest.

[14]  Guang-Wen Pan,et al.  Full wave analysis of microstrip floating line structures by wavelet expansion method , 1995 .

[15]  Weng Cho Chew,et al.  A study of wavelets for the solution of electromagnetic integral equations , 1995 .

[16]  Hao Ling,et al.  A fast moment method algorithm using spectral domain wavelet concepts , 1996 .

[17]  Yehuda Leviatan,et al.  Impedance matrix compression using adaptively constructed basis functions , 1996 .

[18]  Hao Ling,et al.  Moment matrix sparsification using adaptive wavelet packet transform , 1997 .

[19]  W. L. Golik.  Wavelet packets for fast solution of electromagnetic integral equations , 1998 .