吹管音樂濾波去噪基于漢寧創(chuàng)的fir濾波器(參考版)
2025-06-30 23:09本頁面
【正文】 。 grd(2:1:501) 0]。% grd = [0 grd(1:1:500)。% grd = diff(pha)。% pha = unwrap(angle(H))。db = 20*log10((mag+eps)/max(mag))。 w = (w(1:1:501))39。)。% 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。附錄3：function [db,mag,pha,grd,w] = freqz_m(b,a)。m = n alpha + eps。% 理想低通濾波器計算% % [hd] = ideal_lp(wc,M)% hd = 0 to M1之間的理想脈沖響應(yīng)% wc = 截止頻率(弧度) % M = 理想濾波器的長度alpha = (M1)/2。sound (y_fil,fs,bits)。)。ylabel(39。頻率（f）39。)。title(39。plot(f,Y_fil)。grid on。幅度39。)。xlabel(39。濾波后音樂信號時間y39。axis([0,2,])。subplot(3,2,5)。)。ylabel(39。頻率（f）39。)。title(39。plot(f,Y)。grid on。幅度39。)。xlabel(39。加干擾音樂信號時間x139。axis([0,2,])。subplot(3,2,3)。)。ylabel(39。頻率（f）39。)。title(39。plot(f,X)。grid on。幅度39。)。xlabel(39。原始音樂信號時間x39。 axis([0,2,])。 % 計算頻譜取前一半figure(3)subplot(3,2,1)。 % 用設(shè)計好的濾波器對y進行濾波Y_fil=fft(y_fil)。grid on。h(n)39。)。xlabel(39。濾波器脈沖響應(yīng)圖39。axis([0,800,1])。subplot(2,2,4)。axis([0,3,4,4])。相位（pha）39。)。xlabel(39。濾波器相位響應(yīng)圖39。plot(w,pha)。grid on。mag39。)。xlabel(39。以線性為單位的幅度特性39。axis([0,])。subplot(2,2,2)。)。ylabel(39。w/pi39。)。title(39。plot(w/pi,db)。 % 用窗口法計算實際濾波器脈沖響應(yīng)[db,mag,pha,grd,w]=freqz_m(h_bs,1)。 % 調(diào)用自編函數(shù)計算理想帶阻濾波器的脈沖響應(yīng)h_bs=w_ham39。 % 定義時間范圍w_ham=hanning(M)。M=ceil(*pi/dw)+1。 % 將Hz為單位的模擬頻率換算為rad為單位的數(shù)字頻率wsd=fsd/fs*2*pi。wcu=fcu/fs*2*pi。df=min((fsdfpd),(fpufsu))。 % 帶阻濾波器設(shè)計指標fcd=(fpd+fsd)/2。 % FIR濾波器的上下截止頻率Rp=1。fpu=2350。fpd=1450。axis([0,3000,0,3000])。幅度譜39。)。xlabel(39。加入單頻干擾后的音樂信號頻譜39。plot(f,Y)。 grid on。)。ylabel(39。時間(t)39。)。title(39。 subplot(2,2,3)。axis([0,3000,0,3000])。幅度譜39。)。xlabel(39。原始音樂信號頻譜39。plot(f,X)。 %加上標題和橫坐標名稱grid on。)。ylabel(39。時間(t)39。)。 % 布局為2*2的四個小圖title(39。 % 計算頻譜頻率范圍figure(1)subplot(2,2,1)。 % 截取前半部分deltaf=fs/ N。 % 對原始信號和加噪信號進行fft變換X=X(1:N/2)。 % 可以明顯聽出有尖銳的單頻嘯叫聲X=abs(fft(x))。y=x+sin(fn*2*pi*t)。 % 單頻噪聲頻率t=0:1/fs:(N1)/fs。 % 按指定的采樣率和每樣本編碼位數(shù)回放N=length(x)。 % 輸入?yún)?shù)為文件的全路徑和文件名，輸出的第一個參數(shù)是每個樣本的值，fs是生成該波形文件時的采樣率，bits是波形文件每樣本的編碼位數(shù)。h:2013DSP\39。G普羅克斯. 數(shù)字信號處理[M]. 西安：西安交通大學出版社，[6][M].：西安電子科技出版社，2002． 附錄1：語音信號濾波去噪設(shè)計源程序清
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1