吹管音樂濾波去噪使用雙線性變換法設(shè)計(jì)切比雪夫ii型濾波器(參考版)
2025-06-06 22:30本頁面
【正文】 % grd = median(grd)*500/pi。 grd。 % grd = [grd(1) grd]。 grd = grpdelay(b,a,w)。 pha = angle(H)。 mag = abs(H)。 H = (H(1:1:501))39。whole39。 % freqz 子程序的改進(jìn)版本 % % [db,mag,pha,grd,w] = freqz_m(b,a)。,N) [b,a]=u_chb2ap(N,As,Ws)。 fprintf(39。 g = sqrt(A*A1)/ep。 OmegaC = Wp。) end ep = sqrt(10^(Rp/10)1)。) end if (Rp = 0) | (As 0) error(39。) 鄒磊 《吹管音樂 濾波去噪 — 使用 雙線性變換設(shè)計(jì)切比雪夫 II型濾波器》 第 22頁 共 23頁 end if Ws = Wp error(39。 (As 0) % if Wp = 0 error(39。 Ws Wp 0 % Rp = Passband ripple in +dB。 % b = Numerator coefficients of Ha(s) % a = Denominator coefficients of Ha(s) % Wp = Passband edge frequency in rad/sec。 附錄二： afd_chb2 函數(shù) function [b,a] = afd_chb2(Wp,Ws,Rp,As)。grid on。plot(f,20*log10(abs(Y_fil)/max(abs(X))))。grid on。plot(f,20*log10(abs(Y)/max(abs(X))))。)。title(39。幅度譜 39。)。xlabel(39。grid axis([0,3,1,1]) subplot(3,2,6)。濾波后的吹管音樂 39。)。ylabel(39。時(shí)間 (單位 :s)39。plot(t,y_fil)。)。title(39。幅度譜 39。)。xlabel(39。grid subplot(3,2,4)。加噪后的吹管音樂39。)。ylabel(39。時(shí)間 (單位 :s)39。plot(t,y)。)。title(39。幅度譜 39。)。xlabel(39。grid subplot(3,2,2)。原始吹管音樂 39。)。ylabel(39。時(shí)間 (單位 :s)39。plot(t,x)。Y_fil=Y_fil(1:N/2)。 %對加噪信號進(jìn)行濾波處理并 畫出加噪前后的信號時(shí)域和頻譜： y_fil=filter(b,a,y)。plot(ha)。plot(w/pi,pha)。plot(w/pi,mag)。plot(w/pi,db)。 figure(2)。 %繪制數(shù)字濾波器頻率響應(yīng)幅度圖 delta=[1,zeros(1,99)]。 %選擇濾波器最小階數(shù) *** Chebyshev2 Filter Order = 11 [b,a] = bilinear(cs,ds,Fs)。 Omegas=(2/T)*tan(ws/2)。Fs = 1/T。 ws=fc/fs*2*pi。As=35。fc=fn50。 axis([0,4*10^3,0,6000])。加噪后的吹管音樂幅度譜圖 39。)。ylabel(39。頻率 (單位 :Hz)39。plot(f,Y)。)。title(39。幅度 39。)。xlabel(39。grid subplot(2,2,3)。)。title(39。幅度譜 39。)。xlabel(39。grid subplot(2,2,2)。原始吹管音樂 39。)。ylabel(39。時(shí)間 (單位 :s)39。plot(t,x)。 % 計(jì)算頻譜的譜線間隔 f=0:deltaf:fs/2deltaf。 Y=Y(1:N/2)。 Y=abs(fft(y))。 sound(y,fs,bits)。 % 計(jì)算時(shí)間范圍，樣本數(shù)除以采樣頻率 x=x(:,1)39。 % 計(jì)算信號 x 的長度 fn=2200。 sound(x,fs,bits)。)。 % 程序作者 ：鄒磊 /% 最后修改日期 ： 2021316 [x,fs,bits]=wavread(39。 最后向老師和同學(xué)們表示最真誠的謝意。 通過這次課程設(shè)計(jì)，讓我對這些抽象的概鄒磊 《吹管音樂 濾波去噪 — 使用 雙線性變換設(shè)計(jì)切比雪夫 II型濾波器》 第 16頁 共 23頁 念的理解轉(zhuǎn)化為實(shí)實(shí)在在的運(yùn)用與擴(kuò)展，從中感受到 不同與課本的知識 運(yùn)用，加強(qiáng)了我對本科目的理解和升華。 我們在 學(xué)習(xí)這門課程時(shí) ,普遍感到數(shù)字信號處理的概念抽象 ,對其中的分析方法與基本理論不能很好地理解與掌握。 5 結(jié)束語 本次 數(shù)字信號處理 課程設(shè)計(jì)歷時(shí)兩個(gè)星期， 此次課程設(shè)計(jì) 采用 Matlab 進(jìn)行數(shù)字信號處理課程設(shè)計(jì)， 明顯可以看出 ， 通過這次設(shè)計(jì)，讓我們 加深了對課堂抽象概念的理解，鞏固了課堂上所學(xué)的理論知識，并能很好地理解與掌握數(shù)字信號處理中的基本概念、基本原理、基本 分析方法。 鄒磊 《吹管音樂 濾波去噪 — 使用 雙線性變換設(shè)計(jì)切比雪夫 II型濾波器》 第 15頁 共 23頁 4 出現(xiàn)的問題及解決方法 （ 1）采集的語音信號 需要設(shè)置成 8000Hz， 8 位，單聲道 ， 在得到理想的波形時(shí)，通過多次嘗試，和查找書籍及同學(xué)討論，最后猜得到理想的吹管音樂的時(shí)域圖和頻譜圖； （ 2）最好將角頻率的頻譜圖轉(zhuǎn)化為頻率的頻譜圖，這樣在觀察通帶截止頻率和阻帶截止頻率是更加精確，誤差更小； （ 3） 由于使用的 MATLAB 軟件版本不同，導(dǎo)致在自己機(jī)器上始終不能畫出正確的波形，最后只能在其他機(jī)器上 完成此次的波形圖； （ 4） 我們設(shè)計(jì)的濾波器的阻帶截止頻率，在畫出的頻譜圖上應(yīng)該更改坐標(biāo)，由此觀察設(shè)計(jì)的濾波器是否滿足要求，若不滿足則可以及時(shí)更改。 在 MATLAB 中，經(jīng) sound 函數(shù)，對經(jīng)過切比雪夫 II 型濾波器之后的信號進(jìn)行回放，可以聽出濾波之后的信號比原始信號更清晰一些，清除了環(huán)境噪音。 并比較加噪與去噪的吹管音樂頻譜圖如圖 所示。Y_fil=Y_fil(1:N/2)。 圖 切比雪夫?yàn)V波器的幅度和相位
畢業(yè)設(shè)計(jì)相關(guān)推薦
鄂ICP備17016276號-1