數(shù)字濾波器及在語音信號分析中的作用_數(shù)字信號處理課程設(shè)計(編輯修改稿)
2025-04-09 14:56 本頁面
【文章內(nèi)容簡介】 1012加單頻余弦信號的語音信號時域波形0 0 . 5 1 1 . 5 2 2 . 5 3 3 . 5x 1 04050001000015000加單頻余弦信號的語音信號頻譜 8 圖 4 加單頻噪聲波形 0 1 2 3 4 5 6 7 8 9 10x 1 0421012加多頻余弦信號的語音信號時域波形0 0 . 5 1 1 . 5 2 2 . 5 3 3 . 5x 1 04050001000015000加多頻余弦信號的語音信號頻譜 圖 5 加多頻余弦信號波形 數(shù)字濾波器的設(shè)計及濾波 矩形窗低通濾波器設(shè)計及濾波結(jié)果 %矩形窗 %低通 wp=*pi。ws=*pi。 Bt=wswp。 N=ceil(*pi/Bt)。 wc=(wp+ws)/2/pi。 hn=fir1(N,wc,boxcar(N+1))。 [H,f]=freqz(hn,1,512,8000)。 %采用 8000Hz 的采樣頻率求出頻率響應(yīng) figure(4) subplot(2,1,1),plot(f,20*log10(abs(H))) xlabel(39。頻率 /Hz39。)。ylabel(39。振幅 /dB39。)。grid on。 9 subplot(2,1,2),plot(f,180/pi*unwrap(angle(H))) xlabel(39。頻率 /Hz39。)。ylabel(39。相位 /^o39。)。grid on。 y4=fftfilt(hn,x3)。 %采用 fftfilt 對輸入信號濾波 %sound(y4,fs,bits)。 figure(5) subplot(2,1,1), plot(t,x3),title(39。輸入信號 39。)%繪出輸入信號波形 subplot(2,1,2),plot(t,y4) %繪出輸出信號波形 0 500 1000 1500 2021 2500 3000 3500 4000 1 0 0 5 0050頻率 / H z振幅/dB0 500 1000 1500 2021 2500 3000 3500 4000 1 0 0 0 8 0 0 6 0 0 4 0 0 2 0 00頻率 / H z相位/o 圖 6 矩形窗低通濾波器 10 0 2 4 6 8 10 1221012輸入信號0 2 4 6 8 10 12 0 . 4 0 . 200 . 20 . 4輸出信號時間 /s 圖 7 矩形窗低通濾波結(jié)果 矩形帶通濾波器設(shè)計及濾波結(jié)果 wls=*pi。wlp=*pi。whs=*pi。whp=*pi。 delta_w=whpwhs。 wp=[ ]。N=48。 %通帶邊界頻率（歸一化頻率）和濾波器階數(shù) Fs=8000。 b=fir1(N,wp)。 %設(shè)計 FIR 帶通濾波器 figure(6) [H,f]=freqz(b,1,512,Fs)。 %以 50Hz 為采樣頻率求出濾波器頻率響應(yīng) subplot(2,1,1),plot(f,20*log10(abs(H))) xlabel(39。頻率 /Hz39。)。ylabel(39。振幅 /dB39。)。grid on。 subplot(2,1,2),plot(f,180/pi*unwrap(angle(H))) xlabel(39。頻率 /Hz39。)。ylabel(39。相位 /^o39。)。grid on。 y5=fftfilt(b,x4)。 %采用 fftfilt 對輸入信號濾波 %sound(y5,fs,bits)。 figure(7) 11 subplot(2,1,1), plot(t,x4),title(39。輸入信號 39。)%繪出輸入信號波形 subplot(2,1,2),plot(t,y5) %繪出輸出信號波形 hold on。plot(N/2**ones(1,2),ylim, 39。r39。) %繪制延遲到的時刻 title(39。輸出信號 39。),xlabel(39。時間 /s39。) 0 500 1000 1500 2021 2500 3000 3500 4000 1 5 0 1 0 0 5 0050頻率 / H z振幅/dB0 500 1000 1500 2021 2500 3000 3500 4000 2 0 0 0 1 5 0 0 1 0 0 0 5 0 00500頻率 / H z相位/o 圖 8 矩形窗帶通濾波器設(shè)計 12 0 2 4 6 8 10 1221012輸入信號0 2 4 6 8 10 12 0 . 4 0 . 200 . 20 . 4輸出信號時間 /s 圖 9 矩形窗帶通濾波器濾波結(jié)果 哈明窗低通濾波器設(shè)計及濾波結(jié)果 %哈明窗 %低通 wp=*pi。ws=*pi。 Bt=wswp。 N=ceil(*pi/Bt)。 wc=(wp+ws)/2/pi。 hn=fir1(N,wc,hamming(N+1))。 [H,f]=freqz(hn,1,512,8000)。 %采用 8000Hz 的采樣頻率求出頻率響應(yīng) figure(8) subplot(2,1,1),plot(f,20*log10(abs(H))) xlabel(39。頻率 /Hz39。)。ylabel(39。振幅 /dB39。)。grid on。 subplot(2,1,2),plot(f,180/pi*unwrap(angle(H))) xlabel(39。頻率 /Hz39。)。ylabel(39。相位 /^o39。)。grid o
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號-1