吹管音樂濾波去噪基于漢寧創(chuàng)的fir濾波器(存儲版)
【正文】 次改動單頻噪聲頻率值，所設(shè)計的濾波器的性能指標(biāo)沒有隨之變化，出現(xiàn)了濾波上的錯誤。最后，在團結(jié)合作方面，雖然每個人都單獨分配了課題，但是總能找到與自己課題相關(guān)或者類似的同學(xué)，跟同學(xué)交流經(jīng)驗成為本次課設(shè)中的一個重要環(huán)節(jié)。h:2013DSP\39。y=x+sin(fn*2*pi*t)。 % 計算頻譜頻率范圍figure(1)subplot(2,2,1)。ylabel(39。原始音樂信號頻譜39。axis([0,3000,0,3000])。時間(t)39。plot(f,Y)。幅度譜39。 % FIR濾波器的上下截止頻率Rp=1。 % 將Hz為單位的模擬頻率換算為rad為單位的數(shù)字頻率wsd=fsd/fs*2*pi。 % 用窗口法計算實際濾波器脈沖響應(yīng)[db,mag,pha,grd,w]=freqz_m(h_bs,1)。w/pi39。axis([0,])。mag39。xlabel(39。subplot(2,2,4)。)。 % 計算頻譜取前一半figure(3)subplot(3,2,1)。)。title(39。)。xlabel(39。plot(f,Y)。ylabel(39。濾波后音樂信號時間y39。grid on。頻率（f）39。% 理想低通濾波器計算% % [hd] = ideal_lp(wc,M)% hd = 0 to M1之間的理想脈沖響應(yīng)% wc = 截止頻率(弧度) % M = 理想濾波器的長度alpha = (M1)/2。)。% grd = diff(pha)。% pha = unwrap(angle(H))。% 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。sound (y_fil,fs,bits)。)。幅度39。axis([0,2,])。頻率（f）39。grid on。加干擾音樂信號時間x139。ylabel(39。plot(f,X)。xlabel(39。 % 用設(shè)計好的濾波器對y進行濾波Y_fil=fft(y_fil)。xlabel(39。axis([0,3,4,4])。濾波器相位響應(yīng)圖39。)。subplot(2,2,2)。)。 % 調(diào)用自編函數(shù)計算理想帶阻濾波器的脈沖響應(yīng)h_bs=w_ham39。wcu=fcu/fs*2*pi。fpu=2350。)。 grid on。)。幅度譜39。plot(f,X)。時間(t)39。 % 截取前半部分deltaf=fs/ N。 % 單頻噪聲頻率t=0:1/fs:(N1)/fs。G普羅克斯. 數(shù)字信號處理[M]. 西安：西安交通大學(xué)出版社，[6][M].：西安電子科技出版社，2002． 附錄1：語音信號濾波去噪設(shè)計源程序清單% 程序名稱：% 程序功能：采用基于漢寧的窗口設(shè)計法，設(shè)計FIR濾波器對含噪語音進行濾波去噪處理。此間，我得到了以下收獲：首先，在學(xué)習(xí)方面，雖然已經(jīng)學(xué)習(xí)過了DSP課程，但是如何融洽的把理論和實際結(jié)合仍然是我需要面對的問題。4出現(xiàn)的問題及解決方法 在設(shè)計課程設(shè)計時，出現(xiàn)了以下幾個問題：1.、由于沒有控制好所需音樂信號的時間，導(dǎo)致結(jié)果不理想。 % 用設(shè)計好的濾波器對y進行濾波Y_fil=fft(y_fil)。 % 產(chǎn)生M階的漢寧窗 hd_bs=ideal_lp(wcd,M)+ideal_lp(pi,M)ideal_lp(wcu,M)。 % 計算上下邊帶中心頻率，和頻率間隔wcd=fcd/fs*2*pi。fsd=1650。grid on。)。幅度39。plot(t,y)。頻率(f)39。 %改變橫縱坐標(biāo)的范圍axis([0,2,])。原始音樂信號39。 % 對原始信號和加噪信號進行fft變換X=X(1:N/2)。 % 計算信號x的長度fn=2100。首先，輸入原始音樂信號并播放一次。漢寧窗函數(shù)的時域幅度與頻域幅度特性曲線的MATLAB實現(xiàn)的曲線圖如圖21所示。窗函數(shù)的類型可根據(jù)最小阻帶衰減AS獨立選擇。設(shè)FIR濾波器的單位脈沖響應(yīng)的長度為，則其頻率響應(yīng)函數(shù)為 （21）一般將表示成如下形式： （22）式中，是的實函數(shù)（可以去負(fù)值）。2 設(shè)計原理 數(shù)字信號處理數(shù)字信號處理（Digital Signal Processing，簡稱DSP）是一門涉及許多學(xué)科而又廣泛應(yīng)用于許多領(lǐng)域的新興學(xué)科[4]。[1]MATLAB軟件包括五大通用功能：數(shù)值計算功能（Nemeric）；符號運算功能（Symbolic）；數(shù)據(jù)可視化功能（Graphic）；數(shù)據(jù)圖形文字統(tǒng)一處理功能（Notebook）和建模仿真可
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1