【正文】 SP 芯片組成的基本系統(tǒng)的相關電路有了比較深的認識。 在此次課程中主要的要求是用麥克風采集一段語音信號，在 MATLAB 軟件平臺上，利用編程加入噪聲函數(shù)， 繪制 出前后的 波形并觀察其頻譜，給定相應技術指標，用 漢寧窗 設計一個滿足指標的 FIR 濾波器， 自己編寫出符合漢寧窗函數(shù)性質的函數(shù)， 對該語音信號進行濾波去噪處理，比較濾波前后 的波形和頻譜并進行分析，根據結果和學過的理論得出合理的結論 ，寫出符合要求的課程設計報告。 %程序作者： xxx %最后修改日期： 2020313 [x,fs,bits]=wavread(‘ ’ )。 sound(x,fs,bits)。 % 計算時間范圍，樣本數(shù)除以采樣頻率 y=x’ +*sin(fn*2*pi*t)。 % 對原始信號和加噪信號進行 fft變換 X=X(1:length(X)/2)。 % 計算頻譜頻率范圍 figure(1) subplot(2,2,1)。)。ylabel(39。 %加上標題和橫坐標名稱 subplot(2,2,2)。)。ylabel(39。 subplot(2,2,3)。)。ylabel(39。 subplot(2,2,4)。)。ylabel(39。 fpd=1700。 % FIR 濾波器的上下截止頻率 Rp=1。 df=min((fsdfpd),(fpufsu))。 % 將 Hz 為單位的模擬頻率換算為 rad 為單位的數(shù)字頻率 wsd=fsd/fs*2*pi。 % 定義時間范圍 w_ham=hanning(M)。 % 用窗口法計算實際濾波器脈沖響應 [db,mag,pha,grd,w]=freqz_m(h_bs,1)。以 db 為單位的幅度特性 39。)。 grid on。以線性為單位的幅度特性 39。)。 grid on。相位特性 39。)。 grid on。脈沖響應 39。)。 grid on。 % 計算頻譜取前一半 figure(3) subplot(3,2,1)。加噪聲前的時域波形 39。)。 grid on。加噪聲前的頻域波形 39。)。 g rid on。加噪聲后的時域波形 39。)。 grid on。加噪聲后的頻域波形 39。)。 grid on。濾波后的時域波形 39。)。 grid on。濾波后的頻域波形 39。)。 grid on。Y_fil39。 xlabel(39。plot(f,Y_fil)。y_fil39。 xlabel(39。plot(t,y_fil)。Y39。 xlabel(39。plot(f,Y)。y39。 xlabel(39。plot(t,y)。X39。 xlabel(39。plot(f,X)。x39。 xlabel(39。 grid on。 % 用設計好的濾波器對 y進行濾波 Y_fil=abs(fft(y_fil))。h_bs39。 23 xlabel(39。plot(h_bs)。pha39。 xlabel(39。plot(w,pha)。mag39。 xlabel(39。plot(w/pi,mag)。db39。 xlabel(39。plot(w/pi,db)。 % 調用自編函數(shù)計算理想帶阻濾波器的脈沖響應 h_bs=w_ham39。 M=ceil(*pi/dw)+1。 wcu=fcu/fs*2*pi。 % 帶阻濾波器設計指標 fcd=(fpd+fsd)/2。fsu=1850。)。f39。 title(39。)。t39。 title(39。)。f39。 title(39。)。t39。 % 布局為 2*2 的四個小圖 title(39。 % 截取前半部分 21 deltaf=fs/2/length(X)。 % 可以明顯聽出有尖銳的單頻嘯叫聲 X=fft(x)。 % 計算信號 x的長度 fn=1800。 [x,fs,bits]=wavread(‘ e:\’ )。 最后，一個學期的學習結束了，特別感謝老師的教導。 數(shù)字濾波器是 DSP 的典型應用，學會了有助于觸類 旁通，利于進一步的學習研究，能做到理解其他基于 DSP 的系統(tǒng)的功能和工作原理。然后自己查了一些相關書籍得知，繪制流程圖時 一般用橢圓表示 “ 開始 ” 與 “ 結束 ” ， 行動方案普通工作環(huán)節(jié)用矩形表示 ， 問題判斷或判定（審核 /審批 /評審）環(huán)節(jié)用菱形表示 ，箭頭代表工作流方向 ， 流程圖實例圖 ， 輸入輸出為平行四邊形 。 4. 繪制出濾波后的波形，發(fā)現(xiàn) FIR 濾波器沒有濾掉單頻噪聲。 17 4. 出現(xiàn)的問題及解決方法 在設計課程設計時，出現(xiàn) 以下 幾個問題： 1. 調用語音時，用了語句 [x,fs,bits]=wavread(‘ ’ )，出現(xiàn)錯誤，后來知道語句的單引號用了中文輸入法，最后將單引號的輸入改為英問輸入，如[x,fs,bits]=wavread(‘ ’ )，這樣就不會出錯了。 16 圖 36 濾波前后的波形對比圖 結果分析 語音信號經過 FIR 濾波器的濾除噪聲的處理， 在 Matlab 中，函數(shù) sound 可以對聲音進行回放。 在同一張大圖里分別繪制原始信號 x，加噪信號 x’ ，濾波去噪信號 y的時域波形和頻譜，以便比較和分析； 主程序如下： y_fil=fftfilt(h_bs,y)。.*hd_bs。 % 計算漢寧窗設計該濾波器時需要的階數(shù) n=0:M1。 dw=df/fs*2*pi。fcu=(fpu+fsu)/2。fpu=1900。再調用 freqz(b,a,512,fs)即可得到該濾波器的頻率響應。將語音重放時，則聽到一種尖尖刺耳的聲音， 13 說明噪聲完全覆蓋了以前的語音信號。 Y=Y(1:length(Y)/2)。 % 可以明顯聽出有尖銳的單頻嘯叫聲 對語音信號加入單頻噪聲前后的波形進行分析， 首先 畫出語音信號的時域波形；然后對語音號進行快速傅里葉變換，得到信號的頻譜特性 。 % 計算信號 x 的長度 fn=1800。 調用的程序如下： [x,fs,bits]=wavread(‘ ’ )。采樣頻率與聲音頻率之間有一 定的關系，根據奎斯特理論，只有采樣頻率高于聲音信號最高頻