双线性变换法

给定模拟信号x(t),要用数字滤波器对它进行滤波。按照双线性变换法,步骤如下(假设为高通滤波器):

(1)选择一个采样频率F_s,对x(t)进行采样,得到时域离散信号x(n)

(2)确定模拟滤波的指标:通带边界频率f_p,阻带边界频率f_s,通带最大衰减\alpha_p,阻带最小衰减\alpha_s

(3)转换成数字滤波器指标。由于F_s对应于2\pi,因此\omega_p=2\pi f_p/F_s, \omega_s=2\pi f_s/F_s\alpha_p\alpha_s不变。

(4)设计模拟滤波器。

        (4.1)这里又要转换成模拟滤波器指标。但应用的是预畸变校正公式:\Omega_p = \frac{2}{T}tan(\frac{\omega_p}{2}),\Omega_s = \frac{2}{T}tan(\frac{\omega_s}{2})。注意这里的T是任意选取的,并不等于前面的1/F_s

        (4.2)设计模拟滤波器,比如巴特沃斯滤波器,得到其系统函数H_a(s)

(5)将H_a(s)转换成数字滤波器系统函数:


从上面的步骤里可以看到,确定了模拟指标,先转成了数字指标,又转成模拟指标。那能不能直接根据第(2)步的模拟指标来设计模拟滤波器(步骤4.2),然后执行步骤(5)?假设下图蓝色曲线是信号的频谱,现在要把右边的一部分滤出来。按照简化的步骤为:

  1. 确定采样频率F_s、边界频率指标f_pf_s
  2. 设计模拟滤波器,得到橙色的幅度响应曲线。
  3. 用双线性变换法转换成数字滤波器。注意模拟频率和数字频率之间的映射关系为:\omega=2arctan(2\pi fT/2)。所以,f_pf_s分别对应于\omega_p=2arctan(2\pi f_p T/2)\omega_s=2arctan(2\pi f_s T/2),而f=+\infty对应于\omega=\pi
  4. 而我们对模拟信号采样后,模拟边界频率对应的数字频率又是\omega=2\pi f/F_s的关系。所以f_pf_s分别对应于\omega'_p=2\pi f_p/F_s\omega'_s=2\pi f_s/F_s,而f=F_s/2对应于\omega=\pi。可以看到,\omega'_p不等于\omega_p\omega'_s不等于\omega_s边界频率产生偏差

但是,如果令第3步的T和第1步的F_s满足T=1/F_s,那么在f较小的时候,两个转换关系得到的\omega是比较接近的。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值