用模拟退火法进行纯永磁波荡器磁块组合优化.pdfVIP

维普资讯 第 17卷 第 1期 强 激 光 与 粒 子 束 Vo1．17，NO．1 2005年 1月 HIGH POWER LASER AND PARTICLE BEAMS Jan．，2005 文章编号 ： 1001—4322(2005)01—0157—04 用模拟退火法进行纯永磁波荡器磁块组合优化 陈 念，何多慧，李 格，贾启卡，张鹏飞，徐宏亮， 蔡根旺 (中国科学技术大学 国家同步辐射实验室，安徽 合肥 230029) 摘 要： 纯永磁波荡器由多个磁块组成，磁块的剩磁离散性会引起波荡器磁场误差，从而影响储存环工 作状态和 自发辐射谱质量。在波荡器磁块安装之前，使用模拟退火法对磁块进行组合排序优化，可以使峰值场 强误差降低到 1O 量级以下，磁场一次积分降低到 1O T ·m量级，二次积分降低到 1O T ·m。量级，优化 结果不依赖于初始状态的选择。给出优化的详细过程 ，提出了根据磁块剩磁快速计算波荡器峰值场强误差和 积分场的方法。 关键词： 波荡器； 模拟退火法； 磁场误差； 积分场 中图分类号 ： TN544 文献标识码 ： A 波荡器(undulator)是第三代同步辐射光源及 自由电子激光装置中的关键设备，由周期性磁场组成，能够 产生 1em到远紫外甚至软 X射线波段的 自发辐射光，其亮度 比弯铁 的同步辐射光高两个数量级 以上。波荡 器磁场性能影响着储存环工作状态和 自发辐射谱质量，一般采用积分场和多极场分布来衡量波荡器磁场误差 对储存环工作状态的影响，采用峰值场强误差和位相误差来标度 自发辐射谱质量口]。 波荡器安装之前，必须对磁块进行组合优化，文献[2～3]使用遗传算法以峰值场强误差作为 目标函数，文 献-14]使用模拟退火法以位相误差作为 目标函数，优化结果保证了自发辐射谱质量，但没有考虑到对储存环工 作状态的影响。中国原子能科学研究院在研制 国家同步辐射实验室 (NSRL)光学速调管时，使用模拟退火法 以二次积分曲线峰值误差作为 目标函数进行磁块组合优化，计算量过大，对积分场的优化结果不够理想。因 此，在NSRL光学速调管改造中，我们使用 了多 目标函数 (峰值场强误差、磁场一次积分绝对值和磁场二次积 分绝对值)利用模拟退火法进行磁块组合优化，为了提高计算速度，提出了快速计算峰值场强误差和积分场的 方法。最终可以在 30min之 内给出一组近似最优的磁块组合方式，使波荡器峰值场强误差、磁场一次积分和 二次积分都尽可能小。 1 模拟退火法基本理论 统计热力学表明，在平衡状态下，原子能量的概率分布与温度 T满足Boltzmann方程 Pi一 1 p(昔) ㈩ 其 中：E 为状态i的能量 ；P为原子在状态 i的几率；k为Boltzmann常数；Z为概率分布的标准化因子 。在对 固体进行退火处理时，常先将它加温使其粒子可 自由运动。显然，固体处于能量较低微观态的几率较大。在温 度降低时，能量相对最低的微观态最有可能出现。当温度趋于零时，固体只能处于能量最小的基态上，形成低 能态晶格 。 模拟退火法(simulatedannealingalgorithm)源于对固体退火过程的模拟；采用Metropolis接受准则；并用 冷却进度表控制算法进程，使算法在多项式时间里给出一个近似最优解 引。Metropolis接受准则：首先以组合 问题的初始状态 i作为当前状态，目标函数为厂()。然后对组合进行随机微小变化，得到新状态 ，目标函数 为厂()。如果 f(j)≤厂()，则以 取代 成为当前状态；如果厂()厂()，则在Eo，1)区间产生一个随机数 ， 若 e[，“卜 ，则以J取代 成为当前状态(￡为控制参数)。相应的算法称为Metropolis算法。 模拟退火法具体步骤 ：从某个初始解 出发，经过 L 次变换并执行 Metropoils算法后，可以得到控制参数 t—to时组合优化问题的相对最优解。然后减小t的值 ，重复执行Metropoils算法，经过k次迭代后控制参数t