基于matlab的fft算法實(shí)現(xiàn)(論文)-資料下載頁
2025-06-27 18:15本頁面
【正文】 參考文獻(xiàn)[1] 張志勇，[s].北京航空航天大學(xué)出版社，2006，08[2] 鄭君里 ,[s].高等教育出版社, 2008, 12 [3] ：清華大學(xué)出版社出版，2001[4] 高西全 丁玉美等 數(shù)字信號處理. 北京：電子工業(yè)出版社，2009[5] 余成波，陶紅艷。數(shù)字信號處理及MATLAB 實(shí)現(xiàn)，北京：清華大學(xué)出版社，2008 附錄 源程序清單附錄Clear all。[x,fs,bits]=wavread(‘’,2048)。 %讀取聲音x1=reshape(x,1,4096)。sound(x1,fs,bits)。 %播放語音信號y1=fft(x1)N=length(x1)。 %計(jì)數(shù)讀取信號的點(diǎn)數(shù)t=(1:N)/fs。 %信號的時域采樣點(diǎn)figure(1)plot(t, x1)。 %畫出聲音采樣后的時域波形title(39。原聲音信號的時域波形39。)。 %給圖形加注標(biāo)簽說明xlabel(39。時間/t39。)。ylabel(39。振幅/A39。)。grid 。 %添加網(wǎng)格M=nextpow2(x1)。 % 求x的長度對應(yīng)的2的最低冪次m N=2^M。if length(x1)N x1=[x1,zeros(1,Nlength(x1))]。 % 若x的長度不是2的冪，補(bǔ)零到2的整數(shù)冪 end NV2=N/2。NM1=N1。I=0。J=0。while INM1 if IJ T=x1(J+1)。 x1(J+1)=x1(I+1)。 x1(I+1)=T。 end K=NV2。 while K=JJ=JK。 K=K/2。 end J=J+K。 I=I+1。end %x1。 y=x1。 % 將x倒序排列作為y的初始值 WN=exp(i*2*pi/N)。 %蝶形運(yùn)算 for L=1:M B=2^L/2。 %第L級中，每個蝶形的兩個輸入數(shù)據(jù)相距B個點(diǎn)，每級有B個不同的旋轉(zhuǎn)因子 for J=0:B1 % J代表了不同的旋轉(zhuǎn)因子 p=J*2^(ML)。 WNp=WN^p。 for k=J+1:2^L:N % 本次蝶形運(yùn)算的跨越間隔為2^L kp=k+B。 % 蝶形運(yùn)算的兩個因子對應(yīng)單元下標(biāo)的關(guān)系 t=y(kp)*WNp。 % 蝶形運(yùn)算的乘積項(xiàng) y(kp)=y(k)t。 % 蝶形運(yùn)算， 注意必須先進(jìn)行減法運(yùn)算，然后進(jìn)行加法運(yùn)算，否則要使用中間變量來傳遞y(k) y(k)=y(k)+t。 % 蝶形運(yùn)算 end end end %yfigure(2)[x2,w1]=freqz(x1,1) 。 %繪制原始語音信號的頻率圖plot(w1/pi,20*log10(abs(x2)))。title(39。頻率特性圖39。)xlabel(39。歸一化頻率39。)。ylabel(39。幅度/DB39。)。grid。 figure(3)subplot(2,1,1)。 %把畫圖區(qū)域劃分為2行1列，指定第1 個圖plot(abs(y1))。 %繪制原始語音信號的幅頻響應(yīng)圖title(39。直接運(yùn)算 語音信號FFT頻譜特性39。)。 %給圖形加注標(biāo)簽說明xlabel(39。n39。)。ylabel(39。幅值/A39。)。grid。 %添加網(wǎng)格subplot(2,1,2)。 %把畫圖區(qū)域劃分為2行1列，指定第2個圖plot(abs(y))。title(39。蝶形運(yùn)算 語音信號FFT頻譜特性39。)。 xlabel(39。n39。)。ylabel(39。幅值/A39。)。grid。 %添加網(wǎng)格figure(4)subplot(2,1,1)if y~=0 %判斷指數(shù)是否為0 plot(20*log10(abs(y1)))。 %畫信號頻譜的分貝圖endxlabel(39。n39。)。ylabel(39。振幅/db39。)。title(39。直接運(yùn)算 語音信號FFT頻譜特性39。)。grid 。 %添加網(wǎng)格第 24 頁 共 2
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1