简单自相关代数迭代重建算法.pdf

第26卷，第12期 光谱学与光谱分析 oo 2 6年l2月 and December，2006 spectroscopySpectralAnalySis 简单自相关代数迭代重建算法 宋一中1’2，胡国英3，贺安之1 1．南京理工大学理学院信息物理与工程系，江苏南京210094 2．德州学院物理系，山东德州 253000 3．烟台大学光电信息科学与技术学院，山东烟台 264000 self_correlativereconstruction 摘 要 创建了简单自相关代数迭代重建算法(simple algebraic technique，SS— reconstruction cART)，分析了其重建效果。采用数值模拟的方法，对代数迭代法(algebrajctechnique， 建精度均方误差(meansquare sSCART的MsE)／ART的MsE]，峰值相对误差(peakerror，PE)指标在101数量级上降低了99．9％。计 算上仅仅比ART多1个除数因子，所以sSCART具有算法简单、易于实现、重建精度高、稳定性好等优点。 就作者检索的文献看是目前最优秀的代数迭代法。 主题词数值模拟；重建精度；误差分析；代数迭代法；光学层析 中图分类号：()29文献标识码：A 文章编号：1000—0593(2006)12—2364一04 本文用一个非对称、含有正负相的三Gauss峰函数来模 引 言 拟三维场的1个切片， 利用发射光谱诊断等离子体时口’2]，仅能给出估计范 围l““。应用光学计算机层析技术(optical (1) computerized (y+O．4)2)) tomography，()cT)对等离子体诊断是本课题组正在探讨的函数图像如图1所示。 问题。目前，关于()CT的研究主要还集中在重建算法上。已 经提出了许多算法，这些算法大致可以分为2类：变换类算 法(卷积)和级数类算法(迭代)。级数类算法，典型代表是代 reconstruction 数迭代法f5](aIgebraic technique，ART)适于非 完全数据重建，所以在()CT领域得到广泛研究，并提出了 许多种改进算法，如：ART3[“，SARF“，AVART[…， MAR下9]，QMAR可”1等。这些算法，就收敛速度、稳定性、 精度等某一方面有所改进，但计算复杂了许多，有些算法还 不适合双极场的重建[5’6]。 本文针对复杂流场的显示需要出发，研究适合于双极复 杂场的重建方法，在这些算法的基础上提出了简单自相关代 数迭代重建算法。 Fi昏1M使lelfield 1数值模拟 1．2投影 1．1三维场模拟 重建区域划分成50×50的网格，同方向取100条光线 收稿日期：2005一05—12。修订日期：2005一08—16 作者简介：宋一中，1971年生，南京理工大学理学院光学工程博士后 万方数据 第12期 光谱学与光谱分析 2365 均匀穿过重建区，光线间距设置为压d／2，其