【正文】 濾波前信號的波形39。plot(abs(x),39。 %高通到低通濾波器參數(shù)轉(zhuǎn)換ws1=tan(pi*Fs/Ft)。合成信號高通濾波f1=20。plot(abs(m12),39。 [b12,a12]=butter(n12,wn12,39。Fp2=3000。濾波后的信號波形39。)。[n11,wn11]=buttord(wp,ws,1,50,39。 %繪制x(t)的圖形x2=sin(2*pi*t*f2)。)。 %傅里葉變換z13=filter(num13,den13,s)。ws=wp1*wp/ws1。axes()。sound(z12)。[n12,wn12]=buttord(wp,ws,1,50,39。)。axes()。 [y,fs,nbits]=wavread (39。Fp=1000。n=[1:100]。 %繪制x(t)的圖形x2=sin(2*pi*t*f2)。)。title(39。Fs=3500。s39。帶通Ft=8000。ws=2*pi*Fs/Ft。首先要做的工作就是如何把設(shè)計條件寫入到MATLAB中，其次就是對條件進行分析，主要是得出一些參數(shù)。那么，就可以完全利用數(shù)字信號處理的知識來解決語音及加噪處理問題。GUIDE主要是一個界面設(shè)計工具集，MAYLAB將所有GUI支持度控件都集成在這個環(huán)境中，并提供界面外觀、屬性和行為響應(yīng)方式的設(shè)置方法。axes()。sound(z13)。Fp=4000。)。[h,w]=freqz(num12,den12)。wp1=tan(pi*Fp1/Ft)。)。g39。 ws=2*pi*Fs/Ft。plot(w,abs([X(1:256)]))。title(39。 IIR帶通濾波器的幅頻特性讀入信號對信號傅里葉變換進行頻譜分析，利用filter函數(shù)濾波，對得到的信號傅里葉變換進行頻譜分析，與濾波前進行比較。 Wb =2*fs*tan(wb/2)。s39。)。 wc=fc*2*pi/fs。例如：帶通濾波器用作頻譜分析儀中的選頻裝置；低通濾波器用作數(shù)字信號分析系統(tǒng)中的抗頻混濾波；高通濾波器被用于聲發(fā)射檢測儀中剔除低頻干擾噪聲；帶阻濾波器用作電渦流測振儀中的陷波器，等等。唐 山 學(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界面進行控制。MATLAB是由美國Mathworks公司推出的用于數(shù)值計算和圖形處理的科學(xué)計算系統(tǒng)環(huán)境。 Wb=2*fs*tan(wb/2)。 IIR低通濾波器的幅頻特性讀入信號對信號傅里葉變換進行頻譜分析，利用filter函數(shù)濾波，對得到的信號傅里葉變換進行頻譜分析，與濾波前進行比較。)。 Wc =2*fs*tan(wc/2)。并在濾波前后用sound函數(shù)播放語音信號，觀察濾波前后的變化。合成信號39。 xlabel(39。fp=2*Fp*tan(wp/2)。)。axes()。 wp2=tan(pi*Fp2/Ft)。 z12=filter(num12,den12,x)。axes()。Fs=3500。m13=fft(z13)。plot(x)。GUIDE將設(shè)計好的GUI保存在一個FIG文件中，同時生成M文件框架。我們可以像給一般信號做頻譜分析一樣，來對語音信號做頻譜分析，也可以較容易的用數(shù)字濾波器來對語音進行濾波處理。以上的內(nèi)容全都要通過程序的編寫來完成，不過主要的操作是對函數(shù)的調(diào)用，這就要我們熟悉MATLAB的一些函數(shù)的基本功能以及它們的基本語句。fp=2*Fp*tan(wp/2)。Fp1=1200。)。wp1=tan(pi*Fp/Ft)。IIR高通濾波器39。語音信號頻譜[y,fs,nbits]=wavread (39。x3=sin(2*pi*t*f3)。t=n/2000X=fft(x,512)。Fs=1200。C:\program files\matlab71\work\speech\39。plot(abs(m11),39。語音信號帶通濾波Ft=8000。s39。m12=fft(z12)。plot(s)。[n13,wn13]=cheb1ord(wp,ws,1,50,39。sound(z13)。axes()。x3=sin(2*pi*t*f3)。s39。axes()。)。Fs1=1000。s39。r39。f2=200。wp=1。g39。)。 title(39。 %求濾波后的信號axes()。wp1=tan(pi*Fp/Ft)。)。axes()。)。Fp1=1200。title(39。濾波前信號的頻譜39。fs=2*Fs*tan(ws/2)。x1=sin(2*pi*t*f1)。濾波后信號的頻譜39。 S=fft(s)。wp=1。)。 z12=filter(num12,den12,s)。ws=(wp1*wp2w.^2)/(bw*w)。濾波后的信號波形39。)。 [h,w]=freqz(num11,den11)。語音信號低通濾波Ft=8000。x=sin(2*pi*t*f1)+sin(2*pi*t*f2)+sin(2*pi*t*f3)。x1=sin(2*pi*t*f1)。原始信號波形39。plot(w*21000*,abs(h))。Fp=4000。[n12,wn12]=buttord(wp,ws,1,50,39。)。wp=2*pi*Fp/Ft。完成任務(wù)后才發(fā)現(xiàn)這次訓(xùn)練的內(nèi)容并不算難。課題的特色在于它將語音信號看作一個向量，于是就把語音數(shù)字化了。在MATLAB中GUI是一中包含多種對象的圖形窗口，并為GUI開發(fā)提供一個方便高效的集成開發(fā)環(huán)境GUIDE。)。 z13=filter(num13,den13,x)。 Ft=8000。濾波后信號頻譜39。[num12,den12]=bilinear(num2,den2,)。Ap=1。濾波前信號波形39。plot(abs(x),39。 wp=2*pi*Fp/Ft。axes()。plot(x)。)。 wc =[fc1 fc2]*2*pi/fs。,39。IIR低通濾波器幅頻特性39。IIR低通濾波器的設(shè)計程序為：wb=fb*2*pi/fs。模擬濾波器在測試系統(tǒng)或?qū)Ｓ脙x器儀表中是一種常用的變換裝置。設(shè)計要求：根據(jù)題目要求進行數(shù)字濾波器總體設(shè)計。MATLAB是英文MATrix LABoratory（矩形實驗室）的縮寫。Wc=2*fs*tan(wc/2)。并在濾波前后用sound函數(shù)播放語音信號，觀察濾波前后的變化。 [bd,ad]=bilinear(b,a,fs)。 [n,wn]=buttord(Wb,Wc,ap,as,39。 濾波前后信號的波形和頻譜比較分析濾波前后信號的變化，波形圖變窄。)。Hz39。 )。