FFT棰戣氨娴嬭瘯娉昅atlab浠ｇ爜 - 鐧惧害鏂囧簱

鏂板缓

涓婁紶

棣栭〉

鍔╂墜

鏈�杩?/div>

璧勬枡搴?/div>

宸ュ叿

FFT

棰戣氨娴嬭瘯娉?/p>

Matlab

浠ｇ爜濡備笅锛?/p>

 

adc_data=[];%

閲囨牱寰楀嚭鐨勬暟鎹�褰㈡垚鐨勭煩闃�

 

data_num=size(adc_data); 

numpt=data_num(1,1);%

閲囨牱鐐圭殑鎬绘暟

 

no=1;vin=2;d9=3;d8=4;d7=5;d6=6;d5=7;d4=8;d3=9;d2=10;d1=11;d0=12;%

鏍?/p>

绀?/p>

鏍?/p>

鏈�鐭╅樀鐨勫�?/p>

 

vref=2; 

fclk=1e6

;%A/D

杞�鎹㈠櫒鐨勯噰鏍烽�戠巼

 

adc_bit=10; 

for i=1:1:numpt 

dout(i)=(adc_data(i,d9)/2+adc_data(i,d8)/4+adc_data(i,d7)/8+adc_data(i,d6)/16+adc_

data(i,d5)/32+adc_data(i,d4)/64+adc_data(i,d3)/128+adc_data(i,d2)/256+adc_data(i,d

1)/512+adc_data(i,d0)/1024)*vref;%

閫氳繃杞�鎹㈠櫒鍑烘潵鐨勭粨鏋滐紝鎭㈠�嶅師鏉ョ殑娉㈠舰

 

 

 

 

 

doute(i)=dout(i)-1;%

灏嗘�ｅ鸡娉㈢殑鍏辨ā鐢靛钩鍋忕疆涓�

0V 

end 

 

figure(3) 

plot(doute); 

doutw=doute.*hanning(numpt)';%

鍔犳眽瀹佺獥

 

dout_spect=fft(doutw);%

杩涜��

FFT

鍙樻崲锛?/p>

鐒跺悗鍒╃敤涓嬮潰鐨勭畻娉曟眰鍑?/p>

SNR

銆?/p>

SINAD

銆?/p>

SFDR 

dout_dB=20*log10(abs(dout_spect)); 

maxdB=max(dout_dB(1:numpt/2)); 

x=(

[0:numpt/2-1]

.*

fclk/numpt

)/1000;

 

 

 

 

 

% 

fclk/numpt

涓洪�戣氨鍒嗚鲸鐜囷�?/p>

鍥犱负

FFT

棰戣氨鍥惧叧浜庝腑蹇冨�圭О鎵�浠ュ彧瑕佸垎鏋愪竴鍗婄殑棰戣氨鍥惧嵆鍙?/p>

 

y1=dout_dB(1:numpt/2)-maxdB; 

figure(1) 

plot(x,y1); 

axis([0,200,-120,5]); 

grid off; 

title('1024-Samples FFT spectrum'); 

xlabel('Frequency (kHz)'); 

ylabel('Magnitude (dB)');%

鐢诲嚭

FFT

鐨勫箙棰戠壒鎬ф洸绾?/p>