基于傅里葉變換的MEMS地震檢波器設(shè)計
Vm表示輸入電壓信號,Vs表示輸出電壓,Cs1與Cs2分別表示固定臂與可動臂之間的兩個電容,外界加速度與輸出電壓的關(guān)系為:
可見,在加速度計的結(jié)構(gòu)和輸入電壓確定的情況下,輸出電壓與加速度呈正比關(guān)系。
A/D轉(zhuǎn)換
MSP430系列單片機具有處理能力強、運行速度快、資源豐富等優(yōu)點,有很高的性價比,內(nèi)部自帶12位A/D,可以選擇多個通道的模擬輸入,轉(zhuǎn)換內(nèi)核由一個采樣保持器和一個轉(zhuǎn)換器組成,對高速變化的信號進行瞬時采樣時,一旦ADC 開始轉(zhuǎn)換,采樣保持器則進行保持,即使現(xiàn)場輸入的信號的變化比較快,也不會影響到ADC的轉(zhuǎn)換工作。
采樣信號高的時候采樣,低的時候轉(zhuǎn)換,自動將轉(zhuǎn)換的結(jié)果保存到相應(yīng)的存儲器里。ADC12一共有12個轉(zhuǎn)換通道,有16個轉(zhuǎn)換存儲器,存儲的數(shù)據(jù)再經(jīng)過FFT變換,得到相應(yīng)的幅頻特性。
FFT變換
經(jīng)過A/D轉(zhuǎn)換后的數(shù)據(jù),利用傅里葉變換可以把信號從時域轉(zhuǎn)換到頻域,進行頻域分析,可以看到各個頻率下的信號信息,有利于對地震波進行更準(zhǔn)確的分析。
離散傅里葉變換分析如下:一個周期為的函數(shù)可用傅里葉基數(shù)展開為
其中:
將連續(xù)函數(shù)傅里葉基數(shù)展開式(2)離散化,為了離散化式(2),在周期區(qū)間(0,2, 上等間隔的取個點,取樣間隔為,那么,這里要注意。則的散化序列為,,由此式(2)離散化形式為:
且對分子分母同乘以后變?yōu)?img onload="if(this.width>620)this.width=620;" onclick="window.open(this.src)" style="cursor:pointer" height="32" src="http://editerupload.eepw.com.cn/fetch/20130823/162558_2_17.jpg" width="75" alt="" />,由此可得出xk第項為一個正弦和一個余弦周期函數(shù)和,其頻率為:,其中T為所取序列總的時間長度。隨著k的增大,三角函數(shù)的頻率逐漸增加,周期逐漸增加,周期逐漸減小,其周期為:。當(dāng)時,諧波的頻率最大為:該頻率稱為Nyquist頻率,當(dāng)k從取到N時,其結(jié)果與k從0取到是鏡像對稱時,現(xiàn)在將式(3)的各次諧波寫成如下形式:
其中:,為k次諧波振幅:為k次諧波的初相。
通過以上算法進行編程,實現(xiàn)對采集信號進行FFT變換,得到其頻譜特性圖。
軟件設(shè)計
軟件部分主要由主程序、ADC中斷子程序和FFT程序組成,主程序完成系統(tǒng)初始化以及各軟件模塊的初始化,ADC中斷子程序完成通道的選擇和采樣率的設(shè)置,F(xiàn)FT程序完成采集數(shù)據(jù)經(jīng)A/D轉(zhuǎn)換后的相應(yīng)頻率點的幅值運算,系統(tǒng)流程圖如4所示。 加速度計相關(guān)文章:加速度計原理
評論