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