語(yǔ)音信號(hào)處理與濾波(數(shù)字信號(hào)處理課程設(shè)計(jì))

《數(shù)字信號(hào)處理》課程設(shè)計(jì)報(bào)告題目：語(yǔ)音信號(hào)處理與濾波專業(yè)：通信工程學(xué)號(hào)：姓名：指導(dǎo)教師：周老師成績(jī)：二〇一二年一月三日一、課程設(shè)計(jì)目的綜合運(yùn)用數(shù)字信號(hào)處理的理論知識(shí)進(jìn)行頻譜分析和濾波器設(shè)計(jì)，通過(guò)理論推導(dǎo)得出相應(yīng)結(jié)論，再利用MATLAB作為編程工具進(jìn)行計(jì)算機(jī)實(shí)現(xiàn)，從而加深對(duì)所學(xué)知識(shí)的理解。二、課程設(shè)計(jì)基本要求及分組1、對(duì)知識(shí)點(diǎn)的掌握要求：（1）熟悉離散信號(hào)和系統(tǒng)的時(shí)域特性。（2）掌握序列快速傅里葉變換FFT方法。（3）學(xué)會(huì)MATLAB的使用，掌握MATLAB的程序設(shè)計(jì)方法。（4）利用MATLAB對(duì)語(yǔ)音信號(hào)進(jìn)行頻譜分析。（5）掌握MATLAB設(shè)計(jì)FIR和IIR數(shù)字濾波器的方法。三、詳細(xì)設(shè)計(jì)過(guò)程1、問(wèn)題描述：錄制一段自己的語(yǔ)音信號(hào)，取不同的數(shù)據(jù)點(diǎn)對(duì)語(yǔ)音信號(hào)進(jìn)行頻譜 分析；對(duì)所有數(shù)據(jù)進(jìn)行插值和抽取處理，改變抽樣率再對(duì)信號(hào)進(jìn)行頻譜分析； 設(shè)計(jì)FIR和IIR數(shù)字濾波器，并對(duì)被抽樣后的語(yǔ)音信號(hào)進(jìn)行濾波，分析濾波 后信號(hào)的時(shí)域和頻域特征，回放語(yǔ)音信號(hào)。2、詳細(xì)操作步驟和部分運(yùn)行結(jié)果（1）采集語(yǔ)音信號(hào)：利用Windows下的錄音機(jī)，錄制一段自己的話音，時(shí)間控制在4s左右；然 后在MATLAB軟件平臺(tái)下，利用函數(shù)wavread對(duì)語(yǔ)音信號(hào)進(jìn)行采樣，記住 采樣頻率和采樣點(diǎn)數(shù)。通過(guò)使用wavread函數(shù)，理解采樣頻率、采樣位數(shù)等 概念。wavread函數(shù)調(diào)用格式如下：y=wavread(file)，讀取file所規(guī)定的wav文件，返回采樣值放在向量y中。[y,fs,nbits]=wavread(file)，采樣值放在向量y中，fs表示采樣頻率(Hz)， bits表示采樣位數(shù)。y=wavread(file,N)，讀取前N點(diǎn)的采樣值放在向量y中。y=wavread(file,[N1，N2])，讀取從N1點(diǎn)到N2點(diǎn)的采樣值放在向量y中。錄制并存儲(chǔ)語(yǔ)音的程序代碼：(函數(shù)名'luyin'被主函數(shù)‘main’調(diào)用)fs=8000;t=4;a=msgbox('請(qǐng)說(shuō)話...');x=wavrecord(t*fs,fs,1,'Double');%錄制語(yǔ)音信號(hào)close(a);wavwrite(x,fs,'D:\MATLAB\R2010b\work\數(shù)字信號(hào)課程設(shè)計(jì)杜\課程設(shè)計(jì).wav');%存儲(chǔ)錄音信號(hào)原始語(yǔ)音時(shí)域信號(hào)波形圖：程序代碼：[x1,fs,bits]=wavread('D:\MATLAB\R2010b\work\數(shù)字信號(hào)課程設(shè)計(jì)\課程設(shè)計(jì).wav');sound(x1,fs,bits);%播放語(yǔ)音figure(1);subplot(1,1,1);plot(x1);%做原始語(yǔ)音信號(hào)的時(shí)域圖形title('原始語(yǔ)音信號(hào)作品');xlabel('時(shí)間n');ylabel('音量n');運(yùn)行結(jié)果：（2）對(duì)語(yǔ)音信號(hào)進(jìn)行頻譜分析：在MATLAB中，利用函數(shù)fft對(duì)信號(hào)進(jìn)行快速傅里葉變換，得到信號(hào)的頻 譜特性。首先畫出語(yǔ)音信號(hào)的時(shí)域波形，然后對(duì)語(yǔ)音信號(hào)進(jìn)行頻譜分析。程序代碼：[x1,fs,bits]=wavread('D:\MATLAB\R2010b\work\數(shù)字信號(hào)課程設(shè)計(jì) 忠\課程設(shè)計(jì).wav');y1=fft(x1,1024);f=fs*(0:511)/1024;figure(1)y1=fft(x1);%做length(x1)點(diǎn)的FFTy1=fftshift(y1);%平移，是頻率中心為0derta_fs=fs/length(x1);%設(shè)置頻譜的間隔，分辨率subplot(2,1,1);plot([-fs/2:derta_fs:fs/2-derta_fs],abs(y1));%畫出原始語(yǔ)音信號(hào)的頻譜圖title('原始語(yǔ)音信號(hào)的頻譜杜作品');subplot(2,1,2);plot(abs(y1(1:1024)));title('原始語(yǔ)音信號(hào)FFT頻譜 杜作品');xlabel('點(diǎn)數(shù)N');ylabel('幅值原始語(yǔ)音信號(hào)頻譜和FFT轉(zhuǎn)換后的頻譜圖（3）設(shè)計(jì)數(shù)字濾波器和畫出頻率響應(yīng)：根據(jù)語(yǔ)音信號(hào)的特點(diǎn)給出有關(guān)濾波器的性能指標(biāo)：低通濾波器性能指標(biāo)：fp=1000Hz，fc=1200Hz，As=100dB，Ap=1dB；高通濾波器性能指標(biāo)：fc=4800Hz，fp=5000Hz，As=100dB，Ap=1dB；帶通濾波器性能指標(biāo)：fp1=1200Hz，fp2=3000Hz，fc1=1000Hz，fc2=3200Hz，As=100dB，Ap=1dB。首先用窗函數(shù)法（矩形窗(Rectangularwindow)、三角窗(Triangularwindow)、 漢寧窗(Hanningwindow)、海明窗(Hammingwindow)、布拉克曼窗(Blackman window)、切比雪夫窗(Chebyshevwindow)、巴特里特窗(Bartlettwindow)及凱 塞窗(Kaiserwindow)。）設(shè)計(jì)上面要求的三種濾波器。在MATLAB中，利用 函數(shù)fir1設(shè)計(jì)FIR濾波器；然后再用雙線性變換法設(shè)計(jì)上面要求的三種濾波 器，利用函數(shù)butte、cheby1和ellip設(shè)計(jì)IIR濾波器；最后，利用MATLAB 中的函數(shù)freqz畫出各種濾波器的頻率響應(yīng)。（4）IIR濾波器的設(shè)計(jì)過(guò)程IIR低通濾波器：程序代碼：()fs=8000;T=1/fs;duration=3;[y,fs,bits]=wavread('D:\MATLAB\R2010b\work\數(shù)字信號(hào)課程設(shè)計(jì) 杜\課程設(shè)計(jì).wav');soundview('D:\MATLAB\R2010b\work\數(shù)字信號(hào)課程設(shè)計(jì)杜\課程設(shè) 計(jì).wav');N=length(y);t=(1:N)/fs;df=fs/N;n1=1:N/2;f=(n1-1)*df;subplot(2,1,1);plot(y);title('原始信號(hào)時(shí)域的圖形')xlabel('時(shí)間/t')ylabel('幅值')subplot(2,1,2);q=fft(y);plot(f,20*log10(abs(q(n1))));title('原始信號(hào)頻譜圖形')xlabel('頻率/Hz')ylabel('幅值/db')fs=8000;fp=1000;fst=1200;ap=1;ast=100;wp=2*fp/fs;wst=2*fst/fs;[N,wc]=cheb1ord(wp,wst,ap,ast);[num,den]=cheby1(N,ap,wc,'low');[h,w]=freqz(num,den);subplot(2,1,1);plot(w/pi,abs(h));grid;xlabel('\omega/\pi');ylabel('振幅（幅值）');title('契比雪夫Ⅰ型低通濾波器的幅頻響應(yīng)');subplot(2,1,2);plot(w/pi,20*log10(abs(h)));grid;xlabel('\omega/\pi');ylabel('振幅(分貝)');title('契比雪夫Ⅰ型低通濾波器的幅頻響應(yīng)');y0=filter(num,den,y);fftwav=fft(y0);subplot(2,2,1);plot(y);title('原始信號(hào)的波形')xlabel('時(shí)間/t')ylabel('幅值')subplot(2,2,2);plot(y0);title('通過(guò)契比雪夫Ⅰ型低通濾波器后的信號(hào)圖形')xlabel('時(shí)間/t')ylabel('幅值')subplot(2,2,3);q=fft(y);plot(f,20*log10(abs(q(n1))));title('濾波前的信號(hào)頻譜圖')xlabel('頻率/Hz')ylabel('幅值/db')subplot(2,2,4);plot(f,20*log10(abs(fftwav(n1))));title('濾波后信號(hào)的頻譜圖')xlabel('頻率/Hz')ylabel('幅值/db')wavwrite(y0,fs,'D:\MATLAB\R2010b\work\數(shù)字信號(hào)課程設(shè)計(jì)\ 課程設(shè)計(jì)IIR_LP.wav');[y0,fs,bits]=wavread('D:\MATLAB\R2010b\work\數(shù)字信號(hào)課程設(shè)計(jì) 杜\課程設(shè)計(jì)IIR_LP.wav');%語(yǔ)音回放soundview(y0,fs,'課程設(shè)計(jì)IIR_LP.wav運(yùn)行結(jié)果如下圖： IIR高通濾波器設(shè)計(jì) 程序代碼：fs=8000;T=1/fs;duration=3;[y,fs,bits]=wavread('D:\MATLAB\R2010b\work\數(shù)字信號(hào)課程設(shè)計(jì) 忠\課程設(shè)計(jì).wav');soundview('D:\MATLAB\R2010b\work\數(shù)字信號(hào)課程設(shè)計(jì)杜藝忠\課程設(shè) 計(jì).wav');N=length(y);t=(1:N)/fs;df=fs/N;n1=1:N/2;f=(n1-1)*df;subplot(2,1,1);plot(y);title('原始信號(hào)時(shí)域的圖形')xlabel('時(shí)間/t')ylabel('幅值')subplot(2,1,2);q=fft(y);plot(f,20*log10(abs(q(n1))));title('原始信號(hào)頻譜圖形')xlabel('頻率/Hz')ylabel('幅值/db')Fs=16000;fp=4000;fst=3800;ap=1;ast=100;wp=2*fp/Fs;wst=2*fst/Fs;[N,wc]=cheb1ord(wp,wst,ap,ast);[num,den]=cheby1(N,ap,wc,'high');[h,w]=freqz(num,den);subplot(2,1,1);plot(w/pi,abs(h));grid;xlabel('\omega/\pi');ylabel('振幅（幅值）');title('契比雪夫Ⅰ型高通濾波器的幅頻響應(yīng)');subplot(2,1,2);plot(w/pi,20*log10(abs(h)));grid;xlabel('\omega/\pi');ylabel('振幅(分貝)');title('契比雪夫Ⅰ型高通濾波器的幅頻響應(yīng)');y0=filter(num,den,y);fftwav=fft(y0);subplot(2,2,1);plot(y);title('原始信號(hào)的波形')xlabel('時(shí)間/t')ylabel('幅值')subplot(2,2,2);plot(y0);title('通過(guò)契比雪夫Ⅰ型高通濾波器后的信號(hào)圖形')xlabel('時(shí)間/t')ylabel('幅值')subplot(2,2,3);q=fft(y);plot(f,20*log10(abs(q(n1))));title('濾波前的信號(hào)頻譜圖')xlabel('頻率/Hz')ylabel('幅值/db')subplot(2,2,4);plot(f,20*log10(abs(fftwav(n1))));title('濾波后信號(hào)的頻譜圖')xlabel('頻率/Hz')ylabel('幅值/db')wavwrite(y0,fs,'D:\MATLAB\R2010b\work\數(shù)字信號(hào)課程設(shè)計(jì)\課程 設(shè)計(jì)IIR_HP.wav');[y0,fs,bits]=wavread('D:\MATLAB\R2010b\work\數(shù)字信號(hào)課程設(shè)計(jì) 忠\課程設(shè)計(jì)IIR_HP.wav');%語(yǔ)音回放soundview(y0,fs,'課程設(shè)計(jì)IIR_HP.wav');程序運(yùn)行結(jié)果如下圖：IIR帶通濾波器的設(shè)計(jì)程序代碼:fs=8000;T=1/fs;duration=3;[y,fs,bits]=wavread('D:\MATLAB\R2010b\work\數(shù)字信號(hào)課程設(shè)計(jì) 忠\課程設(shè)計(jì).wav');soundview('D:\MATLAB\R2010b\work\數(shù)字信號(hào)課程設(shè)計(jì)杜\課程設(shè) 計(jì).wav');N=length(y);t=(1:N)/fs;df=fs/N;n1=1:N/2;f=(n1-1)*df;subplot(2,1,1);plot(y);title('原始信號(hào)時(shí)域的圖形')xlabel('時(shí)間/t')ylabel('幅值')subplot(2,1,2);q=fft(y);plot(f,20*log10(abs(q(n1))));title('原始信號(hào)頻譜圖形')xlabel('頻率/Hz')ylabel('幅值/db')fs=8000;ast=100;ap=1;wp=[12003000]*2/fs;wst=[10003200]*2/fs;[N,wn]=cheb1ord(wp,wst,ap,ast);[num,den]=cheby1(N,ap,wn);[h,w]=freqz(num,den);subplot(2,1,1);plot(w/pi,abs(h));grid;xlabel('\omega/\pi');ylabel('振幅（幅值）');title('契比雪夫Ⅰ型帶通濾波器的幅頻響應(yīng)');subplot(2,1,2);plot(w/pi,20*log10(abs(h)));grid;xlabel('\omega/\pi');ylabel('振幅(分貝)');title('契比雪夫Ⅰ型帶通濾波器的幅頻響應(yīng)');y1=filter(num,den,y);fftwav=fft(y1);subplot(2,2,1);plot(y);title('原始信號(hào)的波形')xlabel('時(shí)間/t')ylabel('幅值')subplot(2,2,2);plot(y1);title('通過(guò)IIR帶通濾波器后的信號(hào)圖形')xlabel('時(shí)間/t')ylabel('幅值')subplot(2,2,3);q=fft(y);plot(f,20*log10(abs(q(n1))));title('濾波前的信號(hào)頻譜圖')xlabel('頻率/Hz')ylabel('幅值/db')subplot(2,2,4);plot(f,20*log10(abs(fftwav(n1))));title('帶通IIR濾波器濾波后信號(hào)的頻譜圖')xlabel('頻率/Hz')ylabel('幅值/db')wavwrite(y1,fs,'D:\MATLAB\R2010b\work\數(shù)字信號(hào)課程設(shè)計(jì)杜\課程 設(shè)計(jì)IIR_BP.wav');[y1,fs,bits]=wavread('D:\MATLAB\R2010b\work\數(shù)字信號(hào)課程設(shè)計(jì) 忠\課程設(shè)計(jì)IIR_BP.wav');%語(yǔ)音回放)soundview(y1,fs,'課程設(shè)計(jì)IIR_BP.wav');程序運(yùn)行結(jié)果：（5）FIR濾波器的設(shè)計(jì)FIR低通濾波器：程序代碼：fs=8000;T=1/fs;duration=3;[y,fs,sbits]=wavread('D:\MATLAB\R2010b\work\數(shù)字信號(hào)課程設(shè)計(jì) 杜藝\課程設(shè)計(jì).wav');soundview('D:\MATLAB\R2010b\work\數(shù)字信號(hào)課程設(shè)計(jì)杜\課程設(shè) 計(jì).wav');N=length(y);t=(1:N)/fs;df=fs/N;n1=1:N/2;f=(n1-1)*df;subplot(2,1,1);plot(y);title('原始信號(hào)時(shí)域的圖形')xlabel('時(shí)間/t')ylabel('幅值')subplot(2,1,2);q=fft(y);plot(f,20*log10(abs(q(n1))));title('原始信號(hào)頻譜圖形')xlabel('頻率/Hz')ylabel('幅值/db')fs=8000;fp=1000;fb=1200;ap=1;ast=100;wp=2*pi*fp/fs;wst=2*pi*fb/fs;deltaW=wst-wp;wc=(wst+wp)/2;N=ceil(2*6.4*pi/deltaW);hdWindow=ideallp(wc,N);%---ideallp（）函數(shù)（非系統(tǒng)自有函數(shù)）在系 統(tǒng)安裝目錄的WORK子目錄中，ideallp.mwdWindow=Kaiser(N,10.056);hr=wdWindow.*hdWindow';n=0:N-1;subplot(2,2,1);stem(n,wdWindow);xlabel('時(shí)間')ylabel('振幅作品')title('凱澤窗')subplot(2,2,2);stem(n,hr);xlabel('時(shí)間')ylabel('振幅')title('凱澤窗沖激響應(yīng)')[H,W]=freqz(hr,1);subplot(2,2,3);plot(W/pi,abs(H));xlabel('\omega/\pi')ylabel('絕對(duì)振幅')title('低通濾波器幅頻特性（1）')subplot(2,2,4);plot(W/pi,20*log10(abs(H)));xlabel('\omega/\pi')ylabel('相對(duì)振幅')title('低通濾波器幅頻特性（2）')y0=fftfilt(hr,y);fftwav=fft(y0);subplot(2,2,1);plot(y);title('原始信號(hào)的波形')xlabel('時(shí)間/t')ylabel('幅值')subplot(2,2,2);plot(y0);title('通過(guò)凱澤低通濾波器后的信號(hào)圖形')xlabel('時(shí)間/t')ylabel('幅值')subplot(2,2,3);q=fft(y);plot(f,20*log10(abs(q(n1))));title('濾波前的信號(hào)頻譜圖')xlabel('頻率/Hz')ylabel('幅值/db')subplot(2,2,4);plot(f,20*log10(abs(fftwav(n1))));title('濾波后信號(hào)的頻譜圖')xlabel('頻率/Hz')ylabel('幅值/db')wavwrite(y0,fs,'D:\MATLAB\R2010b\work\數(shù)字信號(hào)課程設(shè)計(jì)藝\ 課程設(shè)計(jì)FIR_LP.wav');[y0,fs,bits]=wavread('D:\MATLAB\R2010b\work\數(shù)字信號(hào)課程設(shè)計(jì) 杜\課程設(shè)計(jì)FIR_LP.wav');%語(yǔ)音回放soundview(y0,fs,'課程設(shè)計(jì)FIR_LP.wav');程序運(yùn)行結(jié)果：FIR高通濾波器：程序代碼：fs=8000;T=1/fs;duration=3;[y,fs,sbits]=wavread('D:\MATLAB\R2010b\work\數(shù)字信號(hào)課程設(shè)計(jì)忠\課程設(shè)計(jì).wav');soundview('D:\MATLAB\R2010b\work\數(shù)字信號(hào)課程設(shè)計(jì)杜\課程設(shè)計(jì).wav');N=length(y);t=(1:N)/fs;df=fs/N;n1=1:N/2;f=(n1-1)*df;subplot(2,1,1);plot(y);title('原始信號(hào)時(shí)域的圖形')xlabel('時(shí)間/t')ylabel('幅值')subplot(2,1,2);q=fft(y);plot(f,20*log10(abs(q(n1))));title('原始信號(hào)頻譜圖形')xlabel('頻率/Hz')ylabel('幅值/db')Fs=16000;fp=4000;fst=3800;ast=100;ap=1;wp=2*pi*fp/Fs;wst=2*pi*fst/Fs;wc=(wst+wp)/2;deltaW=wp-wst;N0=ceil(2*6.4*pi/deltaW);%---查表7-3（P342）凱澤窗N=N0+mod(N0+1,2);%---確保N為奇數(shù)N=ceil(2*6.4*pi/deltaW);hdWindow=ideallp(pi,N)-ideallp(wc,N);%---ideallp（）函數(shù)（非系統(tǒng)自有函數(shù)）在系統(tǒng)安裝目錄的WORK子目錄中，ideallp.mwdWindow=Kaiser(N,10.056);hr=wdWindow.*hdWindow';n=0:N-1;subplot(2,2,1);stem(n,wdWindow);xlabel('時(shí)間')ylabel('振幅')title('凱澤窗')subplot(2,2,2);stem(n,hr);xlabel('時(shí)間')ylabel('振幅')title('凱澤窗沖激響應(yīng)')[H,W]=freqz(hr,1);subplot(2,2,3);plot(W/pi,abs(H));xlabel('\omega/\pi')ylabel('絕對(duì)振幅')title('凱澤高通濾波器幅頻特性（1）')subplot(2,2,4);plot(W/pi,20*log10(abs(H)));xlabel('\omega/\pi')ylabel('相對(duì)振幅')title('凱澤高通濾波器幅頻特性（2）')y0=fftfilt(hr,y);fftwav=fft(y0);subplot(2,2,1);plot(y);title('原始信號(hào)的波形')xlabel('時(shí)間/t')ylabel('幅值')subplot(2,2,2);plot(y0);title('通過(guò)凱澤高通濾波器后的信號(hào)圖形')xlabel('時(shí)間/t')ylabel('幅值')subplot(2,2,3);q=fft(y);plot(f,20*log10(abs(q(n1))));title('濾波前的信號(hào)頻譜圖')xlabel('頻率/Hz')ylabel('幅值/db')subplot(2,2,4);plot(f,20*log10(abs(fftwav(n1))));title('凱澤高通濾波器濾波后信號(hào)的頻譜圖')xlabel('頻率/Hz')ylabel('幅值/db')wavwrite(y0,fs,'D:\MATLAB\R2010b\work\數(shù)字信號(hào)課程設(shè)計(jì)杜\課程 設(shè)計(jì)FIR_HP.wav');[y0,fs,bits]=wavread('D:\MATLAB\R2010b\work\數(shù)字信號(hào)課程設(shè)計(jì) 忠\課程設(shè)計(jì)FIR_HP.wav');%語(yǔ)音回放soundview(y0,fs,'課程設(shè)計(jì)FIR_HP.wav');程序運(yùn)行結(jié)果：FIR帶通濾波器：程序代碼：fs=8000;T=1/fs;duration=3;[y,fs,sbits]=wavread('D:\MATLAB\R2010b\work\數(shù)字信號(hào)課程設(shè)計(jì)\課程設(shè)計(jì).wav');soundview('D:\MATLAB\R2010b\work\數(shù)字信號(hào)課程設(shè)計(jì)杜\課程設(shè)計(jì).wav');N=length(y);t=(1:N)/fs;df=fs/N;n1=1:N/2;f=(n1-1)*df;subplot(2,1,1);plot(y);title('原始信號(hào)時(shí)域的圖形')xlabel('時(shí)間/t')ylabel('幅值')subplot(2,1,2);q=fft(y);plot(f,20*log10(abs(q(n1))));title('原始信號(hào)頻譜圖形')xlabel('頻率/Hz')ylabel('幅值/db')fs=8000;fp1=1200;fp2=3000;fst1=1000;fst2=3200;ast=100;ap=1;wp1=2*pi*fp1/fs;wp2=2*pi*fp2/fs;wst1=2*pi*fst1/fs;wst2=2*pi*fst2/fs;wc1=(wst1+wp1)/2;wc2=(wst2+wp2)/2;deltaW=min((wp1-wst1),(wst2-wp2));N0=ceil(2*6.4*pi/deltaW);%---查表7-3（P342）凱澤窗N=N0+mod(N0+1,2);%---確保N為奇數(shù)N=ceil(2*6.4*pi/deltaW);hdWindow=ideallp(wc2,N)-ideallp(wc1,N);%---ideallp（）函數(shù)（非系統(tǒng) 自有函數(shù)）在系統(tǒng)安裝目錄的WORK子目錄中，ideallp.mwdWindow=Kaiser(N,10.056);hr=wdWindow.*hdWindow';n=0:N-1;subplot(2,2,1);stem(n,wdWindow);xlabel('時(shí)間')ylabel('振幅')title('凱澤窗')subplot(2,2,2);stem(n,hr);xlabel('時(shí)間')ylabel('振幅')title('凱澤窗沖激響應(yīng)')[H,W]=freqz(hr,1);subplot(2,2,3);plot(W/pi,abs(H));xlabel('\omega/\pi')ylabel('絕對(duì)振幅')title('帶通濾波器幅頻特性（1）')subplot(2,2,4);plot(W/pi,20*log10(abs(H)));xlabel('\omega/\pi')ylabel('相對(duì)振幅')title('帶通濾波器幅頻特性（2）')y0=fftfilt(hr,y);fftwav=fft(y0);subplot(2,2,1);plot(y);title('原始信號(hào)的波形')xlabel('時(shí)間/t')ylabel('幅值')subplot(2,2,2);plot(y0);title('通過(guò)帶通濾波器后的信號(hào)圖形')xlabel('時(shí)間/t')ylabel('幅值')subplot(2,2,3);q=fft(y);plot(f,20*log10(abs(q(n1))));title('濾波前的信號(hào)頻譜圖')xlabel('頻率/Hz')ylabel('幅值/db')subplot(2,2,4);plot(f,20*log10(abs(fftwav(n1))));title('帶通濾波器濾波后信號(hào)的頻譜圖')xlabel('頻率/Hz')ylabel('幅值/db')wavwrite(y0,fs,'D:\MATLAB\R2010b\work\數(shù)字信號(hào)課程設(shè)計(jì)\課程 設(shè)計(jì)FIR_BP.wav');[y0,fs,bits]=wavread('D:\MATLAB\R2010b\work\數(shù)字信號(hào)課程設(shè)計(jì)杜藝 忠\課程設(shè)計(jì)FIR_BP.wav');%語(yǔ)音回放soundview(y0,fs,'課程設(shè)計(jì)FIR_BP.wav');程序運(yùn)行結(jié)果：（6）設(shè)計(jì)系統(tǒng)界面在所設(shè)計(jì)的系統(tǒng)界面上可以選擇濾波器的類型，輸入濾波器的參數(shù)，顯示濾 波器的頻率響應(yīng)，選擇信號(hào)等。主函數(shù)程序（main）：b=menu('忠作品版權(quán)所有，翻版必究','錄音','原始信號(hào)采樣后時(shí)域圖和頻譜 圖','FIR濾波器','IIR濾波器','退出');ifb==5b=0;Endifb==1luyin;main;elseifb==2temp=menu('作品版權(quán)所有，翻版必究','播放原始語(yǔ)音并生成時(shí)域圖',' 原始語(yǔ)音頻譜圖及做FFT變換后頻譜圖','返回');iftemp==1cjyy;main; elseiftemp==2