多道γ能谱分析软件中寻峰算法比较总结
差值算法在γ能谱峰面积计算中的应用
天然放射性元素 Ｋ 、ｕ 、Ｔ ｈ 的 能谱可用各种 能谱仪测得 ，仪器标定后 ，即可将 能谱特 征 峰面积 转换 成 Ｋ、ｕ、Ｔ ｈ含量 。因此 ， 能 谱特 征峰 峰面 积 的计 算 精度 非 常重 要 ， 目前 峰 面积 计算法一般有计数法和函数拟合法等…。计数法主要包括总峰面积法 、科沃尔峰面积法和瓦森 峰 面积法 ；函数拟 合法 包括 扣 除本底 后净 峰拟 合法 以及 峰形 函数 加 基底 函数 拟合 的峰 面 积法 Ｊ 。 在 峰面 积计 算 中 ，左 右边界 往 往是 以道 址 确 定 的 ，由于 道址 的整 数 特 性 ，往 往 与 能谱 的能 量 边界不能完全一致 ，因此给峰面积的计算带来计算误差。作者利用插值算法对此进行 了插值计
ｅ ｒ ｒ ｏ ｒ ．Ｅ ｘ ｐ ｅ ｒ ｉ ｍ ｅ ｎ ｔ ｓ ｓ ｈ ｏ ｗ ｔ ｈ ａ ｔ ：ｔ ｈ ｅ １ ０ ２ ４ Ｎ ａ Ｉ（ Ｔ １ ）ｇ ａ ｍｍ ａ ｓ ｐ ｅ ｃ ｔ ｒ ｏ ｍｅ ｔ ｅ ｒ ，ｃ ｏ ｕ ｎ ｔ ｒ ａ ｔ ｅ ｍａ ｘ ｉ ｍ ｕ ｍ ｃ ｏ ｒ ｅ ｃ ｔ ｉ ｏ ｎ
高 红 赖 万 昌 李思威 曹发 明
（ 成都 理 工大 学核 技 术与 自动 化 工程 学院 ，四川 成都
摘
６ １ ０ ０ ５ ９ ）
要 ：野外 多道 能谱仪在计算 能谱 峰面积时 ，左右边界道址 取整数往往 产生较大 的计 算误
差 ，引入拉格 朗 日插值 法对 非整数道的道址进行插值 ，能有效减低 这种误差 。实验表 明，对 １ ０ ２ ４ 道的 Ｎ ａ Ｉ （ 哦） 能谱仪 ，计数率最大校正值可达 １ ． ０ ８ ｃ ｐ ｓ ，测量误差则可降低 １ ０ ． ６ １ ％。
Ｇａ ｍｍ ａ Ｓ ｐ ｅ ｃ ｔ ｒ ａ ｌ Ｐｅ ａ ｋ Ａｒ ｅ ａ
能谱仪的使用方法与峰位分析技巧
能谱仪的使用方法与峰位分析技巧能谱仪是一种常用的实验仪器，用于分析样品中的元素成分。
它通过测量样品中放射性粒子的能量与强度，从而得到样品的能谱图。
本文将介绍能谱仪的基本使用方法，并分享一些峰位分析技巧，帮助读者更好地利用能谱仪进行实验研究。
一、能谱仪的基本使用方法1. 稳定仪器：在使用能谱仪之前，首先要确保仪器的稳定性。
检查仪器连接是否牢固，各部分仪器的状态是否正常，确保仪器处于可正常工作的状态。
2. 样品制备：根据实验需求，制备好待测样品。
样品的制备方法因实验对象的不同而有所差异，如放射性样品的处理需要特殊注意。
3. 样品装入：将制备好的样品装入能谱仪的样品槽中。
注意确保样品的位置准确且固定，以避免在测量过程中的位置偏差。
4. 调整参数：根据实验要求，调整能谱仪的参数。
这些参数包括放大倍数、灵敏度、测量时间等，应根据实验需求来确定。
5. 开始测量：按下测量按钮，启动能谱仪开始测量。
在测量过程中，要保持实验环境的稳定，以获得准确可靠的测量结果。
6. 储存数据：实验完成后，将能谱仪测得的数据储存起来。
数据可以保存在电脑上或其他存储设备中，以备后续分析使用。
同时也可以通过打印或导出文件的方式进行结果的备份和共享。
二、峰位分析技巧1. 峰位识别：在能谱图中，不同元素的能量峰位会表现为清晰的峰状。
通过观察能谱图，我们可以识别出不同元素的峰位，从而确定样品中的元素成分。
2. 峰位测量：利用能谱仪提供的测量功能，可以精确地测量出不同峰位的位置与强度。
这些数据可以用于后续的峰位分析和元素定量分析。
3. 峰位校准：为了提高测量的准确性，可以进行峰位校准。
峰位校准通过测量一系列已知元素的能量峰位，然后根据这些数据来校正未知样品的峰位。
峰位校准可以提高测量结果的准确性和可靠性。
4. 谱峰分析：在能谱图中，有时会出现多个重叠的峰位。
为了准确地确定每个峰位的能量和强度，可以采用谱峰分析的方法。
谱峰分析通过对峰位进行分段或近似处理，以获取单个峰位的尽可能精确的测量结果。
伽马能谱实验报告doc
伽马能谱实验报告篇一：闪烁伽马能谱测量实验报告实验题目：《闪烁γ能谱测量》一、实验目的1加深对γ射线和物质相互作用的理解。
2.掌握NaI(Tl)γ谱仪的原理及使用方法。
3.学会测量分析γ能谱。
4.学会测定γ谱仪的刻度曲线.。
二、实验仪器Cs放射源 Co放射源 FH1901型NaI闪烁谱仪 SR-28双踪示波器三、实验原理1. γ射线与物质相互作用γ射线与物质相互作用主要有光电效应、康普顿散射及电子对效应。
1）光电效应：在光电效应中，原子吸收光子的全部能量，其中一部分消耗于光电子脱离原子束缚所需的电离能，另一部分就作为光电子的动能。
所以，释放出来的光电子能量和该束缚电子所处的电子壳层的结合能B?之差。
因此，E光电子=E??Bi?E?（需要原子核参加） 2）康普顿散射：康普顿散射是γ光子与原子外层电子相互作用的结果。
反冲电子的动能为：Ee?E?2(1?cos?)m0c2?E?(1?cos?)即使入射γ光子的能量是单一的，反冲电子的能量却是随散射角连续变化的。
3）电子对效应：电子对效应是γ光子从原子核旁经过时，在原子核的库伦场作用下，γ光子转化为一个正电子和一个负电子的过程。
根据能量守恒定律，只有当入射光子的能量hν大于2m0c2，即hν〉1.02MeV时，才能发生电子对效应。
（与光电效应相似，需要原子核参加）2． NaI（Tl）γ能谱仪介绍 1）闪烁谱仪装置示意图2）闪烁谱仪的工作原理Γ射线次级电子荧光Γ放射源与闪烁体发闪烁体受光阴极吸收生三种作用激辐射光电子电脉冲定标器记录分析器分析各打拿极逐级放大3)能谱分析（以137Cs为例）全能峰是γ光子与闪烁体发生光电效应产生的，直接反映了γ射线的能量；康普顿坪是由康普顿效应贡献的；逸出的γ射线与闪烁体周围的物质发生康普顿散射，反散射光子进入闪烁体发生光电效应形成反散射峰。
4）谱仪的能量分(原文来自：小草范文网:伽马能谱实验报告)辨率和能量刻度曲线闪烁单晶γ谱仪最主要的指标是能量分辨率和线性。
HD-2001使用说明书
一、简介HD-2001型低本底多道γ能谱仪符合国标GB6566-2001《建筑材料放射性核素限量》、国标GB50325-2001《民用建筑工程室内环境污染控制规范》和行标《放射性矿产地质分析测试实验室质量保证规范》的要求。
HD-2001型低本底多道γ能谱仪（以下简称γ能谱仪）是国内目前最新一代微机多道谱仪系统，技术独有，国内领先。
本仪器采用核工业北京地质研究院具有自主知识产权的新一代数字化全波高速采集系统，最大数据通过率达700Kcps，避免了以往环境γ谱仪在测量中、高含量样品时对采集信号的漏记，不需死时间修正即可实现了仪器量程的宽范围。
仪器的软件包括测量部分和分析部分，其特点在于采用全波采集、全谱测量，全谱参加计算机数据处理，这是目前国内外低本底碘化钠多道γ能谱仪中最先进的谱数据处理技术。
仪器采用中英文界面进行人机对话，操作方便，测量数据准确可靠。
仪器的主要用途及特点：1、常规测量：检测粉末样品中天然放射性核素铀（238U）、镭（226Ra）、钍（232Th）、钾（40K）的比活度和含量。
2、无损检测：检测石材样品中天然放射性核素镭（226Ra）、钍（232Th）、钾（40K）的比活度，同时给出被测样品的内、外照射指数I Ra、Ir。
能充分发挥核辐射测量的优势——非破坏性、快速、方便和成本低等优点的现场无损检测。
3、室内环境氡浓度检测：采用活性炭盒自然吸附法测量室内空气中的氡浓度。
4、活性炭盒法测量土壤中气体的氡浓度：采用活性炭盒主动吸附取样分析法测量土壤中气体的氡浓度。
此外还可以广泛用于核物理、核化学、冶金、地质、海洋、石油、医疗卫生、环保、天文学、考古学、土壤学、生物学等领域。
二、系统功能及主要技术指标1、ADC（模数转换器）（1）道数： 4096、2048、1024、512 道，本系统固定采用2048道；（2）积分非线性：±0.1%；（3）微分非线性：±0.05%（MAX）；（4）输入极性：负。
低放射性钍的多道伽马能谱分析方法探讨
。 实 验研 究
中国 资源 综合 利用 Ｃｈｉｎａ Ｒｅｓｏｕｒｃｅｓ Ｃｏｍｐｒｅｈｅｎｓｉｖｅ Ｕｔｉｌｉｚａｔｉ 性钍 的多道伽 马能谱分析方法探讨
王瑞 芬 ，凤吾 利
（安徽 省核 工 业勘 查技 术 总院 ，安徽 芜 湖 ２４１０００）
伽马能谱分析是一种常规 的测量土壤 、岩石 中的 放射性核 素的物理方法 ，其特点是一次测量 可得到多 个放 射性 核素的含量 ，样 品制备 相对简 单 ，对于土壤 和岩石只需烘干 、粉碎后 ，装入标准样 品容器中 ，填实 、 定量密封 即可 ［１－２］ｏ钍元素的分析结果 ，会受到样 品粒 度 、含水 率 、样 品质量变化 、样 品高度 和谱线 测量时 间以及其他 因素 的影 响。为了减少各种 物理因素对室 内伽马能谱测量结 果 的影 响 】，人们必须对物 理 因素 进行试验分析 ，从而得到稳定准确的结果 。
１ 主 要 仪 器
主要仪器有鄂式破碎 机 、对辊式破碎机 、盘磨机 、 ＧＪ密封 式制 样机 、套 筛 、电子天平 、鼓风 电热 恒温 干燥箱 、烧杯 、空调 、ＡＲＬ多道 能谱仪 。
２ 试验部分
本文 主要考察样 品 的颗粒 度 、测量 时间 、含水 率 、 样 品质 量变 化 、样 品盒 的形状 、温 度对 测量 结 果 的 影 响 。 ２．１ 样 品粒度对测量结果 的影 响
康普顿散射虚拟仿真实验记录数据处理报告
康普顿散射虚拟仿真实验记录数据处理报告电子对效应是高能γ射线与物质相互作用的一种过程。
当γ射线入射至物质时，其能量足够高，能够转化成正负电子对。
这些电子对在物质中相互作用，产生电离作用，并在物质中形成电子对径迹。
电子对径迹在物质中的长度与能量有关，能量越高，径迹越短。
2.康普顿散射实验原理康普顿散射实验是利用康普顿效应测量γ光子能量及微分截面与散射角的关系。
实验装置主要包括放射源、闪烁体探测器、多道分析器和电子学系统等。
放射源发出γ光子，射线与物质相互作用后发生康普顿散射，散射光子被闪烁体探测器探测，多道分析器对探测到的信号进行处理，得到γ能谱。
通过测量γ能谱中康普顿边缘的位置和形状，可以计算出散射光子的能量和微分截面与散射角的关系。
三、实验步骤1.实验前准备：检查实验装置是否正常，调整探测器位置，调节放射源距离探测器的距离，确保实验安全。
2.测量γ能谱：打开实验装置电源，打开多道分析器软件，进行能谱测量。
记录康普顿边缘的位置和形状，计算出散射光子的能量和微分截面与散射角的关系。
3.测量吸收系数：更换不同物质，测量不同能量γ射线在典型物质中的吸收系数，记录实验数据。
4.实验结束：关闭实验装置电源，整理实验数据和记录。
四、注意事项1.实验过程中要注意辐射安全，避免直接接触放射源。
2.实验装置应调整好位置，确保测量精度和安全性。
3.实验数据应认真记录和整理，避免误差产生。
4.实验结束后应及时清理实验装置，保持实验室环境整洁。
当高于1.022MeV的γ光子穿过原子核时，它会在原子核的库仑场作用下转变成一个电子和一个正电子。
其中一部分光子的能量会转变成正负电子的静止能量，而其余部分则会成为它们的动能。
被释放出的电子还能与介质产生激发、电离等作用。
而正电子在失去能量后，会与物质中的负电子相遇并相互湮灭，产生γ射线。
探测这种湮灭辐射是可靠地确定正电子产生的实验方法之一。
闪烁体探测器是一种广泛应用的电离辐射探测器，利用电离辐射在某些物质中产生的闪光来进行探测。
《核数据处理》4 峰位与峰区的确定
i
原始能谱数据 方波函数变换后的数据 类高斯函数变换结果
4.6 各方法比较
表4.1 各种寻峰方法寻峰效果的比较
寻峰方法
弱峰检测能力 重叠峰分辨能力
简单比较法
很差
很差
一阶导数法
差
差
二阶导数法
差
好
三阶导数法
很差
最好
高斯乘积函数法 好
中等
协方差法
中等
好
对称零面积法 中等
好
核数据处理
第四章 峰位与峰区的确定
确定核素指纹，首先必须精确确定能谱中含有的特征射线 峰的位置。
寻峰算法应具备以下性能： ■ 能识别弱峰，特别是高康普顿平台上的弱峰 ■ 能有效地剔除假峰 ■ 比较高的重峰分辨能力
4.1 简单比较法
Counts
2000 1500
Fitted values of counts at characteristic peaks ,which have been discovered by the peak-searching method
图4-1 单峰导数图
图4-2 重叠的组合峰导数图
4.4 协方差寻峰法
设 yi j fi C j bi j ,i m, m
m
m
m
m
f g gg C gy C (ggC C ) g y ffii
m m w C j CCj jyiymij jm j mi
1000
4.2 高斯乘积函数找峰法
假设特征峰为
yi
A
2
imp
e 2 2
2
取能谱中四个点i-m-j, i-m, i, i+j，那么
j jm
伽玛能谱光滑方法综述
伽玛能谱光滑方法综述摘要：伽玛能谱仪探测到的伽玛能谱数据，因统计涨落而产生误差。
为了降低这一误差人们提出了多种伽玛能谱光滑处理方法。
本文概述了这些方法的基本原理与数学模型，比较它们的优缺点，为实际应用提供指导。
关键词：伽玛能谱；光滑；滤波；噪声引言由于核衰变和探测器中固有的统计涨落、电子学系统的噪声影响,所以在伽玛能谱的测量过程中，测得的谱数据不可避免带有很大的统计涨落和干扰噪声。
这对伽玛能谱的定性定量分析产生误差。
在伽玛能谱的分析中,为了减少能谱测量数据的统计涨落,又保留谱峰的全部重要的特征,以便可靠地定性和定量分析伽玛能谱,必须首先对实测伽玛能谱原始数据进行光滑或去噪处理。
本文对常用的光滑方法进行综述，为伽玛能谱分析软件中的谱光滑功能的实现提供指导。
基本思想由于能谱数据是按整数道址离散存储的，所以谱光滑处理是逐道进行的；以待处理道为中心，用其左右m 道的测量数据，对该道数据作修正，消除统计涨落的影响。
传统的伽玛谱光滑方法有：平均移动法；重心法；多项式最小二乘拟合法；离散函数褶积滑动变换法；傅立叶变换法。
近年来，小波变换法开始被应用到伽玛能谱的数据光滑中。
光滑方法及原理算术滑动平均法该方法思想如下:设i y 为待光滑的第i 道数据，左右各取m 道，则共有2m +1个点，用所有2m +1个点的算术平均值作为这道的修正值。
公式为：121m i i j j my y m +=-=+∑ (1) 式中i y 为原始谱数据，i y为光滑后的谱数据。
此方法两端各有m 个点得不到平滑，称为边沿损失。
重心法重心法就是选取加权因子和归一化因子，使光滑后的数据成为原来数据的重心。
由于道数是整数，没有半道的情况存在，若用2道的数据取重心，则第i 道计数的重心(平均值)为()11124i i i i y y y y -+=++ (2) 上式即为第i 道计数i y 的3点重心法光滑公式。
按照此推理的公式可以导出常用的5点、7点重心法等公式。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
自动寻峰由于谱结构的复杂和统计涨落的影响，从谱中正确地找到全部存在的峰是比较困难的。
尤其是找到位于很高本底上的弱峰，分辨出相互靠得很近的重峰更为困难。
谱分析对寻峰方法的基本要求如下：(1)比较高的重峰分辨能力。
能确定相互距离很近的峰的峰位。
(2)能识别弱峰，特别是位于高本底上的弱峰。
(3)假峰出现的几率要小。
(4)不仅能计算出峰位的整数道址，还能计算出峰位的精确值，某些情况下要求峰位的误差小于 0.2 道。
很多作者对寻峰方法进行了研究，提出了很多有效的寻峰方法。
目的：判断有没有峰存在确定峰位（高斯分布的数学期望），以便把峰位对应的道址，转换成能量确定峰边界——为计算峰面积服务（峰边界道的确定，直接影响峰面积的计算）分为两个步骤：谱变换和峰判定要求：支持手动 /自动寻峰，参数输入，同时计算并显示峰半高宽、精确峰位、峰宽等信息，能够区分康普顿边沿和假峰感兴区内寻峰人工设置感兴趣大小，然后在感兴区内采用简单方法寻峰重点研究：对感兴区内的弱峰寻峰、重峰的分解对于一个单峰区，当峰形在峰位两侧比较对称时，可以由峰的FWHM 计算峰区的左、右边界道址。
峰区的宽度取为3FWHM ，FWHM 的值可以根据峰位m p由测量系统的 FWHM刻度公式计算。
由于峰形对称，左、右边界道和峰位的距离都是 1.5FWHNM 。
m L INT (m p 1.5FWHM0.5)m R INT (m p 1.5FWHM0.5)式中 m p是峰位， INT 的含义是取整数。
对于存在有低能尾部的峰，其峰形函数描述（参见图）。
y m HEXP[ (m m p ) 2 / 22]，m≥mp－Jy m HEXP [ J (2m2m p J ) / 22 ] ，m≤mp－J式中 H 为峰高， mp 为峰位，是高斯函数的标准偏差，J 为接点的道址和峰位之间的距离。
在峰位的左侧，有一个接点，其道址为mp－J。
在接点的右侧，峰函数是高斯函数。
在接点的左侧，峰函数用指数曲线来描述。
这时峰区的左、右边界道址为m INT (mp 1.12FWHM 2/ J 0.5J 0.5)Lm R INT (m p 1.5FWHM0.5)带有低能尾部的峰函数的图形全谱自动寻峰基于核素库法：能量刻度完成后，根据核素库中的能量计算对应的道址，在各个道址附近（左右 10 道附近）采用简单的寻峰方法（导数法）方法：根据仪器选择开发IF 函数法 /简单比较法 (适于寻找强单峰，速度快)满足条件： data i m data i k data i data i m可认为有峰存在然后在 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)高斯乘积函数找峰法（可靠性差，不建议采用）描述谱峰形状的函数主要是高斯函数G (i )A exp(i i 0 ) 2 / 2 2则由相邻的数据点定义2一个新的函数（第一高斯乘积函数，只与FWHM 2.3556有关）：P m (i )G(i )G (i m 1)exp( 11.092m ) m2G (i 2)G (i m)H 2m 是步长（用道表示），是高斯乘积函数的阶数，则Pm(i)称为第 m 阶高斯乘积函数。
找峰的灵敏度与 m 有关，随 m 的增加灵敏度提高。
为避免基线参数的影响，最好扣除本底后，再应用高斯乘积函数找峰。
考虑统计涨落的影响，把判断无峰存在的 1 变为一个“单位带”。
即峰的判断为：1k /y i无峰P m (i )k /y i k 3(1) 有峰峰位的确定：由Pm(i)过 1 的两点求平均来确定；峰边界的确定：“单位带”下限的两个最端点；半高宽的确定：函数Pm(i)在“1上”的截距；组合峰的确定：在乘积函数的两个峰之间没有处于“带内”的乘积函数值导数法（一阶、二阶、三阶）y i'1mCjyi j N m j mNm 为规范化常数， Cj 平滑的变换系数。
3 次多项式 5 点光滑一阶导数公式：（可以采用）y i'1 ( y i 28 y i 18y i 1y i 2)峰位确定：一阶导数值由正变负=0处；峰边界确定：一阶导数12由负变正 =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~20peakposition[p++]=j+0.5;// 保持峰位对应的道址}5点光滑二阶导数公式（软件中推荐采用）''1y i(2 y i 2y i 1 2 y i y i 1 2 y i 2 )//7 点二阶导数(5*(countsdata[j-3]+countsdata[j+3])-3*(countsdata[j-1]+countsdata[j+1])-4*countsdata[j])/42;y i'1(22.0 y i 3 67.0 y i 2 58.0 y i 1 58.0 y i 1 67.0 y i 2 22.0 y i 3 )252.0软件中推荐采用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~20peakposition[p++]=j+0.5;// 保持峰位对应的道址 }1c 0j n iD ini 12n i n ijj1j kk1/ 2D ic j n i j ;c 2j n ijjkjkSignifican ce of 2nd Derivative :SD i /c j 100 p 2 j 2 exp 21 j 2 p 2p 2where p is the assumedpeak width.k : Go upto c j10- 6Peak ‘ found ’whenS > Threshold试验：系列 1 为处理后的原始能谱，系列 2 为 5 点一阶导数，系列 3 为 5 点二阶导数，系列 4 为对称零面积法寻峰只要选择好合适的寻峰阈值，足以满足准确寻找到全能峰，并剔除假峰（如康普顿边沿，反散射峰）5 点光滑三阶导数公式判定各感兴区是单峰还是重峰y i'''1( y i 2 2 y i 1 2 y i 1 y i 2 ) 2峰位确定：三阶导数由负变正 =0 处；峰边界确定：三阶导数由正变负=0 处判定峰是否有意义 0.8FWHM ≤ N ≤ 3FWHM峰高判定条件| y m |max TRH ym p e 0.5 /这个公式就是在一阶导数法寻峰程序中实际应用的峰高判定条件。
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;}}return n1;if ((nmin-nmax)>0.8*fwhm && (nmin-nmax)<3*fwhm)else return (0);}对称零面积法（推荐自动寻峰中采用，可探测弱峰和重峰）面积为零的“窗”函数与实验谱数据进行褶积变换，且要求“窗”函数为对称函数。
对线性基底的褶积变换将为零，只有存在峰的地方不为零。
y i 'mmC j yi jC j 0 C jC jjmjm匹配滤波器法 （类峰形函数） C jj 2 1mk 2 2 ] exp[2 ]exp[2m 221 k mm~C j data i jy ij m1f峰判定准则 R i ~y im2C 2j data i jjm2m+1 为变换宽度，FWHM 2.355 6 为峰宽参数， 若变换后的 y'和其均方根误差的比值超过预先给定的寻峰阈值（ f ），则认为找到了一个峰。