【導讀】和頻域波形的程序設(shè)計。在課程設(shè)計中，系統(tǒng)操作平臺為WindowsXP，程序設(shè)。編函數(shù)，設(shè)計出FIR濾波器。此課程設(shè)計首先是用麥克風采集一段語音信號，通。圖進行對比分析。程序經(jīng)過多次的調(diào)試運行，基本實現(xiàn)了設(shè)計的目標，并且經(jīng)過。學科而又廣泛應(yīng)用于許多領(lǐng)域的新興學科。信息技術(shù)的飛速發(fā)展，數(shù)字信號處理技術(shù)應(yīng)運而生并得到迅速的發(fā)展。二十多年時間里，數(shù)字信號處理已經(jīng)在通信等領(lǐng)域得到極為廣泛的應(yīng)用。估值、增強、壓縮、識別等處理，以得到符合人們需要的信號形式。率采樣法、和等波紋最佳逼近法；第二類是最優(yōu)設(shè)計法。數(shù)字信號處理的原理進行驗證的一個過程。和交互式環(huán)境，主要包括MATLAB和Simulink兩大部分。Maple并稱為三大數(shù)學軟件。MATLAB可以進行矩陣運算、繪制函數(shù)和數(shù)據(jù)、實現(xiàn)算法、創(chuàng)建用戶界面、圖像處理、信號檢測、金融建模設(shè)計與分析等領(lǐng)域。在新的版本中也加入了對C，F(xiàn)ORTRAN，C++，JAVA的支持。

　　

【正文】 =wavread(‘ e:\’ )。 % 輸入?yún)?shù)為文件的全路徑和文件名，輸出的第一個參數(shù)是每個樣本的值， fs是生成該波形文件時的采樣率， bits 是波形文件每樣本的編碼位數(shù)。 sound(x,fs,bits)。 % 按指定的采樣率和每樣本編碼位數(shù)回放 N=length(x)。 % 計算信號 x的長度 fn=1800。 % 單頻噪 聲頻率設(shè)為 1800Hz t=0:1/fs:(length(x)1)/fs。 % 計算時間范圍，樣本數(shù)除以采樣頻率 y=x’ +*sin(fn*2*pi*t)。 % 加入一個單頻噪聲 sound(y,fs,bits)。 % 可以明顯聽出有尖銳的單頻嘯叫聲 X=fft(x)。 Y=fft(y)。 % 對原始信號和加噪信號進行 fft變換 X=X(1:length(X)/2)。 Y=Y(1:length(Y)/2)。 % 截取前半部分 21 deltaf=fs/2/length(X)。 % 計算頻譜的譜線間隔 f=0:deltaf:fs/2deltaf。 % 計算頻譜頻率范圍 figure(1) subplot(2,2,1)。plot(t,x)。 % 布局為 2*2 的四個小圖 title(39。加噪聲前的時域波形圖 39。)。 xlabel(39。t39。)。ylabel(39。x39。)。 grid on。 %加上標題和橫坐標名稱 subplot(2,2,2)。plot(f,X)。 title(39。加噪聲前的頻域波形圖 39。)。 xlabel(39。f39。)。ylabel(39。X39。)。 grid on。 subplot(2,2,3)。plot(t,y)。 title(39。加噪聲后的時域波形圖 39。)。 xlabel(39。t39。)。ylabel(39。y39。)。 grid on。 subplot(2,2,4)。plot(f,Y)。 title(39。加噪聲后的頻域波形圖 39。)。 xlabel(39。f39。)。ylabel(39。Y39。)。 grid on。 fpd=1700。fsd=1750。fsu=1850。fpu=1900。 % FIR 濾波器的上下截止頻率 Rp=1。As=36。 % 帶阻濾波器設(shè)計指標 fcd=(fpd+fsd)/2。fcu=(fpu+fsu)/2。 df=min((fsdfpd),(fpufsu))。 22 % 計算上下邊帶中心頻率，和頻率間隔 wcd=fcd/fs*2*pi。 wcu=fcu/fs*2*pi。 dw=df/fs*2*pi。 % 將 Hz 為單位的模擬頻率換算為 rad 為單位的數(shù)字頻率 wsd=fsd/fs*2*pi。 wsu=fsu/fs*2*pi。 M=ceil(*pi/dw)+1。 % 計算漢寧窗設(shè)計該濾波器時需要的階數(shù) n=0:M1。 % 定義時間范圍 w_ham=hanning(M)。 % 產(chǎn)生 M 階的漢寧窗 hd_bs=ideal_lp(wcd,M)+ideal_lp(pi,M)ideal_lp(wcu,M)。 % 調(diào)用自編函數(shù)計算理想帶阻濾波器的脈沖響應(yīng) h_bs=w_ham39。.*hd_bs。 % 用窗口法計算實際濾波器脈沖響應(yīng) [db,mag,pha,grd,w]=freqz_m(h_bs,1)。 % 調(diào)用自編函數(shù)計算濾波器的頻率特性 figure(2) subplot(2,2,1)。plot(w/pi,db)。 title(39。以 db 為單位的幅度特性 39。)。 xlabel(39。w/pi39。)。ylabel(39。db39。)。 grid on。 subplot(2,2,2)。plot(w/pi,mag)。 title(39。以線性為單位的幅度特性 39。)。 xlabel(39。w/pi39。)。ylabel(39。mag39。)。 grid on。 subplot(2,2,3)。plot(w,pha)。 title(39。相位特性 39。)。 xlabel(39。w39。)。ylabel(39。pha39。)。 grid on。 subplot(2,2,4)。plot(h_bs)。 title(39。脈沖響應(yīng) 39。)。 23 xlabel(39。t39。)。ylabel(39。h_bs39。)。 grid on。 % 加上網(wǎng)格 y_fil=fftfilt(h_bs,y)。 % 用設(shè)計好的濾波器對 y進行濾波 Y_fil=abs(fft(y_fil))。Y_fil=Y_fil(1:length(Y_fil)/2)。 % 計算頻譜取前一半 figure(3) subplot(3,2,1)。plot(t,x)。 grid on。 title(39。加噪聲前的時域波形 39。)。 xlabel(39。t39。)。ylabel(39。x39。)。 grid on。 subplot(3,2,2)。plot(f,X)。 title(39。加噪聲前的頻域波形 39。)。 xlabel(39。f39。)。ylabel(39。X39。)。 g rid on。 subplot(3,2,3)。plot(t,y)。 title(39。加噪聲后的時域波形 39。)。 xlabel(39。t39。)。ylabel(39。y39。)。 grid on。 subplot(3,2,4)。plot(f,Y)。 title(39。加噪聲后的頻域波形 39。)。 xlabel(39。f39。)。ylabel(39。Y39。)。 grid on。 subplot(3,2,5)。plot(t,y_fil)。 title(39。濾波后的時域波形 39。)。 xlabel(39。t39。)。ylabel(39。y_fil39。) 。 grid on。 subplot(3,2,6)。plot(f,Y_fil)。 title(39。濾波后的頻域波形 39。)。 xlabel(39。f39。)。ylabel(39。Y_fil39。)。 grid on。 24