數(shù)字信號課程設(shè)計(jì)-基于matlab的噪聲語音信號的濾波設(shè)計(jì)與噪聲去除(已改無錯字)
2022-12-23 01:30:26 本頁面
【正文】 tan(Wp/2)。 %將模擬指標(biāo)轉(zhuǎn)換成數(shù)字指標(biāo) Ws1=2/Ts*tan(Ws/2)。 [N,Wn]=cheb2ord(Wp1,Ws1,Rp,Rs,39。s39。)。 %選擇濾波器的最小階數(shù) [Z,P,K]=cheb2ap(N,Rs)。 %創(chuàng)建切比雪夫模擬濾波器 [Bap,Aap]=zp2tf(Z,P,K)。 [b,a]=lp2hp(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)告 10 %=============================IIR 帶 通濾波器 =================== global fs y f z Ts=1/fs。 fb1=2200。fb2=3000。fc1=2020。fc2=3200。 W1=2*fb1*pi/fs。W2=2*fc1*pi/fs。 W3=2*fb2*pi/fs。W4=2*fc2*pi/fs。 Wp=[W1,W3]。Ws=[W2,W4]。 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]=cheb2ord(Wp1,Ws1,Rp,Rs,39。s39。)。 %選擇濾波器的最小階數(shù) [Z,P,K]=cheb2ap(N,Rs)。 %創(chuàng)建切比雪夫模擬濾波器 [Bap,Aap]=zp2tf(Z,P,K)。 [b,a]=lp2bp(Bap,Aap,2200*2*pi,2020*2*pi)。 [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 哈爾濱工業(yè)大學(xué) （威海） 課程 設(shè)計(jì) 報(bào)告 11 title(39。IIR 帶通濾波器 39。) xlabel(39。頻率／ Hz39。)。ylabel(39。頻率響應(yīng)幅度 39。) 對語音信號進(jìn)行濾波處理 用自己設(shè)計(jì)的各濾波器分別對加噪的語音信號進(jìn)行濾波，在 Matlab 中， IIR 濾波器利用函數(shù) filter 對信號進(jìn)行濾波。函數(shù) filter 的調(diào)用格式：yn=filter(B， A， xn)它是按照直線型結(jié)構(gòu)實(shí)現(xiàn)對 xn 的濾波。其中 xn 是輸入信號向量， yn 輸出信號向量。 濾波的主要程序 及濾波后信號的頻域波形 z=filter(bz,az,y)。 F=fft(z,1024)。 axes()。plot(f,abs(F(1:512)))。 title(39。IIR 濾波器濾波后的頻譜 39。) xlabel(39。頻率 /Hz39。)。ylabel(39。幅值 39。)。 %畫出濾波后的頻譜圖 哈爾濱工業(yè)大學(xué) （威海） 課程 設(shè)計(jì) 報(bào)告 12 設(shè)計(jì)一個信號處理系統(tǒng)界面 圖形用戶界面（ graphical user interfaces ， GUI）則是由窗口、光標(biāo)、按
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1