曾 真，王元慶

(南京大學(xué)電子科學(xué)與工程學(xué)院，江蘇南京210093)

1 引言

現(xiàn)有測(cè)相方法主要可以分為自動(dòng)數(shù)字測(cè)相法、數(shù)字同步解調(diào)法、頻域數(shù)字測(cè)相法。而頻域數(shù)字測(cè)相方法與其他測(cè)相方法相比有以下的優(yōu)點(diǎn):首先，頻域數(shù)字測(cè)相方法的測(cè)距精度很高，同等條件下頻域數(shù)字測(cè)相的測(cè)相誤差為數(shù)字同步解調(diào)方法的二分之一;且由DSP對(duì)數(shù)字序列做N點(diǎn)FFT變換，速度快，電路實(shí)現(xiàn)簡(jiǎn)單［1］;其次，由于設(shè)計(jì)為相位敏感系統(tǒng)，當(dāng)信號(hào)的幅度有所變化時(shí)，若使用DSP做信號(hào)處理，測(cè)相器有固定的相位分辨率，可以克服由于模擬信號(hào)的幅度變化引起的測(cè)相誤差;再次，頻域數(shù)字測(cè)相方法的測(cè)相速度快，計(jì)算量比較小，特別是在位數(shù)較多的時(shí)候，這一優(yōu)點(diǎn)更加明顯，在實(shí)際應(yīng)用中更能體現(xiàn)出實(shí)時(shí)的特性［2］?，F(xiàn)普遍使用的頻域數(shù)字測(cè)相方法為FFT譜分析法［2－3］，由于FFT和譜分析只能在有限區(qū)間內(nèi)進(jìn)行，這就不可避免地存在由于時(shí)域截?cái)喈a(chǎn)生的能量泄漏，而譜線的離散化也不可避免地產(chǎn)生柵欄效應(yīng)，使譜峰值變小，相位精度降低。為了提高相位測(cè)距的精度，本文選用了密集頻譜細(xì)化及校正的方法，并使用了傳統(tǒng)方法與本文方法對(duì)仿真的激光回波信號(hào)進(jìn)行對(duì)比研究。

2 相位式激光雷達(dá)測(cè)距原理

激光相位測(cè)距的基本原理是通過(guò)發(fā)射強(qiáng)度調(diào)制的激光照射目標(biāo)，使用高敏感度的光電器件，如APD，采集反射回來(lái)的信號(hào)。通過(guò)計(jì)算發(fā)射信號(hào)與回波信號(hào)之間的相位差Δφ，即可以得到目標(biāo)的距離值 d，計(jì)算公式如下［4］:

其中，c為光速;f0為激光的調(diào)制頻率;N為發(fā)射信號(hào)與回波信號(hào)相差的整數(shù)個(gè)波數(shù)。從式(1)可以看出，發(fā)射信號(hào)與回波信號(hào)之間的相位差Δφ的精度決定了測(cè)距的精度［5］。

加Hanning窗時(shí)，只進(jìn)行單諧波頻率成分的幅值恢復(fù)時(shí)的最大誤差仍高達(dá)15．3%，相位誤差更高達(dá)90°。而對(duì)于一個(gè)頻率為20 kHz的理想正弦信號(hào)，使用200 kHz的采樣頻率采256個(gè)點(diǎn)，此時(shí)頻率分辨率為頻率分辨率即 781．25 Hz，頻率結(jié)果誤差高達(dá)4%，相位的測(cè)量精度不高［7］。其中，相位的測(cè)量精度又是由譜線查找和相位計(jì)算精度決定，因而本文使用復(fù)調(diào)制細(xì)化選帶頻譜分析方法提高譜線查找精度，相位差校正法提高相位計(jì)算精度。

3 相位精確估計(jì)方法的選用

3．1 復(fù)調(diào)制細(xì)化選帶頻譜分析算法

復(fù)調(diào)制細(xì)化選帶頻譜分析方法使用復(fù)解析帶通濾波器，復(fù)調(diào)制移頻點(diǎn)數(shù)為傳統(tǒng)方法的1/D，較傳統(tǒng)復(fù)調(diào)制細(xì)化選帶頻譜分析方法計(jì)算量小。它不但降低了FFT的點(diǎn)數(shù)［8］，計(jì)算量小，且所需的內(nèi)存空間也少得多，而且也不需要進(jìn)行復(fù)雜的頻率調(diào)整［9－12］。由于激光相位測(cè)距中數(shù)據(jù)量大，實(shí)時(shí)處理要求高，因而采用一級(jí)FIR復(fù)解析帶通濾波器，能保證分析的精度誤差在1%以內(nèi)。具體的作法如下:先確定中心頻率及細(xì)化倍數(shù)，在頻率(f1－f2)范圍內(nèi)進(jìn)行頻率細(xì)化分析，則再構(gòu)造一個(gè)寬度為1/D的復(fù)解析帶通濾波器，對(duì)樣本信號(hào)作選抽濾波。選抽后的信號(hào)進(jìn)行復(fù)調(diào)制移頻，將細(xì)化的起始頻率移到零頻點(diǎn)，最后做N點(diǎn)FFT和譜分析，不需要進(jìn)行頻率調(diào)整就可以得到具有N條獨(dú)立譜線的細(xì)化頻譜。其數(shù)據(jù)處理的流程如圖1所示。最后得到的細(xì)化頻譜譜線間隔為普通FFT的1/D，即在幅度譜上找到相應(yīng)發(fā)射頻率的頻譜精度提高D倍［8］。復(fù)調(diào)制細(xì)化選帶頻譜分析方法可以細(xì)化譜線間隔，從而提高譜線查找的精度。

圖1 復(fù)調(diào)制譜細(xì)化算法的數(shù)據(jù)處理流程圖Fig．1 data processing flow chart ofmultiple modulation zoom spectrum analysis algorithm

3．2 相位差校正法

在常用的四種校正方法中，比值校正法和能量重心校正法是一種近似的校正方法，F(xiàn)FT+FT連續(xù)細(xì)化傅里葉變換法計(jì)算量大，而相位差法得到的是真實(shí)解且相位的校正不依賴于窗譜函數(shù)，因而采用相位差校正法。

先對(duì)原始單頻率成分信號(hào)采(L+N)點(diǎn)樣本，然后對(duì)這1到N和L到(L+N)兩段信號(hào)分別進(jìn)行作傅里葉變換，利用其對(duì)應(yīng)離散譜線譜峰處的相位差進(jìn)行校正。第一段的峰值譜線號(hào)i(i=0，1，2，…，N－1)，設(shè)頻率校正量為Δf=difs/N，其中di為歸一化的譜線號(hào)修正量［9］，其值范圍為［－ 0．5，0．5］。第二段進(jìn)FFT(或DFT)后與之對(duì)應(yīng)的峰值譜線號(hào)為i，T=N/fs，有相位差［7］:

設(shè)信號(hào)FFT的實(shí)部為Ri，虛部為Ii，由對(duì)稱窗函數(shù)相位特點(diǎn)，歸一化頻率校正量為di時(shí)，則校正后的相位為:

4 激光相位測(cè)距的仿真計(jì)算與比較

設(shè)定發(fā)射信號(hào)頻率f0=3 MHz，采樣頻率為fs=32 MHz，f1為信號(hào)歸一化頻率。加Hanning窗，進(jìn)行整型FFT變換，進(jìn)行選抽濾波之后做N=4096的FFT變換。利用計(jì)算機(jī)生成仿真的激光發(fā)射和回波信號(hào):

對(duì)于激光相位式測(cè)距，距離信息是由發(fā)射信號(hào)與接收信號(hào)的相位差計(jì)算得到的，因而對(duì)于測(cè)距有意義的是兩信號(hào)間的相位差值。使用密集頻譜細(xì)化及相位差校正法和傳統(tǒng)方法得到兩相位差值的分析結(jié)果(其中細(xì)化倍數(shù)為150倍，中心頻率為3 MHz，采用一級(jí)FIR濾波器，濾波半階數(shù)為500)。本文方法與傳統(tǒng)方法的仿真計(jì)算角度值如圖1所示，誤差結(jié)果如圖2所示，測(cè)距結(jié)果如表1表示。

圖2 本文方法與傳統(tǒng)方法分析結(jié)果Fig．2 results of the proposed method and the traditional one

圖3 本文方法與傳統(tǒng)方法分析誤差Fig．3 errors of the proposed method and the traditional one

表1 測(cè)距結(jié)果對(duì)比表Tab．1 comparing table of the ranging result

從以上數(shù)據(jù)可以看出，相位誤差高達(dá)10°，即對(duì)于3 MHz的發(fā)射頻率，測(cè)距精度為3 m，只能通過(guò)增加FFT的點(diǎn)數(shù)來(lái)提高精度。本文采用先進(jìn)行頻譜細(xì)化，再進(jìn)行頻譜校正的方法，誤差最大僅為1°即對(duì)于3 MHz的發(fā)射頻率，測(cè)距精度為0．3 m?？朔祟l域數(shù)字測(cè)相法分辨率不高的缺陷。

4 結(jié)論

在傳統(tǒng)的相位測(cè)距方法中，只能通過(guò)增加FFT的點(diǎn)數(shù)來(lái)增加相位計(jì)算的精度，使得運(yùn)算量增加、速度降低。相比較而言，頻譜細(xì)化算法中基于復(fù)解析帶通濾波器的復(fù)調(diào)制譜細(xì)化算法所需的運(yùn)算量少，計(jì)算速度最快;頻譜校正算法中相位差校正算法校正是一種精確地校正方式，校正精度很高;頻譜細(xì)化和頻譜校正技術(shù)能大大提高頻譜分辨率。在工程應(yīng)用中，本文方法運(yùn)用于激光相位測(cè)距中能大幅度提高頻域數(shù)字測(cè)相的精度。

［1］ Zhang Ting，et al．Study of phase-shift laser measuring based on matlab［J］．Laser ＆ Infrared，2010，40(1):22 －27．(in Chinese)張婷，等．基于Matlab的相位式激光測(cè)距研究［J］．激光與紅外，2010，40(1):22 －27．

［2］ Hui Zhang．The realization of laser distance-measuring equipment on embedded system［D］．Shanghai:Tongji U-niversity，2007．(in Chinese)張慧．相位法激光測(cè)距算法的嵌入式系統(tǒng)實(shí)現(xiàn)［D］．上海:同濟(jì)大學(xué)，2007．

［3］ Xie Lei，Li Ji．Realization of digital phase detector based on FFT algorithm in laser distance measurement system［J］．Chinese Journal of Quantum Electronics，2003，20(1):85 －88．(in Chinese)謝蕾，李季．基于FFT的激光測(cè)距數(shù)字相位計(jì)的實(shí)現(xiàn)［J］．量子電子學(xué)報(bào)，2003，20(1):85 －88．

［4］ Wang Tao．Research of the phase laser ranging techonology［J］．Laser ＆ Infrared，2007，37(1):29 － 31．(in Chinese)汪濤．相位激光測(cè)距技術(shù)的研究［J］．激光與紅外，2007，37(1):29 －31．

［5］ Gao Song．The research on some key questions about digital phase laser range measurement system［D］．Dalian:Dalian Maritime University，2007．(in Chinese)高篙．?dāng)?shù)字鑒相式激光測(cè)距系統(tǒng)幾個(gè)關(guān)鍵問(wèn)題的研究［D］．大連:大連海事大學(xué)，2007．

［6］ Jin Yanhua．Study of a novel zoom spectrum analysis，international conference on communications and mobile computing［J］．Communications and Mobile Computing，2009，(1):367－371．

［7］ Ding Kang，et al．Synthesized correctingmethod of phase difference on discrete spectrum［J］．Journal of Vibration Engineering，2002，15(1):114 －119．(in Chinese)丁康，等．離散頻譜綜合相位差校正法［J］．振動(dòng)工程學(xué)報(bào)，2002，15(1):114 －119．

［8］ Jing Gai．Frequency estimation of uncooperative coherent pulse radars ［J］．Military Communications Conference，2007:1－7．

［9］ QiGuoqing．Analysis of the influence of observation starting time on FFT-based phase estimation of sinusoid［J］．Signal Processing(ICSP)，2010:26 －29．

［10］ Jiang Bo，etal．Research and realization ofmultiplemodulation ZoomFFT algorithm in partial Zoom spectrum analysis［J］．，2010，(7):48 －49．(in Chinese)江波，等．基于復(fù)調(diào)制的ZoomFFT算法在局部頻譜細(xì)化中的研究與實(shí)現(xiàn)［J］．大眾科技，2010，(7):48 －49．

［11］ Ding Kang，Xie Ming，et al．Principle and method ofmultiplemodulation zoom spectrum analysis based onmultiple analytical band-pass filter［J］．Journal of Vibration Engineering，2001，14(1):30 －35．(in Chinese)丁康，謝明，等．基于復(fù)解析帶通濾波器的復(fù)調(diào)制細(xì)化譜分析原理和方法［J］．振動(dòng)工程學(xué)報(bào)，2001，14(1):30－35．

［12］ Shang Haitao，Xue Hongxi．Study on the ZoomFFT algorithm based on band-pass filter［J］．Process Automation Instrumentation，2011，32(4):10 －12．(in Chinese)尚海濤，薛紅喜．基于帶通濾波器的ZoomFFT算法研究［J］．自動(dòng)化儀表，2011，32(4):10 －12．