【正文】 wp2=tan(pi*Fp2/Ft)。Fs2=3200。語音信號帶通濾波Ft=8000。plot(z11)。 title(39。濾波后信號的頻譜39。plot(abs(m11),39。title(39。 axes()。 S=fft(s)。C:\program files\matlab71\work\speech\39。)。s39。ws=2*pi*Fs/Ft。Fs=1200。合成信號頻譜圖39。ylabel(39。plot(w,abs([X(1:256)]))。t=n/2000X=fft(x,512)。 %繪制x(t)的圖形x2=sin(2*pi*t*f2)。f2=200。title(39。x3=sin(2*pi*t*f3)。f3=500。title(39。 %回放語音信號N= length (y) 。語音信號頻譜[y,fs,nbits]=wavread (39。plot(y)。sound(y,fs,nbits)。)。IIR高通濾波器39。 %利用雙線性變換實現(xiàn)S域到Z域轉(zhuǎn)換[h,w]=freqz(num13,den13)。s39。[n13,wn13]=cheb1ord(wp,ws,1,50,39。wp1=tan(pi*Fp/Ft)。)。plot(w*8000*,abs(h))。 [num2,den2]=lp2bp(b12,a12,sqrt(wp1*wp2),bw)。)。wp=1。ws1=tan(pi*Fs1/Ft)。As=100。Fp1=1200。legend(39。 [num11,den11]=bilinear(b11,a11,)。)。fp=2*Fp*tan(wp/2)。As=100 。因為你們的幫助才使得我的這次課程設(shè)計得意圓滿完成，并從中學(xué)到了很多知識，不僅鞏固了已有知識,而且強化了一些關(guān)于數(shù)學(xué)MATLAB軟件應(yīng)用方面的知識，非常感謝你們的幫助。當(dāng)然，經(jīng)過了課程設(shè)計，我也發(fā)現(xiàn)了自己的很多不足。以上的內(nèi)容全都要通過程序的編寫來完成，不過主要的操作是對函數(shù)的調(diào)用，這就要我們熟悉MATLAB的一些函數(shù)的基本功能以及它們的基本語句。由于MATLAB具有強大的接口功能，仿真后的結(jié)果可以很方便的移植到 DSP、CPLD。6 致謝通過這次課設(shè)，我認識到了MATLAB功能非常的強大，使得我們在使用的時候用戶直接調(diào)用這些庫函數(shù)并賦予實際參數(shù)就能解決實際問題，具有極高的變成效率。通過濾波前后的對比，低通濾波后效果最好，高通濾波后的效果最差。我們可以像給一般信號做頻譜分析一樣，來對語音信號做頻譜分析，也可以較容易的用數(shù)字濾波器來對語音進行濾波處理。這一過程的實現(xiàn)，用到了處理數(shù)字信號的強有力工具MATLAB。 系統(tǒng)運行界面5 結(jié)論語音信號處理是語音學(xué)與數(shù)字信號處理技術(shù)相結(jié)合的交叉學(xué)科，課題在這里不討論語音學(xué)，而是將語音當(dāng)做一種特殊的信號，即一種“復(fù)雜向量”來看待。第一步：通過設(shè)置GUIDE應(yīng)用程序的選項來運行GUIDE；第二步：使用界面設(shè)計編輯器進行面設(shè)計；第三步：編寫控件行為響應(yīng)控制（即回調(diào)函數(shù)）代碼。GUIDE將設(shè)計好的GUI保存在一個FIG文件中，同時生成M文件框架。4軟件簡介—GUI概述圖形用戶界面(graphical user interfaces)，GUI則是由窗口、光標、按鍵、菜單、文字說明等對象(objects)構(gòu)成的一個用戶界面。)。axes()。plot(x)。title(39。axes()。)。m13=fft(z13)。 [num13,den13]=bilinear(num,den,)。 [b13,a13]=cheby1(n13,1,wn13,39。ws=wp1*wp/ws1。Fs=3500。頻譜變化表明：信號的低頻率段和高頻率段被濾除，中間頻率段被保留。濾波后信號波形39。)。axes()。)。)。g39。 z12=filter(num12,den12,x)。)。s39。bw=wp2wp1。 wp2=tan(pi*Fp2/Ft)。Fs2=3200。Ft=8000。)。axes()。plot(x)。title(39。axes()。)。m11=fft(z11)。)。s39。fp=2*Fp*tan(wp/2)。 As=100 。)。頻率響應(yīng)幅度39。 xlabel(39。w=(0:255)/256*1000。x3=sin(2*pi*t*f3)。f3=500。合成信號39。 x=sin(2*pi*t*f1)+sin(2*pi*t*f2)+sin(2*pi*t*f3)。 t=(1:100)/2000。 對合成信號進行濾波的濾波器設(shè)計產(chǎn)生一個連續(xù)信號，包含低頻、中頻、高頻分量，對其進行采樣，進行頻譜分析，分別設(shè)計三種高通、低通、帶通濾波器對信號進行濾波處理，觀察濾波后信號的頻譜。并在濾波前后用sound函數(shù)播放語音信號，觀察濾波前后的變化。 title(39。 [h,w]=freqz(bd,ad)。 [b,a]=butter(n,wn,39。 Wc =2*fs*tan(wc/2)。分析濾波前后的聲音變化，濾波后聲音明顯變高而尖銳，這說明濾波器設(shè)計基本符合指標要求。IIR高通濾波器幅頻特性39。subplot(1,1,1)。)。 [b,a]=butter(n,wn,39。 Wc = 2*fs*tan(wc/2)。將指標轉(zhuǎn)換成歸一化模擬低通濾波器的指標，通過歸一化的模擬低通濾波器階數(shù)N和3dB截止頻率的計算，將模擬域頻率變換成模擬低濾波器H(s),并用雙線性變換法將H(s)轉(zhuǎn)換成數(shù)字低通濾波器H(z)，由此得到低通數(shù)字濾波器。 IIR低通濾波器的幅頻特性讀入信號對信號傅里葉變換進行頻譜分析，利用filter函數(shù)濾波，對得到的信號傅里葉變換進行頻譜分析，與濾波前進行比較。plot(w*fs/(2*pi),20*log10(abs(h)))。 [bd,ad]=bilinear(b,a,fs)。)。 Wb=2*fs*tan(wb/2)。本課程設(shè)計設(shè)計思想：首先利用MATLAB分別產(chǎn)生低頻中頻高頻信號，然后進行疊加得到連續(xù)時間信號；對所產(chǎn)生的連續(xù)時間信號進行采樣，得到數(shù)字信號；對信號進行FFT頻譜分析，繪制其頻譜圖；根據(jù)信號頻譜分析的結(jié)果，分別設(shè)計高通，低通，帶通濾波器，得到濾波器的幅頻及相頻特性。 利用雙線性變換法設(shè)計IIR數(shù)字濾波器變換原理：雙線性變換法是使數(shù)字濾波器的頻率響應(yīng)與模擬濾波器的頻率響應(yīng)相似的一種變換方法。目前，Mathworks公司推出了18種工具箱。MATLAB是由美國Mathworks公司推出的用于數(shù)值計算和圖形處理的科學(xué)計算系統(tǒng)環(huán)境?；蛘哒f，把輸入信號變成一定的輸出信號，從而達到改變信號頻譜的目的。4．書寫設(shè)計說明書。（2）確定設(shè)計方法、設(shè)計指標。唐 山 學(xué) 院 數(shù)字信號處理 課 程 設(shè) 計 題 目 基于MATLAB的數(shù)字濾波器設(shè)計 系 (部) 信息工程系 班 級 10通信本2 姓 名 學(xué) 號 指導(dǎo)教師 2013 年 1 月 7 日至 1 月 11 日 共 1 周 數(shù)字信號處理 課程設(shè)計任務(wù)書一、設(shè)計題目、內(nèi)容及要求設(shè)計題目：基于MATLAB的數(shù)字濾波器設(shè)計設(shè)計內(nèi)容：所設(shè)計的數(shù)字濾波器應(yīng)完成以下功能：（1）設(shè)計低通、帶通、高通數(shù)字濾波器；（2）可以對合成信號（含低頻、中頻、高頻分量）、語音信號進行濾波；（3）通過GUI界面進行控制。（1）輸入信號的選定。（2）運行結(jié)果圖。它的作用是利用離散時間系統(tǒng)的特性對輸入信號波形或頻率進行加工處理。例如：帶通濾波器用作頻譜分析儀中的選頻裝置；低通濾波器用作數(shù)字信號分析系統(tǒng)中的