To save finite-difference time-domain（FDTD） computing time, several methods are proposed to convert the time domain FDTD output into frequency domain. The Pad6 approximation with Baker＇s algorithm and the program are introduced to simulate photonic crystal structures. For a simple pole system with frequency 160THz and quality factor of 5000, the intensity spectrum obtained by the Padé approximation from a 2^8-item sequence output is more exact than that obtained by fast Fourier transformation from a 2^20-item sequence output. The mode frequencies and quality factors are calculated at different wave vectors for the photonic crystal slab from a much shorter FDTD output than that required by the FFT method, and then the band diagrams are obatined. In addition, mode frequencies and Q-factors are calculated for photonic crystal microcavity.
