数字滤波器(三)--模拟滤波器的设计

•数字滤波器 专栏收录该内容
6 篇文章 0 订阅

数字滤波器(一)–IIR与FIR的基本结构与MATLAB实现
数字滤波器(二)–最小相位延时系统和全通系统

前言

IIR滤波也就是无限长单位冲击响应滤波器,基本结构和MATLAB实现可以参考本专栏第一篇博客:数字滤波器(一)–IIR与FIR的基本结构与MATLAB实现。 IIR滤波器的设计步骤主要由下图所示的三个大步骤组成。首先我们要设计出模拟原型滤波器,然后将滤波器数字化为相应的数字低通滤波器,然后在记性数字域的转换称为我们希望得到的选频滤波器(低通、高通、带通、带阻)。
在这里插入图片描述
再来说一下IIR滤波器的优缺点:

  • IIR滤波器的优点

    • 可以利用模拟滤波器设计,而模拟滤波器的设计有大量图表可查,方便简单。
  • IIR滤波器的缺点

    • 相位是非线性的
    • 冲击响应无限长,所以不能做快速卷积运算
    • 存在稳定性问题

    本篇博客主要讲述了怎么设计以及实现常用的模拟滤波器,巴特沃斯滤波器与切比雪夫Ⅰ型滤波器。

1. 模拟滤波器的系统函数

模拟滤波器的系统函数通常用 H a ( s ) H_a(s) Ha(s)表示,而模拟滤波器的响应可以用 H a ( j Ω ) H_a(jΩ) Ha(jΩ)表示,所以模拟滤波器的幅度平方函数的表达式为:
在这里插入图片描述
对于实滤波器,零极点在s平面是以共轭复数对的形式出现,如下图所示:
在这里插入图片描述
我们来看一下s平面和z平面的对应关系:

在这里插入图片描述
s平面中 σ < 0 σ<0 σ<0的区域对应的就是z平面上的单位圆,为什么要讨论这个关系呢?因为滤波器肯定都要求是因果稳定的系统,其系统函数在z平面的极点都在单位圆内部,所以对应在s平面,极点应该分布在左半平面。如果要求是最小相位延时系统,那么系统函数的零点也都应该分布在s平面的左半平面。

2 巴特沃斯滤波器

2.1 设计指标

首先我们了解一下巴特沃斯滤波器的形状,如下图所示:
在这里插入图片描述
从上图可以看出,巴特沃斯模拟滤波器具有最大平坦振幅特性,阻带单调变化。阶数N越大,曲线越抖,滤波器特性越好。

设计巴特沃斯滤波器需要考虑一下指标(指标对应下图):

  • N: 滤波器阶数
  • Ω p Ω_p Ωp:通带截止频率
  • Ω c Ω_c Ωc:滤波器截止频率,又称3dB带宽,其对应的赋值为 1 2 \frac{1}{\sqrt{2}} 2 1
  • Ω s t Ω_st Ωst: 阻带截止频率
  • δ p δ_p δp: 通带衰减(波纹)(dB)
  • δ s δ_s δs: 阻带衰减(波纹)(dB)
    在这里插入图片描述
2.2 设计流程
  • 第一步:根据指标要求确定滤波器阶数N(向上取整)与滤波器截止频率

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

  • 第二步: 根据阶数进行查表,得到归一化的系统函数 H a N ( s ) H_aN(s) HaN(s):
    在这里插入图片描述
    系统函数的归一化表达式为:

在这里插入图片描述

  • 第三步:将系统函数进行去归一化,得到最终的系统函数:
    在这里插入图片描述

下面我们通过一个例子来感受一下设计巴特沃斯滤波器的过程:

在这里插入图片描述
在这里插入图片描述

3. 切比雪夫Ⅰ型滤波器

3.1 设计指标

切比雪夫模拟滤波器分为Ⅰ型和Ⅱ型两种,我们来看下一两种滤波器的响应曲线。

  • 切比雪夫Ⅰ型
    振幅特征在通带内是等波纹的,在阻带内单调下降:
    在这里插入图片描述
    对于滤波器的阶数N,为奇数时和偶数时滤波器的通带波纹是不一样的,N为奇数时, Ω = 0 Ω=0 Ω=0时的幅度为1,N为偶数时, Ω = 0 Ω=0 Ω=0时的幅度为 1 1 + ε 2 \frac{1}{\sqrt{1+ε^2}} 1+ε2 1
  • 切比雪夫Ⅱ型
    振幅特性在阻带内是等波纹的,通带内是平稳的:
    在这里插入图片描述
    本篇博文重点讨论切比雪夫Ⅰ型滤波器。
    设计比雪夫Ⅰ型滤波器需要考虑以下指标:
  • N N N:滤波器阶数
  • Ω c Ω_c Ωc:通带宽度
  • δ p δ_p δp: 通带衰减(波纹)(dB)
  • δ s δ_s δs: 阻带衰减(波纹)(dB)
3.2 设计流程
  • 第一步:根据指标要求确定滤波器阶数
    滤波器的阶数可以由通带、阻带衰减确定:
    在这里插入图片描述
    其中双曲余弦函数 a r c c h ( x ) = l n ( x + x 2 − 1 ) arcch(x)=ln(x+\sqrt{x^2-1}) arcch(x)=ln(x+x21 )

  • 第二步:根据阶数和通带波纹进行查表
    此时截止频率是归一化的:
    在这里插入图片描述
    通过查表可以确定模拟滤波器归一化系统函数:
    在这里插入图片描述
    其中的系数 d 0 d_0 d0可以通过 H a ( j 0 ) = H a N ( s ) ∣ s = 0 H_a(j0)=H_{aN}(s)|_{s=0} Ha(j0)=HaN(s)s=0确定, 当阶数N为奇数时, H a N ( s ) ∣ s = 0 = 1 H_{aN}(s)|_{s=0}=1 HaN(s)s=0=1,当阶数N为偶数时, H a N ( s ) ∣ s = 0 = 1 1 + ε 2 H_{aN}(s)|_{s=0}=\frac{1}{\sqrt{1+ε^2}} HaN(s)s=0=1+ε2 1,其中参数 ε ε ε可以通过公式: ε 2 = 1 0 0.1 δ p − 1 ε^2=10^{0.1δ_p}-1 ε2=100.1δp1 来确定。

  • 第三步:将系统函数进行去归一化,得到最终的系统函数
    在这里插入图片描述
    下面我们通过一个例子来感受一下设计切比雪夫Ⅰ型滤波器的过程:
    在这里插入图片描述
    在这里插入图片描述

  • 4
    点赞
  • 0
    评论
  • 8
    收藏
  • 一键三连
    一键三连
  • 扫一扫,分享海报

打赏
文章很值,打赏犒劳作者一下
相关推荐
©️2020 CSDN 皮肤主题: 技术黑板 设计师:CSDN官方博客 返回首页

打赏

爱吃骨头的猫、

你的鼓励将是我创作的最大动力

¥2 ¥4 ¥6 ¥10 ¥20
输入1-500的整数
余额支付 (余额:-- )
扫码支付
扫码支付:¥2
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值