Math Problems On IIR Filter Coefficients Calculation
Math Problems On IIR Filter Coefficients Calculation
** Obtain coefficients of lowpass IIR filter using impulse invariant method based on followings.
Transfer function : 𝟏
𝑯(𝒔) =
𝒔𝟐 + √𝟐 𝒔 + 𝟏
Cutoff frequency : 150 Hz
Sampling frequency : 1.28 kHz
[Hints for Calculator Use: Assume pi = 3.1416 when it is not under any sin or cos function and pi =
180 when it is under any sin or cos function]
𝑠
Here, 𝛼 = 2𝜋𝑓𝑐 = 2 × 3.1416 × 150 = 942.778. So putting 𝑠 = 𝛼 in the transfer function we get,
𝛼2 𝐶1 𝐶2
𝐻(𝑠) = = +
𝑠2 + √2 𝛼𝑠 + 𝛼2 𝑠 − 𝑝1 𝑠 − 𝑝2
Where,
𝑝1 = −666.4324(1 − 𝑗), 𝑝2 = 𝑝1∗
𝐶1 = −𝑗666.4324, 𝐶2 = 𝐶1∗
(Solve the second order equation 𝒔𝟐 + √𝟐 𝜶𝒔 + 𝜶𝟐 to find out 𝒑𝟏 and 𝒑𝟐 . Then calculate 𝑪𝟏 and 𝑪𝟐
using partial fraction method).
Since poles are complex conjugate, constants are also complex conjugate. So,
2𝐶𝑟 − [(𝐶𝑟 cos(𝑝𝑖 𝑇) + 𝐶𝑖 sin ( 𝑝𝑖 𝑇)]2𝑒 𝑝𝑟 𝑇 𝑧 −1
𝐻(𝑧) =
1 − 2𝑒 𝑝𝑟 𝑇 cos(𝑝𝑖 𝑇) 𝑧 −1 + 𝑒 2𝑝𝑟 𝑇 𝑧 −2
Here, 𝐶𝑟 = 0, 𝑝𝑖 𝑇 = 0.5207, 𝐶𝑖 = −666.4324, 𝑝𝑟 𝑇 = −0.5207, 𝑒 𝑝𝑟𝑇 = 0.5941, sin ( 𝑝𝑖 𝑇) =
0.4974, cos ( 𝑝𝑖 𝑇) = 0.8675, 𝑒 2𝑝𝑟𝑇 = 0.3530. Now,
393.9264𝑧 −1
𝐻(𝑧) =
1 − 1.0308𝑧 −1 + 0.3530𝑧 −2
If we substitute 𝑧 = 𝑒 𝑗𝜔𝑇 and 𝜔 = 0, we get 𝐻(𝑧) = 1223 which is close to sampling frequency so dividing
numerator by sampling frequency we get,
0.30784𝑧 −1
𝐻(𝑧) =
1 − 1.0308𝑧 −1 + 0.3530𝑧 −2
Here coefficients are 𝑏0 = 0, 𝑏1 = 0.30784, 𝑎1 = −1.0308, 𝑎2 = 0.3530
IIR Coefficient Calculation: Bilinear z-Transform Method
We will apply this method by a numerical example.
** Obtain coefficients of low pass IIR filter using bilinear z-transform method based on followings.
Transfer function : 𝟏
𝑯(𝒔) =
𝟐
𝒔 + √𝟐 𝒔 + 𝟏
Cutoff frequency : 150 Hz
Sampling frequency : 1.28 kHz
[Hints for Calculator Use: Assume pi = 3.1416 when it is not under any sin or cos function and pi =
180 when it is under any sin or cos function]
Critical frequency is 𝜔𝑝 = 2𝜋𝑓𝑐 = 2 × 3.1416 × 150 = 942.778
𝜔𝑝 𝑇 𝑠
Pre-wrapped frequency is 𝜔𝑝′ = tan ( ) = 0.3857. Putting 𝑠 = ′ in the transfer function we get,
2 𝜔𝑝
2
𝜔𝑝′ 0.1488
𝐻(𝑠) = 2 =
𝑠 2 + √2 𝑠𝜔𝑝′ + 𝜔𝑝′ 𝑠2 + 0.5455𝑠 + 0.1488
𝑧−1
Applying BZT, i.e. 𝑠 = we get,
𝑧+1
0.0878(1 + 2𝑧 −1 + 𝑧 −2 )
1 − 1.0048𝑧 −1 + 0.3561𝑧 −2