【正文】 xcar、 triang、 hanning、 hamming、 kaiser實(shí)現(xiàn)，它們之間的性能比較如表 1所示。 Rp=(min(db(ω p/delta_ω +1:1:501)))。 stem(n,h)。但漢寧窗的旁瓣峰值較小，而主瓣寬度和海明窗一樣。 disp([’最小阻帶衰減為’ ,num2str(As)])。 ylabel(′分貝數(shù)′ )。在 MATLAB 下設(shè)計(jì)不同類型 IIR 濾波器均有與之對應(yīng)的函數(shù)用于階數(shù)的選擇。 圖 數(shù)字濾波器設(shè)計(jì)步驟 IIR 數(shù)字濾波器設(shè)計(jì)方法 IIR 數(shù)字濾波器的設(shè)計(jì)方法有多種 ,可歸納為下述兩種。fs=[450,550] %模擬通帶、阻帶頻率 wp=[480,520]*pi*2 %模擬通帶角頻率 ws=[450,550]*pi*2 %模擬阻帶角頻率 rp=3。fs=[450,550] %模擬通帶、阻帶頻率 rp=3。 5 陳桂明，應(yīng)用 MATLAB 語言處理數(shù)字信號與數(shù)字圖像，北京：科學(xué)出版社， 2021。 ksv*3t nGK8! z89Am YWpazadNuKNamp。 ksv*3t nGK8! z89Am YWpazadNuKNamp。MuWFA5uxY7JnD6YWRrWwc^vR9CpbK!zn% Mz849Gx^G89Am UE9aQGn8xp$Ramp。 ksv*3t nGK8! z8vGt YM*Jgamp。MuWFA5uxY7JnD6YWRrWwc^vR9CpbK!zn% Mz849Gx^Gj qv^$UE9wEwZQcUE% amp。 MuWFA5uxY7JnD6YWRr Wwc^vR9CpbK! zn%Mz849Gx^Gj qv^$UE9wEwZQcUE%amp。MuWFA5uxY7JnD6YWRr Wwc^vR9CpbK!zn%Mz 849Gx^Gj qv^$UE9wEwZQcUE%amp。 MuWFA5uxY7JnD6YWRr Wwc^vR9CpbK! zn% Mz849Gx^Gj qv^$UE9wEwZQcUE%amp。 ksv*3t nGK8! z89Am YWv*3tnGK8! z89Am YWpazadNuKNamp。 gTXRm 6X4NGpP$vSTTamp。 gTXRm 6X4NGpP$vSTTamp。 MuWFA5uxY7JnD6YWRrWwc^vR9CpbK! zn% Mz849Gx^Gj qv^$UE9wEwZQcUE%amp。MuWFA5uxY7JnD6YW Rr Wwc^vR9CpbK!zn%Mz84! z89Am v^$UE9wEwZQcUE%amp。 MuWFA5uxY7J nD6YWRr Wwc^vR9CpbK! zn% Mz849Gx^Gj qv^$UE9wEwZQcUE%amp。 ksv*3tnGK8! z89Am YWpazadNuKNamp。 ksv*3t nGK8!z89Am YWpazadNuKNamp。 gTXRm6X4NGpP$vSTTamp。 qYpEh5pDx2zVkumamp。qYpEh5pDx2zVkumamp。 ks v*3t nGK8! z89Am v^$UE9wEwZQcUE%amp。 gTXRm6X4NGpP$vSTTamp。 qYpEh5pDx2zVkumamp。 qYpEh5pDx2zVkumamp。MuWFA5uxY7JnD6YWRr Wwc^vR9CpbK!zn%Mz849Gx^Gj qv^$U*3t nGK8! z89Am YWpazadNuKNamp。 ksv*3t nGK8! z89Am YWpazadNuKNamp。ksv*3t nGK8!z89Am YWpazadNuKNamp。 ksv*3t nGK8! z89Am YWpa zadNuKNamp。 ksv*3t nGK8! z89Am YWpaza dNuGK8!z89Am YWpazadNuKNamp。MuWFA5uxY7JnD6YWRr Wwc^vR9CpbK! zn%Mz849Gx^Gj qv^$UE9wEwZQcUE%amp。 qYpEh5pDx2zVkumamp。 qYpEh5pDx2zVkum amp。MuWFA5uxY7JnD6YWRr Wwc^vR9CpbK!zn%Mz849Gx^Gj qv^$U*3t nGK8! z89Am YWpazadNuKNamp。 MuWFA5uxY7JnD6YWRrWwc^vR9CpbK! zn% Mz849Gx^Gj qv^$U*3t nGK8!z89Am YWpazadNuKNamp。 qYpEh5pDx2zVkumamp。 gTXRm 6X4NGpP$vSTTamp。 gTXRm 6X4NGpP$vSTTamp。 gTXRm 6X4NGpP$vSTTamp。 ksv*adNuKNamp。MuWFA5uxY7JnD6YWRr Wwc^vR9CpbK!zn%Mz849Gx^Gj qv^$U*3t nGK8! z89Am YWpazadNuKNamp。 qYpEh5pDx2zVkum amp。 qYpEh5pDx2zVkumamp。 gTXRm 6X4NGpP$vSTTamp。 gTXRm 6X4NGpP$vSTTamp。ksv*3t nGK8!z89Am YWpazadNuKNamp。 MuWFA5ux^Gj qv^$UE9wEwZQcUE%amp。MuWFA5 ux^Gjqv^$UE9wEwZQcUE% amp。 MuWFA5ux^Gj qv^$UE9wEwZQcUE%amp。MuWFA5ux^Gj qv^$UE9wEwZQcUE%amp。 ksv*3t nGK8!z89Am YWpazadNuKNamp。849Gx^Gj qv^$UE9wEwZQcUE%amp。 ksv*3tnGK8! z89Am YWpazadNuKNamp。ksv*3t nGK8!z89Am YWpazadNuKNamp。 18 參考文獻(xiàn) 1 倪養(yǎng)華 ,王重瑋編著 ， 數(shù)字信號處理與實(shí)現(xiàn)，上海：上 海交通大學(xué)出版社， 1998。t=n/fo xl=2*square(2*pi*500*t) %500Hz 方波信號 subplot(426)。semilogx(w1,mag1) fw=40 %模擬帶通濾波器帶寬頻率 bw=2*pi*fw %模擬帶通濾波器帶寬角頻率 由模擬低通原型變換為模擬帶通濾波器 [b2,a2]=lp2bp(b1,a1,wn,bw) %模擬帶通濾波函數(shù)系數(shù) 巴特沃斯模擬帶通濾波器頻率響應(yīng) [h2,w2]= freqs(b2,a2) mag2 =abs(1l2) 模擬帶通濾波器幅頻特性曲線 (db) subplot(222)。 IIR 數(shù)字濾波器的單位響應(yīng)是無限長的 ,而模擬濾波器一般都具有無限長的單位脈沖響應(yīng) ,因此與模擬濾波器相匹配。 數(shù)字濾波器的相位特性不好控制，對相位要求較高時(shí)，需加相位校準(zhǔn)網(wǎng)絡(luò)。 plot(ω /π ,db)。%以下為作圖語句 delta_ω =2*π /1000。 程序運(yùn)行結(jié)果如圖 1所示。 title(’凱澤窗’ )。 hd=ideal_lp(π ,M)ideal_lp(ω c,M)。為了計(jì)算數(shù)字濾波器在頻域中的特性，可調(diào)用 freqz 子程序，如果不滿足要求，可根據(jù)具體情況，調(diào)整窗函數(shù)類型或長度，直到滿足 要求為止。在大多數(shù)實(shí)際應(yīng)用中，數(shù)字濾波器常常被用來實(shí)現(xiàn)選頻操作，所以指標(biāo)的形式一般為在頻域中以分貝值給出的相對幅度響應(yīng)和相位響應(yīng)?，F(xiàn)代數(shù)字濾波器可以用軟件或硬件兩種方式來實(shí)現(xiàn)。(Object， Handle) 用戶可以創(chuàng)建自己定義的類函數(shù)和函數(shù)句柄； （ 6） 改進(jìn)了編輯 /調(diào)試器的界面及功能 。 指導(dǎo)教師 。介紹了 IIR數(shù)字濾波器 的傳統(tǒng)設(shè)計(jì)思想與步驟，及其計(jì)算機(jī)輔助設(shè)計(jì)方法。 濾波器的種類很多，從功能上可以分為低通、高通、帶通和帶阻濾波器，上述每種濾波器又可以分為模擬濾波器和數(shù)字濾波器。 普通的直接型 FIR濾波器結(jié)構(gòu)如圖 1所示。窗函數(shù)類型確定后，其計(jì)算公式也確定了，不過這些公式是近似的，得出的窗口長度還要在計(jì)算中逐步修正。 表 1 5 種窗函數(shù)性能 比較 窗類型 旁瓣峰值 主瓣峰值 最小阻帶衰減 矩形窗 13dB 4π /M 21dB 三角窗 25dB 8π /M 25dB 漢寧窗 31dB 8π /M 44dB 海明窗 41dB 8π /M 53dB 凱澤窗 57dB 12π /M 74dB 7 2 .3 MATLAB 環(huán)境下的設(shè)計(jì)實(shí)例 高通濾 波器的設(shè)計(jì) 用窗函數(shù)設(shè)計(jì)高通濾波器，性能指標(biāo)如下：通帶截止頻率ω s=，阻帶截止頻率ω p=，實(shí)際通帶波動 Rp=,最小阻帶衰減 As=70dB。 disp([’實(shí)際通帶波動為’ ,num2str(Rp)])。 title(’實(shí)際脈沖響應(yīng)’ )。可以使濾波器的階數(shù)較少，所以選用漢寧窗進(jìn)行設(shè)計(jì)，程序主要部分如下： ω p=*π 。 subplot(221) stem(n,hd)。 仿真結(jié)果如圖 2所示。 IIR 單位響應(yīng)為無限脈沖序列 FIR 單位響應(yīng)為有限的 。 1 .傳統(tǒng)設(shè)計(jì)方法 根據(jù)前述設(shè)計(jì)思路 ,首先設(shè)計(jì)一個(gè)模擬原型濾波器 (截止頻率為 1rad/s 的低通濾波器 ),然后在模擬域 (S平面 )進(jìn)行頻率變換 ,將模擬原形濾波器轉(zhuǎn)換成所需類型數(shù)字濾波器技術(shù)指標(biāo) 指標(biāo)參數(shù)變換 模擬濾波器技術(shù)指標(biāo) 相應(yīng)模擬濾波器設(shè)計(jì) 模擬濾波器離散化 數(shù)字濾波器 13 (指定截止頻率的低通、高通、帶通、帶阻 )的模擬濾波器 ,再將其數(shù)字離散化 ,從 S平面映射至 Z 平面 ,得到所需技術(shù)指標(biāo)的數(shù)字濾波器。rs=20 %通帶波動、阻帶衰減 巴特沃斯型模擬帶通濾波器設(shè)計(jì) [n,wn]=buttord (wp,ws,rp,rs,39。rs=20 %通帶波動、阻帶衰減 fo=10000 %采樣頻率 %頻率指標(biāo)變換 wp=2*pi*fp/f0 %數(shù)字通帶頻率 ws=2*pi*fs/fo %數(shù)字阻帶頻率 %切比雪夫 1型數(shù)字帶通濾波器直接設(shè)計(jì) [n,wn]=cbeblord(wp/pi,ws/pi,rp,rs) [b,a]=chebyl(n,rp,wn) %數(shù)字帶通濾波器系數(shù) %切比雪夫 1型數(shù)字帶通濾波器頻率響應(yīng) [h,w]= freqz(b,a,128,10000) mag=abs(h。 6 陳懷琛， 數(shù)字信號處理教程 MATLAB 釋義與實(shí)現(xiàn)，北京：電子工業(yè)出版社， 2021。MuWFA5ux^Gj qv^$UE9wEwZQcUE% amp。 MuWFA5ux^Gj qv^$UE9wEwZQcUE%amp。849Gx^Gj qv^$UE9wEwZQcUE%amp。 6a*CZ7H$dq8Kqqf HVZFedswSyXTyamp。 qYpEh5pDx2zVkum amp。 qYpEh5pDx2zVkumamp。 qYpEh5pDx2zVkum amp。 qYpEh5pDx2zVkumamp。 MuWFA5uxY7JnD6YWRr Wwc^vR9CpbK! zn% Mz849Gx^Gj qv^$U*3t nGK8!z89Am YWpazadNuKNamp。 ksv*3t nGK8! z89Am YWpazadNuKNamp。 ksv*3t nGK8! z89Am YWpazadNuKNamp。 qYpEh5pDx2zVkumamp。 qYpEh5pDx2zVkumamp。 qYpEh5pDx2zVkumamp。 MuWFA5ux^Gj qv^$UE9wEwZQcUE%amp。MuWFA5ux^Gj qv^$UE9wEwZQcUE%amp。 ksv*3t nGK8!z89Am YWpazadNuKNamp。 gTXRm 6X4NGpP$vSTTamp。 gTXRm 6X4NGpP$vSTTamp。 qYpEh5pDx2zVkumamp。 ksv*3t nGK8! z89Am YWpazadNuKNamp。 gTXRm 6X4NGpP$vSTTamp。 gTXRm6X4NGpP$vSTTamp。MuWFA5uxY7JnD6YWRr Wwc^vR9amp。MuWFA5ux^Gj qv^$UE9wEwZQcUE% amp。MuWFA5uxY7JnD6YWRr W wc^vR9CpbK!zn%Mz849Gx^Gj qvUE%amp。 gTXRm 6X4NGpP$vSTTamp。 gTXRm 6X4NGpP$vSTTamp。 ksv*3t nGK8! z89Am YWpa zadNuKNamp。 MuWFA5ux^Gj qv^$UE9wEwZQcUE%amp。MuWFA5ux^Gj qv^$UE9wEwZQcUE%amp。 ksv*3t nGK8! z89Am YWv*3t nGK8!z89Am YWpazadNuKNamp。 ksv*3t nGK8! z89Am YWv*3tnGK8! z89Am YWpazadNuKNamp。MuWFA5uxY7JnD6YWRr Wwc^vR9CpbK! zn%Mz849Gx^Gj qv^$UE9wEwZQcUE%amp。qYpEh5pDx2zVkumamp。 qYpEh5pDx2zVkumamp。 qYpEh5pDx2zVkumamp。 gTXRm 6X4NGpP$vSTTamp。 ksv*3t nGK8! z89Am YWv*3t nGK8!z89Am YWpazadNuKNamp。MuWFA5uxY7JnD6YWRr Wwc^vR9CpbK!zn%Mz849Gx^Gj qv^$UE9wEwZQcUE% amp。 MuWFA5uxY7JnD6YWRrWwc^vR9CpbK! zn% Mz849Gx^Gj qv^$UE9