基于小波方法的平滑算法在引力透镜中的应用

摘要

基于离散小波变换(DWT)方法,提出了一种可用于计算三维数值模拟样本面密度的平滑算法.为检验方法的有效性,利用该算法研究了两组不同质量解析度的引力透镜数值模拟样本,样本采用了暗物质晕的等温椭球模型,使用蒙特卡罗方法生成.计算结果表明此算法能够在很高的精度上构建引力透镜模拟样本的面密度分布轮廓,由面密度计算出来的透镜的临界曲线和焦散曲线也能较好地和理论曲线吻合,结果是令人满意的.同时比较了三组不同的小波基的计算结果,包括Daub4,Daub6和B-spline 3th,给出了最优的选择.在不损失平滑效果的同时,此算法具有非常高的速度,非常适合于处理以后更高精度的N体数值模拟.