基于拟合算法的傅里叶变换光谱仪光谱反演

摘要：在傅里叶变换光谱仪（FTS）中，精确获取每个干涉信号采样点处的光程差是获取光谱图的关键。动镜速度的波动会对等时间采样造成采样误差，而当目标光源波长较短时无法直接利用传统方法进行等光程差采样。分析了速度波动误差对光谱的影响，并提出了一种基于拟合算法计算干涉信号每个采样点处光程差的方法，在动镜速度波动较大的情况下对可见光波段的干涉信号进行采样反演。实验结果表明，此方法准确度高，适用于各种傅里叶变换光谱仪。

关键词：傅里叶变换光谱仪； 光谱反演； 拟合算法； 可见光谱

中图分类号： TP 911.73 文献标志码： A doi： 10.3969/j.issn.10055630.2015.06.007

Abstract：In the Fourier transform spectrometer （FTS）， obtaining the optical path difference of each sampling points accurately is the key of data processing. The velocity fluctuation of moving mirror will cause equaltime sampling error， and traditional equal optical path difference sampling method can not be used for visible wavelengths. The influence of speed fluctuation is analyzed and a method to calculate the optical path difference of each sampling points based on fitting algorithm is developed. The method is verified by processing the interference signal of visible light with large velocity fluctuation. The experimental results verify its correctness and effectiveness. This method can also be used in any kind of Fourier transform spectrometer.

Keywords： FTS； spectrum recovery； fitting algorithm； vis spectrum

引 言

傅里叶变换光谱仪（FTS）因其具有高光谱分辨率、高光通量[12]、杂散光影响小、波数精确度高等优点在化学分析、大气探测等领域应用越来越广泛。其中可以获得非常高的光谱分辨率是傅里叶变换光谱仪最重要的优点，而光谱分辨率虽然受限于仪器自身，但是干涉信号采样误差和反演算法带来的误差同样会影响最终的光谱分辨率[34]。在傅里叶变换光谱仪中，随着动镜的移动，两干涉光束产生光程差，从而产生一个随光程差变化的干涉信号。干涉信号和光谱数据是一对傅里叶变换对[4]，对获取到的干涉信号以光程差为变量做傅里叶变换即可得到入射光的光谱数据，因此要进行光谱反演就需要在采样干涉信号的同时能精确地获取每个干涉信号采样点处的光程差。

如果动镜运动速度为绝对匀速，那么对干涉信号进行等时间采样等同于等光程差采样。但在实际情况中，动镜速度不可能是绝对匀速的，如果继续采用等时间采样则不可避免会导致非常大的采样误差，从而在得到的光谱图中产生误差[5]。常用的方法是采用激光作为参考信号，对其进行滤波整形后作为干涉信号采样的触发信号，这种方法只适用于目标光源波长大于参考激光两倍的傅里叶变换光谱仪中，例如中长波红外傅里叶变换光谱仪，当目标光源在近红外和可见光波段时，继续采用上述方法则会因不满足奈奎斯特采样定理而无法从离散的干涉信号中复原出光谱数据。本文基于曲线拟合算法提出了一种精确获取每个干涉信号等时间采样点处光程差的方法，仍然采用激光作为参考信号，对于波长在可见光波段甚至紫外波段的目标光源，即使速度波动较大，只要等时间采样频率满足采样定理，本方法都可以获取每个数据点处的光程差，从而对干涉信号进行精确反演获取光谱数据。

1 傅里叶变换光谱仪干涉信号分析

傅里叶变换光谱仪属于调频的干涉光谱仪，主要用来观测光谱信息，常用的傅里叶变换光谱仪采用的是经典迈克尔逊结构，光程差是动镜行程的两倍，其光学结构如图1所示。动镜采用角镜，这样可以消除动镜倾斜造成的影响[6]，而且在相同的动镜行程下可以获取四倍的光程差，从而得到更高的光谱分辨率。入射光被分束器分成强度相等的两束光，经分束器反射的光束经过左侧的定镜1和定镜2到达角镜左侧，然后反射回来，经分束器反射和投射分别到达探测器2和探测器1；而经分束器透射的光束经右侧的定镜3和定镜4到达角镜右侧，然后反射回来，经分束器反射和投射到达探测器1和探测器2。当动镜沿着光轴来回做直线运动时，在探测器上汇合的两光束的光程差发生周期性变化，形成干涉信号。在傅里叶变换光谱仪中，动镜作匀速扫描的过程，相当于在整个采样过程中以速度u对干涉光信号进行调制，其结果是把频率很高的光波调制成频率很低的电信号[7]。干涉信号被红外探测器转换为电信号，经过放大器和模拟滤波器，进入模数转换器进行等光程差采样，得到数字干涉图。

式（6）中后两项就是所谓的鬼线。如果采用等时间采样采集上述干涉信号，再直接对其做傅里叶变换，在所得到的光谱数据中会有额外的尖峰出现。如果入射光源为连续光谱，则所得到的光谱图会有很多鬼线出现。常用的方法是采用稳频激光作为参考信号，这是因为激光单色性好而且波长λ0已知，动镜每移动λ0/4便产生一个干涉信号周期，即每当激光干涉信号出现过零点时就知道光程差变化了λ0/2，因此只要以参考激光干涉信号过零点作为采样触发信号即可实现等光程差采样[9]。然而从采样点中复原出原谱线必须保证采样频率满足奈奎斯特采样定理，如果目标光源的波长λ全都大于激光波长的两倍，那么当光程差每变化一个λ，激光干涉信号已经至少出现了4个过零点，此时满足奈奎斯特采样定理；但是当目标光源的波长小于激光波长的两倍时，光程差每变化一个λ，激光干涉信号只出现一个过零点，甚至不出现。即对于工作波段为0.4～1.0 μm的傅里叶变换光谱仪，继续采用上述方法进行等光程采样则会造成干涉信息丢失，无法得到完整干涉信号，更无法得到光谱信息，因此如何对可见光干涉信号进行等光程差采样是能否反演出可见光光谱的关键。

2 利用曲线拟合算法获取光程差

在现在智能的控制方案下，虽然动镜往复运动的实际速度是不断变化的，但其速度不会存在突变的现象，因此在傅里叶变换光谱仪中，不可能存在动镜移动0.10～0.25 μm（λ/4）的过程中速度发生很大变化的情况，一般在动镜移动这么短的距离过程中，速度基本可以看作是匀速的。那么在光程差变化一个激光波长的过程中，产生的激光干涉信号应该是一个标准的带有初始相位的单频率正弦信号。

对于采集到的离散信号，很少能直接求得其函数表达式，一般是采用插值和拟合的方法，利用采集到的离散数据点来得出一条近似正确的连续的曲线。如果采样得到的离散信号点没有误差则采用插值方法，如果采样得到的离散信号点与真实值有差距则采用曲线拟合方法。常用的获取曲线的方法有样条插值、多项式插值、多项式拟合、基于遗传算法的拟合算法等，其中多项式拟合是最小二乘拟合的一种常用形式。最小二乘法是应用最广泛的曲线拟合算法，其核心思想是寻找合适的函数参数使得函数与所采样得到的数值之间的误差平方和达到最小[10]，即使下式达到最小值：

式中：N为采样总数。当用标准正弦信号去拟合时，即f（x）=Asin（2πfx+φ）+B，其中幅值A、频率f、相位φ、直流分量B四个参数都是未知的[11]，而且误差平方和是这四个参数的非线性函数，因此无法求出误差平方和的闭合解，只能用迭代法求出每个参数的局部最优解。

本文基于最小二乘法用f（x）=a+bsin（cx）+dcos（cx）对每个激光干涉信号周期进行正弦拟合，求出每个周期对应正弦信号的四个参数和拟合误差，如果能用该方法拟合出结果而且拟合误差非常小，则可近似认为：在光程差变化一个激光波长的过程中，产生的激光干涉信号是一个标准的单频率正弦信号。在傅里叶变换光谱仪中，由于参考激光干涉信号和目标光源干涉信号所经过电路的延迟不同，即使是采用等光程差采样仍然会有采样误差产生，根据文献[1213]的计算，只要速度波动的相对误差小于2%，这种误差对光谱的影响基本可以忽略。因此只要每个干涉信号周期的拟合误差小于2%，则基本可以忽略。对一个完整干涉信号的拟合流程图如图2所示。

3 实验验证

为验证本文所提出的方法，对所用傅里叶变换光谱仪进行速度开环控制，在速度有较大波动的情况下进行干涉信号采样。动镜采用无刷直流电机驱动，电机本身速度较快，经过减速箱减速，动镜运动的平均速度控制在3 cm/s，根据干涉信号频率和速度、波数的关系，可得干涉信号的频率范围为100～400 kHz，因此选择采样频率为2 MS/s，采样位数为16位。

参考激光采用波长为0.685 2 μm的稳频氦氖激光器，因对可见光波段气体的吸收率较低，为方便进行结果对比，目标光源分别采用波长为0.659 8 μm的红色激光和波长为0.532 μm的绿色激光。对动镜移动一个单程所获得的数据进行处理可以得到实际速度波动情况，如图3（a）所示，计算其波动的相对误差值为10.04%。拟合误差如图3（b）所示，最大拟合误差为1.30%，拟合误差的平均值为0.37%，满足误差小于2%的要求。图4和图5分别为用本文提出的算法对波长为0.659 8 μm和波长为0.532 μm的光源进行光谱反演的结果，并给出峰值光谱处的细节，从实验结果中可以看出，当速度变化较大时，采用本文所用的方法所得到的反演结果波数准确度高，而且基本可以达到仪器自身的光谱分辨率0.05 cm-1。

4 结 论

傅里叶变换光谱仪的动镜速度不可避免地会有波动，对等时间采样的干涉信号直接进行傅里叶变换将会产生很大的光谱误差。当目标光源为中长波红外光时，可以用参考激光干涉信号触发采样以实现等光程差采样，但是当目标光源的波长较短时，无法直接进行等光程差采样。本文提出了一种基于曲线拟合的光程差获取算法，能在速度波动较大时仍然能准确地反演出光谱图，而且适用于任何波段的目标光源，同时对速度均匀性无严格要求，降低了对控制系统精度的要求。

参考文献：

[1]夏翔.干涉光谱仪多速度模式扫描技术研究[D].上海：中国科学院研究生院（上海技术物理研究所），2014：1115.

[2]吴航行，华建文，王模昌.新型红外空间遥感用傅里叶变换光谱仪[J].红外与激光工程，2004，33（4）：397400.

[3]盛灏.傅里叶变换光谱仪干涉信号数据获取研究[D].上海：中国科学院研究生院（上海技术物理研究所），2014：1421.

[4]范世福.光谱技术和仪器的新发展[J].光学仪器，2000，22（4）：3540.

[5]王明艳.基于DSP的弹光调制非均匀干涉信号获取及反演算法研究[D].太原：中北大学，2013：3135.

[6]王文桂.干涉光谱仪[M].北京：宇航出版社，1988：180190.

[7]杨隆梓，段星辉，魏焕东，等.星载傅立叶光谱仪数据采集技术研究[J].科学技术与工程，2008，8（8）：20342037.

[8]胡盛雯，殷德奎.高光谱大气探测仪干涉信号的获取电路设计[J].红外与激光工程，2006，35（S1）：528533.

[9]胡玲，潘征宇，洪扁，等.双拼激光干涉仪中的数字相位计设计[J].光学仪器，2014，36（3）：258262.

[10]齐国清，吕健.正弦曲线拟合若干问题探讨[J].计算机工程与设计，2008，29（14）：36773680.

[11]梁志国，朱济杰，孟晓风.四参数正弦曲线拟合的一种收敛算法[J].仪器与仪表学报，2006，27（11）：15131519.

[12]WILLIAMS C S，Mirror misalignment in fourier spectroscopy using a Michelson interferometer with circular aperture[J].Applied Optics，1966，5（6）：10841085.

[13]COHEN D L，Performance degradation of a Michelson interferometer when its misalignment angle is a rapidly varying，random time series[J].Applied Optics，1997，36（18）：40344042.

[14]姚涛，吕群波，相里斌，等.非均匀采样干涉数据光谱反演技术研究[J].光谱学与光谱分析，2010，30（5）：14301433.

[15]吕群波，姚涛，相里斌，等.干涉数据光谱反演方法研究[J].光谱学与光谱分析，2010，30（1）：114117.