吹管音樂濾波去噪使用雙線性變換法設計切比雪夫ii型濾波器(參考版)
2025-01-21 15:51本頁面
【正文】 。 grd(2:1:501) 0]。% grd = [0 grd(1:1:500)。% grd = diff(pha)。% pha = unwrap(angle(H))。 db = 20*log10((mag+eps)/max(mag))。 w = (w(1:1:501))39。)。% db = [0 到pi弧度]區(qū)間內(nèi)的相對振幅(db)% mag = [0 到pi弧度]區(qū)間內(nèi)的絕對振幅% pha = [0 到pi弧度]區(qū)間內(nèi)的相位響應% 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。附錄三：function [db,mag,pha,grd,w] = freqz_m(b,a)。\n*** Chebyshev2 Filter Order = % \n39。N = ceil(log10(g+sqrt(g*g1))/log10(OmegaR+sqrt(OmegaR*OmegaR1)))。OmegaR = Ws/Wp。A = 10^(As/20)。PB ripple and/or SB attenuation ust be larger than 039。Stopband edge must be larger than Passband edge39。Passband edge must be larger than 039。 (Rp 0)% As = Stopband attenuation in +dB。 Wp 0% Ws = Stopband edge frequency in rad/sec。% Analog Lowpass Filter Design: Chebyshev2% % [b,a] = afd_chb2(Wp,Ws,Rp,As)。% 繪制去噪吹管音樂頻譜圖sound(y_fil,fs,bits)。axis tight。% 繪制加噪吹管音樂頻譜圖subplot(212)。axis tight。gridaxis([0,4*10^3,0,2000])figuresubplot(211)。濾波后的吹管音樂幅度譜圖39。)。ylabel(39。頻率(單位:Hz)39。plot(f,Y_fil)。)。title(39。幅度39。)。xlabel(39。gridaxis([0,4*10^3,0,2000])subplot(3,2,5)。加噪后的吹管音樂幅度譜圖39。)。ylabel(39。頻率(單位:Hz)39。plot(f,Y)。)。title(39。幅度39。)。xlabel(39。gridaxis([0,4*10^3,0,2000])subplot(3,2,3)。原始吹管音樂幅度譜圖39。)。ylabel(39。頻率(單位:Hz)39。plot(f,X)。)。title(39。幅度39。)。xlabel(39。 % 計算頻譜取前一半subplot(3,2,1)。 % 用設計好的濾波器對y進行濾波Y_fil=abs(fft(y_fil))。figure。Subplot(2,2,4)。Subplot(2,2,3)。Subplot(2,2,2)。Subplot(2,2,1)。ha=filter(b,a,delta)。[C,B,A] = dir2cas(b,a) %雙線性變換法實現(xiàn)模擬濾波器到數(shù)字濾波器的轉換C = B = 0A = 0[db,mag,pha,grd,w]=freqz_m(b,a)。 %截止頻率預畸變[cs,ds] = afd_chb2(Omegap,Omegas,Rp,As)。 %定義采樣間隔Omegap=(2/T)*tan(wp/2)。 %計算對應的數(shù)字頻率T=1。 % 定義通帶波紋和阻帶衰減wp=fp/fs*2*pi。 %定義通帶和阻帶截止頻率Rp=1。grid使用雙線性變換法設計切比雪夫II型濾波器：fp=fn200。)。title(39。幅度譜39。)。xlabel(39。gridsubplot(2,2,4)。加噪后的吹管音樂39。)。ylabel(39。時間(單位:s)39。plot(t,y)。 axis([0,4*10^3,0,6000])。原
環(huán)評公示相關推薦
鄂ICP備17016276號-1