【正文】 。ws1=tan(pi*Fs1/Ft)。Fp1=1200。 [num11,den11]=bilinear(b11,a11,)。fp=2*Fp*tan(wp/2)。因為你們的幫助才使得我的這次課程設計得意圓滿完成，并從中學到了很多知識，不僅鞏固了已有知識,而且強化了一些關(guān)于數(shù)學MATLAB軟件應用方面的知識，非常感謝你們的幫助。以上的內(nèi)容全都要通過程序的編寫來完成，不過主要的操作是對函數(shù)的調(diào)用，這就要我們熟悉MATLAB的一些函數(shù)的基本功能以及它們的基本語句。6 致謝通過這次課設，我認識到了MATLAB功能非常的強大，使得我們在使用的時候用戶直接調(diào)用這些庫函數(shù)并賦予實際參數(shù)就能解決實際問題，具有極高的變成效率。我們可以像給一般信號做頻譜分析一樣，來對語音信號做頻譜分析，也可以較容易的用數(shù)字濾波器來對語音進行濾波處理。 系統(tǒng)運行界面5 結(jié)論語音信號處理是語音學與數(shù)字信號處理技術(shù)相結(jié)合的交叉學科，課題在這里不討論語音學，而是將語音當做一種特殊的信號，即一種“復雜向量”來看待。GUIDE將設計好的GUI保存在一個FIG文件中，同時生成M文件框架。)。plot(x)。axes()。m13=fft(z13)。 [b13,a13]=cheby1(n13,1,wn13,39。Fs=3500。濾波后信號波形39。axes()。)。 z12=filter(num12,den12,x)。s39。 wp2=tan(pi*Fp2/Ft)。Ft=8000。axes()。title(39。)。)。fp=2*Fp*tan(wp/2)。)。 xlabel(39。x3=sin(2*pi*t*f3)。合成信號39。 t=(1:100)/2000。并在濾波前后用sound函數(shù)播放語音信號，觀察濾波前后的變化。 [h,w]=freqz(bd,ad)。 Wc =2*fs*tan(wc/2)。IIR高通濾波器幅頻特性39。)。 Wc = 2*fs*tan(wc/2)。 IIR低通濾波器的幅頻特性讀入信號對信號傅里葉變換進行頻譜分析，利用filter函數(shù)濾波，對得到的信號傅里葉變換進行頻譜分析，與濾波前進行比較。 [bd,ad]=bilinear(b,a,fs)。 Wb=2*fs*tan(wb/2)。 利用雙線性變換法設計IIR數(shù)字濾波器變換原理：雙線性變換法是使數(shù)字濾波器的頻率響應與模擬濾波器的頻率響應相似的一種變換方法。MATLAB是由美國Mathworks公司推出的用于數(shù)值計算和圖形處理的科學計算系統(tǒng)環(huán)境。4．書寫設計說明書。唐 山 學 院 數(shù)字信號處理 課 程 設 計 題 目 基于MATLAB的數(shù)字濾波器設計 系 (部) 信息工程系 班 級 10通信本2 姓 名 學 號 指導教師 2013 年 1 月 7 日至 1 月 11 日 共 1 周 數(shù)字信號處理 課程設計任務書一、設計題目、內(nèi)容及要求設計題目：基于MATLAB的數(shù)字濾波器設計設計內(nèi)容：所設計的數(shù)字濾波器應完成以下功能：（1）設計低通、帶通、高通數(shù)字濾波器；（2）可以對合成信號（含低頻、中頻、高頻分量）、語音信號進行濾波；（3）通過GUI界面進行控制。（2）運行結(jié)果圖。例如：帶通濾波器用作頻譜分析儀中的選頻裝置；低通濾波器用作數(shù)字信號分析系統(tǒng)中的抗頻混濾波；高通濾波器被用于聲發(fā)射檢測儀中剔除低頻干擾噪聲；帶阻濾波器用作電渦流測振儀中的陷波器，等等。它們是用單位采樣響應h(n)的特性來區(qū)分的，IIR濾波器的h(n)是無限長序列，而FIR濾波器的h(n)是有限長序列。 wc=fc*2*pi/fs。)。)。Wb = 2*fs*tan(wb/2)。s39。title(39。 Wb =2*fs*tan(wb/2)。 [bd,ad]=bilinear(b,a,fs)。 IIR帶通濾波器的幅頻特性讀入信號對信號傅里葉變換進行頻譜分析，利用filter函數(shù)濾波，對得到的信號傅里葉變換進行頻譜分析，與濾波前進行比較。f3=500。title(39。 x2=sin(2*pi*t*f2)。plot(w,abs([X(1:256)]))。合成信號頻譜39。 ws=2*pi*Fs/Ft。s39。g39。)。)。這說明濾波器設計基本符合指標要求。wp1=tan(pi*Fp1/Ft)。[n12,wn12]=buttord(wp,ws,1,50,39。[h,w]=freqz(num12,den12)。濾波前信號頻譜39。)。title(39。Fp=4000。)。sound(z13)。)。axes()。濾波后信號波形39。GUIDE主要是一個界面設計工具集，MAYLAB將所有GUI支持度控件都集成在這個環(huán)境中，并提供界面外觀、屬性和行為響應方式的設置方法。最后，單擊相應的按鈕，便跳出對應的界面，進行操作。那么，就可以完全利用數(shù)字信號處理的知識來解決語音及加噪處理問題。設計出了一個簡易的用戶應用界面，可以讓人實現(xiàn)界面操作，更加方便的進行語音的頻譜分析與加噪濾波處理。首先要做的工作就是如何把設計條件寫入到MATLAB中，其次就是對條件進行分析，主要是得出一些參數(shù)。最后,我要在這里對給我這次課程設計提供幫助的老師還有同學們說聲誠摯的：謝謝。ws=2*pi*Fs/Ft。)。帶通Ft=8000。 wp2=tan(pi*Fp2/Ft)。s39。axes()。Fs=3500。 %求模擬的低通濾波器階數(shù)和截止頻率[b13,a13]=cheby1(n13,1,wn13,39。title(39。)。)。plot(abs(Y))。 %繪制x(t)的圖形x2=sin(2*pi*t*f2)。合成信號頻譜f1=20。n=[1:100]。)。Fp=1000。[n11,wn11]=buttord(wp,ws,1,50,39。 [y,fs,nbits]=wavread (39。m11=fft(z11)。axes()。plot(s)。)。wp1=tan(pi*Fp1/Ft)。[n12,wn12]=buttord(wp,ws,1,50,39。[h,w]=freqz(num12,den12)。sound(z12)。)。axes()。濾波后的信號波形39。ws=wp1*wp/ws1。 %將S域低通參數(shù)轉(zhuǎn)為高通的[num13,den13]=bilinear(num,den,)。 %傅里葉變換z13=filter(num13,den13,s)。濾波前信號的頻譜39。)。title(39。 %繪制x(t)的圖形x2=sin(2*pi*t*f2)。Fp=1000。[n11,wn11]=buttord(wp,ws,1,50,39。 z11=filter(num11,den11,x)。)。axes()。濾波后的信號波形39。x3=sin(2*pi*t*f3)。Fp2=3000。ws2=tan(pi*Fs2/Ft)。 [b12,a12]=butter(n12,wn12,39。 axes()。plot(abs(m12),39。 title(39。合成信號高通濾波f1=20。n=[1:100]。 %高通到低通濾波器參數(shù)轉(zhuǎn)換ws1=tan(pi*Fs/Ft)。)。plot(abs(x),39。r39。濾波前信號的波