彈撥音樂濾波去噪使用三角窗設(shè)計fir濾波器要點-資料下載頁
2025-06-29 12:50本頁面
【正文】 =max(round(db(ws1/dw+1:1:501)))As = 22n=0:M1。 % 定義時間范圍 w_tri=(triang(M))39。 % 產(chǎn)生M階的三角窗 hd=ideal_lp(wc1,M)+ideal_lp(pi,M)ideal_lp(wc2,M)。 % 調(diào)用自編函數(shù)計算理想帶阻濾波器的脈沖響應(yīng) h=hd.*w_tri。 % 用窗口法計算實際濾波器脈沖響應(yīng) [db,mag,pha,grd,w]=freqz_m(h,1)。 % 調(diào)用自編函數(shù)計算濾波器的頻率特性 subplot(2,2,1)。plot(w/pi,db)。xlabel(39。w39。)。ylabel(39。db39。)。 axis([0 ,300 100])。 subplot(2,2,2)。plot(w/pi,mag)。xlabel(39。w39。)。ylabel(39。mag39。)。 axis([0 ,0 2])。 subplot(2,2,3)。plot(w/pi,pha)。xlabel(39。w39。)。ylabel(39。pha39。)。 axis([0 ,4 4])。 subplot(2,2,4)。stem(n,w_tri)。xlabel(39。n39。)。ylabel(39。w(n)39。)。 y_fil=filter(h,1,y)。 % 用設(shè)計好的濾波器對加噪信號y進(jìn)行濾波 Y_fil=fft(y_fil)。Y_fil=Y_fil(1:N/2)。 % 計算頻譜取前一半 subplot(3,2,1)。plot(t,x)。xlabel(39。時間（單位：s）39。)。ylabel(39。幅度39。)。title(39。原始語音信號39。) subplot(3,2,2)。plot(f,X)。xlabel(39。頻率（單位：Hz）39。)。ylabel(39。幅度譜39。)。title(39。原始語音信號幅度譜39。) 。axis([0 6000,0 100]) subplot(3,2,3)。plot(t,y)。xlabel(39。時間（單位：s）39。)。ylabel(39。幅度39。)。title(39。加入單頻干擾后的語音信號39。) subplot(3,2,4)。plot(f,Y)。xlabel(39。頻率（單位：Hz）39。)。ylabel(39。幅度譜39。)。title(39。加入干擾后的語音信號幅度譜39。)。axis([0 6000,0 100]) subplot(3,2,5)。plot(t,y_fil)。xlabel(39。時間（單位：s）39。)。ylabel(39。幅度39。)。title(39。濾波后的語音信號39。) subplot(3,2,6)。 plot(f,Y_fil)。xlabel(39。頻率（單位：Hz）39。)。ylabel(39。幅度譜39。)。title(39。濾波后的語音信號幅度譜39。)。axis([0 6000,0 100])附錄2：ideal_lp函數(shù)源程序清單function hd = ideal_lp(wc,M)。% 理想低通濾波器計算% % [hd] = ideal_lp(wc,M)% hd = 0 to M1之間的理想脈沖響應(yīng)% wc = 截止頻率(弧度) % M = 理想濾波器的長度%alpha = (M1)/2。n = [0:1:(M1)]。m = n alpha + eps。hd = sin(wc*m) ./ (pi*m)。附錄3： freqz_m函數(shù)源程序清單function [db,mag,pha,grd,w] = freqz_m(b,a)。% freqz 子程序的改進(jìn)版本% % [db,mag,pha,grd,w] = freqz_m(b,a)。% db = [0 到pi弧度]區(qū)間內(nèi)的相對振幅(db)% mag = [0 到pi弧度]區(qū)間內(nèi)的絕對振幅% pha = [0 到pi弧度]區(qū)間內(nèi)的相位響應(yīng)% grd = [0 到pi弧度]區(qū)間內(nèi)的群遲延% w = [0 到pi弧度]區(qū)間內(nèi)的501個頻率樣本向量% b = Ha(z)的分子多項式系數(shù)(對FIR b=h)% a = Ha(z)的分母多項式系數(shù)(對 FIR: a=[1])%[H,w] = freqz(b,a,1000,39。whole39。)。H = (H(1:1:501))39。 w = (w(1:1:501))39。mag = abs(H)。db = 20*log10((mag+eps)/max(mag))。pha = angle(H)。% pha = unwrap(angle(H))。grd = grpdelay(b,a,w)。% grd = diff(pha)。% grd = [grd(1) grd]。% grd = [0 grd(1:1:500)。 grd。 grd(2:1:501) 0]。% grd = median(grd)*500/pi。
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1