无限冲激响应
無限脈衝響應濾波器，簡稱IIR數位濾波器（），是數位濾波器的一種。由於無限脈衝響應濾波器中存在反饋迴路，因此對於脈衝輸入信號的響應是無限延續的。濾波器的任務是通過一定運算關係，改變輸入信號的頻譜。數字濾波器是利用計算機程序、專用芯片等軟、硬體改變數字信號頻譜。如果要處理的是模擬信號可以通過A/D在信號形式上進行轉換，在利用數字濾波器處理後經過D/A恢復為模擬信號。
利用選頻濾波器的特性可以提取有用的信號頻率分量，這類濾波器適用輸入信號中有用信號頻帶與干擾信號頻帶不同的情況。 選頻數字濾波器設計過程一般歸納為以下三個步驟:
無限脈衝響應濾波器的脈衝響應無限長，不像有限脈衝響應濾波器，所以有以下問題： 能使用最小相位濾波器（minimum phase filter）來解決以上問題，使一無限脈衝響應濾波器為穩定且能量集中在formula_1附近。無限脈衝響應濾波器（IIR filter）的優點： 無限脈衝響應濾波器（IIR filter）的缺點：
IIR數字濾波器設計有以下幾種基本方法:
脈冲響應不變法也稱冲激響應不變法，是滿足s平面到z平面的映射兩個基本要求的常用方法之一。實現方法是已知模擬濾波器的衝冲響應h(t)，讓數字濾波器h(n)的包絡為h(t)，即 故此得名脈衝不變法，這種方法的實質上是時域採樣法，基本設計步驟是先得到模擬濾波器的系統函數H(s)；然後對沖激響應h(t)採樣，得到數字濾波器的單位脈冲響應h(n)；最後，h(n)所對應的z變換，正是所要求的數字系統的系統函數 H(z)，即 脈衝不變法又稱標準z變換法。根據理想採樣序列z變換與拉普拉斯變換的關係式，我們可以得到 其映射關係為formula_27 或 formula_28 由(1)式當s = jΩ 時，得到的數字濾波器的響應為 由(2)式可見，數字濾波器的頻率響應不是簡單的重現模擬濾波器的頻率響應，而是模擬濾波器頻率響應的週期展開，其中頻率映射關係為formula_30 若模擬濾波器是帶限的，即formula_31 ,則 這時，數字濾波器可在摺疊頻率內不失真的重現模擬濾波器的頻率響應。但是實際模擬濾波器頻率響應不可能是真正帶限的，所以脈冲響應不辨法總會有混疊失真。 下面以H均為單極點為例，討論脈衝不變法設計IIR濾波器的具體方式。H可展開為部分分式 上面式子中formula_34是H_a(s)的極點，對應的模擬系統單位沖激響應為 對h(t)取樣，所有連續時間便量t由離散時間便量nT代替，則 對H(n)取z變換 上式中formula_38是z平面的極點，當極點s在s平面的左半平面時，formula_39，則 formula_40 ，formula_41 ，H(z)的極點在單位圓內。由此可知，s平面與z平面的極點一一對應，並且穩定的模擬系統變換後仍為穩定的數字系統。 在(3)式formula_42，考慮到若取樣頻率很高時，會使得數字系統的福頻特性的增益很高。為防止溢出，通常在取樣時做修正，令h(n)=Th(n)，這樣實際數字濾波器的系統函數及頻率響應特性分別為 由以上分析可見，按照formula_45變換過程，從模擬濾波器H(s)得到數字濾波器H(z)，設計還是相當冗繁，不過由於s平面與z平面的極點一一對應，我們可以簡化變換步驟。 設H(s)均為一階極點，省略中間步驟，從formula_46設計過程分為兩個步驟: 如果H(s)還有m階的重極點，我們給出對應的z變換的有理分式為 除了單極點與重極點外，對H(s)的共軛極點，我們也給出對應的z變換的一般有理分式 這樣，當H(s)部分分式展開式中如果有重極點或共軛極點，在上述設計過程第二步做相應修改即可。 脈衝不變法的頻率變換formula_52是線性關係，所以在幅度及相位上數字濾波器都保留著模擬濾波器的頻率響應特點。但是由於脈沖不變法本身具有頻譜周期延展拓寬的特性，所以失真較大。尤其是對衰減特性差的濾波器頻率響應會有嚴重的混疊現象，適用範圍受到限制。 下面給出用脈沖不變法設計數字濾波器的一般步驟 在設計過程中，除了第一步求數字臨界頻率formula_53時，要用到取樣間隔T或取樣頻率f以外，最後的結果與其他各步驟中T或f的取值無關。所以為了簡化運算，在實際運算時，除了第一步外，通常取T=1。 這是由解出和類比濾波器在同樣的取樣時間下有一樣輸出值的T(z)所獲得,且當這兩者的輸入接為一個脈衝時才適用。 有一點需要注意的是,藉由此種方法產生的數位濾波器的所有輸入皆是一個近似值，除了只對於脈衝輸入非常精確。這是一種最簡單的IIR濾波器設計方法，它在低頻是最精確的，所以通常被用在低通濾波器上。
對於拉普拉斯轉換(Laplace transform)或是z-transform，轉換過後的輸出只是輸入乘上相對應的轉換函數，T(s)或T(z)。Y(s)和Y(z)分別是輸入X(s)和輸入X(z)轉換過後的輸出。 然而，當拉普拉斯轉換和z-transform對單位脈衝作轉換時為1，轉換後的輸出結果分別為 現在類比濾波器的輸出在時域上就只是拉普拉斯逆轉換(inverse Laplace transform) 如果我們用nT代替t，可以得到在取樣時間下源於脈衝的輸出，y(nT)，又可以為了便利性表為y(n) 這個離散時間信號可以被z-transform得到T(z) 最後這條等式指出在數學領域描述數位IIR濾波器就是對取樣完經拉普拉斯轉換為T(s)的類比信號作z-transform，通常會簡化為 雖然在數學領域上不是完全正確的。 需特別注意有一個乘法數T出現在式子中，這是因為就算拉普拉斯轉換和z-transform對單位脈衝為1，脈衝本身不一定是一樣的。對類比信號而言，脈衝在t=0時有無限大的值但是面積為1，但在離散時間的脈衝t=0時為1，所以為了這個情況，需要乘法數T的存在。
這通常相較於impulse invariant昰一個更好的設計方法，數位濾波器正是在取樣時有幾段不同常數的輸入，也就是由離散步階(step)組成。step invariant IIR 濾波器和同樣輸入步階信號至ADC相比是較低精確度的，然而，和impulse invariant相比這對任何輸入而言是一個較好的近似法。 當T(z)和T(s)兩者都是步階(step)輸入時，step invariant解決一樣取樣值的困擾，對數位濾波器的輸入為u(s)，對類比濾波器的輸入是u(t)，z-transform和拉普拉斯轉對這兩個輸入作轉換以得到轉換後的輸出信號。 z-transform的轉換函數為formula_68 經z-transform的輸出formula_69 Laplace transform的轉換函數為formula_70 經Laplace transform的輸出formula_71 類比濾波器的輸出是y(t)，就是Y(s)的拉普拉斯逆轉換。如果每經過T秒就被取樣一次，它就是y(n)，也就是Y(z)的逆轉換。這些信號被用來解出數位濾波器的和類比濾波器在取樣時間下有著相同輸出 以下的等式指出T(z)的解法，也就是對類比濾波器取近似的公式
BLT公式利用拉普拉斯轉換(Laplace transform)的特性：對於一個x(t)經Laplace轉換成的X(s)而言，它的積分在沒有初始條件下相當於X(s)/s。所以得出將拉普拉斯轉換的信號乘上1/s相當於對此信號在時域下做積分。 然而，接下來這個不同的等式一樣也是用近似積分x(t)信號的特性達成：利用x(nT)=x(n)。 在沒有初始條件下積分信號相當於計算信號下的面積 此等式告訴我們再經過一個新的數入信號x(n)後計算新的面積y(n)，就只是將先前運算的面積y(n-1) 加上新輸入和先前輸入的平均乘上取樣時間。藉由z轉換先前的等式變為以下的式子 當我們得出X(z)和Y(z)後，我們將解出T(z) 先前的等式是一個數位訊號處理系統中的轉換函數，此函數將一個系統中的類比轉換函數趨近成作積分，並在分子分母同乘上z後開始解s 這個關係被用在任何類比濾波器的拉普拉斯轉換函數，以或取近似於類比濾波器的數位無限脈衝響應(IIR)濾波器T(z)。這個結果被標示在下列公式中，此公式用來運算BLT IIR數位濾波器，並從類比濾波器的拉普拉斯轉換函數開始運算起 這條公式指出近似類比濾波器的bilinear transform數位IIR濾波器可以從把類比轉換函數T(s)中的s替換成和z的關係得出