【正文】 %繪制x(t)的圖形x2=sin(2*pi*t*f2)。Fp=4000。)。sound(z13)。)。axes()。濾波后的信號波形39。plot(z13)。濾波后信號的頻譜39。title(39。 %利用雙線性變換實現(xiàn)S域到Z域轉(zhuǎn)換[h,w]=freqz(num13,den13)。[n13,wn13]=cheb1ord(wp,ws,1,50,39。x=sin(2*pi*t*f1)+sin(2*pi*t*f2)+sin(2*pi*t*f3)。t=(1:100)/2000。axes()。title(39。)。 [num2,den2]=lp2bp(b12,a12,sqrt(wp1*wp2),bw)。wp=1。As=100。t=n/2000X=fft(x,512)。f2=200。濾波前信號的波形39。r39。plot(abs(x),39。 [b11,a11]=butter(n11,wn11,39。Ap=1。n=[1:100]。合成信號低通濾波f1=20。 title(39。plot(abs(m13),39。 %求濾波后的信號axes()。C:\program files\matlab71\work\speech\39。)。Fp=4000。濾波前信號的波形39。r39。plot(abs(S),39。)。 [b12,a12]=butter(n12,wn12,39。ws2=tan(pi*Fs2/Ft)。Fp2=3000。)。)。g39。n = length (y) 。 [b11,a11]=butter(n11,wn11,39。Ap=1。)。x=sin(2*pi*t*f1)+sin(2*pi*t*f2)+sin(2*pi*t*f3)。t=(1:100)/2000。axes()。)合成信號f1=20。)。 %求出語音信號的長度Y=fft(y,N)。legend(39。 %求S域的頻率響應(yīng)的參數(shù)[num,den]=lp2hp(b13,a13,wn13)。wp=1。legend(39。s39。w=wp1*wp2/ws2。Fs1=1000。 axes()。[n11,wn11]=buttord(wp,ws,1,50,39。Fp=1000。通過設(shè)計IIR和FIR的各種濾波器，并對聲音信號進(jìn)行濾波處理，加深對DSP數(shù)字信號處理的理解和對MATLAB的進(jìn)一步學(xué)習(xí)，將理論轉(zhuǎn)化成實踐，并在實踐中檢驗理論。本文次課設(shè)通過一個設(shè)計實例，利用MATLAB實現(xiàn) IIR數(shù)字濾波器設(shè)計與濾波的三種方法，從仿真結(jié)果可以看出它們均可以達(dá)到技術(shù)指標(biāo)要求，而且方法簡單、快捷，大大減輕了工作量。從含噪語音信號的頻譜圖中可以看出含噪聲的語音信號頻譜,在整個頻域范圍內(nèi)分是布均勻。從課題的中心來看，課題“基于MATLAB的有噪聲語音信號處理”是希望將數(shù)字信號處理技術(shù)應(yīng)用于某一實際領(lǐng)域，這里就是指對語音及加噪處理。M文件包括GUI設(shè)計、控件函數(shù)以及定義為子函數(shù)的用戶控件回調(diào)函數(shù)，主要用于控制GUI展開時的各種特征。：分析濾波前后信號的變化，波形圖變窄，頻譜變化表明：信號的高頻率段被保留，低頻率段被濾除。濾波前信號波形39。r39。plot(abs(x),39。)。 ws1=tan(pi*Fs/Ft)。讀入信號對信號傅里葉變換進(jìn)行頻譜分析，利用filter函數(shù)濾波，對得到的信號傅里葉變換進(jìn)行頻譜分析，與濾波前進(jìn)行比較。 title(39。plot(abs(m12),39。 axes()。 [b12,a12]=butter(n12,wn12,39。ws2=tan(pi*Fs2/Ft)。Fp2=3000。title(39。)。濾波前信號頻譜39。 [h,w]=freqz(num11,den11)。 fs=2*Fs*tan(ws/2)。 Fp=1000。)。n=[1:100]。 合成信號波形 合成信號的頻譜程序設(shè)計f1=20。 x2=sin(2*pi*t*f2)。頻譜變化表明：信號的低頻率段和高頻率段被濾除，中間頻率段被保留；分析濾波前后的聲音變化，濾波后聲音一定程度上變的尖銳，與高通濾波器濾波后的聲音相比較低。plot(w*fs/(2*pi),20*log10(abs(h)))。s39。并在濾波前后用sound函數(shù)播放語音信號，觀察濾波前后的變化。 [h,w]=freqz(bd,ad,1024)。s39。 分析濾波前后信號的變化，波形圖變窄，頻譜變化表明:信號的低頻率段被保留，高頻率段被濾除。 figure(1)。[n,wn]=buttord(Wb,Wc,Ap,As,39。IIR濾波器的設(shè)計原理IIR數(shù)字濾波器的設(shè)計一般是利用目前已經(jīng)很成熟的模擬濾波器的設(shè)計方法來進(jìn)行設(shè)計，通常采用模擬濾波器原型有butterworth函數(shù)、chebyshev函數(shù)、bessel函數(shù)、橢圓濾波器函數(shù)等。在MATLAB環(huán)境下，用戶可以集成地進(jìn)行程序設(shè)計，數(shù)值計算，圖形繪制，輸入輸出，文件管理等各項操作。三、要求的設(shè)計成果（課程設(shè)計說明書、設(shè)計實物、圖紙等） 設(shè)計結(jié)果能正確仿真演示設(shè)計說明書一份（包括總體設(shè)計、算法原理圖及說明、系統(tǒng)GUI演示、源程序清單等）四、進(jìn)程安排 周一：資料收集周二：利用MATLAB完成GUI界面繪制周三：程序調(diào)試 周四：書寫課程設(shè)計說明書周五：答辯五、主要參考資料樓順天,—子科技大學(xué)出版社，1998，2000宗孔德,，1997 指導(dǎo)教師（簽名）：教研室主任（簽名）：課程設(shè)計成績評定表出勤情況出勤天數(shù) 缺勤天數(shù)成績評定出勤情況及設(shè)計過程表現(xiàn)（20分）課設(shè)答辯（20分）設(shè)計成果（60分）總成績（100分）提問（答辯）問題情況綜合評定 指導(dǎo)教師簽名： 年 月 日目錄1 緒論 12數(shù)字濾波器的設(shè)計 2 數(shù)字濾波器的基本概念 2 利用雙線性變換法設(shè)計IIR數(shù)字濾波器 23 課程設(shè)計的具體實現(xiàn) 3 3 IIR數(shù)字濾波器的流程圖 3 IIR低通濾波器設(shè)計 3 IIR高通濾波器設(shè)計 5 IIR帶通濾波器設(shè)計 7 對合成信號進(jìn)行濾波的濾波器設(shè)計 8 8 8 9 10 IIR帶通濾波器設(shè)計 11 IIR高通濾波器濾波 134軟件簡介 15—GUI概述 15 155 結(jié)論 176 致謝 18參考文獻(xiàn) 19附錄 20課程設(shè)計說明書1 緒論與模擬濾波器相對應(yīng)，在離散系統(tǒng)中廣泛應(yīng)用數(shù)字濾波器。 2. 完成數(shù)字濾波器具體設(shè)計。3. 濾波器程序的設(shè)計。數(shù)字濾波器一般可以用兩種方法來實現(xiàn)：一種方法是用數(shù)字硬件裝配成一臺專門的設(shè)備，這種設(shè)備稱為數(shù)字信號處理機；另一種方法就是直接利用通用計算機，將所需要的運算編成程序讓通用計算機來完成，即利用計算機軟件來實現(xiàn)。用戶可以根據(jù)自己的工作任務(wù)，開發(fā)自己的工具箱。3 課程設(shè)計的具體實現(xiàn) IIR數(shù)字濾波器的流程圖開始讀入數(shù)字濾波器技術(shù)指標(biāo)將指標(biāo)轉(zhuǎn)換成歸一化模擬濾波器的指標(biāo)計算歸一化的模擬濾波器階數(shù)N和3db截止頻率模擬域頻率變換，將信號變換成模擬濾波器H（s）用雙線性變換法將H(s)轉(zhuǎn)換成數(shù)字濾波器H(z)輸入信號后顯示相應(yīng)結(jié)果結(jié)束 IIR數(shù)字濾波器流程圖 IIR低通濾波器設(shè)計根據(jù)數(shù)字濾波器的設(shè)計原理，首先將數(shù)字域的指標(biāo)轉(zhuǎn)化為模擬域的指標(biāo)設(shè)計模擬低通濾波器，然后應(yīng)用雙線性變換法將模擬濾波器轉(zhuǎn)化為數(shù)字濾波器將指標(biāo)轉(zhuǎn)換成歸一化模擬低通濾波器的指標(biāo)，通過歸一化的模擬低通濾波器階數(shù)N和3dB截止頻率的計算，將模擬域頻率變換成模擬低濾波器H(s),并