湘潭大學物理與光電工程學院 湖南省微光電與系統(tǒng)集成工程實驗室 李 凡 金湘亮
MEMS數(shù)字檢波器用高精度抽取濾波器的設計與實現(xiàn)
湘潭大學物理與光電工程學院 湖南省微光電與系統(tǒng)集成工程實驗室 李 凡 金湘亮
針對類似石油勘探應用中采用一位數(shù)據(jù)輸出的高精度MEMS數(shù)字檢波器,設計數(shù)字抽取濾波器對輸出數(shù)據(jù)做進一步處理。首先基于Matlab的simulink建模平臺搭建完整的MEMS數(shù)字檢波器系統(tǒng)模型,使用Filterbuilder工具箱進行數(shù)字抽取濾波器設計與建模,結果表明經(jīng)過數(shù)字抽取濾波器濾波后,信噪比僅下降1dB。其次使用Verilog硬件語言實現(xiàn)該濾波器,使用雙時鐘實現(xiàn)四通道數(shù)據(jù)分配,降低資源消耗。通過Modelsim與Quarts2進行聯(lián)合仿真,結果顯示通帶波紋小于0.02dB,阻帶衰減小于-130dB,該濾波器效果達到設計要求。
MEMS;數(shù)字抽取濾波器;Verilog
MEMS數(shù)字檢波器廣泛應用于航空航天,石油探測,生物學等各種領域。在諸如石油探測的高精度應用中,MEMS數(shù)字檢波器中的加速度傳感器大多使用sigma-delta調(diào)制器來獲得較好的性能[1][7];文獻[2][3]研究了MEMS加速度閉環(huán)系統(tǒng)仿真參數(shù);文獻[4]中提到搭配一位量化器來獲得最高的線性度,以解決量化臺階之間的非線性失配問題,并構成一個穩(wěn)定的回路。但一位輸出數(shù)據(jù)無法直觀地還原原始信號波形,而過采樣ADC的噪聲整形特性使量化噪聲被排至高頻率區(qū)間。
針對這個問題,本文以簡便的濾波器設計工具設計數(shù)字抽取濾波器,建模仿真并用Verilog硬件語言實現(xiàn),用以濾除高頻段的量化噪聲,同時將一比特數(shù)據(jù)流還原為原始波形,將輸出速率降至奈奎斯特頻率。
傳統(tǒng)MEMS數(shù)字檢波器系統(tǒng)結構如圖1所示,輸入加速度信號通過一個敏感MEMS器件引起電容變化,通過一個C/V讀出電路將電容變化轉為電信號,繼而通過過采樣ADC進行采樣轉為數(shù)字信號,一位量化器將輸出表達為1bit數(shù)據(jù)流,這種方法解決了量化臺階之間的非線性失配問題。1bit數(shù)據(jù)流輸出以靜電力反饋形式反饋至敏感MEMS器件以維持器件的平衡,構成伺服系統(tǒng)。數(shù)字抽取濾波器用以處理量化器后輸出的1bit數(shù)據(jù),輸出一個原始信號波形,并濾除高頻段量化噪聲[12]。
圖1 系統(tǒng)結構圖
MEMS數(shù)字檢波器的過采樣ADC基帶帶寬為500Hz,采樣率為128,故輸出信號速率為128kHz。
針對輸出1bit數(shù)據(jù)流速率為128kHz的高精度MEMS數(shù)字檢波器,對其進行128倍的抽取,輸出速率為1kHz的信號。數(shù)字抽取濾波器的設計指標如表1所示。
表1 數(shù)字抽取濾波器的設計指標
圖2 降采樣示意圖
傳統(tǒng)數(shù)字抽取濾波器設計采用手工計算FIR濾波器,過程繁瑣[5][6]?,F(xiàn)今各種濾波器設計軟件已經(jīng)成熟,本文利用Matlab中Filterbuilder工具箱對以上濾波器進行設計并導出simulink模型,輸入具體設計參數(shù)如表2所示。
作為整個數(shù)字抽取濾波器的第一級,CIC濾波器的主要作用是實現(xiàn)較大的抽取數(shù),為后級濾波器節(jié)省資源,同時CIC濾波器可以增大阻帶衰減,接近設計指標[9]。在Filterbuilder工具箱輸入各項指標后可生成CIC幅頻響應圖。
圖3 CIC濾波器幅頻圖
如圖3所示,依據(jù)我們的設計指標生成了一個6階16抽取的CIC濾波器。
CIC濾波器在實現(xiàn)較大幅度的阻帶衰減的同時,通帶中也會有相應較大衰減,這時需要一個補償濾波器來補償其通帶中的衰減,這里需要注意的是,補償濾波器的阻帶起始頻率如果越接近最終半帶輸出的阻帶起始頻率,濾波效果會更好但是濾波器階數(shù)會變得非常大,這樣做會消耗大量資源,這里經(jīng)過權衡后將阻帶起始頻率設置在1500Hz,生成43階的補償濾波器,使用fvtool可視化工具可以查看補償效果如圖4所示。
圖4 補償濾波器補償效果圖
可以看到經(jīng)過補償后,通帶截止頻率點400Hz處的幅度為-0.0034dB,滿足設計指標。
圖5 濾波器系統(tǒng)整體級聯(lián)圖
使用半帶濾波器作為整個濾波器最后一級,與一般的FIR濾波器相比,半帶濾波器的通帶特性比較平坦,過渡帶的帶寬也比較容易控制,并且濾波器的系數(shù)有一半為零,因此減少了設計的硬件開支,減少了濾波器整體的面積和功耗。
通過Fvtool工具將濾波器整體級聯(lián)圖作出如圖5,可以看到通帶內(nèi)波紋遠小于0.02dB,阻帶衰減也達到小于-130dB要求。
在simulink平臺中搭建MEMS數(shù)字檢波器的模型如圖6所示。輸入加速度經(jīng)由MEMS加速度傳感器前端電路,輸出一位數(shù)據(jù)流,一位數(shù)據(jù)流經(jīng)一個反饋系數(shù)轉為靜電力反饋至MEMS數(shù)字檢波器的敏感器件,用以將MEMS敏感器件回歸原位,構成伺服系統(tǒng)[4]。
圖6 MEMS數(shù)字檢波器系統(tǒng)simulink模型
輸入信號為一個低頻的幅度為1g的加速度信號與一個幅度為0.4g的高頻噪聲信號相疊加,圖7分別為輸入信號,輸出1bit數(shù)據(jù)與最終輸出信號的示波器示意圖。
圖7 示波器示意圖
可以看到,輸入信號通過MEMS數(shù)字檢波器前段電路后輸出極高速率的1bit數(shù)據(jù)流,表現(xiàn)為示波器中密集的1與-1,再通過所設計生成的數(shù)字抽取濾波器組進行信號處理,輸出的最終信號完美地保留并還原了原始的加速度信號波形,濾除了高頻噪聲信號。對輸出1bit數(shù)據(jù)流與最終輸出數(shù)據(jù)進行頻譜分析如圖8,9所示[11]。
從圖8中可以看到輸出的1bit數(shù)據(jù)流包含了頻率固定的低頻的加速度信號,高頻噪聲信號以及高頻段大量的量化噪聲信號,帶寬內(nèi)信噪比為135.7dB。圖9中可以看到高頻段的噪聲信號與量化噪聲信號都已完全被濾除,只保留了帶寬內(nèi)的加速度信號,同時帶寬內(nèi)信噪比只下降了約1個dB。
圖8 1bit輸出數(shù)據(jù)頻譜分析圖
圖9 濾波器輸出數(shù)據(jù)頻譜分析圖
圖10 系統(tǒng)框架圖
在Matlab中對濾波器進行建模仿真后,確定simulink模型的正確性后,可以進行編寫Verilog代碼進行濾波器的硬件實現(xiàn)。Filterbuilder支持直接生成Verilog代碼,但是自動生成的Verilog代碼冗長繁雜[11]。在simulink模型中可以查看底層電路和濾波器系數(shù),依據(jù)底層電路與濾波器系數(shù)可以編寫Verilog代碼,整個硬件系統(tǒng)框架圖如圖10所示。
系統(tǒng)采用計數(shù)器的方式進行時鐘分頻,分別進行16倍,64倍,128倍分頻,并向補償濾波器產(chǎn)生一個降采樣用分頻時鐘,這個降采樣分頻時鐘頻率與64倍分頻后的時鐘頻率一樣,相位相差半個周期,目的是以兩個時鐘的上升沿與下降沿實現(xiàn)補償濾波器中四通道數(shù)據(jù)分配,以滿足補償濾波器降4倍降采樣的需要,實現(xiàn)波形如圖11所示。
圖11 雙時鐘實現(xiàn)四通道數(shù)據(jù)分配
在實現(xiàn)濾波器中的系數(shù)時,采用CSD編碼可以大大減少資源消耗[6]。其原理是將系數(shù)轉化為由0,1和-1組成的編碼,進行優(yōu)化后產(chǎn)生最少非零元素的編碼,保證用最少的加減運算和移位運算實現(xiàn)系數(shù)進行乘法運算后的結果,在公式(1)中通過一個例子演示了CSD編碼的過程。
可以看出原本需要7次移位和加法運算在經(jīng)過CSD編碼后只需4次移位和加減法運算便可實現(xiàn)原系數(shù)乘法運算的結果。
圖12 聯(lián)合仿真波形
在通過Modelsim軟件功能仿真后,在Quaters2軟件中配置Altera公司的cycloneIII-EP3C120F484C8型號FPGA板,編寫Testbeach測試文件,通過編譯并進行仿真[7]。將Matlab工作區(qū)的MEMS數(shù)字檢波器量化器后1bit數(shù)據(jù)輸出導出,利用Testbeach輸入至Quaters2軟件中,與Modelsim進行聯(lián)合仿真得到下圖波形結果。
圖12中由上至下依次為輸入1bit數(shù)據(jù),CIC濾波器輸出數(shù)據(jù)波形,補償濾波器輸出數(shù)據(jù)波形與半帶濾波器即最終輸出數(shù)據(jù)波形??梢钥吹?,數(shù)字抽取濾波器成功實現(xiàn)了抽取與濾波的功能。
本文針對一位輸出的高精度MEMS數(shù)字檢波器利用Matlab設計了一套完整的數(shù)字抽取濾波器系統(tǒng),在simulink平臺中建模并仿真,結果顯示濾波后的信噪比只下降了1dB,表明濾波器性能達到期望指標。其次使用Verilog代碼實現(xiàn)整個數(shù)字抽取濾波器系統(tǒng),通過優(yōu)化代碼減少資源消耗。最后通過Quarts2軟件配置FPGA開發(fā)板與Modelsim軟件聯(lián)合仿真,結果表明該濾波器成功實現(xiàn)了抽取與濾波的功能,滿足設計要求,有較好的實際應用價值。
[1]Liang W,Yang C,Qiao D.On the design of feed-forward interface for MEMS based accelerometers[J].Analog Integrated Circuits and Signal Processing,2012,72(1):3-10.
[2]Petkov V P,Boser B E.High-order electromechanical sigma-delta modulation in micro-machined inertial sensors[J].Circuits and Systems I:IEEE Transactions on,2006,53(5):1016-1022.
[3]ZHANG F,JIN X,WANG L.Modeling and test result of closedloop MEMS accelerometer with wide dynamic range[J].Microsystem Techn ologies,2016,22(3):653-657.
[4]Wang L,Xie L,Jin X.A design of order-adjustable sigma-delta modulator for high performance micro-machined accelerometers[J].Analog Integrated Circuits and Signal Processing,2014,80(1):77-84.
[5]LIU X.A High Speed Digital Decimation Filter with Parallel Cascaded Integrator-Comb Pre-Filters[C]//International Congress on Image and Signal Processing.IEEE,2009:1-4.
[6]LI L,XIAO yue L I,ZENG J,et al.Design of 24-bit Delta-Sigma A/D digital decimation filter[J].Journal of Beijing Jiaotong University,2016.
[7]李宗敏,張法全,王國富,等.基于FPGA高精度微震信號采集系統(tǒng)設計[J].儀表技術與傳感器,2016(4):108-110.
[8]李秀花,馬雯.基于FPGA的IIR數(shù)字濾波器設計與實現(xiàn)[J].電子技術,2016(5).
[9]楊靜.寬帶∑-△A/D轉換器中數(shù)字抽取濾波器的設計與驗證[J].電子設計工程,2013,21(22):168-170.
[10]劉忠超,張長春,李衛(wèi),等.0.18μm CMOS Σ-ΔADC用數(shù)字抽取濾波器設計[J].微電子學與計算機,2014(6):44-47.
[11]楊翠娥.CIC抽取濾波器的MATLAB設計及FPGA實現(xiàn)[J].山西電子技術,2017(1):46-48.
[12]駱麗,李曉玥,曾俊琦,等.一種24位Delta-SigmaA/D數(shù)字抽取濾波器設計[J].北京交通大學學報,2016,40(5):45-49.
李凡(1993—),碩士研究生,現(xiàn)就讀于湘潭大學,主要研究方向為信號處理與集成電路設計。
金湘亮【通訊作者】(1974—),教授,主要研究方向為傳感與信息獲取、集成電路與系統(tǒng)。
Design and realization of high precision decimation filter for MEMS digital geophone
Li Fan,Jin Xiang-liang
(Hunan Engineering Laboratory for Microelectronics ,Optoelectronics and System on a chip,F(xiàn)aculty of Physics and Optoelectronics Xiangtan University,Xiangtan 411105 China)
In view of the application of oil exploration,a high-precision MEMS digital geophone with one bit output data is used,the digital decimation filter is designed to further process the output data.A MEMS digital geophone system model is complete based on Simulink platform of Matlab,the Filterbuilder toolbox is used to design and modeling the digital decimation filter,the results show that after the digital decimation filter,the SNR decreased by 1dB .Then use the Verilog HDL to realize the filter,using dual clock to achieve four channel data allocation,reduce resource consumption.Through the joint simulation simulation of Modelsim and Quarts2,the simulation result shows that the passband ripple is less than 0.02dB,and the stopband attenuation is smaller than -130dB,the filter effect meets the design requirements.
MEMS;Digital decimation filter;Verilog
國家自然科學基金重點項目(61233010);國家自然科學基金面上項目(61774129);國家自然科學基金青年科學基金項目(61704145);湖南省自然科學杰出青年基金項目(2015JJ1014)。