李濤　華建文　孫曉杰　陳仁

摘要：在傅里葉變換光譜儀（FTS）中，精確獲取每個(gè)干涉信號(hào)采樣點(diǎn)處的光程差是獲取光譜圖的關(guān)鍵。動(dòng)鏡速度的波動(dòng)會(huì)對(duì)等時(shí)間采樣造成采樣誤差，而當(dāng)目標(biāo)光源波長(zhǎng)較短時(shí)無(wú)法直接利用傳統(tǒng)方法進(jìn)行等光程差采樣。分析了速度波動(dòng)誤差對(duì)光譜的影響，并提出了一種基于擬合算法計(jì)算干涉信號(hào)每個(gè)采樣點(diǎn)處光程差的方法，在動(dòng)鏡速度波動(dòng)較大的情況下對(duì)可見(jiàn)光波段的干涉信號(hào)進(jìn)行采樣反演。實(shí)驗(yàn)結(jié)果表明，此方法準(zhǔn)確度高，適用于各種傅里葉變換光譜儀。

關(guān)鍵詞：傅里葉變換光譜儀； 光譜反演； 擬合算法； 可見(jiàn)光譜

中圖分類號(hào)： TP 911.73 文獻(xiàn)標(biāo)志碼： 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]、雜散光影響小、波數(shù)精確度高等優(yōu)點(diǎn)在化學(xué)分析、大氣探測(cè)等領(lǐng)域應(yīng)用越來(lái)越廣泛。其中可以獲得非常高的光譜分辨率是傅里葉變換光譜儀最重要的優(yōu)點(diǎn)，而光譜分辨率雖然受限于儀器自身，但是干涉信號(hào)采樣誤差和反演算法帶來(lái)的誤差同樣會(huì)影響最終的光譜分辨率[34]。在傅里葉變換光譜儀中，隨著動(dòng)鏡的移動(dòng)，兩干涉光束產(chǎn)生光程差，從而產(chǎn)生一個(gè)隨光程差變化的干涉信號(hào)。干涉信號(hào)和光譜數(shù)據(jù)是一對(duì)傅里葉變換對(duì)[4]，對(duì)獲取到的干涉信號(hào)以光程差為變量做傅里葉變換即可得到入射光的光譜數(shù)據(jù)，因此要進(jìn)行光譜反演就需要在采樣干涉信號(hào)的同時(shí)能精確地獲取每個(gè)干涉信號(hào)采樣點(diǎn)處的光程差。

如果動(dòng)鏡運(yùn)動(dòng)速度為絕對(duì)勻速，那么對(duì)干涉信號(hào)進(jìn)行等時(shí)間采樣等同于等光程差采樣。但在實(shí)際情況中，動(dòng)鏡速度不可能是絕對(duì)勻速的，如果繼續(xù)采用等時(shí)間采樣則不可避免會(huì)導(dǎo)致非常大的采樣誤差，從而在得到的光譜圖中產(chǎn)生誤差[5]。常用的方法是采用激光作為參考信號(hào)，對(duì)其進(jìn)行濾波整形后作為干涉信號(hào)采樣的觸發(fā)信號(hào)，這種方法只適用于目標(biāo)光源波長(zhǎng)大于參考激光兩倍的傅里葉變換光譜儀中，例如中長(zhǎng)波紅外傅里葉變換光譜儀，當(dāng)目標(biāo)光源在近紅外和可見(jiàn)光波段時(shí)，繼續(xù)采用上述方法則會(huì)因不滿足奈奎斯特采樣定理而無(wú)法從離散的干涉信號(hào)中復(fù)原出光譜數(shù)據(jù)。本文基于曲線擬合算法提出了一種精確獲取每個(gè)干涉信號(hào)等時(shí)間采樣點(diǎn)處光程差的方法，仍然采用激光作為參考信號(hào)，對(duì)于波長(zhǎng)在可見(jiàn)光波段甚至紫外波段的目標(biāo)光源，即使速度波動(dòng)較大，只要等時(shí)間采樣頻率滿足采樣定理，本方法都可以獲取每個(gè)數(shù)據(jù)點(diǎn)處的光程差，從而對(duì)干涉信號(hào)進(jìn)行精確反演獲取光譜數(shù)據(jù)。

1 傅里葉變換光譜儀干涉信號(hào)分析

傅里葉變換光譜儀屬于調(diào)頻的干涉光譜儀，主要用來(lái)觀測(cè)光譜信息，常用的傅里葉變換光譜儀采用的是經(jīng)典邁克爾遜結(jié)構(gòu)，光程差是動(dòng)鏡行程的兩倍，其光學(xué)結(jié)構(gòu)如圖1所示。動(dòng)鏡采用角鏡，這樣可以消除動(dòng)鏡傾斜造成的影響[6]，而且在相同的動(dòng)鏡行程下可以獲取四倍的光程差，從而得到更高的光譜分辨率。入射光被分束器分成強(qiáng)度相等的兩束光，經(jīng)分束器反射的光束經(jīng)過(guò)左側(cè)的定鏡1和定鏡2到達(dá)角鏡左側(cè)，然后反射回來(lái)，經(jīng)分束器反射和投射分別到達(dá)探測(cè)器2和探測(cè)器1；而經(jīng)分束器透射的光束經(jīng)右側(cè)的定鏡3和定鏡4到達(dá)角鏡右側(cè)，然后反射回來(lái)，經(jīng)分束器反射和投射到達(dá)探測(cè)器1和探測(cè)器2。當(dāng)動(dòng)鏡沿著光軸來(lái)回做直線運(yùn)動(dòng)時(shí)，在探測(cè)器上匯合的兩光束的光程差發(fā)生周期性變化，形成干涉信號(hào)。在傅里葉變換光譜儀中，動(dòng)鏡作勻速掃描的過(guò)程，相當(dāng)于在整個(gè)采樣過(guò)程中以速度u對(duì)干涉光信號(hào)進(jìn)行調(diào)制，其結(jié)果是把頻率很高的光波調(diào)制成頻率很低的電信號(hào)[7]。干涉信號(hào)被紅外探測(cè)器轉(zhuǎn)換為電信號(hào)，經(jīng)過(guò)放大器和模擬濾波器，進(jìn)入模數(shù)轉(zhuǎn)換器進(jìn)行等光程差采樣，得到數(shù)字干涉圖。

式（6）中后兩項(xiàng)就是所謂的鬼線。如果采用等時(shí)間采樣采集上述干涉信號(hào)，再直接對(duì)其做傅里葉變換，在所得到的光譜數(shù)據(jù)中會(huì)有額外的尖峰出現(xiàn)。如果入射光源為連續(xù)光譜，則所得到的光譜圖會(huì)有很多鬼線出現(xiàn)。常用的方法是采用穩(wěn)頻激光作為參考信號(hào)，這是因?yàn)榧す鈫紊院枚也ㄩL(zhǎng)λ0已知，動(dòng)鏡每移動(dòng)λ0/4便產(chǎn)生一個(gè)干涉信號(hào)周期，即每當(dāng)激光干涉信號(hào)出現(xiàn)過(guò)零點(diǎn)時(shí)就知道光程差變化了λ0/2，因此只要以參考激光干涉信號(hào)過(guò)零點(diǎn)作為采樣觸發(fā)信號(hào)即可實(shí)現(xiàn)等光程差采樣[9]。然而從采樣點(diǎn)中復(fù)原出原譜線必須保證采樣頻率滿足奈奎斯特采樣定理，如果目標(biāo)光源的波長(zhǎng)λ全都大于激光波長(zhǎng)的兩倍，那么當(dāng)光程差每變化一個(gè)λ，激光干涉信號(hào)已經(jīng)至少出現(xiàn)了4個(gè)過(guò)零點(diǎn)，此時(shí)滿足奈奎斯特采樣定理；但是當(dāng)目標(biāo)光源的波長(zhǎng)小于激光波長(zhǎng)的兩倍時(shí)，光程差每變化一個(gè)λ，激光干涉信號(hào)只出現(xiàn)一個(gè)過(guò)零點(diǎn)，甚至不出現(xiàn)。即對(duì)于工作波段為0.4～1.0 μm的傅里葉變換光譜儀，繼續(xù)采用上述方法進(jìn)行等光程采樣則會(huì)造成干涉信息丟失，無(wú)法得到完整干涉信號(hào)，更無(wú)法得到光譜信息，因此如何對(duì)可見(jiàn)光干涉信號(hào)進(jìn)行等光程差采樣是能否反演出可見(jiàn)光光譜的關(guān)鍵。

2 利用曲線擬合算法獲取光程差

在現(xiàn)在智能的控制方案下，雖然動(dòng)鏡往復(fù)運(yùn)動(dòng)的實(shí)際速度是不斷變化的，但其速度不會(huì)存在突變的現(xiàn)象，因此在傅里葉變換光譜儀中，不可能存在動(dòng)鏡移動(dòng)0.10～0.25 μm（λ/4）的過(guò)程中速度發(fā)生很大變化的情況，一般在動(dòng)鏡移動(dòng)這么短的距離過(guò)程中，速度基本可以看作是勻速的。那么在光程差變化一個(gè)激光波長(zhǎng)的過(guò)程中，產(chǎn)生的激光干涉信號(hào)應(yīng)該是一個(gè)標(biāo)準(zhǔn)的帶有初始相位的單頻率正弦信號(hào)。

對(duì)于采集到的離散信號(hào)，很少能直接求得其函數(shù)表達(dá)式，一般是采用插值和擬合的方法，利用采集到的離散數(shù)據(jù)點(diǎn)來(lái)得出一條近似正確的連續(xù)的曲線。如果采樣得到的離散信號(hào)點(diǎn)沒(méi)有誤差則采用插值方法，如果采樣得到的離散信號(hào)點(diǎn)與真實(shí)值有差距則采用曲線擬合方法。常用的獲取曲線的方法有樣條插值、多項(xiàng)式插值、多項(xiàng)式擬合、基于遺傳算法的擬合算法等，其中多項(xiàng)式擬合是最小二乘擬合的一種常用形式。最小二乘法是應(yīng)用最廣泛的曲線擬合算法，其核心思想是尋找合適的函數(shù)參數(shù)使得函數(shù)與所采樣得到的數(shù)值之間的誤差平方和達(dá)到最小[10]，即使下式達(dá)到最小值：

式中：N為采樣總數(shù)。當(dāng)用標(biāo)準(zhǔn)正弦信號(hào)去擬合時(shí)，即f（x）=Asin（2πfx+φ）+B，其中幅值A(chǔ)、頻率f、相位φ、直流分量B四個(gè)參數(shù)都是未知的[11]，而且誤差平方和是這四個(gè)參數(shù)的非線性函數(shù)，因此無(wú)法求出誤差平方和的閉合解，只能用迭代法求出每個(gè)參數(shù)的局部最優(yōu)解。

本文基于最小二乘法用f（x）=a+bsin（cx）+dcos（cx）對(duì)每個(gè)激光干涉信號(hào)周期進(jìn)行正弦擬合，求出每個(gè)周期對(duì)應(yīng)正弦信號(hào)的四個(gè)參數(shù)和擬合誤差，如果能用該方法擬合出結(jié)果而且擬合誤差非常小，則可近似認(rèn)為：在光程差變化一個(gè)激光波長(zhǎng)的過(guò)程中，產(chǎn)生的激光干涉信號(hào)是一個(gè)標(biāo)準(zhǔn)的單頻率正弦信號(hào)。在傅里葉變換光譜儀中，由于參考激光干涉信號(hào)和目標(biāo)光源干涉信號(hào)所經(jīng)過(guò)電路的延遲不同，即使是采用等光程差采樣仍然會(huì)有采樣誤差產(chǎn)生，根據(jù)文獻(xiàn)[1213]的計(jì)算，只要速度波動(dòng)的相對(duì)誤差小于2%，這種誤差對(duì)光譜的影響基本可以忽略。因此只要每個(gè)干涉信號(hào)周期的擬合誤差小于2%，則基本可以忽略。對(duì)一個(gè)完整干涉信號(hào)的擬合流程圖如圖2所示。

3 實(shí)驗(yàn)驗(yàn)證

為驗(yàn)證本文所提出的方法，對(duì)所用傅里葉變換光譜儀進(jìn)行速度開(kāi)環(huán)控制，在速度有較大波動(dòng)的情況下進(jìn)行干涉信號(hào)采樣。動(dòng)鏡采用無(wú)刷直流電機(jī)驅(qū)動(dòng)，電機(jī)本身速度較快，經(jīng)過(guò)減速箱減速，動(dòng)鏡運(yùn)動(dòng)的平均速度控制在3 cm/s，根據(jù)干涉信號(hào)頻率和速度、波數(shù)的關(guān)系，可得干涉信號(hào)的頻率范圍為100～400 kHz，因此選擇采樣頻率為2 MS/s，采樣位數(shù)為16位。

參考激光采用波長(zhǎng)為0.685 2 μm的穩(wěn)頻氦氖激光器，因?qū)梢?jiàn)光波段氣體的吸收率較低，為方便進(jìn)行結(jié)果對(duì)比，目標(biāo)光源分別采用波長(zhǎng)為0.659 8 μm的紅色激光和波長(zhǎng)為0.532 μm的綠色激光。對(duì)動(dòng)鏡移動(dòng)一個(gè)單程所獲得的數(shù)據(jù)進(jìn)行處理可以得到實(shí)際速度波動(dòng)情況，如圖3（a）所示，計(jì)算其波動(dòng)的相對(duì)誤差值為10.04%。擬合誤差如圖3（b）所示，最大擬合誤差為1.30%，擬合誤差的平均值為0.37%，滿足誤差小于2%的要求。圖4和圖5分別為用本文提出的算法對(duì)波長(zhǎng)為0.659 8 μm和波長(zhǎng)為0.532 μm的光源進(jìn)行光譜反演的結(jié)果，并給出峰值光譜處的細(xì)節(jié)，從實(shí)驗(yàn)結(jié)果中可以看出，當(dāng)速度變化較大時(shí)，采用本文所用的方法所得到的反演結(jié)果波數(shù)準(zhǔn)確度高，而且基本可以達(dá)到儀器自身的光譜分辨率0.05 cm-1。

4 結(jié) 論

傅里葉變換光譜儀的動(dòng)鏡速度不可避免地會(huì)有波動(dòng)，對(duì)等時(shí)間采樣的干涉信號(hào)直接進(jìn)行傅里葉變換將會(huì)產(chǎn)生很大的光譜誤差。當(dāng)目標(biāo)光源為中長(zhǎng)波紅外光時(shí)，可以用參考激光干涉信號(hào)觸發(fā)采樣以實(shí)現(xiàn)等光程差采樣，但是當(dāng)目標(biāo)光源的波長(zhǎng)較短時(shí)，無(wú)法直接進(jìn)行等光程差采樣。本文提出了一種基于曲線擬合的光程差獲取算法，能在速度波動(dòng)較大時(shí)仍然能準(zhǔn)確地反演出光譜圖，而且適用于任何波段的目標(biāo)光源，同時(shí)對(duì)速度均勻性無(wú)嚴(yán)格要求，降低了對(duì)控制系統(tǒng)精度的要求。

參考文獻(xiàn)：

[1]夏翔.干涉光譜儀多速度模式掃描技術(shù)研究[D].上海：中國(guó)科學(xué)院研究生院（上海技術(shù)物理研究所），2014：1115.

[2]吳航行，華建文，王模昌.新型紅外空間遙感用傅里葉變換光譜儀[J].紅外與激光工程，2004，33（4）：397400.

[3]盛灝.傅里葉變換光譜儀干涉信號(hào)數(shù)據(jù)獲取研究[D].上海：中國(guó)科學(xué)院研究生院（上海技術(shù)物理研究所），2014：1421.

[4]范世福.光譜技術(shù)和儀器的新發(fā)展[J].光學(xué)儀器，2000，22（4）：3540.

[5]王明艷.基于DSP的彈光調(diào)制非均勻干涉信號(hào)獲取及反演算法研究[D].太原：中北大學(xué)，2013：3135.

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

[7]楊隆梓，段星輝，魏煥東，等.星載傅立葉光譜儀數(shù)據(jù)采集技術(shù)研究[J].科學(xué)技術(shù)與工程，2008，8（8）：20342037.

[8]胡盛雯，殷德奎.高光譜大氣探測(cè)儀干涉信號(hào)的獲取電路設(shè)計(jì)[J].紅外與激光工程，2006，35（S1）：528533.

[9]胡玲，潘征宇，洪扁，等.雙拼激光干涉儀中的數(shù)字相位計(jì)設(shè)計(jì)[J].光學(xué)儀器，2014，36（3）：258262.

[10]齊國(guó)清，呂健.正弦曲線擬合若干問(wèn)題探討[J].計(jì)算機(jī)工程與設(shè)計(jì)，2008，29（14）：36773680.

[11]梁志國(guó)，朱濟(jì)杰，孟曉風(fēng).四參數(shù)正弦曲線擬合的一種收斂算法[J].儀器與儀表學(xué)報(bào)，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]姚濤，呂群波，相里斌，等.非均勻采樣干涉數(shù)據(jù)光譜反演技術(shù)研究[J].光譜學(xué)與光譜分析，2010，30（5）：14301433.

[15]呂群波，姚濤，相里斌，等.干涉數(shù)據(jù)光譜反演方法研究[J].光譜學(xué)與光譜分析，2010，30（1）：114117.

（編輯：劉鐵英）