Discretization of the Lippmann-Schwinger integral equation with complex conjugate gradient method and fast Fourier transform (CCGM-FFT) is solved, which can reduce the memory storage and the CPU time compared with the traditional method, MOM. Thus objects with large size and multiple scattering objects could be simulated with CCGM-FFT. The total intensity and the distribution of each field component of the dielectric and metallic objects under the excitation of the TE/TM-polarized wave are calculated with photon scanning tunnelling microscopy (PSTM) at the constant height. The simulating results are analysed and explained reasonably. The results show that the polarization plays an important role for imaging of PSTM.