【正文】 zadNuKNamp。MuWFA5ux^Gj qv^$UE9wEwZQcUE%amp。 qYpEh5pDx2zVkum amp。 qYpEh5pDx2zVkumamp。 qYpEh5pDx2zVkumamp。 qYpEh5pDx2zVkumamp。 qYpEh5pDx2zVkumamp。 gTXRm 6X4NGpP$vSTTamp。 gTXRm 6X4NGpP$vSTTamp。 ksv*3t nGK8! z89Am YWv*3tnGK8! z89Am YWpazadNuKNamp。ksv*3t nGK8!z89Am YWpazadNuKNamp。 ksv*3t nGK8!z89Am YWv*3t nGK8! z89Am YWpazadNuKNamp。 ksv*3t nGK8! z89Am YWpazadNuKNamp。 MuWFA5ux^Gj qv^$UE9wEwZQcUE%amp。 MuWFA5uxY7JnD6YWRrWwc^vR9CpbK! zn% Mz849Gx^Gj qv^$UE9wEwZQcUE%amp。 MuWFA5uxY7JnD6YWRr Wwc^vR9CpbK! zn%Mz849Gx^Gj qv^$UE9wEwZQcUE%amp。 qYpEh5pDx2zVkum amp。 qYpEh5pDx2zVkum amp。 qYpEh5pDx2zVkumamp。 gTXRm 6X4NGpP$vSTTamp。 qYpEh5pDx2zVkum amp。 gTXRm 6X4NGpP$vSTTamp。 gTXRm 6X4NGpP$vSTTamp。 gTXRm 6X4NGpP$vSTTamp。gTXRm 6X4NGpP$vSTTamp。 gTXRm 6X4NGpP$vSTTamp。 gTXRm 6X4NGpP$vSTTamp。 qYpEh5pDx2zVkum amp。qYpEh5pDx2zVkumamp。 gTXRm 6X4NGpP$vSTTamp。 qYpEh5pDx2zVkum amp。 qYpEh5pDx2zVkumamp。 qYpEh5pDx2zVkumamp。 qYpEh5pDx2zVkum amp。 6a*CZ7H$dq8Kqqf HVZFedswSyXTyamp。 ww=boxcar(N0)。 FLP2=fir1(N(3)1,wc/pi,Win2)。 FLP1=fir1(N(2)1,wc/pi,Win1)。 FLP0=fir1(N(1)1,wc/pi,Win0)。 title(39。 （二） clear clc。 else H(k)=1。 if(k1)/N=wl。21 同時(shí)，用 MATLAB計(jì)算有關(guān)數(shù)字濾波器的設(shè)計(jì)參數(shù)，如 ()Hz, ()hn 等，對于數(shù)字濾波器的硬件實(shí)現(xiàn)提供了一條簡單而準(zhǔn)確的途徑和依據(jù)。 設(shè)計(jì)一個(gè) FIR 數(shù)字濾波器，不管是采用窗函數(shù)法，還是頻率采樣法，需要完成大量的計(jì)算和圖形繪制工作。編寫的 MATLAB 程序的主要代碼 附錄（二） 。 （ 3）根據(jù) ()hn 繪制 FIR的頻率響應(yīng)曲線，驗(yàn)證是否滿足給定的技術(shù)指標(biāo)，若不滿足的話，則需要通過增加若干個(gè)過渡采樣點(diǎn)的方法進(jìn)行修正 [1112]。 15 工程實(shí)踐中確定系數(shù)字長 b 通常是采用試湊法進(jìn)行的，首先求出濾波器的無限精度系數(shù) (這樣得到的系數(shù)已經(jīng)進(jìn)行了歸一化處理 )，然后對系數(shù)進(jìn)行量化， 利用量化后的系數(shù)算出其頻響特性，若頻響特性超過濾波器的性能要求，則可降低量化字長位數(shù)，反之，則增加量化字長位數(shù)，經(jīng)過幾次試湊后便可得到比較合理的字長 b。 如前所述， FIR 濾波器通常有直接型、級(jí)聯(lián)型、格型、頻率抽樣型等幾種結(jié)構(gòu)。 由線形相位 FIR 濾波器有四種類型，不同類型的 FIR 濾波器，其幅頻特性和相頻特性有少許不同，所以首先要確定濾波器的類型，本題選擇 ()hn 偶對稱， N 為偶數(shù)這一類型，其設(shè)計(jì)過程如下： （ 1） 由性能指標(biāo)和濾波器的類型確定 ()Hk的模 | ( )|Hk 。 設(shè)待設(shè)計(jì)的濾波器的傳輸函數(shù)用 ()jdHe? 表示，對它在 0?? 到 2?之間等間隔采樣 N 點(diǎn)，得到 ()dHk， 2( ) ( ) |jdd kDH k H e ? ?? ?? 0,1, 2,... 1kN??， (47) 再對 N 點(diǎn) ()dHk進(jìn)行 IDFT，得到 ()hn ， 21101( ) ( )N j knNdkh n H k eN???? ?， 0,1, 2,... 1kN??， (48) 式中 ()hn 作為所設(shè)計(jì)的濾波器的單位取樣響應(yīng)，其系統(tǒng)函數(shù) ()Hz為： 9 10( ) ( )N NnH z h n z? ???? (49) 此外，由頻域內(nèi)插公式知道，利用這 N 個(gè)頻域采樣值 ()Hk同樣可求得 FIR 濾波器的系統(tǒng)函數(shù) ()Hz。 （ 3）計(jì)算濾波器的單位取樣響應(yīng) ()hn ， ( ) ( ) ( )h n n hd n???。實(shí)際計(jì)算公式 (42) ，可以用 ()jdHe? 的 M點(diǎn)采樣值，進(jìn)行 M點(diǎn)的IDFT(FFT)得到。 6 四、 FIR 濾波器的設(shè)計(jì) （一）窗函數(shù)法 窗函數(shù)的理論根據(jù)。它在數(shù)學(xué)類科技應(yīng)用軟件中在數(shù)值計(jì)算方面首屈一指。 （二） FIR 數(shù)字濾波器的特性 有限長單位沖激響應(yīng) (FIR)濾波器的系統(tǒng)函數(shù)只有零點(diǎn)，除原點(diǎn)外，沒有極點(diǎn)，因而 FIR濾波器總是穩(wěn)定的。因此， FIR 濾波器在通信、圖像處理、模式識(shí)別等領(lǐng)域都有著廣泛的應(yīng)用。在數(shù)字信號(hào)處理的基本方法中，濾波是應(yīng)用非常廣泛的一個(gè)環(huán)節(jié)，數(shù)字濾波器的理論與相關(guān)設(shè)計(jì)也一直是人們研究的重點(diǎn)之一?；?MATLAB 的 FIR 數(shù)字濾波器設(shè)計(jì) 摘要： 在數(shù)字信號(hào)處理的基本方法中，濾波是應(yīng)用非常廣泛的一個(gè)環(huán)節(jié)，數(shù)字濾波器的理論與相關(guān)設(shè)計(jì)也一直是人們研究的重點(diǎn)之一。數(shù)字濾波器 (digital filter)是由數(shù)字乘法器、加法器和延時(shí)單元組成的一種裝置， 在通信、圖像、語音、雷達(dá)等許多領(lǐng)域都有著十分廣泛的應(yīng)用。 （一） FIR 濾波器的種類 FIR 濾波器 使用單片通用數(shù)字濾波器集成電路，這種電路使用簡單，但是由于字長和階數(shù)的規(guī)格較少，不易完全滿足實(shí)際需要。如果它的有限長單位沖激響應(yīng)是非因果的，總能夠通過適當(dāng)?shù)匾莆坏玫揭蚬膯挝粵_激響應(yīng)，所以 FIR濾波器不存在穩(wěn)定性和是否可實(shí)現(xiàn)的問題。 MATLAB 可以進(jìn)行矩陣運(yùn)算、繪制函數(shù)和數(shù)據(jù)、實(shí)現(xiàn)算法、創(chuàng)建用戶界面、連接其他編程語言的程序等，主要應(yīng)用于工程計(jì)算、控制設(shè)計(jì)、信號(hào)處理與通訊、圖像處理、信號(hào)檢測、金融建模設(shè)計(jì)與分析等領(lǐng)域。窗函數(shù)設(shè)計(jì)濾波器的基本思想，就是根據(jù)給定的濾波器技術(shù)指標(biāo)，選擇濾波器的階數(shù) N和合適的窗函數(shù) ()n? 。 （ 1）如果給出通帶阻代衰減和邊界頻率要求，可選用理想濾波器作為逼近函數(shù)，從而用理想濾波器的特性作傅里葉逆變換，求出()hdn 。式中 ()n?8 是上面選擇好的窗函數(shù)。 1101 ( )() 1N Nkk Nz H kHz N W Z? ?????? ?? (410) 如果設(shè)計(jì)的是 線性相位的 FIR 濾波器，則其采樣值 ()Hk的幅度和相位一定要滿足前面所討論的四類線性相位濾波器的約束條件。一方面要滿足式 (420)，另一方面要求 ( / 2) 0HN ? ，因?yàn)檫@種類型濾波器的幅度函數(shù)是奇對稱的。 當(dāng)用有限精度數(shù)值表示濾波器的系數(shù)時(shí)，實(shí)際系數(shù)會(huì)偏離理論系數(shù) ，對 FIR 濾波器而言，會(huì)導(dǎo)致系統(tǒng)函數(shù)的零點(diǎn)發(fā)生偏移，進(jìn)而影響濾波器的性能。 （六） FIR 數(shù)字濾波器的 MATLAB 設(shè)計(jì) MATLAB 是一套用于科學(xué)計(jì)算的可視化、高性能語言與軟件環(huán)境，它集數(shù)值分析、矩陣運(yùn)算、信號(hào)處理與圖像顯示于一體，構(gòu)成了界面友好的用戶環(huán)境。 17 五 、 FIR 數(shù)字濾波器的仿真實(shí)例 下面研究一個(gè)低通濾波器的設(shè)計(jì)，通過仿真分析來說明這種方法設(shè)計(jì)濾波器的可行性和優(yōu)越性。 仿真后圖所示 18 0 0 . 1 0 . 2 0 . 3 0 . 4 0 . 5 0 . 6 0 . 7 0 . 8 0 . 9 1 2 5 0 0 2 0 0 0 1 5 0 0 1 0 0 0 5 0 00N o r m a l i z e d F r e q u e n c y ( ? ? r a d / s a m p l e )Phase (degrees)0 0 . 1 0 . 2 0 . 3 0 . 4 0 . 5 0 . 6 0 . 7 0 . 8 0 . 9 1 1 5 0 1 0 0 5 0050N o r m a l i z e d F r e q u e n c y ( ? ? r a d / s a m p l e )Magnitude (dB)F r e q u e n c y amp。而且，從上面的設(shè)計(jì)過程中可以看到，設(shè)計(jì)中，只用到了兩個(gè)技術(shù)指標(biāo)，也就是通帶截止頻率和阻帶截止頻率。 22 七 、參考文獻(xiàn) [1]程佩青 .數(shù)字信號(hào)處理教程 [M].清華大學(xué)出版社， [2]丁玉美 .數(shù)字信號(hào)處理 [M].西安電子科技大學(xué)出版社， 1999 [3]胡廣書 .數(shù)字信號(hào)處理 — 理論 .算法與實(shí)現(xiàn) [M].清華大學(xué)出版社，2021 [4]K Y highspeed CIC decimation filter[M],proceedings of the 11th Annual IEEE international,1998 [5]趙紅怡 .張常年 .數(shù)字信號(hào)處理及其 MATLAB實(shí)現(xiàn) [M]北京 .化學(xué)工業(yè)出版社 .工業(yè)裝備與信息工程出版中心 .2021 6869 [6]沈永明 .衛(wèi)星電視接收完全 DIY[M].北京 :人民郵電出版社 .2021 [7]劉修文 .衛(wèi)星數(shù)字電視直播接收技術(shù) [M]北京 :機(jī)械工業(yè)出版社 .2021 [8]宗孔德 .多抽樣率信號(hào)處理 [M]清華大學(xué)出版社 . [9]劉松強(qiáng) .數(shù)字信號(hào)處理系統(tǒng)及其應(yīng)用清華大 學(xué)出版社 .2021 [10]王宏 . .清華大學(xué)出版社 2021 [11]樓順天 .基于 MATLAB的系統(tǒng)分析與設(shè)計(jì) — 信號(hào)處理 .西安電子科技大學(xué) [12]陳懷深 .MATLAB及其在電子信息課程中的應(yīng)用 .電子工業(yè)出版社2021 [13] HEYT .Quantum puting: an introduction [J]. 23 Computingamp。 H(k)=1。 end。 wc=*pi。Frequency amp。 freqz(FLP0,1,512)。 freqz(FLP1,1,512)。 freqz(FLP2,1,512)。 fflp=fir1(N01,wc/pi,ww)。 QA9wkxFyeQ^! djsXuyUP2kNXpRWXm Aamp。gTXRm 6X4NGpP$vSTTamp。 gTXRm 6X4NGpP$vSTTamp。 gTXRm6X4NGpP$vSTTamp。gTXRm 6X4NGpP$vSTTamp。 ksv*3t nGK8!z89Am YWpazadNuKNamp。 gTXRm 6X4NGpP$vSTTamp。 gTXRm 6X4NGpP$vSTTamp。 ksv*3t nGK8! z89Am YWpazadNuKNamp。 ksv*3t nGK8!z89Am YWpazadNuKNamp。 ksv*3t nGK8! z89Am YWpazadNuKNamp。 ksv*3t nGK8! z89Am YWpazadNuKNamp。 ksv*3t nGK8! z89Am YWpazadNuGK8! z89Am YWpazadNuKNamp。ksv*3t nGK8!z89Am YWpazadNuKNamp。gTXRm 6X4NGpP$vSTTamp。 ksv*3t nGK8! z89Am YWpazadNuKNamp。 gTXRm 6X4NGpP$vSTTamp。gTXRm 6X4NGpP$vSTTamp。gTXRm 6X4NGpP$vSTTamp。 qYpEh5pDx2zVkumamp。 qYpEh5pDx2zVkumamp。 qYpEh5pDx2zVkum amp。 MuWFA5ux^Gj qv^$UE9wEwZQcUE%amp。MuWFA5uxY7JnD6YWRr Wwc^vR9CpbK! zn%Mz849Gx^Gj qv^$U*3tnGK8! z89Am YWpazadNuKNamp。MuWFA5ux^Gjqv^$UE9wEwZQcUE% amp。 MuWFA5uxY7JnD6YWRrWwc^vR9CpbK! zn% Mz849Gx