数字滤波器(二)--最小相位延时系统和全通系统

数字滤波器(一)–IIR与FIR的基本结构与MATLAB实现

1. 系统相位

当输入信号经过一个滤波器系统的时候,其相位可能会发生一定的移动,那么一个系统的相位到底跟什么有关系呢?对于给定的系统:
在这里插入图片描述
其系统相位的延迟可以这么表达:
在这里插入图片描述

2. 最小相位延时系统

2.1 最小相位延时系统

什么是最小相位延时系统呢?顾名思义,就是对输入信号的相位移动延时最小的系统。
对于第一节中的系统 H ( z ) H(z) H(z),我们假设系统有M个零点,N个极点,我们可以定义一下参数:
m i m_i mi: 单位圆内的零点数
m o m_o mo:单位圆外的零点数
p i p_i pi: 单位圆内的极点数
p o p_o po: 单位圆外的极点数
所以我们可以得到 m i + m o = M m_i+m_o=M mi+mo=M, p i + p o = N p_i+p_o=N pi+po=N

对于因果稳定系统,其系统函数的全部极点都在单位圆内(划重点,要考的!),所以 p i = N , p o = 0 p_i=N, p_o=0 pi=N,po=0

我们来观察一下系统零极点对相位的影响:
当数字频率ω从0变化到2π时,相位的变化可以表示为:
在这里插入图片描述
此时幅角变化为负,因此相位产生了延时,此时系统称为相位延时系统。相位角的变化大小与单位圆外的零点数目有关

  • 当零点全部位于单位圆之内时 m i = M , m o = 0 m_i=M,m_o=0 mi=M,mo=0,此时相位变化最小:
    在这里插入图片描述
    此系统为最小相位延迟系统
  • 当零点全部位于单位圆之外时 m i = 0 , m o = M m_i=0,m_o=M mi=0,mo=M,此时相位变化最大:
    在这里插入图片描述
    称此系统为最大相位延迟系统。

对于延时系统的定义,我们可以参照此表:(我们可以通过判断系统的零点和极点相对于单位圆的位置来对延迟系统进行分类)
在这里插入图片描述

2.2 最小相位延时系统的性质

最小相位延迟系统有怎么样的性质呢?

  • 在幅频响应相同的所有系统中,最小相位延迟系统的单位序列响应的能量集中在n=0附近。
  • 在幅频响应相同的所有系统中,最小相位延迟系统具有最小的群延时。
  • 任何因果稳定的非最小相位延时系统 H ( z ) H(z) H(z)都可以表示为全通系统与最小相位延时系统的级联。(这个性质在后面讲解滤波器设计的时候很重要)

3. 全通系统

3.1 定义与系统函数

系统频率响应的幅度在所有频率下均为1或常数,这样的系统称为全通系统,也就是说全通系统是指不改变信号幅度谱,只改变相位谱的信号处理系统。

全通系统的系统函数可以表示为:
在这里插入图片描述
其中 D ( z ) = 1 + a 1 z − 1 + a 2 z − 2 + . . . + a N z − N D(z)=1+a_1z^{-1}+a_2z^{-2}+...+a_Nz^{-N} D(z)=1+a1z1+a2z2+...+aNzN,是关于 z − 1 的 N 次 多 项 式 z^{-1}的N次多项式 z1N,N是系统的阶数。

一阶全通系统的系统函数可以表示为:
在这里插入图片描述
二阶全通系统的系统函数可以表示为:
在这里插入图片描述

3.2 全通系统的零极点关系

对于全通系统:
在这里插入图片描述
我们可以知道,如果系统函数存在极点 z = z 1 = r e j θ z=z_1=re^{jθ} z=z1=rejθ,则系统函数一定存在零点 z = 1 z 1 = 1 r e − j θ z=\frac{1}{z_1}=\frac{1}{r}e^{-jθ} z=z11=r1ejθ
因为系统函数分子分母的系数均为实数,因此零极点均以共轭对的形式出现,所以也会存在对应的共轭极点 z = z 1 ∗ = r e − j θ z=z_1^*=re^{-jθ} z=z1=rejθ和共轭零点 z = 1 z 1 ∗ = 1 r e j θ z=\frac{1}{z_1^*}=\frac{1}{r}e^{jθ} z=z11=r1ejθ.

总结一下,全通系统的零极点可以表示为:
在这里插入图片描述
将这两对零极点作为零极点分布图,可以表示为:
在这里插入图片描述
可以发现,全通系统的零极点出现在共轭反演的位置上。

我们可以通过一个例子来熟悉一下全通系统零极点之间的关系:

在这里插入图片描述
首先,我们将系统的差分方程转化为系统传递函数为:
在这里插入图片描述
从系统函数可以得到该系统的零点为 z 0 = − b z_0=-b z0=b, 极点 z p = a z_p=a zp=a, 根据全通系统零极点互为共轭反演的关系,即 z 0 = 1 z p z_0=\frac{1}{z_p} z0=zp1,可以得到b需要满足以下条件:
在这里插入图片描述

3.3 全通系统的应用

我们可以通过零极点对消的方式,实现系统之间的转化
任何因果稳定的最小相位延时系统,都可以表示为全通系统与最小相位延时系统的级联,即:
在这里插入图片描述
所得系统函数与最小相位系统的差别,只是频率响应的相位不同。我们可以通过下面这个零极点图来表示该转化关系:
在这里插入图片描述
从上图我们可以看出,全通系统中的两个极点(红色叉标注)与最小相位延时系统中的两个对应的零点(红色圈标注)互相抵消了

可以使用级联全通系统的方法,是的一个非稳定的滤波器转化为稳定的系统
我们通过零极点图来反应这一关系:
在这里插入图片描述
因为稳定系统的极点必须都在单位圆内,上图中非稳定系统的蓝色标记的极点显然在圆外面,我们可以通过级联一个全通系统来抵消这个极点。因为全通系统的零极点是共轭反演的,所以两个红色圈标记的零点对应两个单位圆内的零点。通过级联之后,全通系统的单位圆外的零点与非稳定系统的单位圆外的极点相互抵消,从而使系统稳定。

我们可以通过一个例子来更加理解全通系统的应用:

在这里插入图片描述
根据系统传递函数,我们可以得到该系统的零点为:
在这里插入图片描述
该系统的极点为:
在这里插入图片描述
表示在单位圆上为:
在这里插入图片描述
它可以表示为一个最小相位延时系统和一个全通系统的级联:
在这里插入图片描述
最小相位延时系统中的两个红色标记的零点和全通系统中两个红色标记的极点通过级联相互抵消。用公式表示上述关系:
在这里插入图片描述
其中红色框表示最小相位延时系统的传递函数,蓝色框表示全通系统的传递函数。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

爱吃骨头的猫、

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

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

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

打赏作者

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

抵扣说明:

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

余额充值