彈撥音樂濾波去噪使用三角窗設(shè)計fir濾波器要點(文件)
【正文】 間（單位：s）39。)。 axis([0 ,4 4])。頻率（單位：Hz）39。)。axis([0 4000,0 35])。fs2=fn+50。 % 計算上下邊帶中心頻率，和頻率間隔 dw=df/fs*2*pi。ws2=fs2/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。 % 用窗口法計算實際濾波器脈沖響應(yīng) [db,mag,pha,grd,w]=freqz_m(h,1)。w39。)。xlabel(39。mag39。plot(w/pi,pha)。ylabel(39。 subplot(2,2,4)。)。 y_fil=filter(h,1,y)。plot(t,x)。ylabel(39。原始語音信號39。頻率（單位：Hz）39。)。axis([0 6000,0 100]) subplot(3,2,3)。)。title(39。xlabel(39。幅度譜39。)。時間（單位：s）39。)。 plot(f,Y_fil)。ylabel(39。濾波后的語音信號幅度譜39。n = [0:1:(M1)]。% freqz 子程序的改進版本% % [db,mag,pha,grd,w] = freqz_m(b,a)。H = (H(1:1:501))39。pha = angle(H)。% grd = [grd(1) grd]。% grd = median(grd)*500/pi。 grd。grd = grpdelay(b,a,w)。mag = abs(H)。whole39。hd = sin(wc*m) ./ (pi*m)。axis([0 6000,0 100])附錄2：ideal_lp函數(shù)源程序清單function hd = ideal_lp(wc,M)。)。頻率（單位：Hz）39。濾波后的語音信號39。ylabel(39。plot(t,y_fil)。title(39。)。) subplot(3,2,4)。幅度39。xlabel(39。原始語音信號幅度譜39。ylabel(39。plot(f,X)。)。時間（單位：s）39。Y_fil=Y_fil(1:N/2)。w(n)39。xlabel(39。)。w39。 axis([0 ,0 2])。)。 subplot(2,2,2)。ylabel(39。plot(w/pi,db)。 % 產(chǎn)生M階的三角窗 hd=ideal_lp(wc1,M)+ideal_lp(pi,M)ideal_lp(wc2,M)。wp2=fp2/fs*2*pi。wc2=(fs2+fp2)/2。As=15。fs1=fn50。加入干擾后的語音信號幅度譜39。ylabel(39。plot(f,Y)。加入單頻干擾后的語音信號39。ylabel(39。plot(t,y)。原始語音信號幅度譜39。ylabel(39。plot(f,X)。原始語音信號39。ylabel(39。plot(t,x)。 Y=Y(1:N/2)。 %在原信號上添加噪聲sound(y,fs,bits)。 % 單頻噪聲頻率，此參數(shù)可自己設(shè)計t=0:1/fs:(N1)/fs。)。并可綜合運用這些知識解決一定的實際問題，使我在所學知識的綜合運用能力上以及分析問題、解決問題能力上得到一定的提高。通過語音信號濾波去噪——使用三角窗設(shè)計的FIR濾波器課程設(shè)計，我對數(shù)字信號方面的知識得到了更深一步的了解。 語音文件添加好之后就是對語音添加噪聲，進行頻譜分析了，這些基本上都是以前上Matlab課程中接觸到的一些基本知識，只是很多函數(shù)是老師給的，也給我們省下了不少時間。4出現(xiàn)的問題及解決方法在課程設(shè)計的兩周時間里，剛開始的前一周是在機房上機，老師跟我們綜述本次課程設(shè)計的具體流程，在老師的引導下我們可以請教老師不懂的問題及出現(xiàn)錯誤等等。在將三個信號的時域波形和頻譜圖比較之后，我們還要通過回放去濾波去噪音樂信號，來跟原信號相比，以檢驗濾波器的效果。 % 用設(shè)計好的濾波器對加噪信號y進行濾波 Y_fil=fft(y_fil)。 % 用窗口法計算實際濾波器脈沖響應(yīng) [db,mag,pha,grd,w]=freqz_m(h,1)。 % 計算三角窗設(shè)計該濾波器時需要的階數(shù) Rp=min(db(wp1/dw+1:wp2/dw))Rp = As=max(round(db(ws1/dw+1:1:501)))As = 22n=0:M1。ws2=fs2/fs*2*pi。 % 計算上下邊帶中心頻率，和頻率間隔 dw=df/fs*2*pi。fs2=fn+50。所以，一定要將通帶截止db值設(shè)置的小于25db，所以，我將其設(shè)置為15db。 % 計算頻譜頻率范圍 音樂信號加入噪聲前后的時域頻域圖 設(shè)計數(shù)字濾波器和畫出其頻率響應(yīng)截止頻率也可以任意自選，在單頻噪聲干擾附近即可。 % 對原始信號和加噪信號進行fft變換 X=X(1:N/2)。+(fn*2*pi*t)。 % 計算信號x的長度fn=1700。D:\39。然后再
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1