前言
最近看这两个滤波器设计嘛,就试着来写一下
主要从双线性变换,z变换,然后举一个例子来进行一下滤波器的设计。
例子
我们考虑传递函数如下的low pass filter

IIR滤波器
脉冲响应不变法
考虑这个传递函数H(s)=w02s2+w0Qs+w02H(s)=\frac{w_0^2}{s^2+\frac{w_0}{Q}s+w_0^2}H(s)=s2+Qw0s+w02w02
为了更好地将H(s)拆分开来,我们考虑分母的因式分解
假设H(s)H(s)H(s)存在两个根,解一元二次方程,有:
p1,2=−w02Q±1−14Q2 p_{1,2}=-\frac{w_0}{2Q} \pm \sqrt{1-\frac{1}{4Q^2}} p1,2=−2Qw0±1−4Q21
为了方便计算,我们记
α=−w02Q \alpha =-\frac{w_0}{2Q} α=−2Qw0
β=w01−14Q2 \beta = w_0\sqrt{1-\frac{1}{4Q^2}} β=w01−4Q21
因此我们知道
p1=α+jβ p_1=\alpha + j\beta p1=α+jβ
p2=α−jβ p_2=\alpha-j\beta p2=α−jβ
那我们的H(s)此时变成了:
H(s)=w02(s−p1)(s−p2) H(s)=\frac{w_0^2}{(s-p_1)(s-p_2)} H(s)=(s−p1)(s−p2)w02
我们假设这个式子可以展开,那么我们可以得到H(s)如下:
H(s)=As−p1+Bs−p2 H(s)=\frac{A}{s-p1} + \frac{B}{s-p_2} H(s)=s−p1A+s−p2B
这里A和B是待定的,我们有
w02(s−p1)(s−p2)=As−p1+Bs−p2 \frac{w_0^2}{(s-p_1)(s-p_2)}=\frac{A}{s-p_1}+\frac{B}{s-p_2} (s−p1)(s−p2)w02=s−p1A+s−p2B
这里我们其实可以计算出来,两边同时乘以分母,就有
A(s−p2)+B(s−p1)=w02 A(s-p_2)+B(s-p_1)=w_0^2 A(s−p2)+B(s−p1)=w02
求得A和B如下:
A=w02p1−p2B=w02p2−p1 A = \frac{w_0^2}{p_1-p_2} \\ B=\frac{w_0^2}{p_2-p_1} A=p1−p2w02B=

最低0.47元/天 解锁文章
1138

被折叠的 条评论
为什么被折叠?



