橢圓帶通濾波器的設(shè)計(jì)(matlab)
《橢圓帶通濾波器的設(shè)計(jì)(matlab)》由會(huì)員分享,可在線閱讀,更多相關(guān)《橢圓帶通濾波器的設(shè)計(jì)(matlab)(17頁(yè)珍藏版)》請(qǐng)?jiān)谘b配圖網(wǎng)上搜索。
1、 燕山大學(xué) 課 程 設(shè) 計(jì) 說(shuō) 明 書(shū) 題目:橢圓帶通濾波器的設(shè)計(jì) 學(xué)院(系):電氣工程學(xué)院 年級(jí)專(zhuān)業(yè): 12級(jí) 學(xué) 號(hào): 學(xué)生姓名: 指導(dǎo)教師: 教師職稱(chēng): 電氣工程學(xué)院《課程設(shè)計(jì)》任務(wù)書(shū) 課程名稱(chēng): 數(shù)字信號(hào)處理課程設(shè)計(jì) 基層教學(xué)單位:儀器科學(xué)與工程系 指導(dǎo)教師: 學(xué)號(hào) 學(xué)生姓名 (專(zhuān)業(yè))班級(jí) 設(shè)計(jì)題目 2橢圓帶通濾波器的設(shè)計(jì) 設(shè)
2、 計(jì) 技 術(shù) 參 數(shù) 采樣頻率100Hz,采樣點(diǎn)數(shù)100,低頻、中頻、高頻信號(hào)頻率分別為5Hz、15Hz、30Hz 設(shè) 計(jì) 要 求 產(chǎn)生一個(gè)連續(xù)信號(hào),包含低頻,中頻,高頻分量,對(duì)其進(jìn)行采樣,進(jìn)行頻譜分析。設(shè)計(jì)高通濾波器對(duì)信號(hào)進(jìn)行濾波處理,觀察濾波后信號(hào)的頻譜。 (熟悉函數(shù)freqz,ellip,filter,fft) 參 考 資 料 數(shù)字信號(hào)處理方面資料 MATLAB方面資料 周次 前半周 后半周 應(yīng) 完 成 內(nèi) 容 收集消化資料、學(xué)習(xí)MATLAB軟件,進(jìn)行相關(guān)參數(shù)計(jì)算 編寫(xiě)仿真程序、調(diào)試 指導(dǎo)教 師簽字 基層教學(xué)單位主任
3、簽字 說(shuō)明:1、此表一式四份,系、指導(dǎo)教師、學(xué)生各一份,報(bào)送院教務(wù)科一份。 2、學(xué)生那份任務(wù)書(shū)要求裝訂到課程設(shè)計(jì)報(bào)告前面。 電氣工程學(xué)院 教務(wù)科 摘要 所謂數(shù)字濾波器,是指輸入輸出均為數(shù)字信號(hào),通過(guò)數(shù)值運(yùn)算處理改變輸入信號(hào)所含頻率成分的相對(duì)比例,或者濾除某些頻率成分的數(shù)字器件或程序。數(shù)字濾波器處理精度高、穩(wěn)定、體積小、重量輕、靈活、不存在阻抗匹配問(wèn)題。典型的模擬濾波器有巴特沃斯濾波器、切比雪夫?yàn)V波器和橢圓濾波器和貝塞爾濾波器,其中橢圓濾波器具有較好的性能。 Matlab是一套集數(shù)值計(jì)算、符號(hào)運(yùn)算及圖形處理等強(qiáng)大功能于一體的科學(xué)計(jì)算
4、軟件。作為強(qiáng)大的科學(xué)計(jì)算平臺(tái),它幾乎能夠滿足所有的計(jì)算需求。 本課結(jié)合MATLAB設(shè)計(jì)模擬橢圓濾波器。 目錄 第一章 概論 1 第二章 信號(hào)處理原理 2 2.1橢圓濾波器的基本理論 2 2.2采樣定理 3 第三章 軟件仿真設(shè)計(jì) 4 3.1橢圓濾波器設(shè)計(jì)結(jié)構(gòu)圖 4 3.2設(shè)計(jì)橢圓模擬帶通濾波器的步驟 4 3.3 MATLAB相關(guān)函數(shù)介紹 4 第四章 程序和仿真結(jié)果分析 8 4.1帶通通濾波器設(shè)計(jì)程序 8 4.2仿真結(jié)果KK及分析 9 參考文獻(xiàn) 13 第1章 概論 橢圓濾波器又稱(chēng)考爾濾波器。是在通帶和阻帶等波紋的一
5、種濾波器。橢圓濾波器相比其他類(lèi)型的濾波器在階數(shù)相同的條件下有著最小的通帶和阻帶波動(dòng)。它在通帶和阻帶的波動(dòng)相同,這一點(diǎn)區(qū)別于在通帶和阻帶都平坦的巴特沃斯濾波器以及通帶平坦、阻帶等波紋或是阻帶平坦、通帶等波紋的切比雪夫?yàn)V波器。 橢圓濾波器傳輸函數(shù)是一種較復(fù)雜的逼近函數(shù),利用傳統(tǒng)的設(shè)計(jì)方法進(jìn)行電路網(wǎng)絡(luò)綜合要進(jìn)行繁瑣的計(jì)算,還要根據(jù)計(jì)算結(jié)果進(jìn)行查表,整個(gè)設(shè)計(jì)、調(diào)整都十分困難和繁瑣。有許多方法都是希望能快速簡(jiǎn)便地設(shè)計(jì)并實(shí)現(xiàn)橢圓濾波器從而把電子電路設(shè)計(jì)者從煩瑣的模擬濾波器設(shè)計(jì)中解放出來(lái)。Matlab的信號(hào)處理工具箱提供了設(shè)計(jì)橢圓濾波器的函數(shù):ellipord函數(shù)和ellip函數(shù)。通過(guò)編程可以很容易由濾
6、波器的技術(shù)指標(biāo)得到所需濾波器的階數(shù)實(shí)現(xiàn)各種類(lèi)型的橢圓濾波器,大大簡(jiǎn)化了橢圓濾波器的設(shè)計(jì) 。 本文采用的方法是MATLAB設(shè)計(jì)出濾波器的傳輸函數(shù),然后再用通用的可編程濾波器來(lái)實(shí)現(xiàn)。 第二章 信號(hào)處理原理 2.1橢圓濾波器的基本理論 模擬濾波器是電子設(shè)備中最重要的部分之一。常用的濾波器有巴特沃斯濾波器、切比雪夫?yàn)V波器以及橢圓濾波器;貝塞爾濾波器巴特沃斯和切比雪夫?yàn)V波器的傳輸函數(shù)都是一個(gè)常數(shù)除以一個(gè)多項(xiàng)式,為全極點(diǎn)網(wǎng)絡(luò),僅在無(wú)限大處阻帶衰減為無(wú)限大,而橢圓函數(shù)濾波器在有限頻率上既有零點(diǎn)又有極點(diǎn)。零、極點(diǎn)在通帶內(nèi)產(chǎn)生等紋波,阻帶內(nèi)的有限傳輸零點(diǎn)減少了過(guò)渡區(qū),可
7、獲得極為陡峭的衰減曲線。也就是說(shuō)對(duì)于給定的階數(shù)和波紋要求,橢圓濾波器能獲得較其它濾波器更窄的過(guò)渡帶寬,就這點(diǎn)而言,橢圓濾波器是最優(yōu)的。 橢圓濾波器有如下特點(diǎn): 1、橢圓高通濾波器是一種零、極點(diǎn)型濾波器,它在有限頻率范圍內(nèi)存在傳輸零點(diǎn)和極點(diǎn)。 2、橢圓高通濾波器的通帶和阻帶都具有等波紋特性,因此通帶,阻帶逼近特性良好。 3、對(duì)于同樣的性能要求,它比前兩種濾波器所需用的階數(shù)都低,而且它的過(guò)渡帶比較窄。 它的振幅平方函數(shù)為 其中RN(x)是雅可比(Jacobi) 橢圓函數(shù),ε為與通帶衰減有關(guān)的參數(shù)。 2.2采樣定理 模擬信號(hào)經(jīng)過(guò) (
8、A/D) 變換轉(zhuǎn)換為數(shù)字信號(hào)的過(guò)程稱(chēng)之為采樣,信號(hào)采樣后其頻譜產(chǎn)生了周期延拓,每隔一個(gè)采樣頻率fs,重復(fù)出現(xiàn)一次。理想采樣就是假設(shè)采樣開(kāi)關(guān)閉合時(shí)間無(wú)限短,即t0的極限情況。此時(shí)采樣序列可表示為一個(gè)沖激函數(shù)序列。 采樣定理:要想采樣后能夠不失真地還原出原模擬信號(hào),則采樣頻率必須大于兩倍原模擬信號(hào)頻譜的最高截止頻率(2)。 第三章 軟件仿真設(shè)計(jì) 3.1橢圓濾波器設(shè)計(jì)結(jié)構(gòu)圖 橢圓濾波器設(shè)計(jì)結(jié)構(gòu)圖如圖所示: 采樣 (100HZ) 連續(xù)混合 信號(hào) 帶通濾波器 輸出 圖3.1橢圓濾波器結(jié)構(gòu)框圖 3.2設(shè)計(jì)橢圓模擬帶通濾波器的步驟 1. 確定模擬濾波器的性能指標(biāo)。
9、2.由性能指標(biāo)計(jì)算出濾波器階次。 3.通過(guò)歸一化及去歸一化求出的模擬濾波器。 3.3 MATLAB相關(guān)函數(shù)介紹 1. Matlab的信號(hào)處理軟件提供了設(shè)計(jì)橢圓濾波器的函數(shù):ellipord函數(shù)和ellip函數(shù)。 (1).Ellipord函數(shù)的功能是求濾波器的最小階數(shù)其調(diào)用格式為 其中:N為橢圓濾波器最小階數(shù)、Wp是橢圓濾波器通帶頻率(Hz)、Ws為橢圓濾波器阻帶頻率(Hz)、Rp為通帶波紋最大衰減(dB)、Rs為阻帶最小衰減(dB)。 在本次設(shè)計(jì)中,Wp分別為10Hz和20Hz、Ws分別為8Hz和22Hz,、Rp和Rs各為1dB和40
10、db。具體程序如下: Wp=[10,20]; Ws=[8,22]; Rp=1; Rs=40; [N,Wn]=ellipord(Wp,Ws,Rp,Rs,s); 運(yùn)行可得 Wn =10 20 (2).Ellip函數(shù)的功能是用來(lái)設(shè)計(jì)橢圓濾波器。其調(diào)用格式 其中為高通濾波器、為低通濾波器、為帶阻濾波器。 返回長(zhǎng)度為n+1的濾波器系數(shù)行向量b和a,進(jìn)而求得傳遞函數(shù)H(z) 2.頻譜分析函數(shù):fft函數(shù)、filter函數(shù)和freqz函數(shù)。 1)
11、.fft函數(shù)功能:對(duì)信號(hào)進(jìn)行離散傅里葉變換。 其調(diào)用格式:fft(X) fft(X,N) fft(X,[],DIM)或fft(X,N,DIM) 說(shuō)明:fft(X)是對(duì)輸入信號(hào)X的離散傅里葉變換。 fft(X,N)是N點(diǎn)傅里葉變換,如果X少于N點(diǎn)則補(bǔ)0湊齊位數(shù),長(zhǎng)于N點(diǎn)則截?cái)唷H绻鹸是個(gè)矩
12、陣,列的長(zhǎng)度將會(huì)以同樣的方式調(diào)整,fft會(huì)對(duì)每列進(jìn)行傅里葉變換,并返回一個(gè)相同維數(shù)的矩陣。 fft(X,[],DIM)或fft(X,N,DIM)是離散傅里葉變換在DIM尺度上的應(yīng)用。DIM可適應(yīng)于任意維度的fft運(yùn)算。 2).filter函數(shù)功能:利用IIR濾波器和FIR濾波器對(duì)數(shù)據(jù)進(jìn)行濾波。 其調(diào)用格式: y=filter(b,a,x) [y,zf]=filter(b,a,x) y=filter(b,a,x,zi) 說(shuō)明:filter采用數(shù)字濾波器對(duì)數(shù)據(jù)進(jìn)行濾波,其實(shí)現(xiàn)采用移位直接Ⅱ型結(jié)構(gòu),因而適用于IIR和FIR濾波器。濾波器
13、的系統(tǒng)函數(shù)為 即濾波器系數(shù)a=[a0 a1 a2 ...an],b=[b0 b1 ...bm],輸入序列矢量為x。這里,標(biāo)準(zhǔn)形式為a0=1,如果輸入矢量a時(shí),a0≠1,則MATLAB將自動(dòng)進(jìn)行歸一化系數(shù)的操作;如果a0=0則給出出錯(cuò)信息。 y=filter(b,a,x)利用給定系數(shù)矢量a和b對(duì)x中的數(shù)據(jù)進(jìn)行濾波,結(jié)果放入y矢量中,y的長(zhǎng)度取max(N,M)。 y=filter(b,a,x,zi)可在zi中指定x的初始狀態(tài)。 [y,zf]=filter(b,a,x)除得到矢量y外,還得到x的最終狀態(tài)矢量
14、zf。 3).freqz函數(shù)功能:離散時(shí)間系統(tǒng)的頻率響應(yīng)。 其調(diào)用格式: [h,w]=freqz(b,a,n) [h,f]=freqz(b,a,n,Fs) h=freqz(b,a,w) h=freqz(b,a,f,Fs) 說(shuō)明: freqz 用于計(jì)算數(shù)字濾
15、波器H(Z)的頻率響應(yīng)函數(shù)H(ejω)。 [h,w]=freqz(b,a,n)可得到數(shù)字濾波器的n點(diǎn)幅頻響應(yīng)值,這n個(gè)點(diǎn)均勻地分布在[0,π]上,并將這n個(gè)頻點(diǎn)的頻率記錄在w中,相應(yīng)的頻響值記錄在h中。要求n為大于零的整數(shù),最好為2的整數(shù)次冪,以便采用FFT計(jì)算,提高速度。缺省時(shí)n =512。 [h,f]=freqz(b,a,n,Fs)用于對(duì)H(ejω)在[0,Fs/2]上等間隔采樣n點(diǎn),采樣點(diǎn)頻率及相應(yīng)頻響值分別記錄在f 和h中。由用戶指定FS(以HZ為單位)值。 h=freqz(b,a,w)用于對(duì)H(ejω)在[0,2π]上進(jìn)行采樣,采樣頻率點(diǎn)由矢量w指定。 h=fr
16、eqz(b,a,f,Fs) 用于對(duì)H(ejω)在[0,FS]上采樣,采樣頻率點(diǎn)由矢量f指定。 freqz(b,a,n) 用于在當(dāng)前圖形窗口中繪制幅頻和相頻特性曲線。 第四章 程序和仿真結(jié)果分析 4.1帶通通濾波器設(shè)計(jì)程序 Wp=[10,20]; Ws=[8,22]; Rp=1; Rs=40; [N,Wn]=ellipord(Wp,Ws,Rp,Rs,s); Fs=100; q=(0:255)/256*50; t=(1:100)/Fs; c=sin(2*pi*t*5)+sin(2*pi*t*15)+sin(2*pi*t*30); figure(1);
17、 [b,a]=ellip(N,Rp,Rs,Wn/50); [H,w]=freqz(b,a,512); plot(w*Fs/(2*pi),abs(H)); title(帶通濾波器頻率特性); xlabel(頻率(HZ)); ylabel(頻率響應(yīng)圖); sf=filter(b,a,c); figure(2); subplot(211); plot(t,c); xlabel(時(shí)間(秒)); ylabel(幅值); title(連續(xù)信號(hào)); subplot(212); plot(t,sf); title(通過(guò)濾波器后信號(hào)波形); xlabel(時(shí)間(秒));
18、 ylabel(幅值); cc=fft(c,512); sff=fft(sf,512); figure(3); subplot(211); plot(q,abs(cc(1:256))); title(原連續(xù)信號(hào)頻譜特性); xlabel(頻率(HZ)); ylabel(頻率響應(yīng)); subplot(212); plot(q,abs(sff(1:256))); title(通過(guò)濾波器后信號(hào)頻譜特性); xlabel(頻率(HZ)); ylabel(頻率響應(yīng)); 4.2仿真結(jié)果KK及分析 1. 濾波器的頻率特性圖 圖4.1帶通濾波器頻率特性 結(jié)果
19、分析:從仿真圖可以看出橢圓濾波器的過(guò)渡帶非常窄,且在通帶和阻帶都有波動(dòng),符合橢圓濾波器的特點(diǎn)。 2. 通過(guò)帶通濾波器前后的波形圖 圖4.2通過(guò)濾波器前后波形圖對(duì)比 結(jié)果分析:由圖可以看出,通過(guò)橢圓濾波器后的信號(hào)成為頻率單一的正弦信號(hào),信號(hào)的頻率為15Hz。只通過(guò)了中頻信號(hào),因此實(shí)現(xiàn)了帶通濾波器的功能。 3. 信號(hào)通過(guò)帶通濾波器前后的頻譜特性 圖4.3通過(guò)濾波器前后信號(hào)的頻譜圖 結(jié)果分析:通過(guò)帶通濾波器前,低、中、高頻信號(hào)都有頻率響應(yīng)。而通過(guò)濾波器之后,只保留的中頻信號(hào),因此低
20、、高頻信號(hào)的頻率響應(yīng)為零。除此之外,通過(guò)濾波器后信號(hào)的響應(yīng)幅值有所減小。 第五章 心得體會(huì) 《數(shù)字信號(hào)處理》是一門(mén)重要的專(zhuān)業(yè)基礎(chǔ)課程。通過(guò)這周的課程設(shè)計(jì)我加深了對(duì)數(shù)字信號(hào)處理這門(mén)學(xué)科所學(xué)理論知識(shí)的理解。為了完成這次橢圓帶通濾波器的設(shè)計(jì),我花了許多時(shí)間去學(xué)習(xí)matlab這個(gè)軟件,也進(jìn)一步理解了matlab這一軟件的強(qiáng)大之處。相信這會(huì)對(duì)以后的工作和學(xué)習(xí)有很大的幫助。 這次課程設(shè)計(jì)讓我把理論知識(shí)與實(shí)際應(yīng)用結(jié)合了起來(lái),提高自己的實(shí)際動(dòng)手能力和獨(dú)立思考的能力。通過(guò)不斷查閱大量的設(shè)計(jì)資料讓自己的設(shè)計(jì)更加完善
21、,這加強(qiáng)了我對(duì)信息的獲取能力。 總之,這次課程設(shè)計(jì)是自己的能力得到了一定的鍛煉。在這次課程設(shè)計(jì)中,我也認(rèn)識(shí)到了自己的一些不足,這些都會(huì)對(duì)以后成長(zhǎng)有很大的好處。 參考文獻(xiàn) [1]謝平 王娜 林洪彬.信號(hào)處理及應(yīng)用.機(jī)械工業(yè)出版社 [2]王靖 李永全.橢圓濾波器Matlb設(shè)計(jì)與實(shí)現(xiàn)[J].現(xiàn)代電子技術(shù),2007,(6) [3]高西全丁玉美. 數(shù)字信號(hào)處理. 第三版. 西安西安電子科技大學(xué)出版社2008 [4]彭軍李宏. 信號(hào)與信息處理基礎(chǔ). 北京中國(guó)鐵道出版社2009 [5]A. V. 奧本海姆等. 信號(hào)與系統(tǒng). 劉樹(shù)棠譯. 西安西安交通大學(xué)出版社1985 [6]董長(zhǎng)虹. MATLAB信號(hào)處理與應(yīng)用. 北京國(guó)防工業(yè)出版社2005
- 溫馨提示:
1: 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
2: 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
3.本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
5. 裝配圖網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 2025年防凍教育安全教育班會(huì)全文PPT
- 2025年寒假安全教育班會(huì)全文PPT
- 初中2025年冬季防溺水安全教育全文PPT
- 初中臘八節(jié)2024年專(zhuān)題PPT
- 主播直播培訓(xùn)提升人氣的方法正確的直播方式如何留住游客
- XX地區(qū)機(jī)關(guān)工委2024年度年終黨建工作總結(jié)述職匯報(bào)
- 心肺復(fù)蘇培訓(xùn)(心臟驟停的臨床表現(xiàn)與診斷)
- 我的大學(xué)生活介紹
- XX單位2024年終專(zhuān)題組織生活會(huì)理論學(xué)習(xí)理論學(xué)習(xí)強(qiáng)黨性凝心聚力建新功
- 2024年XX單位個(gè)人述職述廉報(bào)告
- 一文解讀2025中央經(jīng)濟(jì)工作會(huì)議精神(使社會(huì)信心有效提振經(jīng)濟(jì)明顯回升)
- 2025職業(yè)生涯規(guī)劃報(bào)告自我評(píng)估職業(yè)探索目標(biāo)設(shè)定發(fā)展策略
- 2024年度XX縣縣委書(shū)記個(gè)人述職報(bào)告及2025年工作計(jì)劃
- 寒假計(jì)劃中學(xué)生寒假計(jì)劃安排表(規(guī)劃好寒假的每個(gè)階段)
- 中央經(jīng)濟(jì)工作會(huì)議九大看點(diǎn)學(xué)思想強(qiáng)黨性重實(shí)踐建新功