利用matlab結(jié)合頻率取樣法設(shè)計(jì)數(shù)字高通fir濾波器課程設(shè)計(jì)任務(wù)書-資料下載頁
2025-07-20 11:48本頁面
【正文】 de39。,39。manual39。,39。YTickLabels39。,[39。5039。39。 039。])% 調(diào)試結(jié)果 圖5 頻率采樣技術(shù)：高通，最優(yōu)法結(jié)果分析：第一幅圖為要高通濾波器原型，可以看到它在過渡帶添加了兩個(gè)采樣點(diǎn)，以增加阻帶衰減；第二幅圖為系統(tǒng)函數(shù)單位脈沖響應(yīng)圖形，可以看出，它以中點(diǎn)成偶對(duì)稱，由于采樣點(diǎn)數(shù)為奇數(shù)，故在對(duì)稱軸處有取值；第三幅圖（左下）為根據(jù)頻率取樣法設(shè)計(jì)出的濾波器振幅響應(yīng)，可以看出它在采樣點(diǎn)處的取值與原高通濾波器精確一致，在其他點(diǎn)處與原高通濾波器取值逼近有上下波動(dòng)；第四幅圖為用分貝數(shù)表示的幅度響應(yīng)，可以看到采用線性最優(yōu)法設(shè)計(jì)的高通濾波器的阻帶衰減大于50db。設(shè)計(jì)取得了良好的效果。4 心得體會(huì)Matlab的課程設(shè)計(jì)做到現(xiàn)在已經(jīng)基本接近尾聲了，既然學(xué)習(xí)一門課程，簡(jiǎn)單的總結(jié)是必須要有的。以前在《信號(hào)與系統(tǒng)》和《數(shù)字信號(hào)處理》的實(shí)驗(yàn)中已經(jīng)接觸過matlab，所以上手并不是很難，不過在設(shè)計(jì)的時(shí)候還是遇到了不少問題，首先是對(duì)頻率取樣法掌握的不到位，重新學(xué)習(xí)了頻率取樣法后，發(fā)現(xiàn)如何利用程序?qū)崿F(xiàn)頻率取樣法成了一個(gè)問題。通過自己在網(wǎng)上查找資料，看從圖書館借來的書以及對(duì)照著老師的PPT，不斷的調(diào)試，終于做出了成果。課程設(shè)計(jì)雖然做完了，但現(xiàn)在學(xué)的這點(diǎn)知識(shí)還遠(yuǎn)遠(yuǎn)不夠，特別是這個(gè)軟件的函數(shù)非常多，要能夠熟練運(yùn)用我們還有很多要學(xué)習(xí)。不過我覺得Matlab的函數(shù)設(shè)計(jì)都比較合理，她總是從函數(shù)本身的意義出發(fā)命名，這使我們記不會(huì)很難。 總之這次課程設(shè)計(jì)完成的還算順利，雖然也遇到過一些問題，但通過和同學(xué)討論一起學(xué)習(xí)都能解決。當(dāng)然，我們也都明白matlab的確是一個(gè)很實(shí)用的工具，在今后的學(xué)習(xí)中我們會(huì)不斷的邊學(xué)邊運(yùn)用它，而且我們還可以將它用在我們專業(yè)的學(xué)習(xí)中。 5 參考文獻(xiàn)[1] 劉 泉，數(shù)字信號(hào)處理原理與實(shí)現(xiàn)，電子工業(yè)出版社，2009[2] 郭仕劍，人民郵電出版社，2007[3] 陳懷琛，MATLAB及在電子信息課程中的應(yīng)用，電子工業(yè)出版社，2006[4] 高會(huì)生，MATLAB實(shí)用教程(第2版)，電子工業(yè)出版社，2010 [5] 陳懷琛，數(shù)字信號(hào)處理教程——MATLAB釋義與實(shí)現(xiàn)，電子工業(yè)出版社，2004 附錄 輔助函數(shù)1function [Hr,w,a,L] = Hr_Type1(h)。% 計(jì)算第一種低通濾波器設(shè)計(jì)的振幅響應(yīng)Hr(w) % % [Hr,w,a,L] = Hr_Type1(h)% Hr = 振幅響應(yīng)% w = 在[0 pi] 區(qū)間內(nèi)計(jì)算Hr 的500個(gè)頻率點(diǎn)% a = 第一種低通濾波器的系數(shù)% L = Hr的階次% h = 第一種低通濾波器的頻率響應(yīng)% M = length(h)。 L = (M1)/2。 a = [h(L+1) 2*h(L:1:1)]。 % 1乘(L+1) 行向量row vector n = [0:1:L]。 % (L+1)乘1 列向量 w = [0:1:500]39。*pi/500。Hr = cos(w*n)*a39。輔助函數(shù)2function [db,mag,pha,grd,w] = freqz_m(b,a)。% freqz 子程序的改進(jìn)版本% % [db,mag,pha,grd,w] = freqz_m(b,a)。% db = [0 到pi弧度]區(qū)間內(nèi)的相對(duì)振幅(db)% mag = [0 到pi弧度]區(qū)間內(nèi)的絕對(duì)振幅% pha = [0 到pi弧度]區(qū)間內(nèi)的相位響應(yīng)% grd = [0 到pi弧度]區(qū)間內(nèi)的群遲延% w = [0 到pi弧度]區(qū)間內(nèi)的501個(gè)頻率樣本向量% b = Ha(z)的分子多項(xiàng)式系數(shù)(對(duì)FIR b=h)% a = Ha(z)的分母多項(xiàng)式系數(shù)(對(duì) FIR: a=[1])%[H,w] = freqz(b,a,1000,39。whole39。)。 H = (H(1:1:501))39。 w = (w(1:1:501))39。 mag = abs(H)。 db = 20*log10((mag+eps)/max(mag))。 pha = angle(H)。% pha = unwrap(angle(H))。 grd = grpdelay(b,a,w)。% grd = diff(pha)。% grd = [grd(1) grd]。% grd = [0 grd(1:1:500)。 grd。 grd(2:1:501) 0]。% grd = median(grd)*500/pi。
范文總結(jié)相關(guān)推薦
鄂ICP備17016276號(hào)-1