基于改进型模拟退火模型的X荧光光谱寻峰方法的.pdfVIP

第33卷，第9期 光谱学与光谱分析 201 and September，2013 3年9月 SpectroscopySpectralAnalysis 基于改进型模拟退火模型的X荧光光谱寻峰方法的研究 曾国强，葛良全，罗耀耀，张庆贤，程 峰，王广西 成都理工大学，四川成都610059 摘要在能量色散X荧光分析技术分析中，对谱信息的处理一直是研究的重点，谱光滑、寻峰、峰面积处 理都是重中之重，本工作基于模拟退火算法原理，建立一种新的寻峰模型算法，该算法利用模拟退火寻找全 局最优点的收敛特性，以Metropolis准则作为峰谷判断的基础，引入新的判定标准和峰谷数组，同时从谱道 址两端同时收敛，以收敛到同一最优解为终止条件。同时，利用该算法与简单比较法、三阶导数法进行对比 实验，结果证明该算法在x荧光光谱寻峰中有较强能力，在实际生产中具有一定价值。 关键词能量色散X荧光分析；寻峰模型；改进型模拟退火 中图分类号：0657．3文献标识码：AⅨ)I：10．3964／j．issrL1000-0593(2013}09—2583—03 于是本文对模拟退火算法进行了如下改进。 引 言 (1)初始值设置，初始峰值T，初始状态S，每个峰值T 的迭代次数L，预存判断数组D1作为波谷数组，功作为波 在光谱寻峰方法中，常用的方法为简单比较法、一阶导 峰数组。(2)触发一个震荡，产生新的状态S，。(3)计算△T 数法、二阶导数法、三阶导数法、对称零面积变换法等[1‘3]。 =C(S，)一c(S)，C作为评价函数，表示是否达到理想状态。 其中简单比较法只能用于寻找简单的峰，不适合寻找重峰。 一阶导数法同样不适用于重峰的分辨，二阶导数法的灵敏度 接受S，作为新状态，如果△T0，则以概率exp(～△T／D 不够，三阶导数法容易识别出假峰，对称零面积变换法的变 接受S，作为新状态，同时对比上一次的△T状态，如果第一 换宽度不易确认[1]。 次AT0，第二次△T0，那么在D1数组添加一组新数组， 模拟退火算法(simulatedannealing，SA)是一种基于蒙记录当前峰值丁，并利用道址转换公式计算当前道址。反之 特卡洛(Monte-Carlo)迭代求解的随机寻优算法，在工程生产在D2数组添加一组新数组。(5)当模型满足终止条件的时 中已经广泛应用[45]。模拟退火在设计初始模型和初始值的候，结束算法，否则重复2～4。(6)由于模型本身是寻找全局 选择空间都比较大，且不需要再计算雅克比矩阵，在收敛过 最优点，因此，在整个X荧光光谱中，需要同时从首道址和 程中，会跳入并跳出局部最小点，直到寻找到全局最优 尾道址一起向中部收敛，因此称之为双退火模型，并且其最 解is,7]。本文利用模拟退火算法的思想，对模型和迭代过程优解要为同一峰值才终止整个迭代。(7)整理数组D1和功 进行改进，使之适用于X荧光光谱寻峰。 并对应到光谱峰值。 1模拟退火算法及改进 2实验部分 利用模拟退火算法进行寻峰，最为重要的是利用其在收 采用的数据自新疆西天山化探样品，以道址一计数形式 敛过程中，会不断跳出局部最小点这一特性[6’8’9]，添加新应用，单一样品仪器谱见图1，横坐标为道址，纵坐标为计 的触发条件，保存峰值数组。在经典的模拟退火算法模型 数。通过刻度6．403k∈V为209道址，17．218keV为840道 中，初始状态的设置对模型的影响不大[7’9]，然而普通的模址，再利用图2将其他峰进行刻度，在寻峰过程中，就能自 拟退火寻峰模型，在寻峰过程中，无论从任何一处进行退火 动识别其他峰并标记出来。 模拟，都只能进行单向迭代收敛，