基于雙線性變換法的IIR數(shù)字高通濾波器設(shè)計(jì)教材

1、基于雙線性變換法的 IIR 數(shù)字高通濾波器設(shè)計(jì)基于雙線性變換法的 IIR 數(shù)字高通濾波器設(shè)計(jì)摘要隨著信息時(shí)代和數(shù)字世界的到來(lái)， 數(shù)字信號(hào)處理已成為當(dāng)今一門極其重要的 學(xué)科和技術(shù)領(lǐng)域。在數(shù)字信號(hào)處理中起著重要的作用并已獲得廣泛應(yīng)用的是數(shù)字 濾波器( DF，Digital Filter)。數(shù)字濾波器是一種用來(lái)過(guò)濾時(shí)間離散信號(hào)的數(shù)字系統(tǒng)，通過(guò)對(duì)抽樣數(shù)據(jù)進(jìn)行數(shù)學(xué)處理來(lái)達(dá)到頻域?yàn)V波的目的。 實(shí)現(xiàn) IIR 濾波器 的階次較低，所用的存儲(chǔ)單元較少，效率高，精度高，而且能夠保留一些模擬濾 波器的優(yōu)良特性，因此應(yīng)用很廣。 Matlab 軟件以矩陣運(yùn)算為基礎(chǔ)，把計(jì)算、可 視化及程序設(shè)計(jì)有機(jī)融合到交互式工作環(huán)境中，

2、 并且為數(shù)字濾波的研究和應(yīng)用提 供了一個(gè)直觀、高效、便捷的利器。尤其是 Matlab 中的信號(hào)處理工具箱使各個(gè) 領(lǐng)域的研究人員可以直觀方便地進(jìn)行科學(xué)研究與工程應(yīng)用。本次主要利用 matlab 的信號(hào)處理工具箱設(shè)計(jì)一個(gè)數(shù)字高通 IIR 濾波器，并 用此濾波器處理一段音頻信號(hào)。 此報(bào)告重點(diǎn)介紹了用雙線性不變法設(shè)計(jì) IIR 數(shù)字 濾波器的基本流程， 比較了各種設(shè)計(jì)方法的優(yōu)缺點(diǎn)， 總結(jié)了模擬濾波器的性能特 征。最后以雙線性不變法設(shè)計(jì)了一個(gè)高通巴特沃斯 IIR 數(shù)字濾波器， 介紹了設(shè)計(jì) 步驟，然后在 Matlab 環(huán)境下進(jìn)行了仿真與調(diào)試，實(shí)現(xiàn)了設(shè)計(jì)目標(biāo)。關(guān)鍵字 IIR 數(shù)字濾波器 巴特沃斯 高通 雙線性

3、變換 MATLAB基于雙線性變換法的 IIR 數(shù)字高通濾波器設(shè)計(jì)AbstractWith the information era and the advent of the digital world, digitalsignal processing has become one of the important fields of science and technology. In the digital signal processing plays an important role and has been widely applied in the digital filter

4、. Digital filter is a discrete time signal is used to filter the digital system,through the mathematical processing of sampling data to reach the goal of frequency domain filterin. Implementation of IIR filter order is low, less storage unit, with high efficiency, high precision, and can keep some s

5、imulation characteristics of filter , so a very wide application. Matlab software based on matrix computation, the calculation,visualization and program design of organic integration to interactiveenvironment, and providesan intuitive, efficient, convenienttool for the research andapplication of dig

6、ital filter .The main use of signal processing toolbox of MATLAB to design a digital high-passIIR filter , and an audio signal with the filter . This report focuses on the basic process of using the bilinearinvariant IIR digital filter design, compares the advantages and disadvantages of various des

7、ign methods, summarizes thecharacteristics of the analogfilter .Finally, by bilinear invariant to design a high Butterworth IIR digitalfilter, introducesthe design steps, andthen in the Matlab environment for thesimulation anddebugging, achieves the design goal.Keywords IIR Butterworth high-pass bil

8、inear transform MATLABdigital filterII基于雙線性變換法的 IIR 數(shù)字高通濾波器設(shè)計(jì)目錄摘 要 . IAbstract . II1 數(shù)字濾波器 . 11.1數(shù)字濾波器介紹 . 11.2IIR 數(shù)字濾波器設(shè)計(jì)原理 . . 21.2.2 雙線性變換法 . 32數(shù)字濾波器設(shè)計(jì)實(shí)現(xiàn) . 62.1典型模擬濾波器比較 . 62.2設(shè)計(jì)步驟 . 62.3程序流程圖 . 72.4音頻信號(hào)部分程序 . . 82.5仿真結(jié)果 . 82.5.1 濾波器性能仿真 . 92.5.2 濾波器除噪性能 . 103總結(jié) . 10致謝 . 11參考文獻(xiàn) . 12附 錄 . 13基于雙線性變

9、換法的 IIR 數(shù)字高通濾波器設(shè)計(jì)1 數(shù)字濾波器1.1數(shù)字濾波器介紹數(shù)字濾波器是具有一定傳輸選擇特性的數(shù)字信號(hào)處理裝置 , 其輸入、輸出均 為數(shù)字信號(hào) , 實(shí)質(zhì)上是一個(gè)由有限精度算法實(shí)現(xiàn)的線性時(shí)不變離散系統(tǒng)。它的基 本工作原理是利用離散系統(tǒng)特性對(duì)系統(tǒng)輸入信號(hào)進(jìn)行加工和變換 , 改變輸入序列 的頻譜或信號(hào)波形 ,讓有用頻率的信號(hào)分量通過(guò) ,抑制無(wú)用的信號(hào)分量輸出。 數(shù)字 濾波器和模擬濾波器有著相同的濾波概念 , 根據(jù)其頻率響應(yīng)特性可分為低通、高 通、帶通、帶阻等類型 , 與模擬濾波器相比 , 數(shù)字濾波器除了具有數(shù)字信號(hào)處理的 固有優(yōu)點(diǎn)外 ,還有濾波精度高 (與系統(tǒng)字長(zhǎng)有關(guān) ) 、穩(wěn)定性好(僅運(yùn)行

10、在 0與 l 兩個(gè) 電平狀態(tài) ) 、靈活性強(qiáng)等優(yōu)點(diǎn)。時(shí)域離散系統(tǒng)的頻域特性 :Yejw X ejw H ejw ，其中Yejw ，X ejw 分別 是數(shù)字濾波器的輸出序列和輸入序列的頻域特性 (或稱為頻譜特性) ， H ejw 是 數(shù)字濾波器的單位取樣響應(yīng)的頻譜， 又稱為數(shù)字濾波器的頻域響應(yīng)。 輸入序列的 頻譜 X ejw 經(jīng)過(guò)濾波后 X e jw H ejw 。因此，只要按照輸入信號(hào)頻譜的特點(diǎn)和處 理信號(hào)的目的， 適當(dāng)選擇 H ejw ，使得濾波后的 X ejw H ejw 滿足設(shè)計(jì)的要求， 這就是數(shù)字濾波器的濾波原理。數(shù)字濾波器根據(jù)其沖激響應(yīng)函數(shù)的時(shí)域特性， 可分為兩種， 即無(wú)限長(zhǎng)沖激響

11、應(yīng) (IIR) 數(shù)字濾波器和有限長(zhǎng)沖激響應(yīng) (FIR) 數(shù)字濾波器。IIR 數(shù)字濾波器的特征 是，具有無(wú)限持續(xù)時(shí)間沖激響應(yīng)，需要用遞歸模型NN 來(lái)實(shí)現(xiàn)，其差分方程為： y n aix n ibi y n ii 0 i 1Mbrzr系統(tǒng)函數(shù)為： H zr 0N1akz kk1設(shè)計(jì)IIR 濾波器的任務(wù)就是尋求一個(gè)物理上可實(shí)現(xiàn)的系統(tǒng)函數(shù) H(z) ，使其頻 率響應(yīng)H(z) 滿足所希望得到的頻域指標(biāo)， 即符合給定的通帶截止頻率、 阻帶截止 頻率、通帶衰減系數(shù)和阻帶衰減系數(shù)?；陔p線性變換法的 IIR 數(shù)字高通濾波器設(shè)計(jì)1.2 IIR 數(shù)字濾波器設(shè)計(jì)原理IIR 數(shù)字濾波器是一種離散時(shí)間系統(tǒng)，其系統(tǒng)函數(shù)為

12、Mbk z kHzk0Nk1ak zk1YzXz假設(shè)MN，當(dāng)MN時(shí),系統(tǒng)函數(shù)可以看作一個(gè) IIR 的子系統(tǒng)和一個(gè) (M-N)的FIR 子系統(tǒng)的級(jí)聯(lián)。 IIR 數(shù)字濾波器的設(shè)計(jì)實(shí)際上是求解濾波器的系數(shù) ak和bk ，它是 數(shù)學(xué)上的一種逼近問(wèn)題， 即在規(guī)定意義上 (通常采用最小均方誤差準(zhǔn)則) 去逼近 系統(tǒng)的特性。如果在S平面上去逼近，就得到模擬濾波器； 如果在z平面上去逼近， 就得到數(shù)字濾波器。設(shè)計(jì)通數(shù)字濾波器通?？梢詺w納為如圖所示的兩種常用方法。圖 1-1 數(shù)字帶通濾波器設(shè)計(jì)的兩種方法方法 1: 首先設(shè)計(jì)一個(gè)模擬原型低通濾波器，然后通過(guò)頻率變換成所需要的 模擬高通濾波器，最后再使用沖激不變法或雙

13、線性變換成相應(yīng)的數(shù)字高通濾波 器。方法 2: 先設(shè)計(jì)一個(gè)模擬原型低通濾波器，然后采用沖激響應(yīng)不變法或雙線 性變換法將它轉(zhuǎn)換成數(shù)字原型低通濾波器， 最后通過(guò)頻率變換把數(shù)字原型低通濾 波器變換成所需要的數(shù)字高通濾波器。本課程設(shè)計(jì)采用第一種設(shè)計(jì)方法， 先構(gòu)造一個(gè)巴特沃斯模擬低通濾波器， 然 方法 1方法 2后將模擬低通濾波器轉(zhuǎn)換成模擬高通濾波器， 最后利用雙線性變換將模擬高通濾 波器轉(zhuǎn)換成數(shù)字高通濾波器。S平面2tanT21T基于雙線性變換法的 IIR 數(shù)字高通濾波器設(shè)計(jì)1.2.2 雙線性變換法為了克服沖激響應(yīng)法可能產(chǎn)生的頻率響應(yīng)的混疊失真， 這是因?yàn)閺?S 平面到 平面是多值的映射關(guān)系所造成的。

14、為了克服這一缺點(diǎn)， 可以采用非線性頻率壓 縮方法，將整個(gè)頻率軸上的頻率范圍壓縮到 -/T/T 之間，再用 z=esT轉(zhuǎn)換 到 Z 平面上。也就是說(shuō)，第一步先將整個(gè) S平面壓縮映射到 S1平面的 -/T / T 一條橫帶里；第二步再通過(guò)標(biāo)準(zhǔn)變換關(guān)系 z=es1T 將此橫帶變換到整個(gè) Z 平面 上去。這樣就使 S平面與 Z 平面建立了一一對(duì)應(yīng)的單值關(guān)系， 消除了多值變換性， 也就消除了頻譜混疊現(xiàn)象，映射關(guān)系如圖 1-3為了將 S平面的整個(gè)虛軸 j 壓縮到 S1平面 j 1軸上的-/T到/T段上，可以通過(guò)以下的正切變換實(shí)現(xiàn)式中, T仍是采樣間隔。當(dāng)1由-/T經(jīng)過(guò) 0變化到/T時(shí)，由-經(jīng)過(guò) 0 變化到

15、+，也即映射 了整個(gè) j 軸。將式（ 1-9 ）寫成j2 ej 1T/2 ej 1T /2jT ej 1T/2 e j 1T /2將此關(guān)系解析延拓到整個(gè) S平面和 S1平面，令 j =s，j 1=s1，則得2 es1T /2e s1T /22s1T2 1e s1Ts sT /2 sT/2 tanh sTT es1T /2e s1T /2T2T 1e s1T再將 S1平面通過(guò)以下標(biāo)準(zhǔn)變換關(guān)系映射到 Z 平面 z=es1T從而得到 S平面和 Z 平面的單值映射關(guān)系為：j圖 1-3 雙線性變換的映射關(guān)系21 z 1T1 z 11 T s2T1s22sT2sT2j2其次，將 s jw 代入，得2 ta

16、nT基于雙線性變換法的 IIR 數(shù)字高通濾波器設(shè)計(jì)這兩個(gè)關(guān)系式是 S平面與 Z 平面之間的單值映射關(guān)系， 這種變換都是兩個(gè)線 性函數(shù)之比，因此稱為雙線性變換 首先, 把z ejw可得21 e j2s j j tanT 1 e jT即 S 平面的虛軸映射到 Z 平面的單位圓。2jz T2jjTj因此由此看出，當(dāng) 0 時(shí)，| z|0 時(shí)，| z|1 。也就是說(shuō)， S平面的左 半平面映射到 Z 平面的單位圓內(nèi)， S平面的右半平面映射到 Z 平面的單位圓外， S平面的虛軸映射到 Z 平面的單位圓上。因此，穩(wěn)定的模擬濾波器經(jīng)雙線性變換 后所得的數(shù)字濾波器也一定是穩(wěn)定的。雙線性變換法優(yōu)缺點(diǎn)： 雙線性變換法

17、與脈沖響應(yīng)不變法相比， 其主要的優(yōu)點(diǎn) 是避免了頻率響應(yīng)的混疊現(xiàn)象。 這是因?yàn)?S平面與 Z平面是單值的一一對(duì)應(yīng)關(guān)系。 S平面整個(gè) j 軸單值地對(duì)應(yīng)于 Z平面單位圓一周，即頻率軸是單值變換關(guān)系。 這個(gè)關(guān)系重寫如下：上式表明， S平面上與 Z平面的成非線性的正切關(guān)系，如圖 1-4 所示基于雙線性變換法的 IIR 數(shù)字高通濾波器設(shè)計(jì)由圖 1-4 看出，在零頻率附近， 模擬角頻率 與數(shù)字頻率 之間的變換關(guān)系接近于線性關(guān)系；但當(dāng) 進(jìn)一步增加時(shí)， 增長(zhǎng)得越來(lái)越慢，最后當(dāng) 時(shí)， 終止在折疊頻率 =處，因而雙線性變換就不會(huì)出現(xiàn)由于高頻部分超過(guò)折疊但是雙線性變換的這個(gè)特點(diǎn)是靠頻率的嚴(yán)重非線性關(guān)系而得到的，圖 1

18、-4 所示。由于這種頻率之間的非線性變換關(guān)系，就產(chǎn)生了新的問(wèn)題。首先，一個(gè)線 性相位的模擬濾波器經(jīng)雙線性變換后得到非線性相位的數(shù)字濾波器， 不再保持原 有的線性相位了； 其次，這種非線性關(guān)系要求模擬濾波器的幅頻響應(yīng)必須是分段 常數(shù)型的，即某一頻率段的幅頻響應(yīng)近似等于某一常數(shù) （這正是一般典型的低通、 高通、帶通、帶阻型濾波器的響應(yīng)特性） ，不然變換所產(chǎn)生的數(shù)字濾波器幅頻響 應(yīng)相對(duì)于原模擬濾波器的幅頻響應(yīng)會(huì)有畸變，如圖 1-5 所示。對(duì)于分段常數(shù)的濾波器， 雙線性變換后， 仍得到幅頻特性為分段常數(shù)的濾波 器，但是各個(gè)分段邊緣的臨界頻率點(diǎn)產(chǎn)生了畸變， 這種頻率的畸變， 可以通過(guò)頻 率的預(yù)畸來(lái)加以校

19、正。 也就是將臨界模擬頻率事先加以畸變， 然后經(jīng)變換后正好 映射到所需要的數(shù)字頻率上。基于雙線性變換法的 IIR 數(shù)字高通濾波器設(shè)計(jì)2數(shù)字濾波器設(shè)計(jì)實(shí)現(xiàn)2.1典型模擬濾波器比較1、Butterworth 巴特沃斯濾波器：它具有單調(diào)下降的幅頻特性；即最平幅度。2、Chebyshev切比雪夫?yàn)V波器：在通帶或阻帶等波紋，可提高選擇性。3、Bessel 貝塞爾濾波器：在通帶內(nèi)有較好的線性相位特性。4、Ellipse 橢圓濾波器：其選擇性相對(duì)前三種是最好的。此處選擇巴特沃斯主要是想獲得最平穩(wěn)的幅頻響應(yīng)。而不計(jì)較相位特性，而 用雙線性變換法也會(huì)將線性相位變?yōu)榉蔷€性相位。2.2設(shè)計(jì)步驟根據(jù)以上 IIR 數(shù)字

20、濾波器設(shè)計(jì)方法，下面運(yùn)用雙線性變換法基于 MATLAB設(shè) 計(jì)一個(gè) IIR 高通濾波器。(1)確定性能指標(biāo)在設(shè)計(jì)高通濾波器之前 , 首先根據(jù)工程實(shí)際的需要確定濾波器的技術(shù)指標(biāo) : 通帶 截止頻率 fp=3200Hz；阻帶截止頻率 fs=3000Hz；阻帶最小衰減 As=100dB和通帶 最大衰減 Ap=1dB(2)把頻率轉(zhuǎn)化為數(shù)字角頻率wp=2*pi*fp*T;ws=2*pi*fs*T;(3)頻率預(yù)畸變以及模擬高通性能指標(biāo)轉(zhuǎn)換成模擬低通性能指標(biāo)用=2/T*tan(w/2) 對(duì)高通數(shù)字濾波器 H(z) 的數(shù)字邊界頻率預(yù)畸變 , 得到高 通模擬濾波器 H(s) 的邊界頻率主要是通帶截止頻率 Wp;阻

21、帶截止頻率 Ws的轉(zhuǎn)換。 抽樣頻率 Fs=44100Hz 。通帶截止頻率 Wp2=(2/T)*tan(wp/2)阻帶截止頻率 Ws2=(2/T)*tan(ws/2)(4)歸一化處理Wp2 1 1 ；p 1 ； p1；p Wp2 p p61N 12lgH(z)=H(s) |1zH aL ( p) |1zpc1z1z基于雙線性變換法的 IIR 數(shù)字高通濾波器設(shè)計(jì)Ws2Wp2由式2-3 ，2-4 ，2-5計(jì)算出 N，查表可得模擬低通濾波器的階數(shù)，從而由下式 確定模擬高通濾波器的參數(shù)H z H s 1 z 1 s1 z 1(5)模擬低通濾波器的構(gòu)造 借助巴特沃斯濾波器模型得到模擬低通濾波器的傳輸函數(shù)

22、Ha(s) N,Wn=buttord(wp2,ws2,Ap,As,s);z,p,k=buttap(N);Bap,Aap=zp2tf(z,p,k);G1=tf(Bap,Aap);(6)模擬低通濾波器轉(zhuǎn)換成模擬高通濾波器調(diào)用 lp2bp 函數(shù)將模擬低通濾波器轉(zhuǎn)化為模擬高通濾波器。 Bbs,Abs=lp2hp(Bap,Aap,Wn);(7)模擬高通濾波器轉(zhuǎn)換成數(shù)字高通濾波器 利用雙線性變換法將模擬高通濾波器 Ha(s) 轉(zhuǎn)換成數(shù)字高通濾波器 H(z) Bbz,Abz=bilinear(Bbs,Abs,FS);2.3程序流程圖首先確定性能指標(biāo)， 把頻率轉(zhuǎn)化為數(shù)字角頻率， 進(jìn)而在進(jìn)行頻率預(yù)畸變， 用 =

23、2/T*tan(w/2) 對(duì)高通數(shù)字濾波器 H(z) 的數(shù)字邊界頻率預(yù)畸變 , 得到高通模擬 濾波器 H(s) 的邊界頻率主要是通帶截止頻率 Wp;阻帶截止頻率 Ws的轉(zhuǎn)換。抽樣 頻率 Fs=44100Hz。上述準(zhǔn)備工作做好之后，就先把模擬高通性能指標(biāo)轉(zhuǎn)換成模擬低通性能1sAs10 10 1指 標(biāo)，然后設(shè)計(jì)模擬低通濾波器， 借助巴特沃斯濾波器模型得到模擬低通濾波器的 傳輸函數(shù) Ha(s) 。然后將模擬低通濾波器轉(zhuǎn)化為模擬高通濾波器。最后利用雙線 性變換法將模擬高通濾波器 Ha(s) 轉(zhuǎn)換成數(shù)字高通濾波器 H(z) 。開始基于雙線性變換法的 IIR 數(shù)字高通濾波器設(shè)計(jì)讀入數(shù)字濾波器技術(shù)指標(biāo)將指標(biāo)

24、轉(zhuǎn)換成歸一化模擬低通濾波器的指標(biāo)設(shè)計(jì)歸一化的模擬低通濾波器階數(shù) N和 3db 截止頻率模擬低通濾波器變換成模擬高通濾波器模擬高通濾波器轉(zhuǎn)換成數(shù)字高通濾波器結(jié)束圖 2-1 程序流程圖2.4音頻信號(hào)部分程序1)寫入聲音信號(hào)fnose =10;T=length(x)/fs; t=T/length(x):T/length(x):T;j=1:1:length(x);2)產(chǎn)生噪聲并加到聲音中nose =sin(2*pi*fnose*T*j); y=x+5*nose; b,a=butter(N,Wn,high);Y=fft(y) ;2.5仿真結(jié)果基于雙線性變換法的 IIR 數(shù)字高通濾波器設(shè)計(jì)2.5.1 濾波

25、器性能仿真源程序設(shè)計(jì)了模擬低通濾波器、 模擬高通濾波器與數(shù)字高通濾波器， 對(duì)數(shù)字 高通濾波器的性能仿真如下：92.5.2 濾波器除噪性能為了實(shí)現(xiàn)濾波器的應(yīng)用， 程序中加入了有噪聲的音頻信號(hào)， 通過(guò)對(duì)其濾波處理，來(lái)顯示數(shù)字高通濾波器的功能， 下面顯示未加入噪聲， 加入噪聲和濾波后的 頻譜波形圖。基于雙線性變換法的 IIR 數(shù)字高通濾波器設(shè)計(jì)10基于雙線性變換法的 IIR 數(shù)字高通濾波器設(shè)計(jì)3總結(jié)這次課設(shè)下來(lái)， 對(duì)設(shè)計(jì)高通數(shù)字濾波器的整個(gè)過(guò)程有了很好的掌握， 懂得了 設(shè)計(jì)濾波器的基本方法， 對(duì)雙線性變換法， 巴特沃斯濾波器有了一定了解， 同時(shí) 呢也熟悉了 MATLAB的環(huán)境，鞏固了相關(guān)知識(shí)。 初步

26、掌握了 MATLAB語(yǔ)言在數(shù)字信 號(hào)處理中一些基本庫(kù)函數(shù)的調(diào)用和編寫基本程序等應(yīng)用； 熟悉了濾波器設(shè)計(jì)的一 般原理，對(duì)濾波器有了一個(gè)感性的認(rèn)識(shí)； 學(xué)會(huì)了數(shù)字高通濾波器設(shè)計(jì)的一般步驟； 加深了對(duì)濾波器設(shè)計(jì)中產(chǎn)生誤差的原因以及雙線性變換法優(yōu)缺點(diǎn)的理解和認(rèn)識(shí)。 總之，使理論聯(lián)系了實(shí)際， 濾波前加噪信號(hào)的波形和頻譜濾波后信號(hào)的波形和頻譜鞏固并深化了對(duì)課本基本知識(shí)的認(rèn)識(shí)和理解， 使理論 得以升華。在做本次課程設(shè)計(jì)的過(guò)程中， 我深深地感受到了自己所學(xué)到知識(shí)的有限， 明 白了只學(xué)好課本上的知識(shí)是不夠的， 要通過(guò)圖書館和互聯(lián)網(wǎng)等各種渠道來(lái)擴(kuò)充自 己的知識(shí)。 在實(shí)驗(yàn)過(guò)程中我們?cè)?jīng)遇到過(guò)問(wèn)題。 但是從中我們學(xué)習(xí)到

27、了如何對(duì)待 遇到的困難， 進(jìn)一步培養(yǎng)了我們一絲不茍的科學(xué)態(tài)度和不厭其煩的耐心。 所有的 這些心得會(huì)對(duì)我以后的學(xué)習(xí)和工作有幫助作用， 忠心感謝學(xué)校給我們提供這次實(shí) 驗(yàn)機(jī)會(huì)。致謝本設(shè)計(jì)過(guò)程中得到老師的悉心指導(dǎo)，設(shè)計(jì)過(guò)程中老師多次幫我分析思路， 開拓視角， 在我遇到困難想要放棄的時(shí)候給予我最大的支持和鼓勵(lì)。 老師那嚴(yán)謹(jǐn) 求實(shí)的治學(xué)態(tài)度、 踏實(shí)堅(jiān)韌的工作精神， 將使我終身受益； 另外同學(xué)和朋友的熱 心幫助也使我學(xué)到了不少的知識(shí)。 再多華麗的語(yǔ)言也顯得蒼白， 在此，謹(jǐn)向老師 和同學(xué)們致以誠(chéng)摯的謝意和崇高的敬意。 因?yàn)樗麄兾覍?duì)這次課程設(shè)計(jì)有了深刻的 印象，雖然只有短暫的一周，但在這期間，卻讓我受益匪淺，我

28、深刻地認(rèn)識(shí)到了 “理論聯(lián)系實(shí)際”這句話的重要性與真實(shí)性。11基于雙線性變換法的 IIR 數(shù)字高通濾波器設(shè)計(jì)參考文獻(xiàn)1程佩青，數(shù)字信號(hào)處理 . 北京：清華大學(xué)出版社 ,2007.22景振毅，張澤兵，董霖 . 北京：中國(guó)鐵道出版社， 2008.123從玉良 .數(shù)字信號(hào)處理原理及其 MATLAB 實(shí)現(xiàn)M. 北京：電子工業(yè)出版社.2009.74胡廣書.數(shù)字信號(hào)處理理論、 算法與實(shí)現(xiàn) M. 北京：清華大學(xué)出版社 .2003,85劉波 MATLAB信號(hào)處理北京：電子工業(yè)出版社， 20066施陽(yáng)等 MATLAB語(yǔ)言工具箱西安：西北工業(yè)大學(xué)出版社， 19997丁玉美數(shù)字信號(hào)處理西安電子科技大學(xué)出版社 第二版8倪養(yǎng)華數(shù)字信號(hào)處理與實(shí)現(xiàn)上海：上海交通大學(xué)出版社， 19989張葛祥，李娜 MATLAB仿真技術(shù)與應(yīng)用北京：清華大學(xué)出版