光的衍射是光的波动性的标志之一。对光衍射现象的分析计算涉及到很多光学领域，如：光学系统设计、光学测量、光学全息、光学加工等，所以对衍射的数值计算非常重要。在衍射的数值计算方法上有很多种，包括：快速傅里叶变换法、分数傅立叶变换法、蒙特卡罗方法、复化辛普森法等，这些方法都有各自的特点，适用范围也各异。本文主要用辛普森矩阵法对衍射公式的计算进行了研究。本文分析了复化辛普森法的公式，并且根据复化辛普森计算定积分的方法推导了计算二重定积分的公式，也就是辛普森矩阵法。针对标量衍射理论的基尔霍夫衍射公式及其近似公式和两类瑞利—索末菲衍射公式进行了计算，所用方法为辛普森矩阵法；在计算过程中可知由于衍射公式的表述不同，从而使计算的时间有很大的不同，当衍射公式中的被积函数复杂时计算时间就长；对衍射公式进行相应的简化可以使计算时间大大减少。用复化辛普森法在计算菲涅尔衍射公式上与快速傅立叶变换法进行了比较；计算结果表明快速傅里叶变换法的计算速度高于复化辛普森法，但是计算结果没有复化辛普森法准确。当中分析了快速傅里叶变换法的计算速度比复化辛普森法高的原因，是因为计算中的复指数乘法的计算次数直接影响着计算的速度，而在相同条件下，复化辛普森法在计算菲涅尔衍射中的计算次数比快速傅里叶变换法高很多，所以直接导致计算时间上的差别。