This paper utilizes a kinetic theory model to gain the accurate electromagnetic (EM) wave characteristics in the plasma region, based on the solution of Maxwell-Boltzmann (MB) equation system. The system is solved by finite-difference time-domain (FDTD) algorithm, which gives the results of the electric field intensity and particle velocity distribution function. Furthermore, the validity and effectiveness of the proposed method is verified by comparing the results of reflection coefficient and transmission coefficient of the wave that impinges on plasma plate region with that of the analytic solution.