【正文】 22 評分表 課程設(shè)計評分表 課程名稱： 數(shù)字信號處理 項(xiàng) 目 評 價 設(shè)計方案的合理性與創(chuàng)造性 設(shè)計與調(diào)試結(jié)果 設(shè)計說明書的質(zhì)量 答辯陳述與回答問題情況 課程設(shè)計周表現(xiàn)情況 綜合成績 教師簽名： 日 期： 23 （注： 1．此頁附在課程設(shè)計報告 最后一頁 ； 2．綜合成績按優(yōu)、良、中、及格和不及格五級評定。在今后的學(xué)習(xí)中還要繼續(xù)努力，認(rèn)真聽講，鞏固自己對于數(shù)字信號處理的知識。通過本次課程設(shè)計使我感受到在了解課本知識的前提下，同時還需要查找資料，這樣才能更好的應(yīng)用這個 MATLAB 這個工具； 在設(shè)計的過程中我還是遇到了不少的問題，比如語音信號的采集，還有就是對于濾波器的濾波原理不是很了解，通過翻書和查找相應(yīng)的資料，認(rèn)真看了各種濾波器，這才是我明白了大多數(shù)濾波器時如何工作的，不在單單只是懂理論，理論與實(shí)際相結(jié)合是很重要的 只有理論知識是遠(yuǎn)遠(yuǎn)不夠的，只有把所學(xué)的理論知識與實(shí)踐相結(jié)合，從理論中得出結(jié)論。由此可見，語音信號主要分布在低頻段，而噪聲主要分布在高頻段。其實(shí)，這正是干擾所造成的。 通過比較加噪前后，語音的頻譜和語音回放，能明顯的感覺到 加入噪聲后回放的聲音與原始的語音信號有很大的不同 ,前者隨較尖銳的干擾嘯叫聲。那么，就 可以完全利用數(shù)字信號處理的知識來解決 語音及加噪處理問題 。通過 MATLAB 里幾個命令函數(shù)的調(diào)用，很輕易的在實(shí)際語音與數(shù)字信號的理論之間搭了一座橋。作為存儲于計算機(jī)中的語音信號，其本身就是離散化了的向量，我們只需將這些離散的量提取出來，就可以對其進(jìn)行處理了。也就是說，課題更多的還是體現(xiàn)了數(shù)字信號處理技術(shù)?？梢圆煊X濾波前后的聲音有明顯的變化。 sound(yt)。其調(diào)用格式： sound(y,fs,nbits)。)。ylabel(39。f(Hz)39。) axis([0,Fs/2,0,])。title(39。plot(f,abs(Y)/max(abs(Y)))。)。ylabel(39。f(Hz)39。) axis([0,Fs/2,0,])。title(39。plot(f,abs(y)/max(abs(y)))。)。ylabel(39。t(s)39。(a)去噪后的時域波形 39。grid。 subplot(222)。yt39。)。) xlabel(39。title(39。plot(t,y)。figure(2)。 yt=filter(hn,1,Si)。f=k/Tp。 t=0:T:(M1)*T。T=1/Fs。)。ylabel(39。f(Hz)39。) axis([0,Fs/2,0,])。title(39。plot(f,abs(Y)/max(abs(Y)))。)。ylabel(39。f(Hz)39。) axis([0,Fs/2,0,])。title(39。plot(f,abs(y)/max(abs(y)))。)。ylabel(39。t(s)39。(a)去噪后的時域波形 39。grid。 subplot(222)。yt39。)。) xlabel(39。title(39。plot(t,y)。 sound(yt)。Y=fft(yt,M)。f=k/Tp。 t=0:T:(M1)*T。T=1/Fs。)。 title(39。 subplot(2,2,4)。)。 title(39。 subplot(2,2,3)。)。 title(39。 subplot(2,2,2)。)。 title(39。g39。 subplot(2,2,1)。 sound(z11)。 S=fft(Si)。 n = length (y) 。39。分貝數(shù) 39。)。xlabel(39。(b)幅度響應(yīng) 39。plot(f,20*log10(abs(H)/max(abs(H))))。)。ylabel(39。n39。) axis([0,N1,min(hn),max(hn)])。title(39。 figure(1) subplot(211)。 hn=remez(n,fo,mo,W)。 dev=[(10^(Rp/20)1)/(10^(Rp/20)+1) 10^(Rs/20),]。 fb=[120,140]。Rp=1。 fp=150。39。利用等波紋最佳逼近準(zhǔn)則設(shè)計線性相位 FIR 數(shù)字濾波器數(shù)學(xué)模型的檢錄及其求解算法的推導(dǎo)復(fù)雜，求解計算必須借助計算機(jī)，幸好濾波器設(shè)計專家已經(jīng)開發(fā)出 MATLAB 工具箱函數(shù) remez，只要簡單的調(diào)用這個函數(shù)就可以完成線性相位 FIR 數(shù)字濾波器的 等波紋最佳逼近法。所以 W(w)取值越大的頻段逼近精度越高，開始設(shè)計時應(yīng)該根據(jù)逼近精度要求確定 W(w)，在 Remez多重交換迭代中 W(w)是確知函數(shù)。 W(w)稱為誤差加權(quán)函數(shù)，用來控制不同頻段（一般指通帶和阻帶）的逼近精度， W(w)和 N 由濾波器設(shè)計指標(biāo)（ Rp、 As 以及過渡帶寬度）確定。 設(shè)希望設(shè)計的濾波器的幅頻響應(yīng)是 Hd(w),實(shí)際設(shè)計的濾波器的幅頻響應(yīng)是 Hg(w)。這種設(shè)計法使最大誤差均勻分布，所以設(shè)計的濾波器性價比最高。) 圖 26窗函數(shù)設(shè)計濾波器 (2)最優(yōu)設(shè)計法 等波紋逼近法師一種優(yōu)化設(shè)計法，它克服了窗函數(shù)設(shè)計法和頻率采樣法的缺點(diǎn)， 使最大誤差最小化，并是誤差在整個逼近頻段上均勻分布。ylabel(39。f(Hz)39。)。grid title(39。 subplot(212)。h(n)39。)。xlabel(39。(a) h(n)波形 39。stem(hn)。 hn=fir1(N1,wc,blackman(N))。 Bt=wswp。 wp=*pi。Rp=1。 fp=120。39。當(dāng) hd(n)的對稱中心點(diǎn)取值為（ N1） /2 時，截取的一段（ N1）/2 對稱，這樣就可以保證所涉及的濾波器具有線性相位。設(shè)截取的一段為 h(n)，則： h(n)=hd(n)*wr(n) 其中 wR(n)={1,0=n=N1。但是通常 Hd(e^jw)在邊界頻率處有不連續(xù)點(diǎn)，這似的 hd(n)是無限長的非因果序列，所以實(shí)際是不 能實(shí)現(xiàn)的。一般情況如下。)。 legend(39。 [h,w]=freqz(Bz,Az)。)。 [B,A]=butter(N,wc,39。s39。Rs=60。 ws=2*Fs*tan(wsz*pi/2)。 wsz=。)。 原程序如下 [y,Fs,nbits]=wavread (39。 Ws表示阻帶截止頻率， Rp表示通帶最大衰減， Rs 表示阻帶最小衰減， s表示設(shè)計模擬濾波器。)。該函數(shù)的調(diào)用格式為 [N,wc]=buttord(wp,ws,Rp,Rs,39。 和 IIR 數(shù)字濾波器的主要程序 ① IIR 數(shù)字濾波器： (1)流程圖 13 (2)巴特沃斯低通濾波器 MATLAB 提供了 buttap 函數(shù)來設(shè)計 N 階歸一化巴特沃斯模擬原型，濾波器，調(diào)用格式為