基于dsp語音信號的fir濾波器處理-資料下載頁
2024-11-10 03:51本頁面
【導(dǎo)讀】和頻域波形的程序設(shè)計(jì)。在課程設(shè)計(jì)中，系統(tǒng)操作平臺為WindowsXP，程序設(shè)。編函數(shù)，設(shè)計(jì)出FIR濾波器。此課程設(shè)計(jì)首先是用麥克風(fēng)采集一段語音信號，通。圖進(jìn)行對比分析。程序經(jīng)過多次的調(diào)試運(yùn)行，基本實(shí)現(xiàn)了設(shè)計(jì)的目標(biāo)，并且經(jīng)過。學(xué)科而又廣泛應(yīng)用于許多領(lǐng)域的新興學(xué)科。信息技術(shù)的飛速發(fā)展，數(shù)字信號處理技術(shù)應(yīng)運(yùn)而生并得到迅速的發(fā)展。二十多年時間里，數(shù)字信號處理已經(jīng)在通信等領(lǐng)域得到極為廣泛的應(yīng)用。估值、增強(qiáng)、壓縮、識別等處理，以得到符合人們需要的信號形式。率采樣法、和等波紋最佳逼近法；第二類是最優(yōu)設(shè)計(jì)法。數(shù)字信號處理的原理進(jìn)行驗(yàn)證的一個過程。和交互式環(huán)境，主要包括MATLAB和Simulink兩大部分。Maple并稱為三大數(shù)學(xué)軟件。MATLAB可以進(jìn)行矩陣運(yùn)算、繪制函數(shù)和數(shù)據(jù)、實(shí)現(xiàn)算法、創(chuàng)建用戶界面、圖像處理、信號檢測、金融建模設(shè)計(jì)與分析等領(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)。 % 計(jì)算信號 x的長度 fn=1800。 % 單頻噪 聲頻率設(shè)為 1800Hz t=0:1/fs:(length(x)1)/fs。 % 計(jì)算時間范圍，樣本數(shù)除以采樣頻率 y=x’ +*sin(fn*2*pi*t)。 % 加入一個單頻噪聲 sound(y,fs,bits)。 % 可以明顯聽出有尖銳的單頻嘯叫聲 X=fft(x)。 Y=fft(y)。 % 對原始信號和加噪信號進(jìn)行 fft變換 X=X(1:length(X)/2)。 Y=Y(1:length(Y)/2)。 % 截取前半部分 21 deltaf=fs/2/length(X)。 % 計(jì)算頻譜的譜線間隔 f=0:deltaf:fs/2deltaf。 % 計(jì)算頻譜頻率范圍 figure(1) subplot(2,2,1)。plot(t,x)。 % 布局為 2*2 的四個小圖 title(39。加噪聲前的時域波形圖 39。)。 xlabel(39。t39。)。ylabel(39。x39。)。 grid on。 %加上標(biāo)題和橫坐標(biāo)名稱 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è)計(jì)指標(biāo) fcd=(fpd+fsd)/2。fcu=(fpu+fsu)/2。 df=min((fsdfpd),(fpufsu))。 22 % 計(jì)算上下邊帶中心頻率，和頻率間隔 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。 % 計(jì)算漢寧窗設(shè)計(jì)該濾波器時需要的階數(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ù)計(jì)算理想帶阻濾波器的脈沖響應(yīng) h_bs=w_ham39。.*hd_bs。 % 用窗口法計(jì)算實(shí)際濾波器脈沖響應(yīng) [db,mag,pha,grd,w]=freqz_m(h_bs,1)。 % 調(diào)用自編函數(shù)計(jì)算濾波器的頻率特性 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è)計(jì)好的濾波器對 y進(jìn)行濾波 Y_fil=abs(fft(y_fil))。Y_fil=Y_fil(1:length(Y_fil)/2)。 % 計(jì)算頻譜取前一半 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
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號-1