本发明提供一种动态星图拖尾星像的数值模拟方法。实现步骤如下：(1)确定龙哥库塔积分步长，应用龙格库塔方法解算基于四元数的星敏感器姿态四元数运动学方程，得到曝光时间内的星敏感器姿态四元数序列；(2)基于星敏感器姿态四元数序列求解视场内所有模拟导航星在星敏感器本体系下的星光矢量序列；(3)依据光学系统模型，将星光矢量序列映射到星敏感器成像阵列面上，得到视场中每颗星的映射位置点序列；(4)以各位置点序列为中心建立二维高斯灰度点扩散函数，所有高斯点扩散函数作叠加，即可模拟出任一星的一条动态拖尾星像。本发明方法符合星光拖尾成像的物理过程，解决了星模拟器动态星图模拟的关键问题，用以实现动态星图模拟。