【正文】 課程設(shè)計報告 任健院 系 城南學(xué)院 專 業(yè) 通信工程 班 級 通信1104班 學(xué) 號 201185250426 學(xué)生姓名 任 健 指導(dǎo)教師 熊文杰 課程成績 完成日期 2014年7月4日 課程設(shè)計成績評定院 系 城南學(xué)院 專 業(yè) 通信工程 班 級 通信1104班 學(xué) 號 201185250426 學(xué)生姓名 任 健 指導(dǎo)教師 熊文杰 完成日期 2014 年 7月 4 日 指導(dǎo)教師對學(xué)生在課程設(shè)計中的評價評分項目優(yōu)良中及格不及格課程設(shè)計中的創(chuàng)造性成果學(xué)生掌握課程內(nèi)容的程度課程設(shè)計完成情況課程設(shè)計動手能力文字表達學(xué)習(xí)態(tài)度規(guī)范要求課程設(shè)計論文的質(zhì)量指導(dǎo)教師對課程設(shè)計的評定意見綜合成績 指導(dǎo)教師簽字 2014年7月 4日課程設(shè)計任務(wù)書 城南 學(xué)院 通信工程 專業(yè) 課程名稱數(shù)字信號處理課程設(shè)計時間2013～2014學(xué)年第二學(xué)期18～19周學(xué)生姓名任健指導(dǎo)老師熊文杰題 目彈撥音樂濾波去噪——使用三角窗設(shè)計的FIR濾波器主要內(nèi)容： 從網(wǎng)站上下載一段彈撥音樂，播放直觀感覺聲音特點，繪制波形并觀察其頻譜特點。（3）處理結(jié)果和分析結(jié)論應(yīng)該一致，而且應(yīng)符合理論。數(shù)字信號處理與模擬信號處理是信號處理的子集。從實現(xiàn)方法上考慮，將濾波器分成兩種，一種稱為無限脈沖響應(yīng)濾波器，簡稱IIR（Infinite Impulse Response）濾波器，另一種稱為FIR（Finite Impulse Response）濾波器[2]。MATLAB是矩陣實驗室（Ma trix Laboratory）的簡稱，是美國MathWorks公司出品的商業(yè)數(shù)學(xué)軟件，用于算法開發(fā)，數(shù)學(xué)可視化，數(shù)據(jù)分析以及數(shù)值計算的高級技術(shù)計算語言和交互式環(huán)境，主要包括MATLAB和Simulink兩大部分。而使數(shù)字信號處理從理論走向?qū)嵱玫氖强焖俑盗⑷~變換(FFT)，F(xiàn)FT的出現(xiàn)大大減少了DFT的運算量，使實時的數(shù)字信號處理成為可能、極大促進了該學(xué)科的發(fā)展。二維(2D)信號：兩個自變量的函數(shù)。窗口設(shè)計法就是截取為有限長的一段因果序列，并用合適的窗口函數(shù)進行加權(quán)作為FIR的單位脈沖響應(yīng) 。而且，就是同一公司的不同系統(tǒng)的DSP芯片，其編程指令也會有所不同，開發(fā)周期較長。圖 打開聲音屬性在其中點擊立即轉(zhuǎn)換，將聲音格式設(shè)置為PCM，8位，單聲道，7KB/秒，然后點擊確定。 % 按指定的采樣率和每樣本編碼位數(shù)回放所錄語音文件N=length(x)。 Y=abs(fft(y))。因為我使用的是三角窗法，在三角窗函數(shù)中，最小阻帶衰減為25db。 % 帶阻濾波器設(shè)計指標(biāo) df=min((fs1fp1),(fp2fs2))。M=ceil(*pi/dw)+1。程序如下： y_fil=filter(h,1,y)。用sound(y_fil,fs,bits)語句回放該濾波去噪信號，便可以感覺到濾波后的語音信號與原信號差不多，該設(shè)計基本符合要求。5 結(jié)束語為期兩周的課程設(shè)計終于結(jié)束了，在這次課程設(shè)計的過程中我又學(xué)到了很多，數(shù)字信號處理課程設(shè)計是數(shù)字信號處理課程的重要實踐性環(huán)節(jié)，是我們在校期間一次較全面的能力訓(xùn)練，使得我鞏固加深、綜合運用所學(xué)的專業(yè)理論知識，對通信工程領(lǐng)域各種技術(shù)的DSP實現(xiàn)的設(shè)計有了一定的了解。D：/39。+*sin(fn*2*pi*t)。 % 計算頻譜頻率范圍 subplot(2,2,1)。title(39。)。 subplot(2,2,3)。title(39。)。 fp1=fn200。wc1=(fs1+fp1)/2。 % 定義時間范圍 w_tri=(triang(M))39。)。w39。xlabel(39。stem(n,w_tri)。 % 用設(shè)計好的濾波器對加噪信號y進行濾波 Y_fil=fft(y_fil)。幅度39。)。plot(t,y)。加入單頻干擾后的語音信號39。)。)。xlabel(39。)。% db = [0 到pi弧度]區(qū)間內(nèi)的相對振幅(db)% mag = [0 到pi弧度]區(qū)間內(nèi)的絕對振幅% pha = [0 到pi弧度]區(qū)間內(nèi)的相位響應(yīng)% 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。% pha = unwrap(angle(H))。