基于四阶平均向量场方法和拟谱方法构造了Dirac方程的高阶整体保能量格式,利用构造的高阶整体保能量格式数值模拟方程孤立波的演化行为.数值模拟结果表明构造的高阶整体保能量格式可以很好地模拟Dirac方程孤立波的演化行为,并且可以精确地保持方程的整体能量守恒特性.