數(shù)字信號處理課程設(shè)計-用雙線性變換法設(shè)計高通和帶通數(shù)字濾波器-資料下載頁
2025-06-06 00:49本頁面
【正文】 換成相應的過度模擬濾波 器指標 wph=[wph1,wph2]。wsh=[wsh1,wsh2]。%帶通指標 [N,wc]=buttord(wph,wsh,ap,as,39。s39。)。%求出模擬濾波器階數(shù) N 和截止頻率 [B,A]=butter(N,wc,39。s39。)。%求出模擬帶通系統(tǒng)函數(shù) H（ s）的分子分母系數(shù)向量 B,A [Bz,Az]=bilinear(B,A,)。%雙線性法把模擬高通轉(zhuǎn)化為數(shù)字高通 [H,w]=freqz(Bz,Az)。 figure(2)。 subplot(2,1,1)。plot(w*Fs/2/pi,abs(H))。title(39。帶通濾波器幅頻特性 39。)。grid on。 xlabel(39。Hz39。)。ylabel(39。幅度 39。)。 subplot(2,1,2)。plot(w*Fs/2/pi,angle(H))。xlabel(39。Hz39。)。ylabel(39。相位 39。)。 title(39。帶通濾波器相頻特性 39。)。grid on。 w=1:n/2。 y1=filter(Bz,Az,y)。Y1=fft(y1,n)。Y1=Y1(1:n/2,:)。 figure(1)。 subplot(3,2,2)。plot(t,y1)。xlabel(39。秒 39。)。title(39。濾波后音樂信號波形 39。)。grid on。 subplot(3,2,4)。plot(w*Fs/n,abs(Y1))。xlabel(39。Hz39。)。title(39。濾波后音樂信號頻譜 39。)。grid on。 subplot(3,2,6)。plot(w*Fs/n,angle(Y1))。xlabel(39。Hz39。)。ylabel(39。 相位 39。)。title(39。 濾波后音樂信號相位39。)。grid on。 sound(y1,Fs,bits)。 figure(3)。 freqz(Bz,Az) 14 高 通濾波器 部分： clc。close all。 %提取音樂信號及畫出前一半頻譜 [y,Fs,bits]=wavread(39。C:\Users\pcadmin\Desktop\數(shù)字信號實驗 11\調(diào)好頻率的音頻 .wav39。)。 n=length(y)。%求出語音信號的長度 t=(0:n1)/Fs。 Y=fft(y,n)。%傅里葉變換 k=1:n/2。 Y=Y(1:n/2,:)。%取一半頻譜 figure(1)。 subplot(3,2,1)。plot(t,y)。xlabel(39。秒 39。)。title(39。濾波前音樂信號波形 39。)。grid on。 subplot(3,2,3)。plot(k*Fs/n,abs(Y))。xlabel(39。Hz39。)。title(39。濾波前音樂信號頻譜 39。)。grid on。 subplot(3,2,5)。plot(k*Fs/n,angle(Y))。xlabel(39。Hz39。)。ylabel(39。相位 39。)。title(39。濾波前音樂信號相位 39。)。grid on。 %高通濾波器設(shè)計 Fs=16000。fp=6000。fs=5000。rp=1。as=70。ap=1。%已知條件 wp=2*pi*fp/Fs。ws=2*pi*fs/Fs。 %模擬 頻率變數(shù)字頻率 wph=tan(wp/2)。wsh=tan(ws/2)。%預畸變校正將數(shù)字濾波器的指標轉(zhuǎn)換成相應的過度模擬濾波器指標 [N,wc]=buttord(wph,wsh,ap,as,39。s39。)。%求出模擬濾波器階數(shù) N 和截至頻率 [B,A]=butter(N,wc,39。high39。,39。s39。)。%求出模擬高通系統(tǒng)函數(shù) H（ s）的分子分母系數(shù)向量 B,A [Bz,Az]=bilinear(B,A,)。%雙線性法把模擬高通轉(zhuǎn)化為數(shù)字高通 [H,w]=freqz(Bz,Az)。 figure(2)。 subplot(2,1,1)。plot(w*Fs/2/pi,abs(H))。title(39。高通濾波器幅頻特性 39。)。grid on。 xlabel(39。Hz39。)。ylabel(39。幅度 39。)。 subplot(2,1,2)。plot(w*Fs/2/pi,angle(H))。xlabel(39。Hz39。)。ylabel(39。相位 39。)。 title(39。高通濾波器相頻特性 39。)。grid on。 w=1:n/2。 y1=filter(Bz,Az,y)。Y1=fft(y1,n)。Y1=Y1(1:n/2,:)。 figure(1)。 subplot(3,2,2)。plot(t,y1)。xlabel(39。秒 39。)。title(39。濾波后音樂信號波形 39。)。grid on。 subplot(3,2,4)。plot(w*Fs/n,abs(Y1))。xlabel(39。Hz39。)。title(39。濾波后音樂信號頻譜 39。)。grid on。 subplot(3,2,6)。plot(w*Fs/n,angle(Y1))。xlabel(39。Hz39。)。ylabel(39。 相位 39。)。title(39。 濾波后音樂信號相位39。)。grid on。 sound(y1,Fs,bits)。 figure(3)。 freqz(Bz,Az)
畢業(yè)設(shè)計相關(guān)推薦
鄂ICP備17016276號-1