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