静态迈克尔逊干涉仪是在传统迈克尔逊干涉仪的基础上改进的一种空间调制型光谱获取装置，将传统动镜采用一个倾斜很小角度的静镜代替，从而产生联系变化的光程差，最终获得待测光源的光谱信息。由于其没有机械移动或扫描结构，所以结构简单紧凑，同时具有良好的稳定性、抗震动和抗干扰能力，这为其应用在更广泛的生产生活领域提供了有利条件[1-2]。与此相似的静态型光谱仪还有静态傅里叶变换光谱仪[3]、沃拉斯顿分光光谱仪[4]、静态Sagnac三角光谱仪[5]、光栅光谱仪[6]等。无论哪种静态光谱仪都需要一个能够实时采集干涉条纹并且实时处理分析，得到光谱数据的处理系统，否则将丧失静态光谱仪的一大优势——实时性。如果想实现实时处理的功能，要求在采集数据及处理分析算法都相应提高，这也是本文的主要研究内容。

    FPGA(Field Programmable Gate Array)为现场可编程门阵列器件[7]，是在GAL、PAL、CPLD等可编程芯片的基础上更集成化的产物。由于其本身就是硬件，所以具有数据处理快、灵活性高、集成度好等优点[8]。它克服了定制电路的缺点，又解决了原有可编程器件的门电路数有限的问题，从而被广泛地应用在数字信号处理领域。          目前，对于静态干涉条纹的采集及处理，主要有三种处理手段： (1)干涉条纹灰度数据由数据采集卡采集并传递给PC机[9]，采用MATLAB、LabVIEW等分析软件实现对干涉条纹的光谱分析； (2)采用定制的专用FFT芯片[10]，再从CCD等光电器件输出后直接进入FFT芯片完成数据处理；(3)采用可编程数字器件，如DSP、FPGA等，通过编程实现采集干涉条纹及FFT等算法的数据处理[11]。三种方法各有特点，方法(1)使用简单、在PC机上实现的功能全，但不能脱机工作、系统复杂，便携性差；方法(2)结构紧凑、效率高、实时性好，但灵活性差、可移植性差，并且对于大点数的FFT不适合采用此法；方法(3)从灵活度、结构性能、可移植性等方面看都相当适中，是目前应用开发较多的一种形式。针对本系统研究内容要求较高的实时性，选用硬件编程器件FPGA作为数字信号处理芯片，可以实现更高速、更灵活的系统性能。
1 静态迈克尔逊干涉仪的原理及结构
    静态迈克尔逊干涉仪的结构原理如图1所示。静态迈克尔逊干涉系统是在传统迈克尔逊干涉系统的基础上，将原有动镜变为倾斜一定角度的固定反射镜，如图中反射镜1所示。当光线1入射系统后，由分束镜分为两束光，一部分由反射镜1反射回柱面镜，这束光由于反射镜1存在一定的夹角而非原路返回，其会与光线2经反射镜2反射的光线发生相干，形成干涉条纹。其他光线也以此类推，与其他光线发生干涉，最终在整个柱面镜上形成静态干涉条纹,再由柱面镜汇聚到CCD上，采集得到干涉条纹的灰度数据。CCD采集的数据传入FPGA中，在硬件FPGA芯片中完成干涉条纹的滤波、去噪，再通过FFT算法将干涉条纹中的频谱信息提取出来，最终显示出其结果。

    在实际情况中，干涉条纹通过引入修正因子，可将光源强度表示为:

2 FPGA处理系统
2.1 总体设计
    对于傅里叶变换的光谱获取需要进行实时处理，即当得到光源干涉条纹的信号后，处理系统需要在处理当前数据时同时接收下一帧的干涉条纹数据，这就需要系统具备高速的数据处理能力。可编程逻辑器件FPGA具有高速并行处理能力，可完成系统实时采集及处理的要求。干涉条纹的灰度数据由CCD探测器采集传入FPGA芯片中，在芯片中完成滤波去噪、快速傅里叶变换、相位校正及标定等算法。系统采用AViiVA M1型线阵CCD探测器,包括1 024个像元，8/10 bit数字输出,可探测波长范围为400~1 100 nm。设计了实时光谱采集及片上处理的系统，其原理图如图2所示。

　　模块一为干涉条纹的采集部分，由驱动控制对AViiVA M1型线阵CCD探测器进行时钟脉冲的提供和控制。当系统开始运行时，其提供起始控制信号（st）、时钟信号（clk）、触发信号（trgi及trgd）等。模拟输出信号可通过示波器观察采集得到的静态干涉条纹，数字化后传给FPGA。模块二为光谱分析处理模块，为了使处理速度达到较高水平，设计了串并转换方式，将数据导入双端口RAM缓存中，当FPGA芯片控制数据并行读取时，数据输入FFT模块并处理，实现数据的快速傅里叶变换。变换后的频谱信息经过预先编写的相位校正和光谱标定程序而获得准确的光谱分布。
2.2 切趾处理
    由于傅里叶变换过程是针对无限大数据的，但实际输入的数据量是有限的，所以干涉图样的采样就相当于与一个矩形函数相乘，也就是频域中光谱和矩形函数的SINC卷积。但由于这样做会导致光谱的失真变形，所以通常采用加窗的处理方式，窗函数选取的好可以使光谱尽量少地泄露。
    在MATLAB仿真软件中，通过对一组660 nm激光干涉条纹数据进行加窗处理，分别采用汉明窗、三角窗、加布莱克曼窗与不加窗的情况进行对比，得到如图3所示的频谱相应。
    切趾函数的选择需要注意主瓣宽度尽量窄，旁瓣幅值尽量低。由图可以看出，加三角窗中心处产生隔断，而加布莱克曼窗时峰峰值受到抑制，相比之下，加汉明窗比较理想，所以本系统中采用汉明窗。窗体程序是先由MATLAB计算汉明窗的值并用二进制补码表示，然后在FPGA的某一个ROM存储器中存放初始化的窗函数值，最后通过乘法器完成切趾功能。
2.3 FFT算法模块

    FPGA设计过程中通常采用VHDL或Verilog HDL编写硬件语言实现时序及逻辑电路。但当输入的数字信号相对复杂时,采用该种方法对编程带来的难度也会加大，硬件工作效率会有所降低。所以在面对复杂问题时，通常采用IP核或者System Generator等算法软件。本文采用IP核完成快速傅里叶变换，采用全精度不缩减及缩减定点算法完成1 024点干涉条纹数据的傅里叶变换。本系统采用基2算法结构，16 bit数据输入与输出，通过设置SCH值确定每级运算后的右移位数。基2算法在整个运算过程中只需采用一个蝶形单元，当第一级数据输入运算模块后，下一级的数据将被放入双端口RAM中,采用乒乓缓存操作处理可以提高工作效率，同时，节约了FPGA的片上资源，实际只需3个硬件乘法器加3块双端口RAM即可实现高速运算，以50 MHz为例，进行1 024点的快速傅里叶变换只需要100 μs左右。
3 实验
    按以上要求搭建实验平台分别采用静态迈克尔逊干涉仪、Atmel公司的AViiVA型12 bit深的线阵CCD探测器、Xilinx公司Virtex VP40型芯片,检测光谱范围为300～1 200 nm，采用1 024个单元作为采集单元，分辨率为10 bit,光源选择FLDT3V10A-LD型660 nm红光激光器。
    通过采集干涉条纹，经滤波、去噪、快速傅里叶变换算法、相位校正、光谱标定等程序功能后，实验的5组数据如表1所示。对比光谱分析仪选用布鲁克公司的Q2-ION型光谱仪。

    分析表1可知，采用本系统得到的光谱分布相比Q2-ION型光谱仪的精度略低一些，但其多次检测产生的最大误差也在3%以下，说明可以基本准确地得到干涉条纹的频谱信息。并且相比Q2-ION型光谱仪而言，由于是采用FPGA硬件编程实现，所以在处理速度上要优于该产品，整个系统周期小于1 s，适应高速实时的处理系统。
    为了提高光谱分析仪器的速度,设计了一种基于FPGA硬件编程完成干涉条纹到光谱分布的处理系统。系统采用滤波、去噪、切趾等处理干涉条纹信息，然后通过快速傅里叶变换完成时频变换，再由相位校正及标定最终输出准确的光谱信息。实验采用静态迈克尔逊干涉仪、660 nm半导体激光器、线阵CCD探测器等完成，与Q2-ION型光谱仪的实验数据进行对比。实验结果显示，误差均在要求范围内，同时，本系统还具备处理速度高、实时性好、灵活性强的优点。
参考文献
[1] ACOSTA Eva, CHAMADOIRA Sara, BLENDOWSKE Ralf. Modified point diffraction interferometer for inspection and evaluation of ophthalmic components[J]. OSA, 2006, 23(3): 632-637.
[2] Jiang Meng, Guan Zuguang, He Sailing. Multiplexing scheme for self-interfering long-period fiber gratings using a low-coherence reflectometry[J]. IEEE Sensors Journal,2007, 28(31):1283-1286.
[3] 张敬波. 基于DSP的傅里叶光谱仪测量控制系统[D].长春：吉林大学,2004.
[4] 孙方,代作晓,华建文,等. 一种傅里叶变换光谱仪动镜速度测量系统[J]. 半导体光电，2007,2(4)：244-248.
[5] 李肖廷，李仰军，王高.基于FPGA马赫-泽德干涉仪的光谱数据采集系统[J].激光与红外,2011,4(1):155-161.
[6] YU K， LEE D. Micro-machined fourier transform spectrometer on silicon optical bench platform[J]. Sensors and Actuators, 2006,130(131):523-530.
[7] 贾立哲, 魏利辉. 相关干涉仪测向算法的FPGA设计实现[J]. 无线电工程, 2006,12(1):12-18.
[8] 杨天博, 罗迅, 李达成.一种基于FPGA的外差干涉仪动态数据采集方法[J].光学技术, 2003,6(2):123-126.
[9] 李丞,娄淑琴,郑毅,等. 基于FPGA的实时互相关运算器[J]. 北京交通大学学报, 2010,2(5):46-49.
[10] 殷世民,相里斌,周锦松, 等.基于FPGA的干涉式成像光谱仪实时数据处理系统研究[J]. 红外与毫米波学报, 2007,4(6):890-896.
[11] 闫敬业,吴季,张升伟, 等.合成孔径辐射计多路高阶数复相关器方案研究[J]. 电子学报,2003,9(1):13-19.