模拟回火马尔可夫链蒙特卡罗全波形分析方法-物理学报.PDF

物理学报 Acta Phys. Sin. Vol. 63, No. 16 (2014) 164205 模拟回火马尔可夫链蒙特卡罗全波形分析方法 尹文也 何伟基 顾国华 陈钱 (南京理工大学光电技术系, 南京 210094) ( 2014 年1 月28 日收到; 2014 年3 月26 日收到修改稿) 针对传统的全波形分析方法不能快速自动处理全波形数据的缺点, 提出了一种模拟回火马尔可夫链蒙特 卡罗全波形分析法, 用于求解全波形数据中的波峰数和峰值位置等参量. 该方法采用Metropolis 更新策略求 解波峰数量和噪声两个参量, 以达到快速求解的目的; 而峰值位置和波峰幅值则采用改进的模拟回火策略求 解, 通过添加的主动干预回火步骤实现对参量更新过程的有效探测, 以满足对速度或运算收敛性的要求. 模 拟回火马尔可夫链蒙特卡罗全波形分析方法以马尔可夫算法为基础, 仍保持马氏链的收敛性, 从而保证本方 法具有良好的鲁棒性, 实现对全波形数据的自动化处理. 关键词: 全波形分析, 模拟回火马尔可夫链蒙特卡罗方法, 主动干预回火, Metropolis 策略 PACS: 42.30.Sy, 42.62.Eh, 42.68.Wt, 42.79.Qx DOI: 10.7498/aps.63.164205 不同的运算结果. 为了提高算法稳定性, Biernacki 1 引 言 等6 提出了期望值多次最大化方法, 统计各种初始 化条件下的运算结果, 这虽然能够得到一个相对稳 1 激光雷达 自问世以来一直备受关注, 但是传 定的运算结果, 但是牺牲了运算时间. Pernkopf 和 统的激光雷达通常只能记录激光光斑内的单个距 Bouchaffra7 提出了基于遗传算法的期望值最大化 离特征, 无法记录并处理光斑内多个距离特征目标 方法解决了有关算法稳定性的问题, 此方法具有 的回波信息. 采用具有完整记录接收信号功能的全 良好的鲁棒性, 但是仍需要预先了解光斑内的距 波形激光雷达可保留光斑内所有回波信号以待后 离特征数目. 第三种是基于马尔可夫链蒙特卡罗 续处理, 这可应对特性未知目标2 的回波信号, 在 8 (MCMC) 算法 的随机方法, 全波形参量的求解 军事和民用领域具有广阔的应用前景. 为了提取 过程就是马氏链的更新过程. 理论上, 进行足够多 光斑内的所有距离信息, 已发展了多种全波形分析 9 的更新就一定能够达到所有马氏链的平稳 状态, 方法. 完成参量求解. 2007 年, Hernández-Marín 等 10 目前, 通过建模方式处理全波形数据的方 在MCMC 算法的基础上, 提出了贝叶斯框架下可 法主要有以下三种: 第一种是基于Levenberg- 逆跳跃MCMC (RJMCMC) 算法, 实现了对特性未 3 Marquardt 算法 的非线性最小二乘法波形拟合 知的目标信号进行自动化处理. 通过延迟拒绝转 4 方法. Clauve 等 根据不同特性目标 (草地、建筑 移的方式