摘要：根据Debye模型、Drude模型和Lorentz模型3种常见色散介质模型频域极化率的特点，利用频域到时域的转换关系jω→Э/Эt，将极化矢量P与电场强度E的频域关系转换成时域内关于P的二阶微分方程，其对3种色散介质模型皆适用，具有统一的形式．然后采用相比于中心差分具有更高精度的Newmark两步算法（Newmark-β-γ法）求解该方程，进而得到E→P的递推公式，再结合本构关系得到D→E的时域递推式．实现了色散介质电磁场量的时域有限差分迭代计算．数值计算结果表明该方法是适用于3种色散介质模型的通用算法，并且相比于移位算子时域有限差分方法等以中心差分为基础的离散方案具有更高的计算精度．

关键词：newmark方法 色散介质 电磁散射 时域有限差分方法 

单位：西安电子科技大学物理与光电工程学院 西安710071