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