數(shù)字信號(hào)處理課程設(shè)計(jì)報(bào)告

1、現(xiàn)代信號(hào)處理課程設(shè)計(jì)任務(wù)實(shí)驗(yàn)報(bào)告 班級(jí)：通信080*姓名：*學(xué)號(hào)：*指導(dǎo)老師：*目 錄第一部分：實(shí)驗(yàn)設(shè)計(jì)目的第二部分：實(shí)驗(yàn)設(shè)計(jì)原理第三部分：實(shí)驗(yàn)設(shè)計(jì)方案第四部分：實(shí)驗(yàn)設(shè)計(jì)過程第五部分：個(gè)人小結(jié)參考文獻(xiàn)第一部分：實(shí)驗(yàn)設(shè)計(jì)目的1全面復(fù)習(xí)課程所學(xué)理論知識(shí)，鞏固所學(xué)知識(shí)重點(diǎn)和難點(diǎn)，將理論與實(shí)踐很好地結(jié)合起來(lái)。2提高綜合運(yùn)用所學(xué)知識(shí)獨(dú)立分析和解決問題的能力；3熟練使用一種高級(jí)語(yǔ)言進(jìn)行編程實(shí)現(xiàn)。第二部分：實(shí)驗(yàn)設(shè)計(jì)原理（一）信號(hào)分析1、編制信號(hào)生成程序，產(chǎn)生下述各序列，繪出它們的時(shí)域波形1） 三角波序列2）反三角波序列3） gaussian序列4） 衰減正弦序列 對(duì)連續(xù)信號(hào) 進(jìn)行采樣，可得到測(cè)試序列令a=5

2、0,采樣周期t=1ms,即fs=1000hz,f0=62.5,a=1005）已調(diào)信號(hào) ，其中調(diào)制信號(hào)，載波為。6）無(wú)限長(zhǎng)序列2. 對(duì)上述信號(hào)完成下列信號(hào)分析1）對(duì)三角波序列和反三角波序列，作n=8點(diǎn)的fft，觀察比較它們的幅頻特性，說(shuō)明它們有什么異同？繪出兩序列及其它們的幅頻特性曲線。 在和的尾部補(bǔ)零，作n=16點(diǎn)的fft，觀察它們的幅頻特性發(fā)生了什么變化？分析說(shuō)明原因。2)、觀察高斯序列，固定信號(hào)中的參數(shù)p=8，令q分別等于2，4，8，觀察它們的時(shí)域和幅頻特性，了解當(dāng)q取不同值時(shí)，對(duì)信號(hào)序列的時(shí)域幅頻特性的影響；固定q=8，令p分別等于8，13，14，觀察參數(shù)p變化對(duì)信號(hào)序列的時(shí)域及幅頻特性

3、的影響，觀察p等于多少時(shí)，會(huì)發(fā)生明顯的泄漏現(xiàn)象，混疊是否也隨之出現(xiàn)？記錄實(shí)驗(yàn)中觀察到的現(xiàn)象，繪出相應(yīng)的時(shí)域序列和幅頻特性曲線。3）觀察衰減正弦序列的時(shí)域和幅頻特性，繪出幅頻特性曲線，改變采樣頻率fs,使fs=300hz， 觀察此時(shí)的頻譜的形狀和譜峰出現(xiàn)位置？說(shuō)明產(chǎn)生現(xiàn)象的原因。4）觀察已調(diào)信號(hào)的時(shí)域與頻域特性，繪出其幅頻特性曲線，并對(duì)它進(jìn)行解調(diào)，以恢復(fù)原調(diào)制信號(hào)。5）對(duì)序列 ，截取序列長(zhǎng)度m=10點(diǎn)，計(jì)算其頻譜，畫出幅頻特性曲線；截取序列長(zhǎng)度m=20點(diǎn)，計(jì)算其頻譜，畫出幅頻特性曲線；（二）、數(shù)字濾波器的設(shè)計(jì)與實(shí)現(xiàn)設(shè)計(jì)和實(shí)現(xiàn)iir或fir數(shù)字濾波器：要求：輸入數(shù)字濾波器的濾波指標(biāo),包括通帶截止頻

4、率,通帶最大衰減,阻帶截止頻率,阻帶最小衰減。設(shè)計(jì)出相應(yīng)的數(shù)字濾波器。顯示得到的濾波器的系數(shù)和階數(shù)，求出系統(tǒng)的單位脈沖響應(yīng)h(n)。并畫出設(shè)計(jì)得到的濾波器的幅頻特性圖(要有坐標(biāo)標(biāo)度)和相頻特性。1對(duì)于iir數(shù)字濾波器：使用雙線性變換法?？蛇x用bw型、cb型和橢圓型等濾波器類型來(lái)進(jìn)行設(shè)計(jì)2對(duì)于fir數(shù)字濾波器：要對(duì)多種窗口（三種及以上）、不同階數(shù)設(shè)計(jì)結(jié)果進(jìn)行分析比較。（三）、數(shù)字濾波1）根據(jù)設(shè)計(jì)的濾波器，選擇2到3個(gè)測(cè)試信號(hào)進(jìn)行濾波。展示濾波后信號(hào)的幅頻與相頻特性，分析是否滿足濾波要求。對(duì)同一濾波要求，根據(jù)輸出信號(hào)頻譜，對(duì)比分析各類濾波器的差異。2）將一段語(yǔ)音作為測(cè)試信號(hào)，通過頻譜展示和語(yǔ)音播放

5、，對(duì)比分析濾波前后語(yǔ)音信號(hào)的變化，進(jìn)一步加深對(duì)數(shù)字信號(hào)處理的理解。第三部分：實(shí)驗(yàn)設(shè)計(jì)方案模擬帶通濾波器在信號(hào)的檢測(cè)和傳輸過程中起著很重要的作用，但在傳統(tǒng)模擬帶通濾波器的設(shè)計(jì)過程中需要大量繁瑣的數(shù)值計(jì)算，如果手工計(jì)算則費(fèi)時(shí)費(fèi)力。也有許多公司開發(fā)了相關(guān)的濾波器設(shè)計(jì)軟件，但是它們大多數(shù)使用起來(lái)不容易上手，而且價(jià)格比較昂貴。matlab中的圖形界面編程環(huán)境guide為用戶開發(fā)軟件界面提供了強(qiáng)有力的工具。在圖形界面編程環(huán)境guide下，設(shè)計(jì)好濾波器的設(shè)計(jì)工具界面后，對(duì)界面上各個(gè)控件的回調(diào)函數(shù)進(jìn)行編寫就能設(shè)計(jì)出交互式濾波器的分析設(shè)計(jì)軟件。本文在matlab gui基礎(chǔ)上設(shè)計(jì)了針對(duì)chebyshev型模擬帶

6、通濾波器的設(shè)計(jì)分析軟件。用戶在使用時(shí)只需輸入要設(shè)計(jì)的濾波器設(shè)計(jì)指標(biāo)值，選擇要設(shè)計(jì)的濾波器類型，即可在分析濾波器特性時(shí)只需選擇相應(yīng)菜單就能得出對(duì)應(yīng)的特性。該工具解決了用戶在設(shè)計(jì)chebyshe型模擬帶通濾波器時(shí)大量繁瑣的數(shù)值計(jì)算，不僅可以作為研究設(shè)計(jì)chebyshev型模擬帶通濾波器的輔助工具，還可以作為設(shè)計(jì)濾波器等相關(guān)課程的輔助教學(xué)軟件。1 模擬帶通濾波器的設(shè)計(jì)理論目前，模擬高通濾波器、帶通濾波器和帶阻濾波器的設(shè)計(jì)方法都是先將要設(shè)計(jì)的濾波器技術(shù)指標(biāo)通過某種頻率轉(zhuǎn)換關(guān)系轉(zhuǎn)換成模擬低通濾波器的技術(shù)指標(biāo)，并依據(jù)這些技術(shù)指標(biāo)設(shè)計(jì)出低通濾波器的轉(zhuǎn)移函數(shù)，然后再依據(jù)頻率轉(zhuǎn)換關(guān)系變成所要設(shè)計(jì)的濾波器的轉(zhuǎn)移函

7、數(shù)。工程實(shí)際中，設(shè)計(jì)高通、帶通、帶阻濾波器的常用方法是借助低通原型濾波器，經(jīng)頻率變換和元件變換得到。在matlab中，圖形用戶界面(graphical user interface，gui)是matlab中一個(gè)專用于gui程序設(shè)計(jì)的向?qū)гO(shè)計(jì)器，而gui是由各種圖形對(duì)象，如圖形窗口、圖軸、菜單、按鈕、文本框等構(gòu)建的用戶界面，是人機(jī)交互的有效工具和方法。通過guide可以很方便地設(shè)計(jì)出各種符合要求的圖形用戶界面。用戶通過一定的方法(如鼠標(biāo)或鍵盤)選擇、激活這些圖形對(duì)象，使計(jì)算機(jī)產(chǎn)生某種動(dòng)作或變化，比如實(shí)現(xiàn)計(jì)算、繪圖等。gui設(shè)計(jì)既能以基本的matlab程序設(shè)計(jì)為主，也能以鼠標(biāo)為主，利用guide工

8、具進(jìn)行設(shè)計(jì)，也可綜合以上兩種方法進(jìn)行設(shè)計(jì)。guide主要是一個(gè)界面設(shè)計(jì)工具集。matlab將所有g(shù)ul支持的用戶控件都集成在這個(gè)環(huán)境中并提供界面外觀、屬性和行為相應(yīng)方式的設(shè)置方法。guide將用戶保存設(shè)計(jì)好的圖形用戶界面保存在一個(gè)fig資源文件中，同時(shí)自動(dòng)生成包含圖形用戶界面初始化和組件界面布局控制代碼的m文件，這個(gè)m文件為實(shí)現(xiàn)回調(diào)函數(shù)的編寫提供了一個(gè)參考框架。fig文件是一個(gè)二進(jìn)制文件，包含系列化的圖形窗口對(duì)象。所有對(duì)象的屬性都是用戶創(chuàng)建圖形窗口時(shí)保存的屬性。該文件最主要的功能是對(duì)象句柄的保存。m文件包含gui設(shè)計(jì)、控制函數(shù)及控件的回調(diào)函數(shù)，主要用來(lái)控制gui展開時(shí)的各種特征。該文件基本上可

9、以分為gui初始化和回調(diào)函數(shù)2個(gè)部分，控件的回調(diào)函數(shù)可根據(jù)用戶與gui的具體交互行為分別調(diào)用。第四部分：實(shí)驗(yàn)設(shè)計(jì)過程（一）信號(hào)分析1、利用matlab編制信號(hào)生成程序，產(chǎn)生下述各序列，繪出它們的時(shí)域波形，并對(duì)各信號(hào)進(jìn)行相關(guān)分析：（1）三角波序列、反三角波序列及其fft三角波序列matlab源程序：subplot(1,1,1); n1=0:3;x1=n1+1;n2=4:7;x2=8-n2; stem(n1,n2,x1,x2,.);axis(-1,9,0,7);反三角波序列matlab源程序：subplot(1,1,1); n1=0:3;x1=4-n1;n2=4:7;x2=n2-3; stem(n

10、1,n2,x1,x2,.);axis(-1,9,0,10);它們的fft變換的matlab源程序：n1=8,n2=16;n1=0:n1-1;x1=1,2,3,4,4,3,2,1;subplot(4,2,1);stem(n1,x1,.);axis(0,10,0,5);title(x3(n);y1=fft(x1,n1);subplot(4,2,2);stem(n1,abs(y1),.);axis(0,10,0,8);title(fft,n=8);x2=4,3,2,1,1,2,3,4;subplot(4,2,3);stem(n1,x2,.);axis(0,10,0,5);title(x4(n);y2

11、=fft(x2,n1);subplot(4,2,4);stem(n1,abs(y2),.);axis(0,10,0,8);title(fft,n=8);n2=0:n2-1;x3=x1,zeros(1,8);subplot(4,2,5);stem(n2,x3,.);axis(0,20,0,5);title(x3(n);y3=fft(x3,n2);subplot(4,2,6);stem(n2,abs(y3),.);axis(0,20,0,20);title(fft,n=16);x4=x2,zeros(1,8);subplot(4,2,7);stem(n2,x4,.);axis(0,20,0,5);

12、title(x4(n);y4=fft(x4,n2);subplot(4,2,8);stem(n2,abs(y4),.);axis(0,20,0,20);title(fft,n=16)生成波形為：（2） gaussian序列及其fftp為定值：matlab源程序：n=0:15;x1=exp(-(n-8).2/2);subplot(3,2,1);stem(n,x1,.);axis(0,20,0,1);title(x5(n),p=8,q=2); y1=fft(x1,16);subplot(3,2,2);stem(n,abs(y1),.);axis(0,20,0,3);title(fft,p=8,q=

13、2);x2=exp(-(n-8).2/4);subplot(3,2,3);stem(n,x2,.);axis(0,20,0,1);title(x5(n),p=8,q=4);y2=fft(x2,16);subplot(3,2,4);stem(n,abs(y2),.);axis(0,20,0,4);title(fft,p=8,q=4); x3=exp(-(n-8).2/8);subplot(3,2,5);stem(n,x3,.);axis(0,20,0,1);title(x5(n),p=8,q=8);y3=fft(x3,16);subplot(3,2,6);stem(n,abs(y3),.);ax

14、is(0,20,0,6);title(x5(n),p=8,q=8);生成波形為：q為定值：n=0:15;x1=exp(-(n-8).2/8);subplot(4,2,1);stem(n,x1,.);axis(0,20,0,1);title(x5(n),p=8,q=8);y1=fft(x1,16);subplot(4,2,2);stem(n,abs(y1),.);axis(0,20,0,6);title(fft,p=8,q=8);x2=exp(-(n-13).2/8);subplot(4,2,3);stem(n,x2,.);axis(0,20,0,1);title(x5(n),p=13,q=8)

15、;y2=fft(x2,16);subplot(4,2,4);stem(n,abs(y2),.);axis(0,20,0,6);title(fft,p=13,q=8);x3=exp(-(n-14).2/8);subplot(4,2,5);stem(n,x3,.);axis(0,20,0,1);title(x5(n),p=14,q=8);y3=fft(x3,16);subplot(4,2,6);stem(n,abs(y3),.);axis(0,20,0,4);title(x5(n),p=14,q=8);x4=exp(-(n-15).2/8);subplot(4,2,7);stem(n,x4,.);

16、axis(0,20,0,1);title(x5(n),p=15,q=8);y4=fft(x4,16);subplot(4,2,8);stem(n,abs(y4),.);axis(0,20,0,4);title(fft,p=15,q=8);生成波形為：（3）衰減正弦序列及其不同采樣頻率時(shí)的fft：n=0:50;x1=50*exp(-0.1.*n).*sin(62.5*pi/500*n);subplot(2,2,1);stem(n,x1,.);axis(0,50,-20,40);title(x7(n),fs=1000hz);y1=fft(x1,51);subplot(2,2,2);stem(n,a

17、bs(y1),.);axis(0,50,0,250);title(fft,fs=1000hz);x2=50*exp(-1/3.*n).*sin(62.5*pi/150*n);subplot(2,2,3);stem(n,x2,.);axis(0,50,-20,40);title(x7(n),fs=300hz);y2=fft(x2,51);subplot(2,2,4);stem(n,abs(y2),.);axis(0,50,0,150);title(fft,fs=300hz生成波形為：（4）載波調(diào)制解調(diào)信號(hào)t=0:0.01:3; x=cos(pi*t).*cos(9*pi*t); plot(t,x

18、); axis(0 3 -2 2); grid; title(x5(t)=cos(pi*t)*cos(9*pi*t) t=0:0.01:3; x=cos(pi*t); plot(t,x); axis(0 3 -1 1); grid; title(解調(diào)信號(hào))（5）無(wú)限長(zhǎng)序列n=0:100; x=0.8.n; subplot(1,1,1);stem(n,x,.);axis(0 100 0 1);title(x6(n)=0.8n)n=0:9 x1=0.8.n; subplot(2,2,1);stem(n,x1,.);axis(0,20,0,1);title(序列x6(n),m=10); y1=fft

19、(x1,10); subplot(2,2,2);stem(n,abs(y1),.);axis(0,20,0,3);title(m=10,fft,幅頻特性);m=0:19 x2=0.8.m; subplot(2,2,3);stem(m,x2,.);axis(0,20,0,1);title(序列x6(n),m=20); y2=fft(x2,20); subplot(2,2,4);stem(m,abs(y2),.);axis(0,20,0,3);title(m=20,fft,幅頻特性)（二）、數(shù)字濾波器的設(shè)計(jì)與實(shí)現(xiàn)采用雙線性變換法設(shè)計(jì)iir數(shù)字濾波器：源程序：低通：global fpglobal f

20、sglobal apglobal asglobal fs wp=2*pi*fp/fs;ws=2*pi*fs/fs;ts=1/fs;wp1=2/ts*tan(wp/2); ws1=2/ts*tan(ws/2); n,wn=buttord(wp1,ws1,ap,as,s); z,p,k=buttap(n); bap,aap=zp2tf(z,p,k);b,a=lp2lp(bap,aap,wn); bz,az=bilinear(b,a,fs); figure(1)freqz(bz,az); y,t=impz(bz,az,50);figure(2)stem(t,y,.);高通：global fpglob

21、al fsglobal apglobal asglobal fs wp=2*pi*fp/fs;ws=2*pi*fs/fs;ts=1/fs;wp1=2/ts*tan(wp/2); ws1=2/ts*tan(ws/2); n,wn=buttord(wp1,ws1,ap,as,s); z,p,k=buttap(n); bap,aap=zp2tf(z,p,k);b,a=lp2hp(bap,aap,wn); bz,az=bilinear(b,a,fs);figure(1)freqz(bz,az); y,t=impz(bz,az,50);figure(2)stem(t,y,.);帶通：global fp1

22、global fp2global fs1global fs2global apglobal asglobal fs wp1=2*pi*fp1/fs;wp2=2*pi*fp2/fs;ws1=2*pi*fs1/fs;ws2=2*pi*fs2/fs;b=ws2-ws1;w0=sqrt(ws1*ws2);wp=max(abs(wp1),abs(wp2);ws=1;n,wc=buttord(wp,ws,ap,as,s);num,den=butter(n,wc,s);numt,dent=lp2bp(num,den,w0,b);numd,dend=bilinear(numt,dent,1);w=linspa

23、ce(0,pi,512);figure(1)freqz(numd,dend,w);y,t=impz(numd,dend,50);figure(2)stem(t,y,.)；帶阻：global fp1global fp2global fs1global fs2global apglobal asglobal fs wp1=2*pi*fp1/fs;wp2=2*pi*fp2/fs;ws1=2*pi*fs1/fs;ws2=2*pi*fs2/fs;b=ws2-ws1;w0=sqrt(ws1*ws2);wp=max(abs(wp1),abs(wp2);ws=1;n,wc=buttord(wp,ws,ap,a

24、s,s);num,den=butter(n,wc,s);numt,dent=lp2bs(num,den,w0,b);numd,dend=bilinear(numt,dent,1);w=linspace(0,pi,512);figure(1)freqz(numd,dend,w);y,t=impz(numd,dend,50);figure(2)stem(t,y,.);生成波形為：低通及其單位脈沖響應(yīng)：高通及其單位脈沖響應(yīng)：帶通及其單位脈沖序列：帶阻及其單位脈沖響應(yīng)：（三）、綜合應(yīng)用語(yǔ)音信號(hào)的分析處理源程序：讀入語(yǔ)音信號(hào)：x1,fs=wavread(c:tada.wav);sound(x1,fs);

25、fs y1=fft(x1,1024); f=fs*(0:511)/1024;figure(1)plot(x1) title(原始語(yǔ)音信號(hào));xlabel(time n);ylabel(fuzhi n);figure(2)freqz(x1) title(頻率響應(yīng)圖)figure(3)subplot(2,1,1);plot(abs(y1(1:512) title(原始語(yǔ)音信號(hào)fft頻譜)subplot(2,1,2);plot(f,abs(y1(1:512);title(原始語(yǔ)音信號(hào)頻譜)xlabel(hz);ylabel(fuzhi);語(yǔ)音信號(hào)低通處理：global fpglobal fsglob

26、al apglobal asx1,fs=wavread(c:tada.wav);t=0:1/fs:(size(x1)-1)/fs;fswp=2*pi*fp/fs;ws=2*pi*fs/fs;ts=1/fs;wp1=2/ts*tan(wp/2); ws1=2/ts*tan(ws/2); n,wn=buttord(wp1,ws1,ap,as,s); z,p,k=buttap(n); bap,aap=zp2tf(z,p,k);b,a=lp2lp(bap,aap,wn); bz,az=bilinear(b,a,fs); h,w=freqz(bz,az); figure(1)plot(w*fs/(2*p

27、i),abs(h)gridxlabel(頻率hz)ylabel(頻率響應(yīng)幅度)title(butterworth)f1=filter(bz,az,x1);figure(2)subplot(2,1,1)plot(t,x1); title(濾波前的時(shí)域波形);subplot(2,1,2)plot(t,f1); title(濾波后的時(shí)域波形);sound(f1,22050); f0=fft(f1,1024);f=fs*(0:511)/1024;figure(3)y2=fft(x1,1024);subplot(2,1,1);plot(f,abs(y2(1:512); title(濾波前的頻譜)xlab

28、el(hz);ylabel(fuzhi);subplot(2,1,2)f1=plot(f,abs(f0(1:512); title(濾波后的頻譜)xlabel(hz);ylabel(fuzhi);語(yǔ)音信號(hào)高通處理：global fpglobal fsglobal apglobal asx1,fs=wavread(c:tada.wav);t=0:1/fs:(size(x1)-1)/fs;fswp=2*pi*fp/fs;ws=2*pi*fs/fs;ts=1/fs;wp1=2/ts*tan(wp/2); ws1=2/ts*tan(ws/2); n,wn=buttord(wp1,ws1,ap,as,s

29、); z,p,k=buttap(n); bap,aap=zp2tf(z,p,k);b,a=lp2hp(bap,aap,wn); bz,az=bilinear(b,a,fs); h,w=freqz(bz,az); figure(1)plot(w*fs/(2*pi),abs(h)gridxlabel(頻率hz)ylabel(頻率響應(yīng)幅度)title(butterworth)f1=filter(bz,az,x1);figure(2)subplot(2,1,1)plot(t,x1); title(濾波前的時(shí)域波形);subplot(2,1,2)plot(t,f1); title(濾波后的時(shí)域波形);s

30、ound(f1,22050); f0=fft(f1,1024);f=fs*(0:511)/1024;figure(3)y2=fft(x1,1024);subplot(2,1,1);plot(f,abs(y2(1:512); title(濾波前的頻譜)xlabel(hz);ylabel(fuzhi);subplot(2,1,2)f1=plot(f,abs(f0(1:512); title(濾波后的頻譜)xlabel(hz);ylabel(fuzhi);語(yǔ)音信號(hào)帶通處理：global fp1global fp2global fs1global fs2global apglobal asx1,fs=

31、wavread(c:tada.wav);t=0:1/fs:(size(x1)-1)/fs;wp1=2*pi*fp1/fs;wp2=2*pi*fp2/fs;ws1=2*pi*fs1/fs;ws2=2*pi*fs2/fs;b=ws2-ws1;w0=sqrt(ws1*ws2);wp=max(abs(wp1),abs(wp2);ws=1;n,wc=buttord(wp,ws,ap,as,s);num,den=butter(n,wc,s);numt,dent=lp2bp(num,den,w0,b);numd,dend=bilinear(numt,dent,1);w=linspace(0,pi,512);

32、figure(3)h,w=freqz(numd,dend,w); figure(1)plot(w*fs/(2*pi),abs(h)gridxlabel(頻率hz)ylabel(頻率響應(yīng)幅度)title(butterworth)f1=filter(numd,dend,x1);figure(2)subplot(2,1,1)plot(t,x1) title(濾波前的時(shí)域波形);subplot(2,1,2)plot(t,f1); title(濾波后的時(shí)域波形);sound(f1,22050); f0=fft(f1,1024);f=fs*(0:511)/1024;figure(3)y2=fft(x1,1

33、024);subplot(2,1,1);plot(f,abs(y2(1:512); title(濾波前的頻譜)xlabel(hz);ylabel(fuzhi);subplot(2,1,2)f1=plot(f,abs(f0(1:512); title(濾波后的頻譜)xlabel(hz);ylabel(fuzhi);語(yǔ)音信號(hào)帶阻處理：global fp1global fp2global fs1global fs2global apglobal asx1,fs=wavread(c:tada.wav);t=0:1/fs:(size(x1)-1)/fs;wp1=2*pi*fp1/fs;wp2=2*pi*

34、fp2/fs;ws1=2*pi*fs1/fs;ws2=2*pi*fs2/fs;b=ws2-ws1;w0=sqrt(ws1*ws2);wp=max(abs(wp1),abs(wp2);ws=1;n,wc=buttord(wp,ws,ap,as,s);num,den=butter(n,wc,s);numt,dent=lp2bs(num,den,w0,b);numd,dend=bilinear(numt,dent,1);w=linspace(0,pi,512);figure(3)h,w=freqz(numd,dend,w); figure(1)plot(w*fs/(2*pi),abs(h)gridx

35、label(頻率hz)ylabel(頻率響應(yīng)幅度)title(butterworth)f1=filter(numd,dend,x1);figure(2)subplot(2,1,1)plot(t,x1) title(濾波前的時(shí)域波形);subplot(2,1,2)plot(t,f1); title(濾波后的時(shí)域波形);sound(f1,22050); f0=fft(f1,1024);f=fs*(0:511)/1024;figure(3)y2=fft(x1,1024);subplot(2,1,1);plot(f,abs(y2(1:512); title(濾波前的頻譜)xlabel(hz);ylabel(fuzhi);subplot(2,1,2)f1=plot(f,abs(f0(1:512); title(濾波后的頻譜)xlabel(hz);ylabel(fuzhi);生成