多道γ能谱分析软件中寻峰算法比较总结
差值算法在γ能谱峰面积计算中的应用
天然放射性元素 Ｋ 、ｕ 、Ｔ ｈ 的 能谱可用各种 能谱仪测得 ，仪器标定后 ，即可将 能谱特 征 峰面积 转换 成 Ｋ、ｕ、Ｔ ｈ含量 。因此 ， 能 谱特 征峰 峰面 积 的计 算 精度 非 常重 要 ， 目前 峰 面积 计算法一般有计数法和函数拟合法等…。计数法主要包括总峰面积法 、科沃尔峰面积法和瓦森 峰 面积法 ；函数拟 合法 包括 扣 除本底 后净 峰拟 合法 以及 峰形 函数 加 基底 函数 拟合 的峰 面 积法 Ｊ 。 在 峰面 积计 算 中 ，左 右边界 往 往是 以道 址 确 定 的 ，由于 道址 的整 数 特 性 ，往 往 与 能谱 的能 量 边界不能完全一致 ，因此给峰面积的计算带来计算误差。作者利用插值算法对此进行 了插值计
ｅ ｒ ｒ ｏ ｒ ．Ｅ ｘ ｐ ｅ ｒ ｉ ｍ ｅ ｎ ｔ ｓ ｓ ｈ ｏ ｗ ｔ ｈ ａ ｔ ：ｔ ｈ ｅ １ ０ ２ ４ Ｎ ａ Ｉ（ Ｔ １ ）ｇ ａ ｍｍ ａ ｓ ｐ ｅ ｃ ｔ ｒ ｏ ｍｅ ｔ ｅ ｒ ，ｃ ｏ ｕ ｎ ｔ ｒ ａ ｔ ｅ ｍａ ｘ ｉ ｍ ｕ ｍ ｃ ｏ ｒ ｅ ｃ ｔ ｉ ｏ ｎ
高 红 赖 万 昌 李思威 曹发 明
（ 成都 理 工大 学核 技 术与 自动 化 工程 学院 ，四川 成都
摘
６ １ ０ ０ ５ ９ ）
要 ：野外 多道 能谱仪在计算 能谱 峰面积时 ，左右边界道址 取整数往往 产生较大 的计 算误
差 ，引入拉格 朗 日插值 法对 非整数道的道址进行插值 ，能有效减低 这种误差 。实验表 明，对 １ ０ ２ ４ 道的 Ｎ ａ Ｉ （ 哦） 能谱仪 ，计数率最大校正值可达 １ ． ０ ８ ｃ ｐ ｓ ，测量误差则可降低 １ ０ ． ６ １ ％。
Ｇａ ｍｍ ａ Ｓ ｐ ｅ ｃ ｔ ｒ ａ ｌ Ｐｅ ａ ｋ Ａｒ ｅ ａ
多道γ能谱分析软件中寻峰算法比较总结
多道γ能谱分析软件中寻峰算法比较总结寻峰算法是检测和识别γ能谱中标准峰的一个过程，是γ能谱的分析过程的关键步骤。
随着γ能谱分析软件的不断发展，出现了多种寻峰算法，优势各有不同，A软件的算法一般用于比较复杂的能谱，而B软件的算法一般用于小范围的能谱分析。
本文将通过对常见的几种γ能谱分析软件中寻峰算法的比较，总结当前寻峰算法的发展趋势和应用前景。
一、Geology 7.0
Geology 7.0是由英国Geochemical Ltd公司开发出的一款专业的γ能谱分析软件。
该软件采用了积分式寻峰算法，即根据信号在每一小段内累积平均值，将信号的变化和峰的确定结合起来。
在寻峰时，Geology 7.0会自动调整峰陡面参数，有效降低对峰形幅度变化的变异性，实现自适应寻峰，从而提高能谱的分析准确度，较好的保持能谱的原始数据。
Geology 7.0的寻峰算法还可以避免背景干扰的出现，以便精确地检测出底峰面，能够有效地抑制低能底峰的异常检测。
此外，它还可以识别出隐藏在较宽能谱峰中的比较窄的峰，让检测结果更加精确。
二、Merck Millipore Type Spectrum
Merck Millipore Type Spectrum是中国德州MILLIPORE公司推出的一款多功能γ能谱分析软件。
软件采用基于峰头的寻峰算法，能够有效识别四种不同形状的γ能谱峰：单峰、双峰、三角峰和脊峰。
,γ能谱测量方法
γ能谱测量是一种用于分析和测量γ射线的能量和强度分布的方法。
以下是几种常见的γ能谱测量方法：
1.闪烁体探测器（Scintillation Detector）：该方法使用具有闪烁性质的物质作为探测器，
当γ射线通过闪烁体时，会产生光闪烁。
这些闪烁信号被转换为电信号，并通过放大和处理后，形成γ能谱。
2.半导体探测器（Semiconductor Detector）：半导体探测器利用半导体材料的特性来测量
γ射线的能量。
当γ射线与半导体相互作用时，会产生电子-空穴对。
通过测量这些电荷对的移动和收集，可以得到γ能谱。
3.多道分析器（Multichannel Analyzer）：多道分析器是一种将不同能量范围内的γ射线
分离并计数的设备。
它通常与闪烁体或半导体探测器一起使用。
多道分析器将接收到的信号根据能量进行离散化，并将其对应到不同的道址上，形成γ能谱。
4.探测器阵列（Detector Array）：此方法使用多个探测器组成的阵列来测量γ射线。
每个
探测器都可以提供关于能量和位置的信息，通过组合分析得到完整的γ能谱。
这些方法在γ能谱测量中具有不同的特点和应用范围，可以根据实验需求选择合适的测量方法。
无论采用哪种方法，γ能谱测量都是研究核物理、放射性衰变以及其他与γ射线相关领域的重要手段。
能谱寻峰比较法
能谱寻峰比较法
能谱寻峰比较法是一种常用的寻峰方法，其基本步骤如下：
1.采集能谱数据：通过实验或检测设备获取能谱数据，得到一组随能量变化的数据点。
2.预处理数据：对采集到的能谱数据进行预处理，包括去除噪声、平滑处理等，以提高寻峰的准确性。
3.确定寻峰范围：根据实际需要，确定寻峰的范围。
通常，寻峰范围是根据实验或检测的目的确定的，例如可能需要寻找某个特定能量范围内的峰。
4.峰检测和识别：在确定的寻峰范围内，使用适当的寻峰算法进行峰检测和识别。
常用的寻峰算法包括对称零面积对合法、比较法寻峰法、导数法寻峰法等。
5.参数调整和优化：根据实际情况，调整和优化寻峰算法的参数。
例如，可以通过调整阈值、平滑参数等来提高寻峰的准确性和稳定性。
6.峰位确定：根据寻峰算法的结果，确定每个峰的峰位。
通常，峰位的确定是根据实验或检测的要求进行的，例如可能需要计算每个峰的精确能量位置。
7.数据分析：对确定峰位的能谱数据进行进一步分析，例如可以计算峰的半高宽、精确峰位、峰宽等信息，以便进行后续的数据处理和分析。
总之，能谱寻峰比较法是一种基于数据驱动的寻峰方法，其关键在于选择合适的寻峰算法和参数调整，以确保寻峰的准确性和稳定性。
在实际应用中，需要根据具体实验或检测条件选择适合的寻峰方法，并对数据进行深入分析和处理，以获得更准确和可靠的结果。
γ能谱滤波方法比对研究
γ能谱滤波方法比对研究
γ能谱滤波方法比对研究是一项重要的研究，旨在分析两种不同的γ能谱滤波方法，以便从中找出更有效和精确的方法。
γ能谱滤波法是用于处理γ射线信号的一种常用方法，它通过对γ射线信号进行滤波来增强信号的信噪比。
相对于传统的信号处理方法，它具有更高的准确性和数据处理能力。
γ能谱滤波方法的比较研究包括对不同滤波方法的相容性、可靠性和可行性进行评估，以及对其他影响因素的考虑。
首先，研究人员将评估不同滤波方法的性能，考虑信号的精度、稳定性和精确性。
其次，将对不同滤波方法的可靠性进行比较，评估其对不同环境条件下的性能。
此外，还将考虑其他因素，例如信号处理的算法、冗余度、信号处理系统的复杂性等。
在研究不同滤波方法之前，研究者需要对γ射线的特性进行分析，找出γ射线的特征参数，以及滤波器的特性参数，并将其作为研究滤波方法的基础。
从而为比较研究奠定基础。
接下来，将采用模拟计算和实验测量等方法，对不同滤波方法进行比较研究。
根据模拟计算和实验测量的结
果，研究者可以得出关于不同滤波方法之间的差异，以及各自的优势和劣势。
最后，研究者将对比较研究结果进行概括和总结，指出不同滤波方法的优势及其在不同应用场景下的适用性。
此外，还可以结合比较研究的结果，分析并建议新的滤波方法或者改进已有的滤波方法，以提高滤波的精度和准确性。
总的来说，γ能谱滤波方法比较研究是一项重要的研究，可以比较不同滤波方法的性能，为选择更有效和精确的滤波方法提供参考。
有效的滤波方法有助于提高γ射线信号处理的准确性，也有助于提高γ射线信号处理系统的性能。
多道γ能谱分析软件中寻峰算法比较总结
自动寻峰由于谱结构的复杂和统计涨落的影响，从谱中正确地找到全部存在的峰是比较困难的。
尤其是找到位于很高本底上的弱峰，分辨出相互靠得很近的重峰更为困难。
谱分析对寻峰方法的基本要求如下：(1) 比较高的重峰分辨能力。
能确定相互距离很近的峰的峰位。
(2) 能识别弱峰，特别是位于高本底上的弱峰。
(3) 假峰出现的几率要小。
(4) 不仅能计算出峰位的整数道址，还能计算出峰位的精确值，某些情况下要求峰位的误差小于0.2道。
很多作者对寻峰方法进行了研究，提出了很多有效的寻峰方法。
目的：判断有没有峰存在确定峰位（高斯分布的数学期望），以便把峰位对应的道址，转换成能量确定峰边界——为计算峰面积服务（峰边界道的确定，直接影响峰面积的计算）分为两个步骤：谱变换和峰判定要求：支持手动/自动寻峰，参数输入，同时计算并显示峰半高宽、精确峰位、峰宽等信息，能够区分康普顿边沿和假峰感兴区内寻峰人工设置感兴趣大小，然后在感兴区内采用简单方法寻峰重点研究：对感兴区内的弱峰寻峰、重峰的分解对于一个单峰区，当峰形在峰位两侧比较对称时，可以由峰的FWHM计算峰区的左、右边界道址。
峰区的宽度取为3FWHM，FWHM的值可以根据峰位m p由测量系统的FWHM 刻度公式计算。
由于峰形对称，左、右边界道和峰位的距离都是1.5FWHNM 。
)5.0FWHM 5.1(INT p L +-=m m )5.0FWHM 5.1(INT p R ++=m m式中m p 是峰位，INT 的含义是取整数。
对于存在有低能尾部的峰，其峰形函数描述（参见图）。
]2/)([H 22p m σ--=m m EXP y ，m ≥mp －J]2/)22([HEXP 2p m σ+-=J m m J y ，m ≤mp －J式中H 为峰高，mp 为峰位，σ是高斯函数的标准偏差，J 为接点的道址和峰位之间的距离。
在峰位的左侧，有一个接点，其道址为mp －J 。
在接点的右侧，峰函数是高斯函数。
最大值法寻峰
最大值法寻峰最大值法是信号处理领域中常用的一种寻峰方法，该方法适用于各种类型的数据集，包括波形、光谱、图像等。
在实际应用中，最大值法寻峰技术已被广泛使用于化学分析、生命科学、物理学、医学等领域。
一、最大值法寻峰原理最大值法寻峰原理基于峰的定义，即峰是局部极大值。
该方法首先寻找数据集中的最大值，并确定该最大值处曲线的方向。
然后，根据局部曲线的方向，对峰进行识别和分类。
对于单峰数据集，最大值法寻峰方法是相对简单的。
在这种情况下，找到最大值并在两个方向上延伸，直到曲线开始下降为止。
此时，峰的边界被定义为曲线开始下降的位置。
对于多峰数据集，最大值法寻峰方法还需要对每个峰进行分类和分离。
这可以通过确定峰的起点和终点来实现。
二、最大值法寻峰优点1. 算法简单直观：使用最大值法寻峰技术可以快速有效地找到数据集中的峰。
2. 对多峰数据集准确性高：最大值法寻峰技术可以精确识别多峰数据集中的各个峰。
3. 可自动识别峰的个数和峰的位置：使用最大值法寻峰技术可以自动识别数据集中峰的数量和位置，无需手动干预。
三、最大值法寻峰缺点1. 对背景噪音敏感：在数据集中存在较大的背景噪音时，最大值法寻峰技术容易出现误判。
2. 不能确定峰的宽度：最大值法寻峰技术不能精确确定峰的宽度，从而限制了其在某些应用领域中的应用。
四、最大值法寻峰实际应用1. 化学分析：在化学分析中，最大值法寻峰技术被广泛用于分析复杂的光谱数据。
2. 生命科学：在生命科学中，最大值法寻峰技术可以用于识别某些生物分子的光谱峰值。
3. 物理学和医学：在物理学和医学领域中，最大值法寻峰技术可以用于分析X射线荧光光谱和MRI等影像数据。
五、总结最大值法寻峰技术是一种简单、直观、高效的寻峰方法，适用于各种数据集类型。
然而，在实际应用中还需要考虑背景噪声和峰的宽度等因素。
尽管如此，最大值法寻峰仍被广泛应用于信号处理领域，包括化学分析、生命科学、物理学和医学等领域。
γ射线能谱分析试验报告
γ射线能谱分析试验一、预习报告实验名称:γ射线能谱分析试验.实验内容:1.学会Ｎa(TI)单晶Y闪烁谱仪整套装置的操作.调整和使用.2.了解多道脉冲分析器在Ｎa(TI)单晶Y闪烁谱仪测量中的数据采集及其基本功能.3.测量Ｎa(TI)单晶Y闪烁谱仪的能量和线性.4.分析137C S单能Y射线谱仪.5. 测量Na(TI)单晶Y闪烁探测器的计数率随工作电压变化的关系曲线.实验目的:1.了解闪烁探测器的结构.原理.2.掌握Ｎa(TI)单晶Y闪烁谱仪的几个性能指标和测量方法.3.了解核电子学仪器的数据采集.记录方法和数据处理原理.4.测量Na(TI)单晶Y闪烁探测器的坪曲线,确定合适的工作电压;实验仪器:Ｎa(TI)单晶Y闪烁探头;微机多道Y(X)谱仪装置;Y放射源137C S一个.实验原理：1. 射线与物质的相互作用γ射线是原子核从激发态跃迁到低能态或基态时发射的波长很短的电磁辐射，研究γ射线的能谱对原子核的能级结构和放射性核素的应用等方面具有重要的意义。
γ射线与物质相互作用，可以有许多方式。
当γ射线的能量在30MeV 以下时，在所有相互作用方式中，最主要的三种，如图 1-1所示。
图1-1 γ射线与物质相互作用示意图(1)光电效应：入射γ粒子把能量全部转移给原子中的束缚电子，光子本身消失而把束缚电子打出来形成光电子这个过程称为光电效应。
由于束缚电子的电离能E i 一般远小于入射γ射线的能量E γ，所以光电子的动能近似等于入射γ射线的能量，即：E 光电 = E γ - E i ≈E γ （1）(2)康普顿散射：核外自由电子与入射γ射线发生康普顿散射。
根据动量守恒的要求，散射与入射只能发生在一个平面内。
设入射γ光子能量为hv ，散射光子能量为hv′，康普顿散射后散射光子能量与散射角θ的关系为：)cos 1(1θ-+='a hv v h（2）式中2c m hva e =，即为入射γ射线能量与电子静止质量m e所对应的能量之比。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
自动寻峰由于谱结构的复杂和统计涨落的影响，从谱中正确地找到全部存在的峰是比较困难的。
尤其是找到位于很高本底上的弱峰，分辨出相互靠得很近的重峰更为困难。
谱分析对寻峰方法的基本要求如下：(1) 比较高的重峰分辨能力。
能确定相互距离很近的峰的峰位。
(2) 能识别弱峰，特别是位于高本底上的弱峰。
(3) 假峰出现的几率要小。
(4) 不仅能计算出峰位的整数道址，还能计算出峰位的精确值，某些情况下要求峰位的误差小于0.2道。
很多作者对寻峰方法进行了研究，提出了很多有效的寻峰方法。
目的：判断有没有峰存在确定峰位（高斯分布的数学期望），以便把峰位对应的道址，转换成能量确定峰边界——为计算峰面积服务（峰边界道的确定，直接影响峰面积的计算）分为两个步骤：谱变换和峰判定要求：支持手动/自动寻峰，参数输入，同时计算并显示峰半高宽、精确峰位、峰宽等信息，能够区分康普顿边沿和假峰感兴区内寻峰人工设置感兴趣大小，然后在感兴区内采用简单方法寻峰重点研究：对感兴区内的弱峰寻峰、重峰的分解对于一个单峰区，当峰形在峰位两侧比较对称时，可以由峰的FWHM计算峰区的左、右边界道址。
峰区的宽度取为3FWHM，FWHM的值可以根据峰位m p由测量系统的FWHM刻度公式计算。
由于峰形对称，左、右边界道和峰位的距离都是1.5FWHNM 。
)5.0FWHM 5.1(INT p L +-=m m )5.0FWHM 5.1(INT p R ++=m m式中m p 是峰位，INT 的含义是取整数。
对于存在有低能尾部的峰，其峰形函数描述（参见图）。
]2/)([H 22p m σ--=m m EXP y ，m ≥mp －J]2/)22([HEXP 2p m σ+-=J m m J y ，m ≤mp －J式中H 为峰高，mp 为峰位，σ是高斯函数的标准偏差，J 为接点的道址和峰位之间的距离。
在峰位的左侧，有一个接点，其道址为mp －J 。
在接点的右侧，峰函数是高斯函数。
在接点的左侧，峰函数用指数曲线来描述。
这时峰区的左、右边界道址为)5.05.0/FWHM 12.1(INT 2p L +--=J J m m)5.0FWHM 5.1(INT p R ++=m m带有低能尾部的峰函数的图形全谱自动寻峰基于核素库法：能量刻度完成后，根据核素库中的能量计算对应的道址，在各个道址附近（左右10道附近）采用简单的寻峰方法（导数法） 方法：根据仪器选择开发IF 函数法/简单比较法(适于寻找强单峰，速度快)满足条件：m i i i m i data data k data data +->-< 可认为有峰存在 然后在data i-m 至data i+m 中找最大值，对应的道值即为峰位 k ：找峰阈值，根据高斯分布，一般k 取值1—1.5 常用5点、7点极大值法（m 取2，3） 判定峰是否有意义一般，用R=N0 / Nb ≥ R0确定峰是否有意义 R 为峰谷比， R0为设定值 （经验值） N0为净峰幅度与基底之和 Nb 为基底计数int CMmcaView::SearPeakCompare(int Beginch, int Endch, int m, float k) 高斯乘积函数找峰法（可靠性差，不建议采用） 描述谱峰形状的函数主要是高斯函数[]2202/)(ex p 2)(σσπi i Ai G --=则由相邻的数据点定义一个新的函数（第一高斯乘积函数，只与σ2.3556FWHM =有关）：2)092.11ex p()()2()1()()(2≥=+--+=m Hmm i G i G m i G i G i P m m 是步长（用道表示），是高斯乘积函数的阶数，则Pm(i)称为第m 阶高斯乘积函数。
找峰的灵敏度与m 有关，随m 的增加灵敏度提高。
为避免基线参数的影响，最好扣除本底后，再应用高斯乘积函数找峰。
考虑统计涨落的影响，把判断无峰存在的1变为一个“单位带”。
即峰的判断为：3)/1(/1)(≥⎪⎩⎪⎨⎧+>±≤=k y k y k i P i im 有峰无峰峰位的确定：由Pm(i)过1的两点求平均来确定；峰边界的确定：“单位带”下限的两个最端点；半高宽的确定：函数Pm(i)在“1”上的截距；组合峰的确定：在乘积函数的两个峰之间没有处于“带内”的乘积函数值 导数法（一阶、二阶、三阶）∑-=+=mmj j i jmiy CN y 1'Nm 为规范化常数，Cj 平滑的变换系数。
3次多项式5点光滑一阶导数公式：（可以采用）)88(1212112'++---+-=i i i i i y y y y y 峰位确定：一阶导数值由正变负=0处；峰边界确定：一阶导数由负变正=0处CalculateDifferential(0, size, m, differ); for (int j = m; j <= size-m; j++) {for(int i=1;i<=m;i++) {if(differ[j-i])>0&&differ[j-i]>maxtemp) {maxtemp=differ[j-i]; nmax=j-i;} if(differ[j+i])<0&&differ[j+i]<mintemp) {mintemp=differ[j+i]; nmin=j+i;} }if ((nmin-nmax)>0.8*fwhm && (nmin-nmax)<3*fwhm) //FWHM 参数根据仪器能量分辨率可人工确定，fwhm~20 peakposition[p++]=j+0.5;//保持峰位对应的道址 }5点光滑二阶导数公式（软件中推荐采用）)222(712112''++--+---=i i i i i i y y y y y y//7点二阶导数(5*(countsdata[j-3]+countsdata[j+3])-3*(countsdata[j-1]+countsdata[j+1])-4*c ountsdata[j])/42;)0.220.670.580.580.670.22(0.2521321123'+++----++--=i i i i i i i y y y y y y y软件中推荐采用11点以上的公式峰位确定：二阶导数最小值对应的道址；峰边界确定：二阶导数正极大值点 for (int j = m; j <= size-m; j++)//m~30 {int maxtemp=-0.5,mintemp=-0.5; If(differ[j]<-0.05) for(int i=1;i<=m;i++) {if(differ[j-i]>maxtemp) {maxtemp=differ[j-i]; nmax=j-i;} if(differ[j+i]>mintemp) {mintemp=differ[j+i]; nmin=j+i;} }if ((nmin-nmax)>0.8*fwhm && (nmin-nmax)<3*fwhm) //FWHM 参数根据仪器能量分辨率可人工确定，fwhm~20 peakposition[p++]=j+0.5;//保持峰位对应的道址 }σσ/:Derivative 2 of ce Significan ;2nd 2/12111i k k j j i j kj k j j i j i j ji j i i i i D S n c n c D nc n n n D =⎥⎦⎤⎢⎣⎡===-+-=∑∑∑-=+=-=+-=+-()6-j 222122210 c upto Go :k .peak width assumed the is p where exp 100=-⎪⎪⎭⎫ ⎝⎛-=p j p j p c j试验：系列1为处理后的原始能谱，系列2为5点一阶导数，系列3为5点二阶导数，系列4为对称零面积法寻峰Peak ‘found’when S > Threshold只要选择好合适的寻峰阈值，足以满足准确寻找到全能峰，并剔除假峰（如康普顿边沿，反散射峰）5点光滑三阶导数公式判定各感兴区是单峰还是重峰)22(212112'''++--+-+-=i i i i i y y y y y峰位确定：三阶导数由负变正=0处；峰边界确定：三阶导数由正变负=0处 判定峰是否有意义0.8FWHM ≤ N ≤ 3FWHM 峰高判定条件σ≥'-/||5.0m ax e ym TRH y p m这个公式就是在一阶导数法寻峰程序中实际应用的峰高判定条件。
CalculateDifferential(Beginch, Endch, m, differ);int CMmcaView::SearPeakDifferential(int Beginch, int Endch, int fwhm, int differ[], int m){int n1=0, differ[Endch-Beginch+1], nmax=0, nmin=0, maxtemp, mintemp,temp; maxtemp=differ[0]; mintemp=differ[0];for (int j = 1; j <= Endch-Beginch; j++){temp=differ[j-1];if(_copysign(temp,differ[j])!=differ[j-1] && differ[j]<0) n1=j+Beginch;if(differ[j]<mintemp) {mintemp=differ[j]; nmin=j+Beginch;}if(differ[j]>maxtemp) {maxtemp=differ[j]; nmax=j+Beginch;}}if ((nmin-nmax)>0.8*fwhm && (nmin-nmax)<3*fwhm) return n1; else return (0); }对称零面积法（推荐自动寻峰中采用，可探测弱峰和重峰）面积为零的“窗”函数与实验谱数据进行褶积变换，且要求“窗”函数为对称函数。
对线性基底的褶积变换将为零，只有存在峰的地方不为零。
j j mmj jmmj ji jiC C Cy Cy --=-=+===∑∑0'匹配滤波器法（类峰形函数）∑-=-+--=mm k j k m j C ]2ex p[121]2ex p[2222σσ 峰判定准则f data C dataC yy R m m j j i j mmj ji jiii >⎪⎪⎭⎫ ⎝⎛=∆=∑∑-=+-=+212~~2m+1为变换宽度，6355.2FWHM =σ为峰宽参数，若变换后的y'和其均方根误差的比值超过预先给定的寻峰阈值（f ），则认为找到了一个峰。