推导了三维抛物型方程(3DPE)方法的相关计算公式,在此基础上,利用CNFD(Crank-Nicolson Finite-Difference)技术和传统非局部吸收边界条件,采用了矩阵预处理技术处理大规模矩阵,大大改善了计算机存储量的限制状况,同时实现了场值的快速递推。将数值计算结果与双线干涉理论值进行比较,两者吻合很好。