中科院现代信号处理的理论和方法Chapter+课件

第四章多抽样率信号处理

4.1信号的抽取与插值4.2滤波器组基础

4.3两通道滤波器组14.1信号的抽取与插值引言信号的抽取信号的插值抽取与插值相结合的抽样率转换信号的多相表示2引言“多抽样率数字信号处理”的核心内容是信号抽样率的转换及滤波器组。抽取与插值滤波器组3信号的抽取减少抽样率以去掉过多数据的过程称为信号的“抽取”。时域表示：4抽取的框图

注：见胡广书《现代信号处理教程》图5.2.15频域表示：67将信号作M倍的抽取后，所得信号

的频谱等于原信号

的频谱先作M倍的扩展，再在轴上作的移位，幅度降为原来的1/M后再叠加。8图5.2.2信号抽取后频谱的变化注：见胡广书《现代信号处理教程》图5.2.2M=39注：见杨绿溪《现代数字信号处理》图8.3M=210抽取滤波器的频响特性：由于M是可变的，所以很难要求在不同的M下都能保证。防止抽取后在中出现混叠的方法是在对抽取前先进行低通滤波，压缩其频带。1112图5.2.4先滤波再抽取后的频谱的变化注：见胡广书《现代信号处理教程》图5.2.413信号的插值增加抽样率以增加数据的过程称为信号的“插值”。时域表示：14图5.3.1信号的插值注：见胡广书《现代信号处理教程》图5.3.115频域表示：16图5.3.2插值后对频域的影响

注：见胡广书《现代信号处理教程》图5.3.217插值滤波器的频响特性：

实现有效插值的方法是将再通过一个低通滤波器，即图5.3.3插值后的滤波

注：见胡广书《现代信号处理教程》图5.3.3181920抽取与插值相结合的抽样率转换图5.4.1插值后的滤波

注：见胡广书《现代信号处理教程》图5.4.12122信号的多相表示使用多相表示可在抽样率转换的过程中去掉许多不必要的计算，因而大大提高运算的速度。

23图5.2.1信号抽取示意图，M=3,横坐标为抽样点数注：见胡广书《现代信号处理教程》图5.2.1

2425例如：2627恒等关系

两个信号分别定标以后再相加后的抽取等于它们各自抽取后再定标和相加。28信号延迟M个样本后作M倍抽取和先抽取再延迟一个样本是等效的。

在M倍抽取器的前后，滤波器的z的幂相差M倍。

29注：见胡广书《现代信号处理教程》图5.6.1-图5.6.630基于多相分解的抽取和插值滤波器实现抽取的滤波器实现：31

注：见胡广书《现代信号处理教程》图5.7.132基于多相结构的信号抽取：33注：见胡广书《现代信号处理教程》图5.7.3344.2滤波器组基础滤波器组的基本概念滤波器组的种类及有关滤波器35滤波器组的基本概念一个滤波器组是指一组滤波器，它们有着共同的输入，或有着共同的输出。

注：见胡广书《现代信号处理教程》图6.1.1363738失真来源：混叠失真；幅度和相位失真；抽取后处理产生的误差。39滤波器组的种类及有关滤波器最大均匀抽取滤波器组正交镜像滤波器组第M带滤波器半带滤波器互补型滤波器401、最大均匀抽取滤波器组均匀滤波器组：设某一滤波器有K个分析滤波器H0(z),…,HK-1(z)，这K个滤波器有如下关系：即最大均匀抽取滤波器组：如果做M倍的抽取，且M=K，则称该滤波器组为最大均匀抽取滤波器组。41DFT调制滤波器组（复系数的M通道均匀滤波器组）：42余弦调制滤波器组（实系数的M通道均匀滤波器组）：43442、正交镜像滤波器组45463、第M带滤波器将分析滤波器组写成多相形式，如果其第0相，也即恒为一常数c，即

那么，其单位抽样响应必有满足上式的滤波器称为“第M带滤波器（Mthfilter）”。

47某一Mth滤波器的单位抽样响应，M=34849定理6.1

若H(z)是一个Mth滤波器，则

如果一个Mth滤波器h(n)，那么将其依次移位后，所得到的M个滤波器的频率响应之和等于1，效果为一全通滤波器。Mth滤波器的这一特性被用于设计M通道滤波器组。504、半带滤波器Mth滤波器中，若令M=2，则所得的滤波器称为半带滤波器，此时51某一半带滤波器的h(n)

52若假定c=1/2，则5354半带滤波器的特点：H(z)的通带与阻带的纹波必须相等，即频率响应关于半带频率是对称的，通带和阻带截止频率相对是等距的；除n=0外，h(n)的所有偶序号项全为0.半带滤波器的长度总是奇数，即可有效的减少计算量。555、互补型滤波器严格互补滤波器：56功率互补滤波器：57多抽样率系统的应用简介语音和图像的子带编码：

在对语音和图像编码时，一个常用的方法是利用滤波器组将信号分解成一个个的子带信号，并令这些子带信号占据不同的频带。由于这些子带信号所具有的能量是不会相同的，因此可按其能量的大小（也即重要性）给出不同的字长，从而达到高效编码的目的。585960时分复用和频分复用：61624.3两通道滤波器组注：见胡广书《现代信号处理教程》图7.1.163各信号间的关系：对图中的分析滤波器组

64对综合滤波器组

6566讨论：若T(z)是全通系统，则整个滤波器组将不会发生幅度失真；若T(z)具有线性相位，则该滤波器组将不会发生相位失真；最简单的情况是令T(z)为纯延迟，即T(z)=cz-k，那么从而实现了PR。存在问题：如何保证F(z)=0?如何保证T(z)为全通系统？如何保证T(z)具有线性相位？能否保证T(z)=cz-k？6768无混叠条件：69707172G0(z)和G1(z)的选取准则：73标准正交镜像滤波器组为了保证G0(z)和G1(z)是FIR的及实现准确重建，需要选择合适的H0(z)和H1(z)。7475多相结构描述的两通道滤波器组中的基本关系：7677QMFB中的幅度失真与相位失真：78存在问题：在要求H0(z)和H1(z)都是FIR的情况下，若想既要保持PR，又要使H0(z)和H1(z)具有实际意义是不可能的。在H0(z)和H1(z)有实际意义的情况下，不可能完全消除幅度失真，因此也做不到PR。解决方法：在去除相位失真的前提下，尽可能地减小幅度失真，从而做到近似PR；－－FIRQMF滤波器组去除幅度失真，不考虑相位失真，近似实现PR。－－IIRQMF滤波器组；放弃H1(z)=H0(-z)的简单形式，取更为合理的形式，从而实现PR。－－CQMFB共轭正交镜像滤波器组79FIR标准正交镜像滤波器组在去除相位失真的前提下，尽可能地减小幅度失真，从而做到近似PR。FIR滤波器容易设计成线性相位的，例如，假定H0(z)是N点FIR低通滤波器，即808182IIR标准正交镜像滤波器组去除幅度失真，不考虑相位失真，近似实现PR。8384共轭正交镜像滤波器组放弃H1(z)=H0(-z)的简单形式，取更为合理的形式，从而实现PR。其中，N为偶数。858687

谱分解定理如果功率谱Pxx(ejω)是平稳随机序列x(n)的有理谱，那么一定存在一个零极点均在单位圆内的有理函数H(z),满足式中，ak,bk都是实数，a0=b0=1,且|αk|＜1,|βk|＜1。88我们总可以用单位圆内部的零极点组成一个系统H(z)（该系统自然是最小相位系统），又因为系统系数是实数，圆外的零极点必定与圆内的零极点共轭对称。这样除了单位圆内部零极点外，用其它零极点组成的系统函数必定是H(z-1)。按照谱分解定理分解出H(z)一定是最小相位系统，它保证了模型的可逆性，即逆系统存在。89

分解方法：我们知道功率谱是cosω的函数，为了对功率谱进行谱分解，下面介绍一种分解方法：（1）用φ代替cosω，得到有理函数V(φ)；（2）求出V(φ)分子、分母的全部根φi;

（3）构造对每个φi的方程:该方程有两个根：Zi和1/Zi

，其中Zi是单位圆内的根；90（4）用单位圆内部极零点构成H(z)，零点是分子多项式的根Zi，极点是分母多项式的根Zj，常数C由功率谱Pxx(ejω)确定。91例1.6.2

已知x(n)的功率谱求其模型的系统函数。解：（1）令φ=cosω，则(2)92(3)(4)93设σ2w=1，对比给定的功率谱，得到：C=2/3，模型系统函数为也可以假定σ2w=4/9，此时模型系统函数为这样得到的模型系统函数的常数因子不同，但常数因子仅影响其幅度大小，不影响问题实质。94由谱分解求H0(z)和H1(z)：9596两通道滤波器组中的制约关系：979899讨论：对一个两通道滤波器组，希望其具有PR性能，四个滤波器具有好的通带和阻带，且都具有线性相位；不论是QMFB还是CQMFB，四个滤波器实际上都来自于同一个滤波器H0(z)，要寻求既具有线性相位又满足PR性能的滤波器组，就不能简单地由H0(z)得到H1(z)。放弃即放弃H0(z)和H1(z)之间的正交性，而是满足双正交关系。100树状滤波器组规则树状滤波器组101非均匀树状滤波器组非均匀综合滤波器组102等效结构103信号的多分辨率分析对频