一种基于高斯拟合的光纤型表面等离子共振信号峰值定位检测方法
【技术领域】
[0001] 本发明属于SPR的信号处理领域，尤其涉及一种基于高斯拟合的光纤型表面等离 子共振信号峰值定位检测方法。
【背景技术】
[0002] 表面等离子共振（surface plasmon resonance, SPR)检测技术，是一项基于物理 光学现象的高灵敏检测技术，具有灵敏度高、样品无需标记，可实现在线检测等优点，已经 成为生物反应、药物筛选、医疗仪器、食品安全、环境监测等领域的一项重要工具。从最初的 棱镜耦合发展到现在的光栅耦合、波导耦合和光纤耦合，多样化的耦合方式为SPR检测技 术的发展助力。其中，光纤耦合凭借体积小，成本低，可以实现远距离检测、分布式检测等优 点，成为近年来SPR技术发展的一个热点。
[0003] 棱镜型SPR检测系统采集到的光谱曲线具有半峰宽较窄、共振峰较尖锐等特点， 而光纤型SPR的光谱数据恰恰不具备这些特点。在对棱镜型SPR信号处理上已有诸多方 法，常见的有一阶导数零点法、局部相似性匹配算法，以及由质心法所衍生的其他算法，比 如Johansen提出的加权质心法；Nenninger提出的追踪质心法；王晓萍提出的无基线质心 法等。一阶导数零点法、局部相似性匹配算法对噪声敏感；质心法及其改进算法对波形的 非对称性敏感，计算的共振波长会因为波形的不对称导致偏移，虽然研究的是共振波长的 偏移量，但每次计算共振波长时的信号不对称性并不能保持一致。因此，很多针对棱镜型的 数据处理方法并不适用于光纤型SPR光谱数据。在实验中，系统还会受到光源波动、环境温 度、机械振动等因素的影响，不可避免的会将噪声混入到理想信号中。选择合适的数据处理 方法不仅可以有效降低随机噪声所造成的影响，还可以提高数据的精度。
【发明内容】
[0004] 光纤型SPR光谱具有半峰宽（full width of half maximum, FWHM)较宽，共振峰 不尖锐等特征，传统的针对棱镜型SPR光谱的峰值定位检测算法无法准确定位此类光谱的 共振峰位置。为了准确计算光纤型SPR光谱的峰值位置（共振波长），本发明提出了基于信 赖域算法（Trust Region Algorithm)的高斯拟合方法，通过对SPR光谱的两次高斯拟合， 准确的定位共振波长的位置，本发明提供一种基于高斯拟合的光纤型表面等离子共振信号 峰值定位检测方法。
[0005] 本发明的技术方案：一种基于高斯拟合的光纤型表面等离子共振信号峰值定位检 测方法，其特征在于包括：
[0006] 步骤1 :选取合适参数的Savitzky-Golay滤波器对采集的SPR信号进行滤波降 噪；
[0007] 步骤2 :选取用于分析的SPR信号数据段，降低运算量，同时可以提高拟合的精 度；
[0008] 步骤3 :确定高斯拟合的峰数，在拟合精度和运算量之间选择一个合适的拟合峰 数；
[0009] 步骤4 :在确定峰数下，基于信赖域算法（Trust region Algorithm)对确定峰数 的高斯拟合系数的求解；
[0010] 步骤5 :在拟合得到的最优高斯拟合曲线上进行一维极值搜索，初步确定共振波 长的位置；
[0011] 步骤6 :以共振波长的位置为中心点，向前后各取一段数据再重新做一次多峰的 高斯拟合，结合信赖域算法，确定拟合函数；
[0012] 步骤7 :在第二次拟合的函数上，利用一维极值搜索来确定最终的共振峰位置。
[0013] 步骤1中的Savitzky-Golay滤波器是光谱分析领域的一种滤波方法，可最大限度 的保留相对极值和宽度等分布特性，调整其参数（拟合阶次、窗宽大小）便可以处理不同的 数据情况。
[0014] 步骤1中的选取合适参数的标准：滤波器参数为拟合阶数和窗宽，根据实时SPR谱 的特点优化滤波器参数，使滤波后光谱的均方误差（mean squared error, MSE)的Stein's 无偏风险估计（Stein' s unbiased risk estimate, SURE)最小。
[0015] 步骤2中的选取标准：检测时光源采用的是卤钨灯光源，光源强度主要集中于一 段波长范围内，信号采集的模块是CCD线阵，CCD的采集范围可能超过光强集中的范围。因 此主要选取光源强度集中的波长范围内的SPR曲线数据。
[0016] 步骤3中的选取标准：以均方根误差（root mean squared error, RMSE)和拟合优 度（R-square)为选取标准，在兼顾拟合效果的前提下，尽量选择较低的峰数。这样可以减 少运算的时间，同时可以为准确定位峰值奠定基础。
[0017] 步骤4中的信赖域算法是求解非线性优化问题的一种数值方法，具有整体收敛 性。在一维搜索中，从Xk点移动到下一点的过程可以表示为xk+akd k。akdk表示为dk方向 上的位移，记为sk。而信赖域算法是根据一定的原则，直接确定位移sk，同时，与一维搜索不 同的是，它并没有先确定搜索方向dk，原则如下：利用二次模型模拟目标函数f(x)，再用二 次模型计算出位移s，根据位移s可以确定下一点x+s，从而可以计算出目标函数的下降量 (下降是最优化的目标），再根据下降量来决定扩大信赖域或缩小信赖域，
[0019] 其中，自变量8是我们要求的位移，應:=样_公-个对称矩阵，且是 Hessian矩阵f/T.v)域其近似，Ak> 〇为置信半径。
[0020] 通过衡量二次模型与目标函数的近似程度，可以做出判定：
[0021] 第k此迭代的实际下降量为：Afk= f k-f(Xk+Sk)
[0022] 第k次迭代的预测下降量为：Δ mk=fk_m(sk)
[0023] 定义比值
[0024] 它用来衡量模型函数与目标函数f的一致性程度。如果根据"某种原则"确定的 位移能使目标函数值充分下降，则扩大信赖域；若不能使目标函数值充分下降，则缩小信赖 域。如此迭代下去，直到收敛。
[0025] 步骤6中的向前后各取一段数据再重新做一次多峰的高斯拟合，向前后各取一段 数据的选取标准：在计算精度和运算量之间选择一个合适的数据点数，再重新做一次多峰 的高斯拟合，多峰为步骤3中确定高斯拟合的峰数。
[0026] 本发明具有的有益效果是：
[0027] 1、本算法采用基于信赖域算法的高斯拟合法对采集到的SPR信号进行处理，算法 本身稳定性强、抑制噪声。
[0028] 2、本算法适应光纤型SPR信号的特点。
[0029] 3、本算法可对高浓度溶液确定的共振峰位。
[0030] 4、本算法具有较快的运算速度，有利于工程应用。
[0031] 5、本算法具有较高的灵敏度
【附图说明】
[0032] 图1是SPR光谱数据滤波前数据
[0033] 图2是SPR光谱数据滤波后数据
[0034] 图3是不同浓度甘油的标准SPR曲线
[0035] 图4是二种不同算法的灵敏度对比
【具体实施方式】
[0036] 下面结合附图对本发明的一种【具体实施方式】做出说明。
[0037] 本发明涉及一种基于高斯拟合的光纤型表面等离子共振信号峰值定位检测方法， 包括：
[0038] 步骤1选取合适参数的Savitzky-Golay滤波器对采集的SPR信号进行滤波降噪。
[0039] Savitzky-Golay滤波器是光谱分析领域的一种滤波方法，可最大限度的保留相对 极值和宽度等分布特性，调整其参数（拟合阶次、窗宽大小）便可以处理不同的数据情况。
[0040] 参数选择标准：滤波器参数为拟合阶数和窗宽，根据实时SPR谱的特点优化滤波 器参数，使滤波后光谱的均方误差（mean squared error, MSE)的Stein' s无偏风险估计 (Stein' s unbiased risk estimate, SURE)最小。
[0041] 光谱数据滤波前后如图I和图2所示，采集的SPR信号数据中存在光源噪声、散粒 噪声等噪声信号，选取窗宽为33、阶次为5的Savitzky-Golay滤波器来对采集的信号进行 预处理。通过滤波，可以有效降低噪声，SPR曲线更加平滑。
[0042] 步骤2合理选取用于分析的SPR信号数据段，降低运算量，同时可以提高拟合的精 度。选取标准：由于光源采用的是卤钨灯光源，光源强度主要集中于一段波长范围内。而信 号采集的模块是CCD线阵，CCD的采集范围可能超过光强集中的范围。因此主要选取光源 强度集中的波长范围内的SPR曲线数据。
[0043] 图1和图2中数据，SPR曲线在400nm和IlOOnm附近抖动非常严重，这是由于光源 输出功率主要集中在450nm-1000nm的范围，超出这个范围的CCD线阵检测到的信号受随机 噪声的影响。因此，本发明选择波长位于450nm-1000nm范围的光谱数据作为待拟合对象。
[0044] 步骤3合理确定高斯拟合的峰数，在拟合精度和运算量之间选择一个合适的拟合 峰数。选取标准：以均方根误差（root mean squared error, RMSE)和拟合优度（R-square) 为选取标准，在兼顾拟合效果的前提下，尽量选择较低的峰数。这样可以提高运算的时间， 同时可以为准确定位峰值奠定基础。
[0045] 理想的SPR光谱假设是由若干个单峰光谱相互叠加形成的，以高斯函数系作为光 谱曲线的基本函数形式，即将y = f (X)设定为高斯函数系，其中每一个高斯函数均由3个 参数决定，峰高A、峰位B和峰宽C。整个高斯函数系写作
[0047] 实验中获得的光谱曲线是由离散点构成的，对这些离散数据进行非线性拟合，也 就是将这些离散数据拟合成高斯