Matlab和雙線性變換的濾波器設(shè)計

1、武漢理工大學(xué)MATLAB課程設(shè)計報告書題目: MATLAB課程設(shè)計 基于 Matlab 和雙線性變換的濾波器設(shè)計初始條件：MATLAB仿真軟件數(shù)字信號處理與圖像處理基礎(chǔ)知識雙線性變換法的原理和算法巴特沃斯高通濾波器的性能指標(biāo)要求完成的主要任務(wù) :利用 MATLAB仿真軟件系統(tǒng)結(jié)合雙線性變換法設(shè)計一個數(shù)字巴特沃斯高通、低通和帶通 IIR 濾波器。說明：自己確定濾波器的不同參數(shù)得到不同的結(jié)果，最后自己分析原因。課程設(shè)計題目由指導(dǎo)教師提供，每人隸屬一組完成任務(wù)，每組完成的內(nèi)容不能雷同（按學(xué)號分組）時間安排：序階段內(nèi)容所需時間號1方案設(shè)計1 天2軟件設(shè)計2 天3系統(tǒng)調(diào)試1 天4答辯1 天合計5 天指導(dǎo)

2、教師簽名：年月日系主任（或責(zé)任教師）簽名：年月日武漢理工大學(xué)MATLAB課程設(shè)計報告書目 錄武漢理工大學(xué)MATLAB課程設(shè)計報告書摘 要數(shù)字濾波器是數(shù)字信號處理的基礎(chǔ)，用來對信號進(jìn)行過濾、檢測與參數(shù)估計等處理，在通信、圖像、語音、雷達(dá)等許多領(lǐng)域都有著十分廣泛的應(yīng)用。尤其在圖像處理、數(shù)據(jù)壓縮等方面取得了令人矚目的進(jìn)展和成就。鑒于此，數(shù)字濾波器的設(shè)計就顯得尤為重要。此報告重點介紹了用雙線性不變法設(shè)計IIR 數(shù)字濾波器的基本流程，比較了各種設(shè)計方法的優(yōu)缺點，總結(jié)了模擬濾波器的性能特征。最后以雙線性不變法設(shè)計了一個高通巴特沃斯 FIR 數(shù)字濾波器，介紹了設(shè)計步驟，然后在Matlab 環(huán)境下進(jìn)行了仿真與

3、調(diào)試，實現(xiàn)了設(shè)計目標(biāo)。關(guān)鍵字：數(shù)字濾波器巴特沃斯MATLAB雙線性變法武漢理工大學(xué)MATLAB課程設(shè)計報告書AbstractMatlab for research and application of digital filter provides an intuitive, efficient and convenient weapon. He matrix operations based on the calculation, visualization, program design integrated into an interactive working environment.

4、 In particular, Matlab toolbox to enable researchers in various fields can be easy and intuitive application of scientific research and engineering.This report introduced with emphasis of the basic flow of designing the IIR digit filter by the bilinear political reform, compared with each kind of de

5、sign method's good and bad points, summarized analog filter's performance characteristic. Finally design one by the bilinear political reform to pass Butterworth high IIR digit filter, introduced the design procedure, then has carried on the simulation and the debugging under the Matlab envi

6、ronment, has achieved the project objective.Keyword: digital filterButterworthMATLAB武漢理工大學(xué)MATLAB課程設(shè)計報告書1 數(shù)字濾波器1.1數(shù)字濾波器的基本概念數(shù)字濾波器是數(shù)字信號處理的重要基礎(chǔ)，是對信號都是過濾檢測與參數(shù)估計等處理過程中，它是使用最為廣泛的一種線性系統(tǒng)。數(shù)字濾波器處理的對象是經(jīng)由采樣期間將模擬信號轉(zhuǎn)換而得到的數(shù)字信號。數(shù)字濾波器是指完成信號濾波處理功能的，用有限精度算法實現(xiàn)的離散時間線性非時變系統(tǒng)。數(shù)字濾波器的輸入是一組數(shù)字量。它本身既可以是用數(shù)字硬件裝配而成的一臺用于完成給定運算的專用數(shù)字

7、計算機(jī)，也可以是將所需的運算編寫的程序通過計算機(jī)來執(zhí)行。數(shù)字濾波器具有穩(wěn)定性高、精度高、靈活性大等突出優(yōu)點。隨著數(shù)字技術(shù)的發(fā)展，用數(shù)字技術(shù)實現(xiàn)濾波器的功能愈來愈受到人們的重視，并得到了廣泛的應(yīng)用。數(shù)字濾波器的數(shù)學(xué)運算通常有兩種實現(xiàn)方式。一種是頻域法，即利用 FFT 快速運算方法對輸入信號進(jìn)行離散傅里葉變換， 分析其頻譜，然后再根據(jù)所希望的頻率特性進(jìn)行濾波，再利用傅里葉反變換得到時域信號。這種方法具有較好的頻域選擇特性和靈活性，并且由于信號頻率與所希望的頻率特性是簡單的相乘關(guān)系，所以它比計算等價的時域卷積要快得多。另一種方法是時域法，這種方法是通過離散的抽樣數(shù)據(jù)做差分?jǐn)?shù)學(xué)運算來達(dá)到濾波目的的。1

8、.2數(shù)字濾波器的分類數(shù)字濾波器按照不同的分類方法，可分為許多種，但總體來講可以分成兩大類。一類稱為經(jīng)典濾波器，即一般濾波器，特點是輸入信號中有用的頻率成分和希望濾除的頻率成分各占有不同的頻帶，通過一個合適的選頻濾波器達(dá)到濾波的目的。但對于一般濾波器，如果信號和干擾的頻帶互不重疊，則不能完成對干擾的有效濾除，這是需要采用另一類所謂的現(xiàn)代濾波器，例如維納濾波器、卡爾曼濾波器、自適應(yīng)濾波器等最佳濾波器，這些濾波器可按照隨機(jī)信號內(nèi)部的一些統(tǒng)計分布規(guī)律，從干擾中最佳地提取信號。從功能上，一般數(shù)字濾波器可以分為低通、高通、帶通、帶阻和全通等，此種分類方法是和模擬濾波器一樣的。武漢理工大學(xué)MATLAB課程設(shè)

9、計報告書1.3數(shù)字濾波器的MATLAB實現(xiàn)MATLAB工具箱未濾波器的設(shè)計應(yīng)用提供了豐富而簡便的方法，如函數(shù)方法和圖形工具方法等，使原來非常繁瑣復(fù)雜的程序設(shè)計變成了簡單的函數(shù)調(diào)用，為濾波器的設(shè)計和實現(xiàn)開辟了廣闊的天地。數(shù)字濾波器的一般設(shè)計過程為：按照實際需要，確定濾波器的性能要求；用一個因果的、穩(wěn)定的離散線性時不變系統(tǒng)，去逼近這一性能指標(biāo)；用有限精度的運算實現(xiàn)所設(shè)計的系統(tǒng)；通過模擬，驗證所設(shè)計的系統(tǒng)是否符合給定性能要求。Dd武漢理工大學(xué)MATLAB課程設(shè)計報告書2 雙線性變換法2.1雙線性變換法知識簡介脈沖響應(yīng)不變法的主要缺點是產(chǎn)生頻率響應(yīng)的混疊失真。這是因為從S 平面到平面是多值的映射關(guān)系所

10、造成的。為了克服這一缺點，可以采用非線性頻率壓縮方法，將整個頻率軸上的頻率范圍壓縮到- / T / T 之間，再用 z=esT 轉(zhuǎn)換到 Z 平面上。也就是說，第一步先將整個 S 平面壓縮映射到S1 平面的 - / T / T 一條橫帶里；第二步再通過標(biāo)準(zhǔn)變換關(guān)系 z=es1T 將此橫帶變換到整個 Z 平面上去。這樣就使 S 平面與 Z 平面建立了一一對應(yīng)的單值關(guān)系，消除了多值變換性，也就消除了頻譜混疊現(xiàn)象。由圖 2-1 看出，在零頻率附近，模擬角頻率 與數(shù)字頻率 之間的變換關(guān)系接近于線性關(guān)系；但當(dāng) 進(jìn)一步增加時， 增長得越來越慢，最后當(dāng) 時， 終止在折疊頻率 =處，因而雙線性變換就不會出現(xiàn)由于

11、高頻部分超過折疊頻率而混淆到低頻部分去的現(xiàn)象，從而消除了頻率混疊現(xiàn)象。2tan 2 T o圖 2-12.2雙線性變換法設(shè)計數(shù)字濾波器原理將 S 平面 j軸壓縮變換到 s1 平面 j軸上的 -/T 到/T 一段，可以采用以下變換關(guān)系：tan1T2這樣變換到 1，0 變換到1 0 ，可將上式寫成T武漢理工大學(xué)MATLAB課程設(shè)計報告書j1Te 2jj 1Te2j1Tj1Te2e2令 js ， j 1s1 ，解析延拓到整個s 平面和 s1 平面，可得s1Ts1Tse 2e 2s1Ts1 Te 2e 2再將 s1 平面通過以下變換關(guān)系映射到z 平面，即z es1T從而得到 s 平面和 z 平面的單值映

12、射關(guān)系為1zs1z111sz1s一般來說，為了使模擬濾波器的某一頻率與數(shù)字濾波器的任一頻率有對應(yīng)關(guān)系，可引入待定常數(shù)c，c tan1T2s c1e s1T1e s1T將 zes1T 代入到上式，可得1zs cz111cszcs在 MATLAB中，雙線性 Z 變換可以通過 bilinear函數(shù)實現(xiàn)，其調(diào)用格式為：Bz ，Az bilinear(B，A，F(xiàn)s) ；，為模擬濾波器傳遞函數(shù)G s的分子分母多項式的系數(shù)向量，而，為數(shù)字濾其中B A( )Bz Az波器的傳遞函數(shù)H z的分子分母多項式的系數(shù)向量。( )武漢理工大學(xué)MATLAB課程設(shè)計報告書3 設(shè)計任務(wù)及方案選擇3.1設(shè)計任務(wù)與要求利用 MA

13、TLAB 仿真軟件系統(tǒng)結(jié)合雙線性變換法設(shè)計一個數(shù)字巴特沃斯高通、低通和帶通 FIR 濾波器。說明：自己確定濾波器的不同參數(shù)得到不同的結(jié)果，最后自己分析原因。假設(shè)一個信號：x ( t )sin 2f 1t0 .5 cos 2f 2 t其中 f1=100Hz，f2=550Hz。試將該信號分別通過低通、高通、帶通濾波器，對輸出信號進(jìn)行分析比較。Step1:Step2:Step3:Step4:3.2方案設(shè)計論證與選擇本次濾波器設(shè)計的主要方法是先設(shè)計低通模擬濾波器，然后轉(zhuǎn)換為高通、帶通或帶阻數(shù)字濾波器。對于其他如高通，帶通，則通過頻率變換轉(zhuǎn)換為設(shè)計相應(yīng)的高通，帶通等。在設(shè)計的全過程的各個步驟，matla

14、b 都提供相應(yīng)的工具箱函數(shù)，使得IIR 數(shù)字濾波器設(shè)計變得非常簡單?？偟膩碚f，我的設(shè)計思路主要有以下兩種：方案一：從歸一化模擬低通原型出發(fā)，先在模擬域內(nèi)經(jīng)頻率變換成為所需類型的模擬濾波器；然后進(jìn)行雙線性變換，由S域變換到 Z域，而得到所需類型的數(shù)字濾波器。歸一化模擬模擬域模擬高 , 帶通雙線性變換法數(shù)字高 ,帶通低通原型或帶阻或帶阻頻率變換圖 3-1 先頻率變換再離散方案二：先進(jìn)行雙線性變換，將模擬低通原型濾波器變換成數(shù)字低通濾波器；然后在Z域內(nèi)經(jīng)數(shù)字頻率變換為所需類型的數(shù)字濾波器。武漢理工大學(xué)MATLAB課程設(shè)計報告書歸一化模擬數(shù)字原型低數(shù)字域數(shù)字高 ,帶通低通原型通頻率變換或帶阻雙線性變換

15、法圖 3-2先離散再頻率變換以上兩種方案都可以，我最后選擇了第一種方案進(jìn)行設(shè)計，即先在模擬域內(nèi)經(jīng)頻率變換成為所需類型的模擬濾波器；然后進(jìn)行雙線性變換，由S域變換到 Z域，而得到所需類型的數(shù)字濾波器。因為高通，帶通濾波器的設(shè)計方法是先將要設(shè)計的濾波器的技術(shù)特性指標(biāo)通過頻率轉(zhuǎn)換成模擬低通濾波器的技術(shù)指標(biāo)，再根據(jù)這些性能指標(biāo)設(shè)計出低通濾波器的傳遞函數(shù)，為了設(shè)計程序的簡潔，故選擇了方案一。武漢理工大學(xué)MATLAB課程設(shè)計報告書4 程序設(shè)計與調(diào)試4.1設(shè)計要求若實驗設(shè)計要求得巴特沃斯高通濾波器的性能指標(biāo)為：通帶邊界頻率 fp： fp=400Hz；通帶波紋Rp： Rp=1dB；阻帶邊界頻率 fs： fs=

16、300Hz；阻帶衰減Rs： Rs=30dB；采樣頻率Fs： Fs=2000Hz；要求巴特沃斯低通濾波器通帶邊界頻率與阻帶邊界頻率與高通的正好相反。模擬帶通濾波器中心頻率為f0=200Hz，帶寬為 Bw=400Hz ，4.2程序設(shè)計與說明4.2.1模擬低通濾波首先應(yīng)用雙線性變換作為數(shù)字低通濾波器到模擬高通濾波器的頻率預(yù)畸變處理：pp2 tanT22T cot wssp222fT利用上式確定的技術(shù)指標(biāo)和已知的參數(shù)Rp、Rs，就可以進(jìn)行模擬原型低通濾波器設(shè)計了。該部分的程序設(shè)計與注釋如下：Fs=2000;Ts=1/Fs;%采樣頻率和采樣周期wp=400*2*pi;ws=300*2*pi;%低通、高通

17、數(shù)字濾波器截止頻率Rp=1;Rs=30;%通帶波紋和阻帶衰減wpl=Fs*2*tan(ws*Ts/2);%按頻率轉(zhuǎn)換公式進(jìn)行轉(zhuǎn)換wsl=Fs*2*tan(wp*Ts/2);N,Wn=buttord(wpl,wsl,Rp,Rs,'s');%計算模擬濾波器的最小階數(shù)z,p,k=buttap(N);%設(shè)計模擬原型低通濾波器武漢理工大學(xué)MATLAB課程設(shè)計報告書b,a=zp2tf(z,p,k);%變零極點增益形式為系統(tǒng)傳遞函數(shù)形式H,w=freqs(b,a);%求模擬低通濾波器的頻率響應(yīng)figure;plot(w,abs(H);grid;%繪制低通頻率特性title(' 模擬低

18、通原型 ');xlabel('w');ylabel(' 幅度 ');4.2.2模擬到數(shù)字的轉(zhuǎn)換（1）從低通向低通的轉(zhuǎn)換 b,a=lp2lp(Bap,Aap,Wn); 功能是把模擬濾波器原型轉(zhuǎn)換成截至頻率為Wn 的低通濾波器。再由得到的模擬低通濾波器通過函數(shù)Bz ，Az bilinear(B ，A ，F(xiàn)s)得到數(shù)字低通濾波器。該部分的程序設(shè)計如下：e1,f1=lp2lp(b,a,300*2*pi) ;%原型低通轉(zhuǎn)換為低通模擬濾波器h1,w1=freqs(e1,f1);%求模擬低通濾波器的頻率響應(yīng)figure;plot(w1/(2*pi),abs(h1);%

19、 繪制高通頻率特性grid;title(' 模擬低通濾波器 ');xlabel('頻率 Hz');ylabel(' 幅度 ');r1,q1=bilinear(e1,f1,1000);%模擬低通濾波器轉(zhuǎn)換為數(shù)字濾波器s1,w2=freqz(r1,q1);%求數(shù)字低通濾波器的頻率響應(yīng)figure;plot(w2/pi, 20*log10(abs(s1);grid;%繪制數(shù)字低通濾波器的頻率響應(yīng)title(' 數(shù)字低通濾波器 ');xlabel('w/pi');ylabel(' 幅度 (dB)');（2）

20、從低通向高通的轉(zhuǎn)換b,a=lp2hp(Bap,Aap,Wn);功能是把模擬濾波器原型轉(zhuǎn)換成截至頻率為 Wn 的模擬高通濾波器，再進(jìn)一步轉(zhuǎn)換成數(shù)字高通濾波器。程序設(shè)計如下：e2,f2=lp2hp(b,a,300*2*pi) ;%原型低通模擬濾波器轉(zhuǎn)換為高通模擬濾波器h2,w3=freqs(e2,f2);%求模擬高通濾波器的頻率響應(yīng)武漢理工大學(xué)MATLAB課程設(shè)計報告書figure;plot(w3/(2*pi),abs(h2);%繪制高通頻率特性grid;title(' 模擬高通濾波器 ');xlabel('頻率 Hz');ylabel(' 幅度 '

21、);r2,q2=bilinear(e2,f2,1000);%模擬高通濾波器轉(zhuǎn)換為數(shù)字濾波器s2,w4=freqz(r2,q2);%求數(shù)字高通濾波器的頻率響應(yīng)figure;plot(w4/pi, 20*log10(abs(s2);grid;%繪制數(shù)字高通濾波器的頻率響應(yīng)title(' 數(shù)字高通濾波器 ');xlabel('w/pi');ylabel(' 幅度 (dB)');（3）從低通向帶通的轉(zhuǎn)換b,a=lp2bp(Bap,Aap,W0， Bw); 功能是把模擬濾波器原型轉(zhuǎn)換成中心頻率為W0、帶寬為 Bw的模擬帶通濾波器，在進(jìn)一步轉(zhuǎn)換成數(shù)字帶通濾波

22、器。程序設(shè)計如下：W0=200*2*pi;Bw=200;e3,f3=lp2bp(b,a,W0,Bw) ;h3,w5=freqs(e3,f3);figure;plot(w5/(2*pi),abs(h3);grid;title(' 模擬帶通濾波器 ');xlabel('頻率 Hz');ylabel(' 幅度 ');r3,q3=bilinear(e3,f3,1000);s3,w6=freqz(r3,q3);figure;plot(w6/pi, 20*log10(abs(s3);grid;%原型低通轉(zhuǎn)換為帶通模擬濾波器%求模擬帶通濾波器的頻率響應(yīng) % 繪

23、制帶通頻率特性%模擬濾波器轉(zhuǎn)換為數(shù)字濾波器%求數(shù)字帶通濾波器的頻率響應(yīng)%繪制數(shù)字帶通濾波器的頻率響應(yīng)title(' 數(shù)字帶通濾波器');xlabel('w/pi');ylabel(' 幅度 (dB)');武漢理工大學(xué)MATLAB課程設(shè)計報告書4.1.3對輸入信號進(jìn)行濾波讀取輸入信號程序：n=0:180;dt=1/Fs;t=n*dt;%時間序列f1=100;f2=550;%輸入信號頻率x=sin(2*pi*f1*t)+cos(2*pi*f2*t)+sin(2*pi*1*t);%輸入信號figure;plot(t,x);grid;%繪制輸入信號tit

24、le(' 輸入信號 ');xlabel('時間 t');ylabel('幅度 ');將輸入信號分別通過得到的低通、高通、帶通濾波器，調(diào)用函數(shù)y=filter(bz,az,x)對輸入信號進(jìn)行濾波。y=filter(r1,q1,x);%對輸入信號進(jìn)行低通濾波figure;plot(t,y);grid;%繪制輸出信號title(' 低通濾波輸出信號 ');xlabel('時間 t');ylabel('幅度 ')再進(jìn)行高通、帶通濾波只需將調(diào)用函數(shù)filter中參量改為 r2 、 q2 和 r3 、 q3 即

25、可。4.3仿真結(jié)果與分析將設(shè)計好的程序在matlab 軟件下進(jìn)行仿真，得到仿真波形。4.3.1仿真結(jié)果記錄武漢理工大學(xué)MATLAB課程設(shè)計報告書圖 4-1模擬低通原型幅頻響應(yīng)（a）低通模擬濾波器的頻率特性曲線（ b）巴特沃斯低通數(shù)字濾波器的頻率特性曲線圖 4-2武漢理工大學(xué)MATLAB課程設(shè)計報告書（ a）高通模擬濾波器的頻率特性曲線（ b）巴特沃斯高通數(shù)字濾波器的頻率特性曲線圖 4-3（ a）帶通模擬濾波器的頻率特性曲線（b）巴特沃斯帶通數(shù)字濾波器的頻率特性曲線圖 4-4武漢理工大學(xué)MATLAB課程設(shè)計報告書圖 4-5 輸入信號波形圖 4-6經(jīng)過巴特沃斯低通濾波器濾波后的輸出波形圖 4-7經(jīng)

26、過巴特沃斯高通濾波器濾波后的輸出波形武漢理工大學(xué)MATLAB課程設(shè)計報告書圖 4-8經(jīng)過巴特沃斯高通濾波器濾波后的輸出波形4.3.2仿真結(jié)果分析由仿真結(jié)果可知，設(shè)計出的數(shù)字低通，高通，帶通的幅頻響應(yīng)均符合各自的特性，從模擬濾波器特性曲線可以看出，所設(shè)計的高通濾波器在大于400Hz 為通帶，其衰減均小于1dB；小于 300Hz 為阻帶，其衰減大于30dB，而設(shè)計的低通濾波器情況正好相反，符合本次設(shè)計的要求。當(dāng)濾波器輸入 100Hz 和 550Hz 兩種信號后，低通濾波器可以濾除頻率為550Hz 的高頻信號，而低頻信號100Hz 的信號可以通過濾波器，驗證了該低通濾波器的功能。而高通濾波器可以濾除

27、頻率為100Hz 的信號，而高頻信號550Hz 的信號可以通過濾波器。4.3.3程序調(diào)試中出現(xiàn)的問題帶通濾波器的仿真結(jié)果雖然符合其帶通特性，但要求得到的性能指標(biāo)存在較大誤差，由于時間緊迫，將在課后進(jìn)一步尋找原因及解決方案。一個線性是不變系統(tǒng)的典型表達(dá)式有狀態(tài)空間型、傳遞函數(shù)型、零極增益型、極點留數(shù)型共四種，在開始的設(shè)計過程中我發(fā)現(xiàn)在把零極點增益型向傳遞函數(shù)型的轉(zhuǎn)換過程中多用了一步，即先把零極點增益型轉(zhuǎn)換為狀態(tài)空間型，進(jìn)行相應(yīng)的頻率變換和雙線性變換后，再由狀態(tài)空間型轉(zhuǎn)換為傳遞函數(shù)型，進(jìn)而繪制出其頻率響應(yīng)圖，但是經(jīng)仿真分析后發(fā)現(xiàn)，性能并不理想，所以并沒有采用。武漢理工大學(xué)MATLAB課程設(shè)計報告書

28、5 心得體會本次課程設(shè)計是基于數(shù)字濾波器的 MATLAB 軟件應(yīng)用，在本學(xué)期中，我在數(shù)字信號處理課程中學(xué)習(xí)了關(guān)于本次課程設(shè)計的一些內(nèi)容， 對于 MATLAB 軟件在數(shù)字濾波器的設(shè)計中的一些函數(shù)也有所認(rèn)識，通過著一周的時間更是讓自己在這個方面的學(xué)習(xí)更進(jìn)一步，使自己在這方面的學(xué)習(xí)深度有所加深。本次用雙線性變換法對巴特沃斯低通，高通及帶通濾波器的設(shè)計，要求我能夠用所學(xué)理論知識對巴特沃斯高通濾波器進(jìn)行設(shè)計，并能夠?qū)υO(shè)計的濾波器的頻率特性進(jìn)行分析。在設(shè)計過程中，首先，我通過查找資料等途徑，對關(guān)于所要設(shè)計的題目進(jìn)行更加全面深入的學(xué)習(xí)，認(rèn)真分析其中的原理，弄明白設(shè)計的思路；其次，根據(jù)自己選擇的性能指標(biāo)要求，

29、結(jié)合自己通過學(xué)習(xí)掌握的知識和資料進(jìn)行程序設(shè)計，根據(jù)自己提出的方案利用軟件實現(xiàn)；最后，運行程序，分析所得結(jié)果是否符合自己的設(shè)計要求，對結(jié)果進(jìn)行分析，實現(xiàn)由理論到實踐，再由實踐結(jié)果驗證理論的過程?？偟膩碚f，這次課程設(shè)計讓我對 MATLAB有了更深刻的了解，對數(shù)字濾波器的設(shè)計流程有了大致的了解，掌握了一些設(shè)計濾波器的基本方法，提高了理論用于實踐的能力，掌握了更多專業(yè)相關(guān)的使用知識與技能。同時，也暴露了我很多的不足，在以后的學(xué)習(xí)中，將進(jìn)一步發(fā)揚有點，克服缺點。武漢理工大學(xué)MATLAB課程設(shè)計報告書參考文獻(xiàn)1 袁潤章，自蔓延高溫合成技術(shù)研究進(jìn)展 . 武漢 : 武漢工業(yè)大學(xué)出版社 ,1994 2Schac

30、ht E.Industrial polysaccharides. Amsterdam:Elsevier Science,19873 郭景坤 . 陶瓷材料的強(qiáng)化與增韌新途徑的探索 . 無機(jī)材料學(xué)報 ,1998,13910;23-264SpriggsG E.A historyof finegrainedhardmetal.IntJ of RefractoryMetaland Hard Material,1995,13:241-2555王為民 . 二硼化鈦陶瓷的自蔓延高溫合成與加工: 博士論文 . 武漢工業(yè)大學(xué) ,19886吳中慶 , 張基田 .SiC 顆粒與液體鋁的界石層. 廣州 :92 秋季中國

31、材料科學(xué)研討會 ,C-MRS,1992.603-609.武漢理工大學(xué)MATLAB課程設(shè)計報告書附錄本次課設(shè)總程序如下：Fs=2000;Ts=1/Fs;%采樣頻率和采樣周期wp=400*2*pi;ws=300*2*pi;%低通、高通數(shù)字濾波器截止頻率Rp=1;Rs=30;%通帶波紋和阻帶衰減wpl=Fs*2*tan(ws*Ts/2);%按頻率轉(zhuǎn)換公式進(jìn)行轉(zhuǎn)換wsl=Fs*2*tan(wp*Ts/2);N,Wn=buttord(wpl,wsl,Rp,Rs,'s');%計算模擬濾波器的最小階數(shù)z,p,k=buttap(N);%設(shè)計模擬原型低通濾波器b,a=zp2tf(z,p,k);%

32、變零極點增益形式為系統(tǒng)傳遞函數(shù)形式H,w=freqs(b,a);%求模擬低通濾波器的頻率響應(yīng)figure;plot(w,abs(H);grid;%繪制低通頻率特性title(' 模擬低通原型 ');xlabel('w');ylabel(' 幅度 ');e1,f1=lp2lp(b,a,300*2*pi) ;%原型低通轉(zhuǎn)換為低通模擬濾波器h1,w1=freqs(e1,f1);%求模擬低通濾波器的頻率響應(yīng)figure;plot(w1/(2*pi),abs(h1);% 繪制高通頻率特性grid;title(' 模擬低通濾波器 ');xla

33、bel('頻率 Hz');ylabel(' 幅度 ');r1,q1=bilinear(e1,f1,1000);%模擬低通濾波器轉(zhuǎn)換為數(shù)字濾波器s1,w2=freqz(r1,q1);%求數(shù)字低通濾波器的頻率響應(yīng)figure;plot(w2/pi, 20*log10(abs(s1);grid;%繪制數(shù)字低通濾波器的頻率響應(yīng)武漢理工大學(xué)MATLAB課程設(shè)計報告書title(' 數(shù)字低通濾波器 ');xlabel('w/pi');ylabel(' 幅度 (dB)');e2,f2=lp2hp(b,a,300*2*pi) ;h

34、2,w3=freqs(e2,f2);figure;plot(w3/(2*pi),abs(h2);grid;title(' 模擬高通濾波器 ');xlabel('頻率 Hz');ylabel(' 幅度 ');r2,q2=bilinear(e2,f2,1000);s2,w4=freqz(r2,q2);figure;plot(w4/pi, 20*log10(abs(s2);grid;%原型低通模擬濾波器轉(zhuǎn)換為高通模擬濾波器%求模擬高通濾波器的頻率響應(yīng)%繪制高通頻率特性%模擬高通濾波器轉(zhuǎn)換為數(shù)字濾波器%求數(shù)字高通濾波器的頻率響應(yīng)%繪制數(shù)字高通濾波器的頻率

35、響應(yīng)title(' 數(shù)字高通濾波器 ');xlabel('w/pi');ylabel(' 幅度 (dB)');W0=200*2*pi;Bw=200;e3,f3=lp2bp(b,a,W0,Bw) ;h3,w5=freqs(e3,f3);figure;plot(w5/(2*pi),abs(h3);grid;title(' 模擬帶通濾波器 ');xlabel('頻率 Hz');ylabel(' 幅度 ');r3,q3=bilinear(e3,f3,1000);s3,w6=freqz(r3,q3);figu

36、re;plot(w6/pi, 20*log10(abs(s3);grid;%原型低通轉(zhuǎn)換為帶通模擬濾波器%求模擬帶通濾波器的頻率響應(yīng)% 繪制帶通頻率特性%模擬濾波器轉(zhuǎn)換為數(shù)字濾波器%求數(shù)字帶通濾波器的頻率響應(yīng)%繪制數(shù)字帶通濾波器的頻率響應(yīng)武漢理工大學(xué)MATLAB課程設(shè)計報告書n=0:180;dt=1/Fs;t=n*dt;%時間序列f1=100;f2=550;%輸入信號頻率x=sin(2*pi*f1*t)+cos(2*pi*f2*t)+sin(2*pi*1*t);%輸入信號figure;plot(t,x);grid;%繪制輸入信號title(' 輸入信號 ');xlabel('時間 t');ylabel('幅度 ');y=filter(r1,q1,x);%對輸入信號進(jìn)行低通濾波figure;plot(t,y);grid;%繪制輸出信號title(' 低通濾波輸出信號 ');xlabel('時間 t');ylabel('幅度 ')y=filter(r2,q2,x);%對輸入信號進(jìn)行高通濾波figure;plot(t,y);grid;%繪制輸出信號title('