一种半导体激光器光谱自适应峰值搜索方法
【技术领域】
[0001] 本发明涉及一种半导体激光器光谱自适应峰值搜索方法。
【背景技术】
[0002] 目前，光纤通信技术由于具有传输容量大、传输速率高、抗电磁干扰、保密性强等 优点，已经在军事、工业、农业、医疗和科学研究等各个领域得到了飞速发展和广泛应用。法 布里-珀罗型半导体激光器（FP-LD)和分布反馈式半导体激光器（DFB-LD)因输出光相干 性和方向性好、谱宽窄、体积小且简单易用，在高速、大容量光纤通信系统中是最常用的光 源，其光谱通常用光谱分析仪测量得到，光谱带宽和纵模数等参数是FP-LD最重要的光谱 参数，边模抑制比和边模偏置等参数反映着DFB-LD的性能，光谱峰值搜索是计算FP-LD和 DFB-LD上述各种光谱参数的基础，实时准确的搜索出峰值对于光谱特性分析和开发高性能 光谱分析仪器都具有十分重要的意义。
[0003] 目前，常用的峰值搜索方法主要有：直接比较法、导数法、对称零面积法、遗传算 法、小波变换等。直接比较法采用直接比较的方式找出比左右相邻点大的数据点作为峰值 点，其原理简单，计算速度快，适合寻找平滑曲线的孤立强峰值，但在复杂的谱线分布中寻 找弱峰时会产生大量虚假峰，抗噪性较差。导数法的基本思想是把谱线看为一条连续曲线， 通过对谱线上各点求导数，依据导数的性质确定谱线峰值的位置，该方法对于光滑曲线具 有较高的搜索准确性，计算速度基本可以满足实时的要求，但对于噪声较大的复杂谱线易 产生虚假峰，因此需要首先对谱线进行平滑去噪，或者设置阈值参数过滤掉虚假峰。对称零 面积法是把面积为零且对称的窗函数与待搜索光谱进行卷积变换，可以证明，面积为零的 对称窗函数，对线性基底的卷积变换结果将为零，只有在存在峰的条件下不为零，且变换值 最大处为峰值所在，该方法峰值搜索效果取决于窗函数选取是否合适，除此之外还需要设 置窗宽、峰的最小半宽、搜索阈值等参数，计算较为复杂，运算速度也较慢。遗传算法能够进 行全局最优搜索，其抗噪能力较强，但该方法需要设置种群规模、初始种群、迭代次数、交叉 概率、变异概率等参数，且需要较长的训练时间来确定最优参数，不适合实时运算。小波变 换法将信号分解成一系列小波函数的叠加，构建一个具有良好时域及频域局部化的时频信 号，从而将时域寻峰转变为小波域对小波系数矩阵的寻峰，该方法搜索准确度较高，并且对 噪声和背景具有较强的抑制能力，但该方法需要选择小波母函数、尺度参数、平移参数等， 计算量较大，也不适合实时运算。
[0004] 由此可见，现有方法对上述半导体激光器光谱进行峰值搜索时存在如下问题：（1) 为得到较高的搜索准确度，大部分方法都需要针对不同类型、不同测量条件下得到的光谱 波形设置不同的最优阈值或参数，人为干预较多，方法无法做到自适应，鲁棒性较差；（2) 为过滤掉虚假峰，需要进行滤波去噪，平滑，背景基线漂移去除等预处理，不仅步骤繁琐，计 算复杂，而且在预处理过程中又不可避免地引入新的阈值或参数；（3)方法本身计算时间 开销比较大，再考虑数据预处理时间花费，无法满足实时性的计算要求。
【发明内容】
[0005] 本发明的目的在于提出一种半导体激光器光谱自适应峰值搜索方法，其采用如下 方案：
[0006] -种半导体激光器光谱自适应峰值搜索方法，包括如下步骤：
[0007] a局部最大矩阵的建立
[0008] 设T代表待搜索半导体激光器的实测光谱数据序列，序列长度为n，T = [tpty h，...tj ;
[0009] 使用窗口宽度N可变的移动窗，比较序列中某一个数据点在N范围内和前向边缘 点及后向边缘点的大小；
[0010] 比较的逻辑关系值作为波动因子，以波动因子作为数据元素建立局部最大矩阵 M(mXn)；
[0011] 其中，m= [length(T)/2]-l = [n/2]-l，符号" □"代表向下取整，窗口宽度 N = 2i(i = 1，2, ...，m);
[0012] 对于每一维的i和j，其中，I < i < m，I < j < n，当i+1 < j < n-i时，波动因子 表示为：
⑴
[0014] 当I < j < i和n-i+1 < j < η时，波动因子表示为：
[0015] Hiij j= rand+1 (2)
[0016] 其中，rand表示在[0, 1]范围内的均匀分布随机数；局部最大矩阵M表示为：
⑶
[0018] b局部最大维度控制
[0019] 统计局部最大矩阵中所有行向量0元素的数目，并找到0元素数目最多的行，记为 P ：
[0020] p = max (Sum (Iiiii j = 0)) (4)
[0021] 其中，max表示最大值，Sum表示求和；
[0022] 重新修正局部最大矩阵M，去除掉行数i大于p的元素，形成新的矩阵M' (pXn):
(S):
[0024] c谱峰定位
[0025] 新矩阵M'的列向量元素值是判别数据序列中某元素是否为局部峰值点的直接判 断标准，当某一列元素全为0时，该列对应的数据序列中的元素即为候选峰值点；
[0026] d谱峰补偿
[0027] 如果数据点对应的新矩阵M'列向量同时满足下列三个条件，则认为该数据点符 合在一定窗口内左侧邻近点上升且右侧邻近点下降这种趋势，加入候选峰值点内：
[0028] (1)首行元素为0 ;
[0029] (2)从首行开始，连续0元素的个数N1S [n/2];
[0030] (3)以该列为中心，窗口为[n/4]范围内的左右各列中0元素的个数N2S [n/2];
[0031] 通过上述步骤即可在无需设置阈值或参数的前提下，准确定位光谱数据峰值点。
[0032] 本发明具有如下优点：
[0033] 本发明方法无需设置阈值或参数，只根据光谱数据序列固有的特征即可较为准确 的搜索出峰值，即具有较强的自适应性；本发明方法对不同条件下测量得到的DFB-LD和 FP-LD光谱均适用，即具有较强的鲁棒性；本发明方法不需要对数据进行滤波去噪，平滑， 背景基线漂移去除等预处理，计算简单易行，检验准确率高；本发明方法计算时间开销小， 完全能够满足实时计算的要求，便于后续谱特性的快速分析和计算。
【附图说明】
[0034] 图1为FP-LD的实测光谱曲线图；
[0035] 图2为DFB-LD的实测光谱曲线图；
[0036] 图3为本发明中一种半导体激光器光谱自适应峰值搜索方法的流程图；
[0037] 图4为使用直接比较法对FP-LD光谱进行峰值搜索的效果图；
[0038] 图5为使用导数法对FP-LD光谱进行峰值搜索的效果图；
[0039] 图6为使用遗传算法对FP-LD光谱进行峰值搜索的效果图；
[0040] 图7为使用本发明自适应峰值搜索方法对FP-LD光谱进行峰值搜索的效果图；
[0041] 图8为使用直接比较法对DFB-LD光谱进行峰值搜索的效果图；
[0042] 图9为使用导数法对DFB-LD光谱进行峰值搜索的效果图；
[0043] 图10为使用遗传算法对DFB-LD光谱进行峰值搜索的效果图；
[0044] 图11为使用本发明自适应峰值搜索方法对DFB-LD光谱进行峰值搜索的效果图；
[0045] 图12为使用直接比较法的检验准确率示意图；
[0046] 图13为使用导数法的检验准确率示意图；
[0047] 图14为使用遗传算法的检验准确率不意图；
[0048] 图15为使用本发明自适应峰值搜索方法的检验准确率示意图；
[0049] 图16为直接比较法、导数法、遗传算法和本发明自适应峰值搜索方法计算时间对 比图。
【具体实施方式】
[0050] 本发明的基本思想是：首先研究了搜索对象--FP-LD和DFB-LD光谱实测数据的 特征；然后提出一种基于局部最大矩阵的半导体激光器光谱自适应峰值搜索方法，最后通 过对比实验进行分析和检验方法效果。具体的，
[0051] 下面结合附图以及【具体实施方式】对本发明作进一步详细说明：
[0052] (一）FP-LD和DFB-LD光谱实测数据特征
[0053] 本发明研究的两种光源光谱实测数据具有不同的特点，FP-LD具有多纵模特性， DFB-LD由于采用波长选择反馈机制，因而其光谱具有单纵模特性。
[0054] 同时，由于光栅分光光谱仪是把由复色光照明的狭缝经过光栅分光系统变为若干 个单色的狭缝像，实际测得的光谱会根据狭缝宽度的不同具有不同的谱线展宽，通过设定 光谱仪不同的分辨带宽，实测光谱会有不同的光谱带宽（FMffH)，其光谱外形或轮廓呈高斯 型。
[0055] 如图1和图2分别为FP-LD和DFB-LD的实测光谱。通过对实测25组DFB-LD和 FP-LD光谱数据分析，待搜索谱峰数据序列具有如下特性：
[0056] (1)数据中电平值较低的部分随机波动大，存在均方根值为常数的随机噪声信号， 即信噪比低，易造成虚假峰的出现。
[0057] (2)存在明显大于噪声的高幅值数据点，即动态范围大。该部分数据包含的大谱峰 较易搜索，但部分光谱数据的大谱峰处曲线不光滑，易造成误判。
[0058] (二）基于局部最大矩阵的半导体激光器光谱自适应峰值搜索方法
[0059] 基于上述数据序列特征，本发明提出了一种基于局部最大矩阵的自适应峰值搜索 方法（简称自适应局部最大法），其方法流程如图3所示：
[0060] 一种半导体激光器光谱自适应峰值搜索方法，包括如下步骤：
[0061] a局部最大矩阵的