彈撥音樂濾波去噪使用三角窗設(shè)計fir濾波器要點(完整版)
【正文】 號處理函數(shù)可調(diào)用，實現(xiàn)FIR濾波器相對簡單，但是由于程序順序執(zhí)行，速度受到限制。錄一段2至3秒的語音，然后左鍵單擊文件/屬性。 sound(x,fs,bits)。 % 應(yīng)該可以明顯聽出有尖銳的單頻嘯叫聲X=abs(fft(x))。這個值一定要根據(jù)我們使用的設(shè)計濾波器的方法來設(shè)定。As=15。wp2=fp2/fs*2*pi。在Matlab中，F(xiàn)IR濾波器利用函數(shù)fftfilt對信號進(jìn)行濾波，IIR濾波器利用函數(shù)filter對信號進(jìn)行濾波。其調(diào)用格式為：sound (x，fs，bits)。 確實，在本次課程設(shè)計過程中遇到了很多的困難，但是同學(xué)們的幫助讓我本次設(shè)計省了不少心，老師在大體上給我們指明方向，同學(xué)們則在我實際動手的過程中給予指點和幫助，讓我受益匪淺。% 程序作者：任健% 最后修改日期：201474 [x,fs,bits]= wavread(39。 y=x39。 % 計算頻譜的譜線間隔f=0:deltaf:fs/2deltaf。)。頻率（單位：Hz）39。 axis([0 ,0 2])。)。頻率（單位：Hz）39。axis([0 4000,0 35])。 % 計算上下邊帶中心頻率，和頻率間隔 dw=df/fs*2*pi。 % 計算三角窗設(shè)計該濾波器時需要的階數(shù) Rp=min(db(wp1/dw+1:wp2/dw))Rp = As=max(round(db(ws1/dw+1:1:501)))As = 22n=0:M1。w39。xlabel(39。plot(w/pi,pha)。 subplot(2,2,4)。 y_fil=filter(h,1,y)。ylabel(39。頻率（單位：Hz）39。axis([0 6000,0 100]) subplot(3,2,3)。title(39。幅度譜39。時間（單位：s）39。 plot(f,Y_fil)。濾波后的語音信號幅度譜39。% freqz 子程序的改進(jìn)版本% % [db,mag,pha,grd,w] = freqz_m(b,a)。pha = angle(H)。% grd = median(grd)*500/pi。grd = grpdelay(b,a,w)。whole39。axis([0 6000,0 100])附錄2：ideal_lp函數(shù)源程序清單function hd = ideal_lp(wc,M)。頻率（單位：Hz）39。ylabel(39。title(39。) subplot(3,2,4)。xlabel(39。ylabel(39。)。Y_fil=Y_fil(1:N/2)。xlabel(39。w39。)。ylabel(39。 % 產(chǎn)生M階的三角窗 hd=ideal_lp(wc1,M)+ideal_lp(pi,M)ideal_lp(wc2,M)。wc2=(fs2+fp2)/2。fs1=fn50。ylabel(39。加入單頻干擾后的語音信號39。plot(t,y)。ylabel(39。原始語音信號39。plot(t,x)。 %在原信號上添加噪聲sound(y,fs,bits)。)。通過語音信號濾波去噪——使用三角窗設(shè)計的FIR濾波器課程設(shè)計，我對數(shù)字信號方面的知識得到了更深一步的了解。4出現(xiàn)的問題及解決方法在課程設(shè)計的兩周時間里，剛開始的前一周是在機房上機，老師跟我們綜述本次課程設(shè)計的具體流程，在老師的引導(dǎo)下我們可以請教老師不懂的問題及出現(xiàn)錯誤等等。 % 用設(shè)計好的濾波器對加噪信號y進(jìn)行濾波 Y_fil=fft(y_fil)。 % 計算三角窗設(shè)計該濾波器時需要的階數(shù) Rp=min(db(wp1/dw+1:wp2/dw))Rp = As=max(round(db(ws1/dw+1:1:501)))As = 22n=0:M1。 % 計算上下邊帶中心頻率，和頻率間隔 dw=df/fs*2*pi。所以，一定要將通帶截止db值設(shè)置的小于25db，所以，我將其設(shè)置為15db。 % 對原始信號和加噪信號進(jìn)行fft變換 X=X(1:N/2)。 % 計算信號x的長度fn=1700。然后再點擊文件/另存為，即可將該語音文件存在MATLAB文件下的work文件夾里。（3）可編程FIR濾波器 還有一種是使用可編程邏輯器件，F(xiàn)PGA／CPLD。常見的窗函數(shù)，可以分為以下主要類型：（1） 冪窗采用時間變量某種冪次的函數(shù)，如矩形、三角形、梯形或其它時間（t）的高次冪； （2）三角函數(shù)窗應(yīng)用三角函數(shù)，即正弦或余弦函數(shù)等組合成復(fù)合函數(shù)，例如漢寧窗、海明窗等；（3） 指數(shù)窗采用指數(shù)時間函數(shù)，例如高斯窗等。多維(MD)信號：多個自變量的函數(shù)。數(shù)字信號處理的特征和分類。是一個可以完成各種精確計算和數(shù)據(jù)處理的、可視化的、強大的計算工具。設(shè)計FIR數(shù)字濾波器的方法有窗函數(shù)法、頻率采樣法和等波紋最佳逼近法等。數(shù)字信號處理的目的是對真實世界的連續(xù)模擬信號進(jìn)行測量或濾波。（4）獨立完成課程設(shè)計并按要求編寫課程設(shè)計報告書。 《彈撥音樂濾波去噪——使用三角窗設(shè)計的FIR濾波器》 第20頁 共20頁 長沙理工大學(xué)城南學(xué)院《數(shù)字信號處理》
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1