移相相關(guān)法計算相位差的研究

劉玉周，趙 斌

（華中科技大學(xué)機械科學(xué)與工程學(xué)院儀器系，武漢430074）

移相相關(guān)法計算相位差的研究

劉玉周，趙 斌*

（華中科技大學(xué)機械科學(xué)與工程學(xué)院儀器系，武漢430074）

為了提高相位式測距儀的測量精度，采用移相相關(guān)方法來估計兩同頻正弦信號的相位差。首先將每路信號移相2π后和原信號做相關(guān)來計算自相關(guān)，以減少噪聲的影響；其次用少許數(shù)據(jù)初步估算相位差，并將一路信號移相，使兩路信號的相位差移到π/2（或3π/2）附近；然后用較多的采樣數(shù)據(jù)計算兩路信號的相位差，將結(jié)果再減去移相量得到最終的相位差。同時分析了頻率誤差對相位差計算精度的影響，進(jìn)行了理論分析和仿真實驗驗證。結(jié)果表明，該方法計算的誤差大大減小。這對提高測距儀的測量精度是有幫助的。

測量與計量；移相相關(guān)法；相位差；頻率誤差

引 言

相位式激光測距在3-D成像［1］、機器人導(dǎo)航［2］、表面檢測［3］等領(lǐng)域有著廣泛的應(yīng)用，它通過測量光波往返的相位差來計算時間延遲從而計算待測距離［4-5］。與模擬測相方法相比，數(shù)字方法測量相位差具有不受環(huán)境干擾、兩路信號間無串?dāng)_等優(yōu)點，但也有不足。過零比較法［6］在信號頻率較高、噪聲較大時誤差較大［7］；當(dāng)實際的傅里葉變換的最大值出現(xiàn)在兩采樣點之間時，快速傅里葉變換［8］會有較大誤差，若頻率擾動，插值快速傅里葉變換法也很難達(dá)到較高的精度［4］；擬合法［9］計算相位差的缺點是計算時間太長并且迭代可能不收斂［10］。在頻率較高時，商業(yè)化相位差測量儀器分辨率很難超過0.05°，而有些文獻(xiàn)中的方法只是達(dá)到了較高的分辨率，而不是較高的精度［4］。相比上面所述的相位差測量方法，相關(guān)法計算相位差［11］被認(rèn)為是最優(yōu)的時間延遲算法［12-13］，采用快速算法［14］后計算速度快，且具有抑制噪聲能力較強并可計算中頻信號的相位差等優(yōu)點。但在采用相關(guān)算法的測距試驗中發(fā)現(xiàn)，距離連續(xù)移動時，非常靠近0°或180°相位差出現(xiàn)的幾率很低，并且在理論相位差處于0°或180°附近時，可能會產(chǎn)生1°以上的相位差誤差。作者對相關(guān)法估計相位差的理論分析后，提出一種移相相關(guān)法，可消除上述缺陷，進(jìn)一步提高相位差測量精度。

1 移相相關(guān)法計算相位差的原理及分析

1.1 普通相關(guān)法估計相位差的誤差分析

設(shè)有頻率為f0的兩個正弦信號x（k）和y（k），幅值分別為A和B，初相位分別為φx和φy，相位差φ=φx-φy?，F(xiàn)對其進(jìn)行采樣，采樣頻率為fs=nf0，整數(shù)n≥3，則數(shù)字頻率f=f0/fs。對x（t）和y（t）均做采樣長度N0同步整周期采樣后的數(shù)據(jù)序列x（k）和y（k）為（k為序號）：

式中，噪聲Nx（k）和Ny（k）一般是互不相關(guān)的加性白噪聲，均值為0，方差為σx2和σy。則普通相關(guān)法按照下式得估計相位差：

這表明普通相關(guān)法對含有噪聲的信號的相位差估計不是無偏估計。

（3）式是關(guān)于隨機變量Nx（k）和Ny（k）的函數(shù)，計算出自相關(guān)的方差為和，互相關(guān)的方差為，按照2維隨機變量非線性函數(shù)方差的計算方法［15］，可得到普通相關(guān)法計算估計相位差的方差：

因此，普通相關(guān)法估計相位差的均方差（mean square error，MSE）E（?φ）為：

在（6）式中取N0=480×103，R=30dB時，均方誤差與實際相位差φ的關(guān)系圖見圖1。

（6）式表明，普通相關(guān)法估計相位差有兩點不足：（1）增加采樣長度N0只能減小（6）式中的v（?φ）項，而偏離量b項不變，精度無法再提高；（2）誤差隨相位差真值φ變化而不同，當(dāng)實際相位差越靠近0°或180°時，估計的偏離量越大。

在激光相位式測距中，采用普通相關(guān)法估計相位差的不足具體表現(xiàn)為：當(dāng)待測距離連續(xù)變化時，測距誤差周期性變化，且很少出現(xiàn)359°～360°，0°～1°和179°～181°內(nèi)的相位差。

1.2 移相相關(guān)法計算相位差的原理

移相是指將采樣序列x（k）或y（k）的前面的幾個數(shù)據(jù)去掉，即采樣序列向左移位，從而改變其初相位的方法。

移相相關(guān)法在計算相位差時，先后有兩次移相：一是用移相360°后的信號與原信號做相關(guān)，代替原來的自相關(guān)運算，簡稱為移相自相關(guān)；二是以較少的數(shù)據(jù)初步估計相位差后，將一路信號移相Δθ，使兩路信號的相位差約為90°或者270°，將此時估計的相位差減去附加移相量Δθ后得到實際相位差，簡稱為移相計算相位差。

1.2.1 移相自相關(guān) 由于不同時間的隨機噪聲是不相關(guān)的，因此，若將信號與向左移位一個周期的信號做相關(guān)，可消除自相關(guān)中的噪聲項：

（9）式是相位差φ的無偏估計。按（9）式估計相位差的均方誤差等于方差v（?φ）：

同樣固定N0=480×103，R=30dB時，（10）式中均方誤差隨相位差真值的變化情況見圖2。對于固定的采樣長度N0和信噪比R，（10）式與（6）式的第2項（方差項）數(shù)量級相當(dāng)，但其比（6）式少了b這一項，因此其真值在90°和270°以外的均方誤差要小得多，在90°和270°附近均方誤差變化也平緩一些，但在0°和180°附近的相位差誤差仍然較大。

1.2.2 移相計算相位差 注意到（10）式在φ=90°和φ=270°時有極小值，這是由于估計相位差（9）式中反余弦的特性導(dǎo)致的。為避免在0°和180°附近較大的計算誤差，首先用較少的數(shù)據(jù)初步估算相位差，然后將y（k）或x（k）移相Δθ（Δθ≤90°），從而使相位差在90°或270°附近，再以較多采樣數(shù)據(jù)重新計算相位差，將結(jié)果減去Δθ得到精確相位差。

信號的一個周期內(nèi)有n個采樣點，移相后的相位差相距90°（或270°）不超過360°/n。因此在不考慮信號頻率誤差時，經(jīng)過上述兩次移相后計算的相位差均方誤差為下式：

1.2.3 頻率誤差引起的相位差估計誤差 晶振的老化和頻率漂移均會引起信號頻率的誤差，頻率誤差將導(dǎo)致兩次移相不準(zhǔn)確，從而產(chǎn)生相位差估計誤差。

若信號數(shù)字頻率的相對誤差（準(zhǔn)確度）為δf，理論數(shù)字頻率為f，實際數(shù)字頻率為f（1+δf），則不考慮噪聲時，x（k）移相360°后的自相關(guān)為：

將上兩式的Rx′（0）和Ry′（0）代入相位差的估計（9）式，得到移相自相關(guān)引起的相位差誤差φf1為：

另外，頻率誤差引起的附加相位差Δθ的誤差為：

式中，ABS表示取絕對值。

因此若在程序中未校正頻率，則頻率誤差引起的總的移相誤差Δφf為：

對于不同數(shù)字頻率準(zhǔn)確度δf和一個信號周期內(nèi)的采樣點數(shù)n，對應(yīng)的最大移相誤差Δφf，max見圖3。在相位式激光測距應(yīng)用中，采用溫補晶振或恒溫晶振驅(qū)動數(shù)字鎖相環(huán)等方法產(chǎn)生的信號的頻率穩(wěn)定度較高，一般其數(shù)字頻率準(zhǔn)確度優(yōu)于10×10-6，當(dāng)n≥6時，由頻率誤差引起的最大移相誤差小于0.002°。即使頻率準(zhǔn)確度為20×10-6，n≥48時，由頻率誤差引起的最大移相誤差小于0.002°，小于預(yù)期精度0.01°的1/5。

在考慮噪聲和頻率誤差時，移相相關(guān)法估計相位差的均方誤差為：

在理論相位差0°至360°整個區(qū)間內(nèi)，移相相關(guān)法估計的均方誤差小于-48dB，即相位差估計誤差小于0.004°。

與普通相關(guān)法相比，移相相關(guān)法估計相位差的誤差在不同相位差真值處的估計誤差基本相同，且誤差大大減小了；在0和180°理論相位差附近，不再有很大的估計偏差；在計算量上只是增加了3次移相和一周期采樣點參與運算的相位差初步估計，因此在精度大大提高的同時，計算時間并沒有明顯增加，當(dāng)N0=480×103時，在普通計算機上的計算時間約為8ms，滿足實時測距的要求。

2 仿真實驗

以Monte Carlo方法生成信噪比R=30dB、采樣長度480×103的兩個正弦信號x（k）和y（k），對于不同的實際相位差φ，按（3）式表述的普通相關(guān)法仿真計算1000次相位差，仿真計算的均方誤差隨實際相位差φ變化的曲線見圖4，帶圈的曲線為按（6）式計算的理論值，帶點的曲線為仿真實驗值，理論值和仿真實驗值一致。圖中均方差的最大仿真值為8.1dB，對應(yīng)的相位差誤差為2.54°。

圖5為移相相關(guān)法估計的相位差估計誤差，仿真條件和圖4一樣，并增加頻率誤差為10×10-6這一條件。其中圖5a為理論最大均方誤差和仿真的均方誤差隨著實際相位差φ變化的曲線，仿真的均方誤差為帶點的曲線，而帶圈的曲線為當(dāng)數(shù)字頻率誤差δf=10×10-6，N0=480×103，n=48時，按（17）式計算的理論最大均方誤差，該圖中最大均方誤差的仿真值為-48dB。而圖5b為仿真的平均誤差和標(biāo)準(zhǔn)差，平均誤差反映的是（16）式所示的頻率誤差引起的相位差誤差。

3 結(jié) 論

移相相關(guān)法計算相位差適用于計算包含噪聲的、頻率相同且已知的兩正弦信號的相位差。該方法不但具有普通相關(guān)法計算速度快和能計算中、低頻信號相位差的優(yōu)點，而且進(jìn)一步抑制了噪聲的影響。在頻率誤差不大于10×10-6、信噪比R= 30dB、采樣長度480×103時，相位差誤差為0.004°。若信號頻率穩(wěn)定度高于20×10-6，且每周期采樣點數(shù)大于48時，頻率誤差引起的移相誤差不超過0.002°，可以滿足0.01°的相位差估計精度。當(dāng)信號頻率準(zhǔn)確度很低且要求很高的相位差估計精度時，可以在計算相位差以前先校正頻率。通過增加采樣長度、提高采樣速率、減小頻率誤差或者多次測量求平均值等辦法，可以進(jìn)一步減小相位差誤差。

［1］ SCH?NER H，BAUER F，DORRINGTON A，et al.Image processing for three-dimensional scans generated by time-of-flight range cameras［J］.Journal of Electronic Imaging，2012，21（2）：023012.

［2］ CHOU Y Sh，LIU J S.A robotic indoor 3-D mapping system using a 2-D laser range finder mounted on a rotating four-bar linkage of a mobile platform［J］.International Journal of Advanced Robotic Systems，2013，10（45）：1-10.

［3］ MAYKOLPINTO N A，ROCHA L F，PAULOMOREIRA A.Object recognition using laser range finder and machine learning techniques［J］.Robotics and Computer-Integrated Manufacturing，2013，29（1）：12-22.

［4］ BAUD C，TAP-B′ETEILLE H，LESCURE M，et al.Analog and digital implementation of an accurate phasemeter for laser range finding［J］.Sensors and Actuators，2006，A132（11）：258-264.

［5］ GAO Y Y，LI Y H，F(xiàn)ENG Q L，et al.Optical design of a laser distance measuring system based on high frequency digital phase detection［J］.Laser Technology，2013，37（3）：353-356（in Chinese）.

［6］ YOON H S，SONG H J，PARK K W.A phase-shift laser scanner based on a time-counting method for high linearity performance［J］.Review of Scientific Instruments，2011，82（7）：075108.

［7］ BERTOTTI F L，HARA M S，ABATTI P J.A simple method to measure phase difference between sinusoidal signals［J］.Review of Scientific Instruments，2010，81（11）：115106.

［8］ ACCATTATIS A，SAGGIO G，GIANNINI F.A real time FFT-based impedance meter with bias compensation［J］.Measurement，2011，44（4）：702-707.

［9］ MARTINO M，LOSITO R，MASI A.Analytical metrological characterization of the three-parameter sine fit algorithm［J］.ISA Transactions，2012，51（2）：262-270.

［10］ RAMOS P M，SERRA A C.A new sine-fitting algorithm for accurate amplitude and phase measurements in two channel acquisition systems［J］.Measurement，2008，41（2）：135-143.

［11］ LIANG Y R，DUAN H Z，YEH H C.Fundamental limits on the digital phase measurement method based on cross-correlation analysis［J］.Review of Scientific Instruments，2012，83（9）：095110.

［12］ WALKER W F，TRAHEY G E.A fundamental limit on delay estimation using partially correlated speckle signals［J］.IEEE Transactions on Ultrasonics，F(xiàn)erroelectrics，and Frequency Control，1995，42（2）：301-308.

［13］ VIOLA F，WALKER W F.A comparison of the performance of time-delay estimators in medical ultrasound［J］.IEEE Transactions on Ultrasonics，F(xiàn)erroelectrics，and Frequency Control，2003，50（4）：392-401.

［14］ LUO J W，KONOFAGOU E E.A fast normalized cross-correlation calculation method for motion estimation［J］.IEEE Transactions on Ultrasonics，F(xiàn)erroelectrics，and Frequency Control，2010，57（6）：1347-1357.

［15］ ZHANG S L，ZHANG K，DUAN H Z.Approximate computation of expectation and variance of nonlinear function of continuous random variable［J］.Journal of Geodesy and Geodynamics，2008，28（4）：107-109（in Chinese）.

Study on phase difference algorithm based on phase-shift correlation analysis

LIU Yuzhou，ZHAO Bin
（Department of Instrumentation，School of Mechanical Science and Engineering，Huazhong University of Science and Technology，Wuhan 430074，China）

In order to improve the accuracy of a phase-shift range finder，a phase-difference algorithm based on phase-shift correlation analysis was proposed to estimate the phase-difference between two sinusoidal signals with same frequency.For reducing the influence of noise，the autocorrelation between the original and 2π shifted signal was calculated firstly.Secondly，the phase difference was estimated approximately with a few sampled data and the initial phase of one signal was shifted by Δθ to make the phase difference between two signals to be near π/2（or 3π/2）.Then，the phasedifference was calculated with whole set of data by correlation method and the final phase difference was obtained by subtracting Δθ.The influence of frequency error was analyzed.Theoretical analysis and simulation shows that the error of this method is greatly reduced.The proposed method can improve the accuracy of a range finder.

measurement and metrology；phase-shift correlation analysis；phase difference；frequency error

TH741

A

10.7510/jgjs.issn.1001-3806.2014.05.013

1001-3806（2014）05-0638-05

國家九七三重點基礎(chǔ)研究發(fā)展計劃資助項目（2013CB035405）

劉玉周（1968-），男，博士研究生，主要研究方向為光電精密測量。

*通訊聯(lián)系人。E-mail：zhaobin63@sohu.com

2013-11-29；

2014-01-13