電子信息工程課程設(shè)計(jì)數(shù)據(jù)采集系統(tǒng)動(dòng)態(tài)有效位測(cè)試

1、電子信息工程課程設(shè)計(jì)實(shí)驗(yàn)報(bào)告數(shù)據(jù)采集系統(tǒng)動(dòng)態(tài)有效位測(cè)試院系：電子信息與光電技術(shù)學(xué)院 專業(yè)：電子信息工程姓名： 指導(dǎo)老師：劉明2013年9月7日一、課程設(shè)計(jì)目的將模擬信號(hào)轉(zhuǎn)變?yōu)閿?shù)字信號(hào)后再進(jìn)行處理，是當(dāng)前信號(hào)處理普遍使用的方法，模數(shù)轉(zhuǎn)換器（ADC）就是將模擬信號(hào)轉(zhuǎn)換為數(shù)字信號(hào)的器件。ADC有8位，10位，12位，16位等，理論上，轉(zhuǎn)換位數(shù)越多，轉(zhuǎn)換精度越高，但由于受到電路噪聲干擾，以及采樣頻率等因素的影響，實(shí)際上所獲得的轉(zhuǎn)換精度（轉(zhuǎn)換位數(shù)）要低于理論值。所以，計(jì)算ADC的有效轉(zhuǎn)換位數(shù)對(duì)系統(tǒng)性能的評(píng)估就顯得尤為重要。本課題就是利用DSP實(shí)驗(yàn)裝置對(duì)模擬信號(hào)進(jìn)行數(shù)字化，并存儲(chǔ)到計(jì)算機(jī)中，利用matla

2、b編程對(duì)信號(hào)進(jìn)行處理，計(jì)算ADC的有效轉(zhuǎn)換位數(shù)（ENOB）,通過調(diào)整信號(hào)幅度、信號(hào)頻率、采樣頻率等參數(shù)，獲得這些參數(shù)與ENOB的關(guān)系曲線，最終計(jì)算出最佳的ENOB值。二、課程設(shè)計(jì)內(nèi)容1、接通信號(hào)發(fā)生器和示波器的電源，調(diào)整信號(hào)發(fā)生器的頻率（小于5KHz），幅度（峰峰值小于1V），利用示波器測(cè)量確認(rèn)后，通過連接電纜將信號(hào)發(fā)生器的輸出連接到DSP實(shí)驗(yàn)箱的INPUT1端口。2、用USB連接電纜連接主機(jī)和DSP實(shí)驗(yàn)箱3、接通DSP實(shí)驗(yàn)箱電源，根據(jù)液晶顯示屏顯示的提示信息進(jìn)行操作。1）上電后，首先選擇 4 （AD），按ENTER鍵確認(rèn)2）通過數(shù)字鍵選擇采樣頻率（符合那奎斯特采樣定理），按ENTER鍵確認(rèn)3

3、）選擇“1”保存，通過主機(jī)上的采集軟件，可將采集的數(shù)據(jù)通過USB線上傳到主機(jī)。選則“2”不保存，可通過DSP試驗(yàn)箱的OUT3接口，通過示波器觀察波形，若系統(tǒng)正常，應(yīng)該能夠看到跟信號(hào)發(fā)生器輸出一致的波形，以此來驗(yàn)證電路系統(tǒng)的正確性。4）若在3）選擇“1”保存后，主機(jī)會(huì)提示安裝USB驅(qū)動(dòng)，正確安裝驅(qū)動(dòng)后，打開主機(jī)上的數(shù)據(jù)采集軟件，如圖所示：若USB驅(qū)動(dòng)安裝正確，會(huì)出現(xiàn) EZ-USB-SX2 提示信息。5）點(diǎn)擊“start”,開始數(shù)據(jù)傳輸，若系統(tǒng)工作正常，Successed Transfers 后會(huì)顯示“5”，表明收到5個(gè)數(shù)據(jù)包，若顯示信息不是5，則將DSP試驗(yàn)箱斷電，重新開始。6）若5）正常，則主

4、機(jī)會(huì)產(chǎn)生一個(gè)數(shù)據(jù)文件 USB.DAT，這就是ADC采集的數(shù)據(jù)，共1024個(gè)采樣點(diǎn)，每個(gè)采樣點(diǎn)為12位有效數(shù)字，表示為2個(gè)字節(jié)，高8位在前（其中高4位為0），低8位在后。4、利用matlab編寫程序，利用公式求出ENOB。三、實(shí)驗(yàn)要求 1）用matlab編寫程序，繪制時(shí)域波形圖、頻譜圖，計(jì)算ENOB2）通過改變信號(hào)頻率、幅度以及采樣頻率，獲取這些參數(shù)與ENOB的對(duì)應(yīng)關(guān)系，繪制出對(duì)應(yīng)曲線圖。 3）設(shè)置最佳參數(shù)，求出最佳的ENOB，要求計(jì)算出的ENOB10 4）通過查找資料自己設(shè)計(jì)最佳的matlab程序算法 5）參考方法：對(duì)采樣數(shù)據(jù)進(jìn)行FFT，求出基波能量和諧波能量，計(jì)算出信噪比（SNR），利用公式

5、求解ENOB 參考公式： ENOB = （SNR-1.79）/6.02 提示：注意采樣頻率、信號(hào)頻率、FFT點(diǎn)數(shù)之間的最佳對(duì)應(yīng)關(guān)系。四、 實(shí)驗(yàn)原理1、 動(dòng)態(tài)有效位ENOB：用來衡量數(shù)據(jù)采集系統(tǒng)實(shí)際工作時(shí)有效的位數(shù)，它是用分辨率來衡量實(shí)際工作時(shí)ADC的噪聲均方值與理想ADC標(biāo)識(shí)分辨率情況下的量化噪聲。ENOB是量化描述模數(shù)轉(zhuǎn)換性能的一種方法，有效位數(shù)越大，表示模數(shù)轉(zhuǎn)換中所記錄的電壓值越精確。2、 測(cè)試動(dòng)態(tài)有效位的原因：對(duì)于理想的模數(shù)轉(zhuǎn)換器，它具有完美的線性特性，僅對(duì)輸入的信號(hào)進(jìn)行簡(jiǎn)單的量化處理。量化過程會(huì)產(chǎn)生量化噪聲，根據(jù)信號(hào)功率與噪聲功率，可以計(jì)算出模數(shù)轉(zhuǎn)換之后的信噪比SNR。若輸入信號(hào)為純正

6、弦波，那么SNR就可以表示為：SNR=1.522B，B為數(shù)模轉(zhuǎn)換器的位數(shù)，將上式表示為db,則B=（SNRdb-1.76）/6.02,其中1.76為理想ADC的量化噪聲，6.02為將log2轉(zhuǎn)化為log10的系數(shù)比。但是現(xiàn)實(shí)生活中理想的模數(shù)轉(zhuǎn)換器并不存在，每種模數(shù)轉(zhuǎn)換器都會(huì)給輸入信號(hào)帶來某種失真。典型的失真有噪聲，非線性輸入特性以及增益與失調(diào)誤差。這些失真會(huì)直接惡化信噪比，所以12bit的模數(shù)轉(zhuǎn)換器的ENOB可能僅為10.5，這意味著，雖然模數(shù)轉(zhuǎn)換器的輸出為12bit,但其達(dá)到的信噪比僅相當(dāng)于理想的10.5bit模數(shù)轉(zhuǎn)換器。所以，測(cè)試模數(shù)轉(zhuǎn)換性能，即動(dòng)態(tài)有效位十分必要。3、 原理公式：五、所需

7、儀器設(shè)備信號(hào)發(fā)生器、示波器、DSP實(shí)驗(yàn)箱、計(jì)算機(jī)、連接電纜6、 實(shí)驗(yàn)結(jié)果與分析采樣頻率fs=8kHz 信號(hào)頻率f=78.125Hz 信號(hào)幅度v=1vFFT點(diǎn)數(shù)N=1024輸出結(jié)果：SNR=38.2214 ENOB=6.0567分析：實(shí)驗(yàn)時(shí)，要注意“頻譜泄漏這個(gè)概念.對(duì)于頻率為fs的正弦序列，它的頻譜應(yīng)該只是在fs處有離散譜。但是，在利用DFT求它的頻譜做了截短，結(jié)果使信號(hào)的頻譜不只是在fs處有離散譜，而是在以fs為中心的頻帶范圍內(nèi)都有譜線出現(xiàn)，它們可以理解為是從fs頻率上“泄露”出去的，這種現(xiàn)象稱 為頻譜“泄露”。周期信號(hào)加窗后做DFT仍然有可能引起頻譜泄露，設(shè)fs為采樣頻率，N為采樣序列長(zhǎng)度

8、，分析頻率為：f0=m*fs/N（m=0，1.)，以cos函數(shù)為例，設(shè)其頻率為f0,如果 f0不=m*fs/N，就會(huì)引起除f0以外的其他m*fs/N點(diǎn)為非零值，即出現(xiàn)了泄露。在理解了頻譜泄漏的概念后，為了達(dá)到ENOB6的實(shí)驗(yàn)要求，且采樣頻率fs10kHz的實(shí)驗(yàn)要求,我們選取了m=100，fs=8kHz，將N=1024代入，得到的正弦信號(hào)頻率為f0=78.125Hz七、實(shí)驗(yàn)截圖ENOB最大時(shí)所對(duì)應(yīng)的時(shí)域波形圖ENOB最大時(shí)所對(duì)應(yīng)的頻譜圖附錄：Matlab程序：clc;clear;close all；%load data.mat data%fid=fopen(E:project214SOCdata

9、f2812sine-1.001khz-4096-73.6040hz.dat,r);% data,N = fscanf(fid,%x);% fclose(fid);FileName,PathName = uigetfile(*.dat,Select the *.dat file);f = fullfile(PathName,filesep,FileName);fid = fopen(f,r);data = fscanf(fid,%x);% data = fread(fid);fclose(fid);data = data(1:2:end)*256 + data(2:2:end);N = 1024;%datsgn = sub_unsign2sign(data,16)/24;datsgn = data;plot(datsgn); %時(shí)域波形datfft = fft(datsgn(1:N);%*4/N/212); %fftnordat = abs(datfft)/max(abs(datfft);figure,stem(abs(nordat); %譜線datfft(1) = 0;maxval,pos = max(abs(datfft(1:N);datfft_hf = nordat(1:N/2);sig_rms=(norm(datfft_hf(pos-5:pos+5); %/sqrt(