摘要：  本音頻信號分析儀由32位MCU為主控制器，通過AD轉(zhuǎn)換，對音頻信號進(jìn)行采樣，把連續(xù)信號離散化，然后通過FFT快速傅氏變換運(yùn)算，在時(shí)域和頻域?qū)σ纛l信號各個(gè)頻率分量以及功率等指標(biāo)進(jìn)行分析和處理，然后通過高分辨率的LCD對信號的頻譜進(jìn)行顯示。該系統(tǒng)能夠精確測量的音頻信號頻率范圍為20Hz-10KHz，其幅度范圍為5mVpp-5Vpp，分辨力分為20Hz和100Hz兩檔。測量功率精確度高達(dá)1%,并且能夠準(zhǔn)確的測量周期信號的周期，是理想的音頻信號分析儀的解決方案。

方案一、用DDS芯片配合FIFO對信號進(jìn)行采集，通過DDS集成芯片產(chǎn)生一個(gè)頻率穩(wěn)定度和精度相當(dāng)高的信號作為FIFO的時(shí)鐘，然后由FIFO對A/D轉(zhuǎn)換的結(jié)果進(jìn)行采集和存儲(chǔ)，最后送MCU處理。

方案二、直接由32位MCU的定時(shí)中斷進(jìn)行信號的采集，然后對信號分析。

由于32位MCU -LPC2148是60M的單指令周期處理器，所以其定時(shí)精確度為16.7ns,已經(jīng)遠(yuǎn)遠(yuǎn)可以實(shí)現(xiàn)我們的40.96KHz的采樣率，而且控制方便成本便宜，所以我們選擇由MCU直接采樣。

由于快速傅立葉變換FFT算法設(shè)計(jì)大量的浮點(diǎn)運(yùn)算，由于一個(gè)浮點(diǎn)占用四個(gè)字節(jié)，所以要占用大量的內(nèi)存，同時(shí)浮點(diǎn)運(yùn)算時(shí)間很慢，所以采用普通的8位MCU一般難以在一定的時(shí)間內(nèi)完成運(yùn)算，所以綜合內(nèi)存的大小以及運(yùn)算速度，我們采用Philips 的32位的單片機(jī)LPC2148，它擁有32K的RAM，并且時(shí)鐘頻率高達(dá)60M，所以對于浮點(diǎn)運(yùn)算不論是在速度上還是在內(nèi)存上都能夠很快的處理。

對于普通的音頻信號，頻率分量一般較多，它不具有周期性。測量周期可以在時(shí)域測量也可以在頻域測量，但是由于頻域測量周期性要求某些頻率點(diǎn)具有由規(guī)律的零點(diǎn)或接近零點(diǎn)出現(xiàn)，所以對于較為復(fù)雜的，頻率分量較多且功率分布較均勻且低信號就無法正確的分析其周期性。

而在時(shí)域分析信號，我們可以先對信號進(jìn)行處理，然后假定具有周期性，然后測出頻率，把采樣的信號進(jìn)行周期均值法和定點(diǎn)分析法的分析后即可以判別出其周期性。

綜上，我們選擇信號在時(shí)域進(jìn)行周期性分析和周期性測量。對于一般的音頻信號，其時(shí)域變化是不規(guī)則的，所以沒有周期性。而對于單頻信號或者由多個(gè)具有最小公倍數(shù)的頻率組合的多頻信號具有周期性。這樣我們可以在頻域?qū)π盘柕念l譜進(jìn)行定量分析，從而得出其周期性。而我們通過先假設(shè)信號是周期的，然后算出頻率值，然后在用此頻率對信號進(jìn)行采樣，采取連續(xù)兩個(gè)周期的信號，對其值進(jìn)行逐次比較和平均比較，若相差太遠(yuǎn)，則認(rèn)為不是周期信號，若相差不遠(yuǎn)（約5%），則可以認(rèn)為是周期信號。

音頻信號經(jīng)過一個(gè)由運(yùn)放和電阻組成的50 Ohm阻抗匹配網(wǎng)絡(luò)后，經(jīng)由量程控制模塊進(jìn)行處理，若是一般的100mV-5V的電壓，我們選擇直通，也就是說信號沒有衰減或者放大，但是若信號太小，12位的A/D轉(zhuǎn)換器在2.5V參考電壓的條件下的最小分辨力為1mV左右，所以如果選擇直通的話其離散化處理的誤差將會(huì)很大，所以若是采集到信號后發(fā)現(xiàn)其值太小，在20mV-250mV之間的話，我們可以將其認(rèn)定為小信號，從而選擇信號經(jīng)過20倍增益的放大器后再進(jìn)行A/D采樣。

file:///C:/Users/chang/AppData/Local/Temp/msohtmlclip1/01/clip_image002.gif經(jīng)過12位A/D轉(zhuǎn)換器ADS7819轉(zhuǎn)換后的數(shù)字信號經(jīng)由32位MCU進(jìn)行FFT變換和處理，分析其頻譜特性和各個(gè)頻率點(diǎn)的功率值，然后將這些值送由Atmega16進(jìn)行顯示。信號由32 位MCU分析后判斷其周期性，然后由Atmegal6進(jìn)行測量，然后進(jìn)行顯示。

信號輸入后通過R5,R6兩個(gè)100Ohm的電阻和一個(gè)高精度儀表運(yùn)放AD620實(shí)現(xiàn)跟隨作用，由于理想運(yùn)放的輸入阻抗為無窮大，所以輸入阻抗即為：R5//R6=50Ohm,阻抗匹配后的通過繼電器控制是對信號直接送給AD轉(zhuǎn)換還是放大20倍后再進(jìn)行AD轉(zhuǎn)換。

在這道題目里，需要檢測各頻率分量及其功率，并且要測量正弦信號的失真度，這就要求在對小信號進(jìn)行放大時(shí)，要盡可能少的引入信號的放大失真。正弦信號的理論計(jì)算失真度為零，對引入的信號失真非常靈敏，所以對信號的放大，運(yùn)放的選擇是個(gè)重點(diǎn)。

我們選擇的運(yùn)放是TI公司的低噪聲、低失真的儀表放大器INA217，其失真度在頻率為1KHz，增益為20dB(100倍放大)時(shí)僅為0.004%，其內(nèi)部原理圖如圖2-2所示。

file:///C:/Users/chang/AppData/Local/Temp/msohtmlclip1/01/clip_image006.jpg

其中放大器A1的輸出電壓計(jì)算公式為

OUT1=1+（R1/RG）*VIN+

R3、R4、R5、R6及A3構(gòu)成減法器，最后得到輸出公式

VOUT=（VIN2-VIN1）*[1+（R1+R2）/RG]

采用12位AD轉(zhuǎn)換器ADS7819進(jìn)行轉(zhuǎn)換，將轉(zhuǎn)換的數(shù)據(jù)送32位控制器進(jìn)行處理。

功率譜測量主要通過對音頻信號進(jìn)行離散化處理，通過FFT運(yùn)算，求出信號各個(gè)離散頻率點(diǎn)的功率值，然后得到離散化的功率譜。

由于題目要求頻率分辨力為100Hz和20Hz兩個(gè)檔，這說明在進(jìn)行FFT運(yùn)算前必須通過調(diào)整采樣頻率(fK)和采樣的點(diǎn)數(shù)(N),使其基波頻率f為100Hz和20Hz。

根據(jù)頻率分辨率與采樣頻率和采樣點(diǎn)數(shù)的關(guān)系：

f=fk/N;

又根據(jù)采樣定理，采樣頻率fk必須不小于信號頻率fm的2倍，即：

                   fk>=2fm;

題目要求的最大頻率為10KHz，所以采樣頻率必須大于20KHz，考慮到FFT運(yùn)算在2的次數(shù)的點(diǎn)數(shù)時(shí)的效率較高，所以我們在20Hz檔時(shí)選擇40.96KHz采樣率，采集2048個(gè)點(diǎn)，而在100檔時(shí)我們選擇51.2KHz采樣率，采集512個(gè)點(diǎn)。

通過FFT 分析出不同的頻率點(diǎn)對應(yīng)的功率后，就可以畫出其功率譜，并可以在頻域計(jì)算其總功率。

主控制芯片為LPC2148,測量周期為Atmega16實(shí)現(xiàn),由于處理器速度較快,所以采用c語言編程方便簡單.軟件流程圖如下：

file:///C:/Users/chang/AppData/Local/Temp/msohtmlclip1/01/clip_image010.jpg

            主流程圖                           周期性分析和測量流程圖

表4-1 總功率測量結(jié)果（室溫條件下）

結(jié)果分析： 由于實(shí)驗(yàn)室提供的能夠模仿音頻信號的且能方便測量的信號只有正弦信號，所以我們用一款比較差點(diǎn)的信號發(fā)生器產(chǎn)生信號，然后進(jìn)行測量，發(fā)現(xiàn)誤差不達(dá)，在+-5%以內(nèi)。我們以音頻信號進(jìn)行測量，由于其實(shí)際值無法測量，所以我們只能根據(jù)時(shí)域和頻域以及估計(jì)其誤差，都在5%以內(nèi)。

表4-2 單個(gè)頻率分量測量結(jié)果（室溫條件下）

結(jié)果分析：我們首先以理論上單一頻率的正弦波為輸入信號，在理想狀況下，其頻譜只在正弦波頻率上有值，而由于有干擾，所以在其他頻點(diǎn)也有很小的功率。

音頻信號由于有多個(gè)頻點(diǎn)，所以沒有一定的規(guī)律性。由于音頻信號波動(dòng)較大，沒有一定的規(guī)律，且實(shí)驗(yàn)室沒有專門配置測量儀器，所以我們只好以正弦波和三角波作為信號進(jìn)行定量分析測量，以及對音頻信號進(jìn)行定性的分析和測量。我們發(fā)現(xiàn)其數(shù)字和用電腦模擬的結(jié)果符合得很近。

由于系統(tǒng)架構(gòu)設(shè)計(jì)合理，功能電路實(shí)現(xiàn)較好，系統(tǒng)性能優(yōu)良、穩(wěn)定，較好地達(dá)到了題目要求的各項(xiàng)指標(biāo)。

[6] 　馬忠梅等． ARM&Linux嵌入式系統(tǒng)教程．北京：北京航空航天大學(xué)出版社，2004年

1、 LPC2148

2、 ATMEGA16

3、 AD620

4、 ADS7819

5、 液晶320*240

1、  低頻信號發(fā)生器      

2、  數(shù)字萬用表            

3、  失真度測量儀

4、  數(shù)字示波器

5、  穩(wěn)壓電源

                     前級放大和AD轉(zhuǎn)換

Atmega16控制板

/*/////////////////////////////////////////////////////////////////////////////////////////////////

      FFT轉(zhuǎn)換函數(shù)，dataR:實(shí)部，datai:虛部，

////////////////////////////////////////////////////////////////////////////////////////////////*/

voidFFT(float *dataR,float *dataI,int n)

{

inti,L,j,k,b,p,xx,qq;

intx[11]={0};

floatTR,TI,temp;

float QQ;

//////////////////////////////////位倒置////////////////////////////////////////////////////

for(i=0;i<count[n];i++)

{ xx=0;

  for(j=0;j<n;j++)

  x[j]=0;

  

  for(j=0;j<n;j++)

  {x[j]=(i/count[j])&0x01;}

for(j=0;j<n;j++)

{xx=xx+x[j]*count[n-j-1];}

dataI[xx]=dataR;