Principal Component Analysis of Results Obtained from Finite-Difference Time-Domain Algorithms

Finite-Differences Time-Domain (FDTD) algorithms are well established tools of computational electromagnetism. Because of their practical implementation as computer codes, they are affected by many numerical artefact and noise. In order to obtain better results we propose using Principal Component Analysis (PCA) based on multivariate statistical techniques. The PCA has been successfully used for the analysis of noise and spatial temporal structure in a sequence of images. It allows a straightforward discrimination between the numerical noise and the actual electromagnetic variables, and the quantitative estimation of their respective contributions. Besides, The GDTD results can be filtered to clean the effect of the noise. In this contribution we will show how the method can be applied to several FDTD simulations: the propagation of a pulse in vacuum, the analysis of two-dimensional photonic crystals. In this last case, PCA has revealed hidden electromagnetic structures related to actual modes of the photonic crystal.

[1]  José Manuel López-Alonso,et al.  Characterization of artifacts in fully digital image-acquisition systems: Application to web cameras , 2004 .

[2]  Javier Alda,et al.  Characterization of dynamic sea scenarios with infrared imagers , 2005 .

[3]  Javier Alda,et al.  Photonic crystal characterization by FDTD and principal component analysis. , 2004, Optics express.

[4]  Javier Alda,et al.  Principal-component characterization of noise for infrared images. , 2002, Applied optics.

[5]  Sacharia Albin,et al.  Numerical techniques for excitation and analysis of defect modes in photonic crystals. , 2003, Optics express.

[6]  D. F. Morrison,et al.  Multivariate Statistical Methods , 1968 .

[7]  Levent Gurel,et al.  Signal-processing techniques to reduce the sinusoidal steady-state error in the FDTD method , 2000 .

[8]  R. Cattell.  The Scree Test For The Number Of Factors. , 1966, Multivariate behavioral research.

[9]  Javier Alda,et al.  Operational parametrization of the 1/f noise of a sequence of frames by means of the principal component analysis in focal plane arrays , 2003 .

[10]  Javier Alda,et al.  Multivariate analysis of photonic crystal microcavities with fabrication defects , 2005, SPIE Microtechnologies.

[11]  Theodoros D. Tsiboukis,et al.  Reduction of numerical dispersion in FDTD method through artificial anisotropy , 2000 .

[12]  Javier Alda,et al.  Bad pixel identification by means of principal components analysis , 2002 .

[13]  J. Alda,et al.  Numerical artifacts in finite-difference time- domain algorithms analyzed by means of principal components , 2005, IEEE Transactions on Antennas and Propagation.

[14]  A. Taflove,et al.  Radar Cross Section of General Three-Dimensional Scatterers , 1983, IEEE Transactions on Electromagnetic Compatibility.

[15]  Javier Alda,et al.  Characterization of photonic crystal microcavities with manufacture imperfections. , 2005, Optics express.

[16]  Orhan Arikan,et al.  An efficient and accurate technique for the incident-wave excitations in the FDTD method , 1998 .