【数据压缩10】完全重建QMF滤波器组的设计

一、Motivation

很多信源在同一时间表现表现出不同的特性（人对不同失真的敏感度不同），不能找到能满足所有情况的一个模型。因此需要将信号分解处理后再合成，采用不同的编码机制，使编码更灵活，从而编码效率更高。信号分解可利用数字滤波器实现。

二、完全重建QMFB滤波器

最常用的滤波器是正交镜像滤波器(Quadrature mirror filters, QMF)。完全重建滤波器QMFB将输入信号经由不同通带的滤波器组滤波之后输出，在接收侧通过综合滤波器将各子带综合为恢复信号。其最大的优点是在对信号进行抽取后，可以根据每个子带的不同特征分别进行处理，而插值和合成环节又能消除信号失真的各种因素。

如果综合滤波器组最后输出的信号与分析滤波器组原来输入的信号有关系：

其中，c和n0均为固定的常数。称是对的完全重建（Perfect Reconstruction），简称PR。也就是说，输出信号是输入信号的延迟样本，在幅度上它们只相差一个固定的倍数，在时间上有一个固定的延迟。对应的频域关系为：

三、两通道正交镜像滤波器理论

一个两通道正交镜像滤波器组如图所示，在分析滤波器组一侧，输入信号（设为宽带信号）被分成K个子频带信号（窄带信号），通过抽取可降低采样率；在综合滤波器一侧，通过零值内插和带通滤波可以重建原来的信号。

对于一个给定的信号，经过分析滤波器后，再进行抽取、编码、传输，可以通过零值内插、综合滤波器滤波、求和运算得到恢复和重建。但是重建后的信号并不能与原始信号完全相同，两者之间存在着误差，主要包括：

1. 混叠失真：由抽取和内插产生的混叠和镜像带来的误差，导致分析滤波器组和综合滤波器组的频带不能完全分开；
2. 幅度失真：由于分析和综合滤波器组的频带在通带内不是全通函数，其幅频特性波纹产生的误差；
3. 相位失真：由滤波器相频特性的非线性所产生的误差；
4. 量化失真：由编、解码产生的误差，与量化噪声相似，这类误差无法完全消除，只能设法减小。

因此，在设计QMF组时，就需要综合考虑如何减小和消除上述的各类误差。消除混叠失真一种简单形式采取：

当两通道无混叠滤波器组的分解滤波器满足：时，该滤波器组为正交镜像滤波器组（quadrature mirror filter bank，QMFB）。且滤波器的幅度特性满足：

四、仿真实现

各滤波器的阶数N、滤波器h0的通带截止频率w、均方误差MSE

1. 确定滤波器的阶数N值。选定一个合适的w值，在不同输入信号的前提下，不断改变N值，并求出此时所对应的MSE，比较找到MSE最小时的N。N=41时MSE最小，重建效果最好。
2. 确定滤波器h0的通带截止频率w值，选定一个合适的N值，在不同输入信号的前提下，不断改变w值，同样是比较找到最小的MSE。w=0.43时，重建效果最好。

滤波器H0和滤波器H1的幅度响应，红色高通滤波器H1，蓝色低通滤波器H2，呈镜像对称：

 滤波器H0和滤波器H1的幅度响应差值：

 h0的幅频响应，达到了通带尽量平，过渡带尽量窄，且阻带尽可能快速衰减的要求：

 输入信号：

 理想输出信号与重建的输出信号：

 理想输出信号与重建输出信号的偏差：