基于高調(diào)制頻率的快速數(shù)字鎖相算法

蔣國亮，楊 暉，李 然

( 上海理工大學(xué)光學(xué)與電子信息工程學(xué)院,上海 200093)

?

基于高調(diào)制頻率的快速數(shù)字鎖相算法

蔣國亮，楊 暉，李 然

( 上海理工大學(xué)光學(xué)與電子信息工程學(xué)院,上海 200093)

為了提高數(shù)字鎖相的檢測速度以及其最高調(diào)制頻率，改善數(shù)字鎖相放大的實時檢測性能，文中提出了一種新型快速數(shù)字鎖相算法。文中根據(jù)方波參考信號傅里葉展開存在奇次諧波，使得調(diào)制頻率等于參考信號頻率的2m-1倍。同時，通過將采樣頻率設(shè)為調(diào)制頻率的4的整數(shù)倍，保證兩路參考信號的正交性。在相敏檢波算法實現(xiàn)中，僅僅存在加減法運算，消除了乘法運算，提高了數(shù)字鎖相算法的運算速度。而且，基于普通微處理器的數(shù)字鎖相放大能達到的最高調(diào)制頻率也得到了有效的提升，改善了鎖相放大的實時性能。實驗表明，該算法既具備較高的速度，又能提高鎖相放大的實時性，使得數(shù)字鎖相對微弱信號的實時檢測在普通微處理器上的實現(xiàn)成為可能。

數(shù)字鎖相算法；調(diào)制頻率；微弱信號；相關(guān)運算；奇次諧波；方波參考

0 引言

微弱信號檢測技術(shù)廣泛地應(yīng)用在光譜學(xué)、光通訊、天文、聲納、雷達以及生物醫(yī)學(xué)工程，用以檢測這些領(lǐng)域中被噪聲背景所淹沒的微弱信號。而基于相關(guān)檢測原理的鎖相放大器(lock-in amplifier，LIA)這一微弱信號檢測儀器是目前準確度較高，應(yīng)用最廣泛的裝置[1-4]。早期的鎖相放大器都是由模擬電路來實現(xiàn)的，但傳統(tǒng)的模擬LIA存在溫度漂移和其他非線性的低頻特性的影響。隨著數(shù)字技術(shù)的不斷發(fā)展，數(shù)字鎖定放大器(DLIA)逐步取代了模擬LIA。DLIA采用了數(shù)字乘法器和數(shù)字濾波器在動態(tài)范圍、噪聲方面遠優(yōu)于模擬器件，大大提高了鎖相放大測量的準確度。

國內(nèi)外眾多科研人員將數(shù)字鎖相檢測應(yīng)用到各自領(lǐng)域，并對提高數(shù)字鎖相檢測的性能做了大量的研究[5-10]。天津大學(xué)的李剛教授，提出了使用過采樣技術(shù)簡化外圍電路和減少量化噪聲[8]。同時，為了提高DLIA的算法速度，把采樣頻率設(shè)為調(diào)制頻率的4的整數(shù)倍，基本去除了鎖相放大算法中的乘法運算并減少了大量的減法運算，降低了運算量和存儲量，極大地提高了數(shù)字鎖相算法的速度[10]。Masciotti等人通過使用矩陣乘法，極大地減少了多調(diào)制頻率數(shù)字鎖相算法的運算量[10]。上述研究極大地提高了數(shù)字鎖相算法的速度，然而，在一些對調(diào)制頻率有很高要求的在線實時系統(tǒng)上，傳統(tǒng)的數(shù)字鎖相算法無法同時滿足系統(tǒng)的高調(diào)制頻率和實時性。

本文針對上述問題，提出了一種基于方波參考信號的快速數(shù)字鎖相算法。首先，把調(diào)制頻率設(shè)為方波參考頻率的2m-1倍，然后再將輸入信號與參考信號相乘。由于方波的傅里葉級數(shù)包含有奇次諧波，上述信號相乘的結(jié)果一定會產(chǎn)生一個含有輸入信號幅值和相位信息的直流分量。在相同的時鐘消耗下，該算法大大提高了數(shù)字鎖相放大的最高調(diào)制頻率，也提高了數(shù)字鎖相系統(tǒng)檢測的實時性。此外，該數(shù)字鎖相算法只有加減法運算，幾乎消除了所有的乘法運算，進一步提高了運算速度并克服了算法實現(xiàn)對處理器的依賴性。

1 數(shù)字鎖相原理

相關(guān)檢測是應(yīng)用信號周期性和噪聲隨機性的特點，通過自相關(guān)或互相關(guān)運算，達到去除噪聲、檢測信號的一種技術(shù)。鎖相放大器就是采用互相關(guān)技術(shù)，通過對待測信號與參考信號的互相關(guān)運算，實現(xiàn)對待測信號中有效信號的提取，達到提高信噪比的目的。傳統(tǒng)的正交數(shù)字鎖相放大器采用了雙相關(guān)器的原理，參考信號由兩路相互正交的方波信號構(gòu)成，其原理框圖如圖1所示。

圖1 正交鎖相放大器原理框圖

假設(shè)輸入信號為

I(t)=Asin(2πf0t+φ)

式中：A為幅值；φ為相位。

以采樣頻率fs=Nfm對輸入信號進行采樣，N≥2，得信號序列為

I[n]=Asin(2πn/N+φ),n=0,1,2，……，M-1

(1)

式中M為總采樣點。

參考信號為與輸入信號頻率相同的方波，且參考信號不用進行采樣，可由微處理器根據(jù)N來產(chǎn)生正交參考序列S[n]和C[n]。

(2)

(3)

輸入信號分別與正交參考序列相乘實現(xiàn)相敏檢波的功能，得到的相關(guān)信號中的直流分量僅與原始信號的幅值和初相位有關(guān)，因此通過數(shù)字低通濾波器取出直流分量即可得到輸入信號的幅值和相位信息。通常采用的低通濾波器為平均濾波器[11]，信號的正交相關(guān)運算和濾波過程如下：

(4)

(5)

式中M為整周期采樣點數(shù)，同時也對應(yīng)著低通濾波器的時間常數(shù)。

輸出的2個分量一個是與cos(φ)成正比的同相分量，一個是與sin(φ)成正比的正交分量，將2個信號經(jīng)如下運算即可得到有用信號的幅度和相位：

(6)

(7)

2 高調(diào)制頻率快速鎖相算法

2.1 算法原理

傳統(tǒng)的鎖相放大器要求調(diào)制頻率與參考頻率一致才能正常工作，而本文提出的鎖相放大器的調(diào)制頻率fm為方波參考頻率fr的2m-1倍。由于方波的傅里葉展開式中包含有奇次諧波的高頻成分，因此當(dāng)調(diào)制信號頻率為參考信號的奇次倍時，其相關(guān)運算結(jié)果仍然能產(chǎn)生一個直流分量僅與原始信號的幅值和初相位有關(guān)。此時，輸入信號的正交相關(guān)運算和濾波過程如下：

(8)

(9)

(10)

(11)

2.2 算法實現(xiàn)

高調(diào)制頻率快速鎖相放大器中，調(diào)制頻率是參考信號的2m-1倍，即采樣1個完整的計算周期M個采樣點，應(yīng)有2m-1個周期輸入信號序列而參考信號僅為1個周期序列。因此，在一個完整的計算周期內(nèi)，由處理器產(chǎn)生的正交參考序列S[n]和C[n]分別為

(12)

(13)

此時，相關(guān)運算的結(jié)果為：

(14)

(15)

圖2(a) 是該算法的相敏檢波原理圖，圖中調(diào)制頻率為參考信號頻率的5倍，圖2(b)為傳統(tǒng)鎖相放大相敏檢波原理圖。圖2(a)和圖2(b)的采樣頻率都設(shè)為參考頻率的20倍，且參考信號與輸入信號相位差為0。

從式(14)和式(15)可以看出，在正交互相關(guān)計算中的乘法運算全部消除，僅僅通過加減法運算就能夠?qū)崿F(xiàn)互相關(guān)運算，大大降低了計算量。對于相同采樣頻率的傳統(tǒng)數(shù)字鎖相算法，若采樣M個點，正交互相關(guān)運算中乘法運算次數(shù)為 2M，加法運算次數(shù)為2(M-1)；而高調(diào)制頻率快速鎖相算法中乘法運算次數(shù)為0，加減法次數(shù)為2(M-1)。顯然，隨著采樣周期的增加，高調(diào)制頻率快速鎖相算法的優(yōu)越性能越能體現(xiàn)。同時，由于參考信號的頻率低于調(diào)制頻率，減少了序列產(chǎn)生占用的時鐘周期，進一步減少了對處理器的時鐘消耗。

(a)高調(diào)制頻率數(shù)字鎖相m=3

(b)傳統(tǒng)數(shù)字鎖相圖2 相敏檢波原理圖

3 實驗和仿真

為了驗證這種高調(diào)制頻率快速數(shù)字鎖相算法的有效性，利用LABVIEW進行仿真，通過改變幅值和調(diào)制頻率與參考信號頻率的倍數(shù)，比較真實值和計算出的幅值和相位。

驗證計算幅值的有效性：通過仿真，產(chǎn)生一系列的頻率為方波參考信號奇次倍的正弦調(diào)制信號，初相位和直流分量為0，采樣頻率設(shè)置為調(diào)制頻率的20倍。改變調(diào)制信號的幅值，并加入白噪聲作為輸入信號，采用該方法計算的幅值如表1所示(保留小數(shù)點后5位)。

表1 不同調(diào)制頻率下幅值測試結(jié)果V

實際幅值計算幅值K=1K=2K=3K=40.001000.001000.001000.001000.001000.010000.010060.010060.010060.010060.100000.100430.100430.100430.100430.500000.502070.502070.502070.502071.000001.004181.004181.004181.00418

驗證計算相位的有效性: 通過仿真，產(chǎn)生一系列的頻率為方波參考信號奇次倍的正弦調(diào)制信號，其中幅值為1 V，直流分量為0，采樣頻率設(shè)置為調(diào)制頻率的20倍。改變調(diào)制信號的初相位，并加入白噪聲作為輸入信號，采用該方法計算的相位如表2所示(保留小數(shù)點后5位)。

表2 不同調(diào)制頻率下相位測試結(jié)果rad

實際相位計算出的相位K=1K=2K=3K=40.349070.349030.349030.349030.349030.523600.523650.523650.523650.523650.785400.785480.785480.785480.785481.047201.047261.047261.047261.047261.570801.570761.570761.570761.57076

在表1和表2中，K=1的情況下可認為是傳統(tǒng)鎖相放大的測量值。由表1可以看出，采用這種高調(diào)制頻率的快速鎖相放大算法測得的幅值與傳統(tǒng)鎖相放大算法相一致不會引入誤差，且測得的準確率大于99.5%。而，由表2可看出，這種高調(diào)制頻率的鎖相放大算法測得的相位與傳統(tǒng)的鎖相放大算法完全一致。

4 結(jié)束語

本文研究了數(shù)字鎖相放大器對微弱信號檢測的相關(guān)理論，推導(dǎo)了參考信號為方波時的相關(guān)運算結(jié)果，利用了方波傅里葉級數(shù)展開包含的奇次諧波高頻成分，提出了使用調(diào)制頻率為參考信號頻率的奇次倍來實現(xiàn)相敏檢波運算，而非傳統(tǒng)的調(diào)制頻率必須等于參考頻率。本文提出的高調(diào)制頻率快速數(shù)字鎖相放大算法，與傳統(tǒng)的數(shù)字鎖相放大相比，去除了幾乎所有有的乘法運算，同時通過使用相對較低的參考頻率減小了處理器的負荷。實驗仿真結(jié)果表明，這種新型的數(shù)字鎖相放大算法并沒有任何的理論誤差。同時，該算法還極大地提高了運算速度，提高了實時鎖相放大檢測的最高調(diào)制頻率，使得基于數(shù)字鎖相放大的微弱信號處理可以在普通的微處理器上實現(xiàn)。

[1] ANDERSSON M，PERSSON L，SVENSSON T,et al.Flexible lock-in detection system based on synchronized computer plug-in boards applied in sensitive gas spectroscopy.Rev.Sci.Instrum,2007,78(11): 113107-113113.

[2] RIPKA P，DRAXLER K，STYBLIKOVA R.Measurement of DC Currents in the Power Grid by Current Transformer.IEEE Trans.Magn,2013,49(1): 73-76.

[3] BENGTSSON L E.A microcontroller-based lock-in amplifier for sub-milliohm resistance measurements.Rev.Sci.Instrum,2012,83(7): 075103-075110.

[4] HSDPSR J，CHEN S F，GORDILLO A,et al.Digital lock in amplifier: study,design and development with a digital signal processor.Microprocessors and microsystems,2004,28(4): 157-162.

[5] BARRAGAN L A,ARTIGAS J I,ALONSO.A modular,low-cost,digital signal processor-based lock-in card for measuring optical attenuation.Rev.Sci.Instrum,2001,72(1): 247-251.

[6] SONNAILLON M,BONETTO F J.Lock-in amplifier error prediction and correction in frequency sweep measurements.Rev.Sci.Instrum,2007,78(1): 014701.

[7] 李剛，張麗君，林凌． 一種新型數(shù)字鎖相放大器的設(shè)計及其優(yōu)化算法．天津大學(xué)學(xué)報,2008,41( 4): 429-432．

[8] SONNAILLON M O,BONETTO F J.A low-cost,high-performance,digital signal processor-based lock-in amplifier capable of measuring multiple frequency sweeps simultaneously.Rev.Sci.Instrum,2005,76(2):024703.

[9] 李剛，周梅，何峰，等．基于數(shù)字鎖相相關(guān)計算結(jié)構(gòu)的優(yōu)化算法．電子與信息學(xué)報,2012,34(3): 744-748．

[10] MASCIOTTI J M，LASKER J M，HIELSCHER A H.Digital Lock-In Detection for Discriminating Multiple Modulation Frequencies With High Accuracy and Computational Efficiency.IEEE Trans.Instrum.Meas,2008,57(1): 182-189.

High Speed Digital Lock-in Algorithm Based on High-modulate Frequency

JIANG Guo-liang，YANG Hui，LI Ran

(School of Optical-Electronical and Computer Engineering,University of Shanghai for Science & Technology,Shanghai 200093,China)

In order to improve the speed and the maximum modulation frequency of the digital lock-in,a high speed digital phase-locked algorithm was presented in this paper.In this paper,a new algorithm that make the modulate frequency equal to 2m-1 of square reference frequency was proposed to implement a lock-in,which was based on odd harmonic components in the square wave reference.The sampling frequency was four times as large as the modulate frequency to insure the reference sequences orthogonal.In the implementation of phase sensitive detection,only additions and subtractions were used,which speed up the computation in lock-in.Furthermore,the maximum modulate frequency of a lock-in was enhanced considerably with a general microprocessor.The results of the experiments show that the new algorithm not only has high speed,but also improves the performance of real-time detection,which makes it possible for a general microprocessor to implement weak signal detection by the new algorithm in real time.

digital lock-in algorithm;modulation frequency;weak signal;correlation operation;the odd harmonic;square wave reference

國家自然科學(xué)基金(61007002)；上海市教育委員會科研創(chuàng)新項目(15ZZ072)

2015-03-16 收修改稿日期：2015-07-02

TP206

A

1002-1841(2015)12-0126-03

蔣國亮(1992—)，碩士研究生，研究方向為信號與信息處理。E-mail:jdangg192@163.com