黃利紅 郭唐永 鄒 彤 吳 濤

(1.中國地震局地震研究所(地震大地測量重點實驗室)，武漢 430071；2.中國地震局地殼應力研究所武漢科技創(chuàng)新基地，武漢 430071)

?

差動式電容傳感器數(shù)字解調(diào)技術(shù)*

黃利紅1,2郭唐永1,2鄒 彤1,2吳 濤1,2

(1.中國地震局地震研究所(地震大地測量重點實驗室)，武漢 430071；2.中國地震局地殼應力研究所武漢科技創(chuàng)新基地，武漢 430071)

在大地測量中，我們對如何檢測出地殼微小的變化很感興趣，利用差動式電容傳感器可以檢測出這種微弱的信號。本文對差動式電容傳感器的輸出信號進行數(shù)字化解調(diào)，介紹了三種解調(diào)方法，分別是數(shù)字包絡解調(diào)、數(shù)字正交相干解調(diào)及基于Goertzel算法的解調(diào)。文中給出了各自的原理，并用Matlab對給出的信號進行了仿真比較，結(jié)果分析。

差動式電容傳感器；數(shù)字解調(diào)；數(shù)字包絡解調(diào)；數(shù)字正交相干解調(diào)；基于Goertzel算法的解調(diào)

0 引言

地殼傾斜的變化極其微小，要對其進行觀測，必須有能測量出微小位移的器件，一種方法是將這種微小的變化放大后再測量，另一種方法是通過電容測微傳感器檢測，如利用差動式電容傳感器。差動電容傳感器的輸出電壓很小，易被噪聲淹沒，因而，必須經(jīng)過高增益放大才能檢測出所需信號。經(jīng)高增益放大后，還需要經(jīng)過鎖相放大器濾波。

鎖相放大器抗噪聲性能非常好,將交流信號放大且變成相應的直流信號,即檢測出并放大相應的同步信號。這種模擬解調(diào)的方法能滿足儀器的各項指標，但數(shù)字信號相比于模擬信號有易于存儲，便于處理及可靠性高等優(yōu)點?？紤]將數(shù)字技術(shù)應用到儀器上面，相對于模擬電路，能提高儀器的可靠性，減少調(diào)試工作量。而儀器的數(shù)字化改造首先要將模擬解調(diào)數(shù)字化，即數(shù)字化解調(diào)方法的實現(xiàn)。

1 數(shù)字化解調(diào)方法

從差動式電容傳感器的檢測原理知，可將已調(diào)制信號設為

y[n]=Ax[n]cos(w0n+θ)

(1)

1.1 數(shù)字包絡解調(diào)

調(diào)制信號的信息包含在載波的包絡中，而載波的各個幅值點組成了載波的信號包絡。通過某種算法來得到幅值點以獲得信號包絡從而實現(xiàn)解調(diào)。

由式(1)知

在一個載波周期內(nèi)有fs/f0個采樣點，將每個載波周期中的最大幅值點取出，并在周期末輸出。當載波頻率非常大以及采樣頻率過高時，可多個載波周期輸出一個幅值。數(shù)字包絡解調(diào)流程圖如圖1所示。在圖1中，描述的是某個載波周期處理的過程，第1步先輸入初始值到存儲大的幅值中儲存；第2步將輸入的值放入比較中；第3步將存儲大的幅值與比較的值進行比較，將大的賦值給存儲大的幅值，重復2、3、4步，一直到得到這個載波周期的最大值，最后一步是第5步，將這個最大的值給輸出。

新載波周期圖1 數(shù)字包絡解調(diào)流程圖

1.2 數(shù)字正交解調(diào)

數(shù)字正交解調(diào)框圖如圖2所示。

圖2 數(shù)字正交解調(diào)框圖

ri[n] =y[n]cos(w0n)

=Ax[n]cos(w0n+θ)cos(w0n)

rq[n] =y[n]sin(w0n)

=Ax[n]cos(w0n+θ)sin(w0n)

經(jīng)低通濾波器濾波后

1.3 基于Goertzel算法解調(diào)

m[n] =y[n]w[n]

做傅里葉變換得

因而，在解調(diào)過程中僅需求載波頻率處的DFT。

期望輸出X(k)=yk(n)n=N，其中k=0,1，2，…，N-1。

式中，WN=e-j2p /N

相應的系統(tǒng)函數(shù)為

(2)

為了避免式(2)所蘊含的復數(shù)乘法和加法，可以合并擁有復共軛極點的諧振器對。這時，系統(tǒng)函數(shù)變?yōu)?/p>

(3)

差分方程描述為

(4)

(5)

2cos(2pk/N)vk(N)vk(N-1)

(6)

初始條件為vk(-1)=vk(-2)=0。

2 對比分析

設差動式電容傳感器載波頻率為100kHz，在Matlab仿真中，選取調(diào)制信號頻率f為0～20Hz，載波頻率f0=100kHz的信號進行解調(diào)，采樣頻率fs=800kHz，分別加入信噪比SNR為0～30dB的高斯白噪聲。選取f=5Hz，SNR=10dB的仿真圖進行比較，其中包絡解調(diào)以10個載波周期輸出一個幅值。程序是在Matlab2010上運行的，選取的是8個調(diào)制信號周期長度的數(shù)據(jù)量。

圖3(a)是待解調(diào)信號的波形圖及頻譜圖。圖3(b)(c)(d)中的子圖a1、a2、a3分別是用三種方法解調(diào)出的信號的波形圖，可以看出，它們都成功實現(xiàn)了解調(diào)的功能；子圖b1、b2、b3是加入噪聲的波形圖；子圖c1、c2、c3是有噪聲解調(diào)出的信號與無噪聲解調(diào)出的信號之差。子圖c1中信號的大小幾乎是子圖b1中信號大小的1/4，子圖c2中信號的大小也幾乎是子圖b2中信號大小的1/4，子圖c3中信號的大小與子圖b3幾乎一樣。因此，在抗干擾方面，基于Goertzel算法的解調(diào)與正交解調(diào)都有一定的抗干擾能力，包絡解調(diào)則沒有。

圖3 解調(diào)效果圖

基于Goertzel算法的解調(diào)，由式(3)它的系統(tǒng)函數(shù)中可以看出，它的分母部分是一個二階系統(tǒng)，相當于一個線性的二階濾波器，因而有一定的濾波功能。從圖2正交解調(diào)框圖中，可以看到，在正交解調(diào)過程中，兩路信號分別都需經(jīng)過一個低通濾波器，因而解調(diào)出來的信號也有一定的濾波效果。在包絡解調(diào)中，由于包絡解調(diào)只是簡單的將每個載波周期的最大值取出，因而無法對信號濾波，在信噪比低時，干擾嚴重，解調(diào)效果不好。

在基于Goertzel算法的解調(diào)中，由于它的系統(tǒng)函數(shù)式(3)合并了擁有復共軛極點的諧振對，因而在它的遞歸關(guān)系式(4)中，系數(shù)均為實數(shù)，將復數(shù)乘法簡化成了實數(shù)乘法，節(jié)約了運算量。正交解調(diào)由于兩路信號先要經(jīng)過乘法器，再經(jīng)過低通濾波器濾波，最后還要實現(xiàn)平方和再求根這個過程，實現(xiàn)過程較復雜，耗用資源。

從圖3子圖a1、a2、a3及圖3(e)中，還可以看出，基于Goertzel算法的解調(diào)還有一定的放大功能，它將信號放大了4倍。

表1～表3列出了f=1、5、10、20Hz，SNR=1、5、10、15、20、25、30dB的條件下，三種方法輸出的幅度最大點,其中包絡解調(diào)以10個載波周期輸出一個幅值的方式。程序是在Matlab2010上運行的，選取的是8個調(diào)制信號周期長度的數(shù)據(jù)量。

表1 包絡解調(diào)輸出幅值最大點

表2 正交解調(diào)輸出幅值最大點

表3 基于Goertzel算法解調(diào)輸出幅值最大點

從3個表格可以看出，在同頻率的情況下，基于Goertzel算法解調(diào)輸出值幅度最大點是正交解調(diào)輸出幅度最大點的4倍，包絡解調(diào)由于選取的是8個點取1個值，計算頻譜值時點數(shù)是正交解調(diào)的1/8，因此，它的幅度最大點是正交解調(diào)的1/8。

3 結(jié)論

數(shù)字相干解調(diào)實現(xiàn)過程復雜，非常耗用資源；數(shù)字包絡解調(diào)實現(xiàn)原理簡單，但在信噪比低的情況下，解調(diào)性能非常差，抗干擾能力差而無法將被噪聲淹沒的信號取出，而基于Goretzel算法的解調(diào)抗干擾能力好，運算過程為實數(shù)運算因而運算速度較快并且易于實現(xiàn)。因此，在儀器的數(shù)字化改造中，與模擬解調(diào)對應的數(shù)字解調(diào)方法可以采用基于Goertzel算法解調(diào)進行嘗試。

[1] Joyce Van de Vegte.數(shù)字信號處理基礎[M].北京：電子工業(yè)出版社，2003

[2] SanjitK.Mitra.數(shù)字信號處理-基于計算機的方法(第四版)[M].北京：電子工業(yè)出版社，2012

[3] 李瑩.靜電反饋技術(shù)在傾斜儀中的應用研究[D].武漢：中國地震局地震研究所，2009

[4] 王慶河,王慶山.數(shù)據(jù)處理中的幾種常用數(shù)字濾波算法[J].計量技術(shù),2003(04):53-54

[5] 張肅文.高頻電子線路(第四版)[M].北京：高等教育出版社，2004

[6] 梁志國.周期信號的諧波分析述評[J].計量技術(shù),2003(02)

[7] 胡國新.新型傾斜儀機械結(jié)構(gòu)研究及誤差分析[D].武漢：中國地震局地震研究所，2012

[8] 陳懷琛.數(shù)字信號處理教程-MATLAB釋義與實現(xiàn)(第3版)[M].北京：電子工業(yè)出版社，2013

[9] 胡廣書.數(shù)字信號處理理論算法與實現(xiàn)(第三版)[M].北京：清華大學出版社，2012

[10] 江志紅.深入淺出數(shù)字信號處理[M].北京：北京航空航天大學出版社，2012

湖北省科技支撐計劃項目(2014BEC090)

10.3969/j.issn.1000-0771.2015.08.01