數(shù)字濾波器及在語音信號分析中的作用_數(shù)字信號處理課程設(shè)計-資料下載頁
2025-03-04 14:56本頁面
【導讀】文檔中包括設(shè)計思路、設(shè)計仿真程序、仿真結(jié)果及相應(yīng)的分析。課程設(shè)計報告裝訂順序：封面、任務(wù)書、目錄、正文、評分、附件。正文的格式:一級標題用3號黑體,二級標題用四號宋體加粗,正文用小四號宋體;行距為22。試；五、總結(jié)；六、附件；七、評分表。使用MATLAB繪出采樣后的語音信號的時域波形和頻譜圖。根據(jù)頻譜圖求出其帶寬，設(shè)置為低頻、高頻、帶限噪聲，或Chirp干擾、沖激干擾。dB;阻帶最小衰減Rs=15dB。矩形窗和哈明窗設(shè)計合理的FIR濾波器。具體設(shè)計流程圖如圖所示：。其中噪音模塊加的有高斯白噪聲，單頻余弦信號，多頻余弦信號噪音。濾波器濾波模塊經(jīng)過多次調(diào)整，濾波效果已趨完善。的，因而濾波器是穩(wěn)定的系統(tǒng)。
【正文】 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。 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) %繪出輸出信號波形 title(39。輸出信號 39。),xlabel(39。時間 /s39。) %矩形窗 %帶通 wls=*pi。wlp=*pi。whs=*pi。whp=*pi。 19 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) 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。) %哈明窗 %低通 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 on。 y6=fftfilt(hn,x3)。 %采用 fftfilt 對輸入信號濾波 %sound(y6,fs,bits)。 figure(9) 20 subplot(2,1,1), plot(t,x3),title(39。輸入信號 39。)%繪出輸入信號波形 subplot(2,1,2),plot(t,y6) %繪出輸出信號波形 title(39。輸出信號 39。),xlabel(39。時間 /s39。) %哈明窗 %帶通 %Samp7_6 wp=[ ]。N=48。 %通帶邊界頻率（歸一化頻率）和濾波器階數(shù) Fs=8000。 b=fir1(N,wp)。 %設(shè)計 FIR 帶通濾波器 figure(10) [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。 y7=fftfilt(b,x4)。 %采用 fftfilt 對輸入信號濾波 sound(y7,fs,bits)。 figure(11) subplot(2,1,1), plot(t,x4),title(39。輸入信號 39。)%繪出輸入信號波形 subplot(2,1,2),plot(t,y7) %繪出輸出信號波形 hold on。plot(N/2**ones(1,2),ylim, 39。r39。) %繪制延遲到的時刻 title(39。輸出信號 39。),xlabel(39。時間 /s39。) 21 評分表 計算機與通信學院課程設(shè)計評分表 課題名稱： 教師簽名： 日 期 ： 項 目 評 價 設(shè)計方案的合理性與創(chuàng)造性 設(shè)計與調(diào)試結(jié)果 設(shè)計說明書的質(zhì)量 答辯陳述與回答問題情況 課程設(shè)計周表現(xiàn)情況 綜合成績
教學課件相關(guān)推薦
鄂ICP備17016276號-1