數(shù)字信號課程設(shè)計(jì)-基于matlab的噪聲語音信號的濾波設(shè)計(jì)與噪聲去除-文庫吧資料
2024-11-18 01:30本頁面
【正文】 00。fb2=3000。) 哈爾濱工業(yè)大學(xué) （威海） 課程 設(shè)計(jì) 報(bào)告 10 %=============================IIR 帶 通濾波器 =================== global fs y f z Ts=1/fs。ylabel(39。頻率／ Hz39。IIR 高通濾波器 39。 %繪制頻率響應(yīng)曲線 axes()。 [bz,az]=bilinear(b,a,fs)。 %創(chuàng)建切比雪夫模擬濾波器 [Bap,Aap]=zp2tf(Z,P,K)。)。 [N,Wn]=cheb2ord(Wp1,Ws1,Rp,Rs,39。 Wp1=2/Ts*tan(Wp/2)。 Rp=1。 Wp=2*pi*5000/fs。頻率響應(yīng)幅度 39。)。) xlabel(39。plot(W*fs/(2*pi),abs(H)) grid title(39。 %用雙線性變換法實(shí)現(xiàn)模擬濾波器到數(shù)字濾波器的轉(zhuǎn)換 [H,W]=freqz(bz,az)。 [b,a]=lp2lp(Bap,Aap,Wn)。 %選擇濾波器的最小階數(shù) [Z,P,K]=buttap(N)。s39。 %將模擬指標(biāo)轉(zhuǎn)換成數(shù)字指標(biāo) ws1=2/Ts*tan(ws/2)。Rs=30。ws=2*pi*1800/fs。 我們根據(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。 我們知道，脈沖響應(yīng)不變法的主要缺點(diǎn)是會產(chǎn)生頻譜混疊現(xiàn)象，使數(shù)字濾波器的頻響偏離模擬濾波器的頻響特性。 （ 5）通過頻率變換將模擬低通轉(zhuǎn)換成相應(yīng)類型的過渡模擬濾波器。 （ 2）將所需類型數(shù)字濾波器的邊界頻率轉(zhuǎn)換成相應(yīng)的模擬濾波器的邊界頻率，轉(zhuǎn)換公式為Ω =2/T tan( ) （ 3） 將相應(yīng)類型的模擬濾波器技術(shù)指標(biāo)轉(zhuǎn)換成模擬低通濾波器技術(shù)指標(biāo)?？梢越柚谀M濾波器的頻率轉(zhuǎn)換設(shè)計(jì)一個(gè)所需類型的過渡模擬濾波器，再經(jīng)過雙線性變換將其轉(zhuǎn)換策劃那個(gè)所需的數(shù)字濾波器。)。ylabel(39。頻率 /Hz39。)。 title(39。 %傅里葉變換 axes()。)。 title(39。 %加入高斯白噪聲 axes()。)。,39。 msgbox(39。String39。 在本次課程設(shè)計(jì)中，我們是用 MATLAB中的 AWGN 函數(shù) 產(chǎn)生噪聲加入到語音信號中，模仿語音信號被污染，并對其頻譜分析。 WGN用于產(chǎn)生高斯白噪聲， AWGN 則用于在某一信號中加入高斯白噪聲。 %繪制頻域譜 程序結(jié)果如下圖： 哈爾濱工業(yè)大學(xué) （威海） 課程 設(shè)計(jì) 報(bào)告 5 噪聲的產(chǎn)生與合成 MATLAB 中產(chǎn)生高斯白噪聲非常方便，有兩個(gè)產(chǎn)生高斯白噪聲的兩個(gè)函數(shù) 。幅值 39。)。 xlabel(39。原始信號的頻域波形 39。plot(f,abs(x1(1:512)))。 %傅里葉變換 f=fs*(0:511)/1024。 % 繪制時(shí)域譜 n=length(x)。原始信號的時(shí)域波形 39。plot(x)。 FileName])。 [x,fs] = wavread([PathName 39。打開音頻文件 39。*.wav39。 我們利用 fft 對語音信號進(jìn)行快速傅里葉變換，就可以得到信號的頻譜特性。再對其進(jìn)行采樣，記住采樣頻率和
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1