數(shù)字信號課程設(shè)計(jì)-基于matlab的噪聲語音信號的濾波設(shè)計(jì)與噪聲去除(編輯修改稿)
2024-12-16 01:30 本頁面
【文章內(nèi)容簡介】 global x f y SNR if isempty(SNR)==1。 msgbox(39。請輸入信噪比 ~O(∩ _∩ )O~39。,39。輸入信噪比 39。)。 else y = awgn(x,SNR)。 %加入高斯白噪聲 axes()。plot(y)。 title(39。加噪語音信號的 時(shí)域波形 39。)。 %繪制時(shí)域譜 y1=fft(y,1024)。 %傅里葉變換 axes()。plot(f,abs(y1(1:512)))。 title(39。加噪語音信號的頻域波形 39。)。 xlabel(39。頻率 /Hz39。)。ylabel(39。幅值 39。)。 %繪制頻域譜 end 程序結(jié)果如下圖： 哈爾濱工業(yè)大學(xué) （威海） 課程 設(shè)計(jì) 報(bào)告 7 數(shù)字濾波器的設(shè)計(jì) 對于數(shù)字 低通、 高通、帶通濾波器的設(shè)計(jì)，通用方法為雙線性變換法?？梢越柚谀M濾波器的頻率轉(zhuǎn)換設(shè)計(jì)一個(gè)所需類型的過渡模擬濾波器，再經(jīng)過雙線性變換將其轉(zhuǎn)換策劃那個(gè)所需的數(shù)字濾波器。具體設(shè)計(jì)步驟如下： （ 1）確定所需類型數(shù)字濾波器的技術(shù)指標(biāo)。 （ 2）將所需類型數(shù)字濾波器的邊界頻率轉(zhuǎn)換成相應(yīng)的模擬濾波器的邊界頻率，轉(zhuǎn)換公式為Ω =2/T tan( ) （ 3） 將相應(yīng)類型的模擬濾波器技術(shù)指標(biāo)轉(zhuǎn)換成模擬低通濾波器技術(shù)指標(biāo)。 （ 4）設(shè)計(jì)模擬低通濾波器。 （ 5）通過頻率變換將模擬低通轉(zhuǎn)換成相應(yīng)類型的過渡模擬濾波器。 （ 6）采用雙線性 變換法將相應(yīng)類型的過渡模擬濾波器轉(zhuǎn)換成所需類型的數(shù)字濾波器。 我們知道，脈沖響應(yīng)不變法的主要缺點(diǎn)是會產(chǎn)生頻譜混疊現(xiàn)象，使數(shù)字濾波器的頻響偏離模擬濾波器的頻響特性。為了克服之一缺點(diǎn)，可以采用雙線性變換法。 我們根據(jù)語音信號的特點(diǎn)給出有關(guān)濾波器的技術(shù)指標(biāo)： ① 低 通濾波器的性能指標(biāo)： fp=1200Hz， fc=1800Hz， Rp=1， Rs=30 哈爾濱工業(yè)大學(xué) （威海） 課程 設(shè)計(jì) 報(bào)告 8 ② 高 通濾波器的性能指標(biāo)： fp=4800Hz， fc=5000Hz， Rp=1， Rs=30 ③ 帶 通濾波器的性能指標(biāo)： fb1=2200， fb2=3000， fc1=2020， fc2=3200， Rp=1， Rs=30 IIR 數(shù)字濾波器的主要程序 及頻率響應(yīng)幅值 ： %=============================IIR 低通濾波器 =================== global fs y f z Ts=1/fs。 wp=2*pi*1200/fs。ws=2*pi*1800/fs。 Rp=1。Rs=30。 wp1=2/Ts*tan(wp/2)。 %將模擬指標(biāo)轉(zhuǎn)換成數(shù)字指標(biāo) ws1=2/Ts*tan(ws/2)。 [N,Wn]=buttord(wp1,ws1,Rp,Rs,39。s39。)。 %選擇濾波器的最小階數(shù) [Z,P,K]=buttap(N)。 %創(chuàng)建 butterworth 模擬濾波器 [Bap,Aap]=zp2tf(Z,P,K)。 [b,a]=lp2lp(Bap,Aap,Wn)。 [bz,az]=bilinear(b,a,fs)。 %用雙線性變換法實(shí)現(xiàn)模擬濾波器到數(shù)字濾波器的轉(zhuǎn)換 [H,W]=freqz(bz,az)。 %繪制頻率響應(yīng)曲線 axes()。plot(W*fs/(2*pi),abs(H)) grid title(39。IIR 低通濾波器 39。) xlabel(39。頻率／ Hz39。)。ylabel(39。頻率響應(yīng)幅度 39。) 哈爾濱工業(yè)大學(xué) （威海） 課程 設(shè)計(jì) 報(bào)告 9 %=============================IIR 高 通濾波器 =================== global fs y f z Ts=1/fs。 Wp=2*pi*5000/fs。Ws=2*pi*4800/fs。 Rp=1。Rs=30。 Wp1=2/Ts*
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1