【正文】 % pha = unwrap(angle(H))。)。)。加入單頻干擾后的語(yǔ)音信號(hào)39。)。 % 用設(shè)計(jì)好的濾波器對(duì)加噪信號(hào)y進(jìn)行濾波 Y_fil=fft(y_fil)。xlabel(39。)。wc1=(fs1+fp1)/2。)。 subplot(2,2,3)。title(39。+*sin(fn*2*pi*t)。5 結(jié)束語(yǔ)為期兩周的課程設(shè)計(jì)終于結(jié)束了，在這次課程設(shè)計(jì)的過程中我又學(xué)到了很多，數(shù)字信號(hào)處理課程設(shè)計(jì)是數(shù)字信號(hào)處理課程的重要實(shí)踐性環(huán)節(jié)，是我們?cè)谛Ｆ陂g一次較全面的能力訓(xùn)練，使得我鞏固加深、綜合運(yùn)用所學(xué)的專業(yè)理論知識(shí)，對(duì)通信工程領(lǐng)域各種技術(shù)的DSP實(shí)現(xiàn)的設(shè)計(jì)有了一定的了解。程序如下： y_fil=filter(h,1,y)。 % 帶阻濾波器設(shè)計(jì)指標(biāo) df=min((fs1fp1),(fp2fs2))。 Y=abs(fft(y))。圖 打開聲音屬性在其中點(diǎn)擊立即轉(zhuǎn)換，將聲音格式設(shè)置為PCM，8位，單聲道，7KB/秒，然后點(diǎn)擊確定。窗口設(shè)計(jì)法就是截取為有限長(zhǎng)的一段因果序列，并用合適的窗口函數(shù)進(jìn)行加權(quán)作為FIR的單位脈沖響應(yīng) 。而使數(shù)字信號(hào)處理從理論走向?qū)嵱玫氖强焖俑盗⑷~變換(FFT)，F(xiàn)FT的出現(xiàn)大大減少了DFT的運(yùn)算量，使實(shí)時(shí)的數(shù)字信號(hào)處理成為可能、極大促進(jìn)了該學(xué)科的發(fā)展。從實(shí)現(xiàn)方法上考慮，將濾波器分成兩種，一種稱為無(wú)限脈沖響應(yīng)濾波器，簡(jiǎn)稱IIR（Infinite Impulse Response）濾波器，另一種稱為FIR（Finite Impulse Response）濾波器[2]。（3）處理結(jié)果和分析結(jié)論應(yīng)該一致，而且應(yīng)符合理論。（4）獨(dú)立完成課程設(shè)計(jì)并按要求編寫課程設(shè)計(jì)報(bào)告書。設(shè)計(jì)FIR數(shù)字濾波器的方法有窗函數(shù)法、頻率采樣法和等波紋最佳逼近法等。數(shù)字信號(hào)處理的特征和分類。常見的窗函數(shù)，可以分為以下主要類型：（1） 冪窗采用時(shí)間變量某種冪次的函數(shù)，如矩形、三角形、梯形或其它時(shí)間（t）的高次冪；　 （2）三角函數(shù)窗應(yīng)用三角函數(shù)，即正弦或余弦函數(shù)等組合成復(fù)合函數(shù)，例如漢寧窗、海明窗等；（3） 指數(shù)窗采用指數(shù)時(shí)間函數(shù)，例如高斯窗等。然后再點(diǎn)擊文件/另存為，即可將該語(yǔ)音文件存在MATLAB文件下的work文件夾里。 % 對(duì)原始信號(hào)和加噪信號(hào)進(jìn)行fft變換 X=X(1:N/2)。 % 計(jì)算上下邊帶中心頻率，和頻率間隔 dw=df/fs*2*pi。 % 用設(shè)計(jì)好的濾波器對(duì)加噪信號(hào)y進(jìn)行濾波 Y_fil=fft(y_fil)。通過語(yǔ)音信號(hào)濾波去噪——使用三角窗設(shè)計(jì)的FIR濾波器課程設(shè)計(jì)，我對(duì)數(shù)字信號(hào)方面的知識(shí)得到了更深一步的了解。 %在原信號(hào)上添加噪聲sound(y,fs,bits)。原始語(yǔ)音信號(hào)39。plot(t,y)。ylabel(39。wc2=(fs2+fp2)/2。ylabel(39。w39。Y_fil=Y_fil(1:N/2)。ylabel(39。) subplot(3,2,4)。ylabel(39。axis([0 6000,0 100])附錄2：ideal_lp函數(shù)源程序清單function hd = ideal_lp(wc,M)。grd = grpdelay(b,a,w)。pha = angle(H)。濾波后的語(yǔ)音信號(hào)幅度譜39。時(shí)間（單位：s）39。title(39。頻率（單位：Hz）39。 y_fil=filter(h,1,y)。plot(w/pi,pha)。w39。 % 計(jì)算上下邊帶中心頻率，和頻率間隔 dw=df/fs*2*pi。頻率（單位：Hz）39。 axis([0 ,0 2])。)。 y=x39。 確實(shí)，在本次課程設(shè)計(jì)過程中遇到了很多的困難，但是同學(xué)們的幫助讓我本次設(shè)計(jì)省了不少心，老師在大體上給我們指明方向，同學(xué)們則在我實(shí)際動(dòng)手的過程中給予指點(diǎn)和幫助，讓我受益匪淺。在Matlab中，F(xiàn)IR濾波器利用函數(shù)fftfilt對(duì)信號(hào)進(jìn)行濾波，IIR濾波器利用函數(shù)filter對(duì)信號(hào)進(jìn)行濾波。As=15。 % 應(yīng)該可以明顯聽出有尖銳的單頻嘯叫聲X=abs(fft(x))。錄一段2至3秒的語(yǔ)音，然后左鍵單擊文件/屬性。因此是無(wú)限長(zhǎng)非因果序列，不能直接作為FIR的單位脈沖響應(yīng)。 數(shù)字信號(hào)處理的核心算法是離散傅立葉變換(DFT)，是DFT使信號(hào)在數(shù)字域和頻域都實(shí)現(xiàn)了離散化，從而可以用通用計(jì)算機(jī)處理離散信號(hào)。從濾波特性方面考慮，數(shù)字濾波器可分成數(shù)字高通、數(shù)字低通、數(shù)字帶通和數(shù)字帶阻等濾波器。（2）設(shè)計(jì)完后應(yīng)檢查其頻率響應(yīng)曲線是否滿足指標(biāo)。應(yīng)當(dāng)提交的文件：（1）課程設(shè)計(jì)學(xué)年論文。 課程設(shè)計(jì)目的用麥克風(fēng)采集一段語(yǔ)音信號(hào)，繪制波形并觀察其頻譜，給定相應(yīng)技術(shù)指標(biāo)，用三角窗設(shè)計(jì)一個(gè)滿足指標(biāo)的FIR濾波器，對(duì)該語(yǔ)音信號(hào)進(jìn)行濾波去噪處理，比較濾波前后的波形和頻譜并進(jìn)行分析，根據(jù)結(jié)果和學(xué)過的理論得出合理的結(jié)論。信號(hào)(signal)是一種物理體現(xiàn)，或是傳遞信息的函數(shù)。：名稱濾波器過渡帶寬最小阻帶衰減名稱濾波器過渡帶寬最小阻帶衰減矩形21dBPARZENWIN56db巴特利特25dBFLATTOPWIN108db漢寧44dBGAUSSWIN60db漢明51dBBARTHANNWIN40db布萊克曼11π/M74dBB