基于提升结构的两通道双正交图滤波器组设计方法

摘要

本发明公开一种基于提升结构的两通道双正交图滤波器组设计方法，其基于单级提升结构，能有效控制分析和综合频谱核的频率特性。其利用单级提升结构，将分析和综合频谱核表示为关于两个提升滤波器的函数，使其自动满足重构条件，从而将对分析和综合子带滤波器的优化转化为对两个提升滤波器的优化，即分析和综合滤波器的频率特性可描述为关于提升滤波器的函数，通过求解带约束的二次规划问题优化，以获得良好的频率特性。本发明方法在保证完全重构特性的前提下获得良好了的频率特性。

权利要求

1.基于提升结构的两通道双正交图滤波器组设计方法，其特征是，包含如下步骤：步骤1、将两通道图滤波器组的分析子带滤波器Hi和综合子带滤波器Gi分别表示为关于分析频谱核和综合频谱核的函数；即      ①

式中，Hi为第i个分析子带滤波器，Gi为第i个综合子带滤波器，hi(λ)为第i个分析子带滤波器频谱核在特征根λ处的取值，gi(λ)为第i个综合子带滤波器频谱核在特征根λ处的取值，Pλ为特征根λ对应的正交投影矩阵，λ为图G拉普拉斯矩阵的特征根，σ(G)为图G拉普拉斯矩阵的特征根集合，i是子带滤波器的序号；

步骤2、进行变量替换x＝λ-1，并引入单级提升结构，将第0个分析子带滤波器频谱核在变量x处的取值 和第0个综合子带滤波器频谱核在变量x处的取值 表示为关于分析提升滤波器和综合提升滤波器的函数，即：   ②

式中， 为第0个分析子带滤波器频谱核在变量x处的取值， 为第0个综合子带滤波器频谱核在变量x处的取值，s(x2)为分析提升滤波器，t(x2)为综合提升滤波器；

步骤3、优化分析提升滤波器s(x2)，即求解以下优化问题：

          ③

式中，λs是阻带截止频率，ch(x)是分析子带滤波器频谱核对应的频率向量，s是由分析提升滤波器系数构成的向量，Uh是分析子带滤波器频谱核对应的阻带能量矩阵，rh是分析子带滤波器频谱核对应的阻带能量向量；

2 2

步骤4、固定分析提升滤波器s(x)，优化综合提升滤波器t(x)，即求解以下优化问题：            ④

式中，λs是阻带截止频率， 是第0个综合子带滤波器频谱核在变量x处的取值，t是由综合提升滤波器系数构成的向量，Ug是综合子带滤波器频谱核对应的阻带能量矩阵，rg是综合子带滤波器频谱核对应的阻带能量向量，cg(x)是综合子带滤波器频谱核对应的频率向量，ε是容忍因子；

步骤5、将步骤3优化得到的分析提升滤波器s(x2)和步骤4优化得到的综合提升滤波器t(x2)带入公式②中，得到第0个分析子带滤波器频谱核在变量x处的取值 和第0个综合子带滤波器频谱核在变量x处的取值步骤6、进行变量反替换λ＝x+1，得到第0个分析子带滤波器频谱核在特征根λ处的取值h0(λ)和第0个综合子带滤波器频谱核在特征根λ处的g0(λ)；

步骤7、将步骤6所得的第0个分析子带滤波器频谱核在特征根λ处的取值h0(λ)和第0个综合子带滤波器频谱核在特征根λ处的g0(λ)带入公式⑤中，得到第1个分析子带滤波器频谱核在特征根λ处的取值h1(λ)和第1个综合子带滤波器频谱核在特征根λ处的g1(λ)，即h1(λ)＝g0(2-λ)，g1(λ)＝h0(2-λ)        ⑤。

说明书

基于提升结构的两通道双正交图滤波器组设计方法

技术领域

[0001] 本发明涉及多速率信号处理领域，具体涉及一种基于提升结构的两通道双正交图滤波器组设计方法。

背景技术

[0002] 图(graph)是一种有效描述非规则的几何结构的模型。社交网络、无线传感器网络、金融网络等均可用图进行描述，图的节点表示数据所在的位置，数据的大小可以用图节点信号来表示，而图中节点之间的连线可以表示不同数据节点之间的关联性。然而，在实际应用中，图数据量十分庞大，对全图的处理将带来十分巨大的计算复杂度，往往无法实现，从而图信号的处理成为了关键，特别是图信号的多分辨分析。

[0003] 目前，有学者提出了适用于图信号处理的小波变换/滤波器组。典型的例子有：适用于交通网络图的类小波变换，适用于无线传感器网络图的两通道可逆滤波器组。在这些工作中，两通道双正交图滤波器组是最为优秀的结构之一，其具备临界采样、紧支撑、完全重构等优点。然而，目前却只有一种方法用于该类结构的设计，但由于该方法基于多项式分解，因而设计中无法兼顾子带滤波器的频率特性。

发明内容

[0004] 本发明所要解决的技术问题是提供一种基于提升结构的两通道双正交图滤波器组设计方法，其所设计的滤波器能在保证重构特性的前提下获得良好的频率特性。

[0005] 为解决上述问题，本发明是通过以下技术方案实现的：

[0006] 基于提升结构的两通道双正交图滤波器组设计方法，包含如下步骤：

[0007] 步骤1、将两通道图滤波器组的分析子带滤波器Hi和综合子带滤波器Gi分别表示为关于分析频谱核和综合频谱核的函数；即

[0008]    ①

[0009] 式中，Hi为第i个分析子带滤波器，Gi为第i个综合子带滤波器，hi(λ)为第i个分析子带滤波器频谱核在特征根λ处的取值，gi(λ)为第i个综合子带滤波器频谱核在特征根λ处的取值，Pλ为特征根λ对应的正交投影矩阵，λ为图G拉普拉斯矩阵的特征根，σ(G)为图G拉普拉斯矩阵的特征根集合，i是子带滤波器的序号；

[0010] 步骤2、进行变量替换x＝λ-1，并引入单级提升结构，将第0个分析子带滤波器频谱核在变量x处的取值 和第0个综合子带滤波器频谱核在变量x处的取值 表示为关于分析提升滤波器和综合提升滤波器的函数，即：

[0011]    ②

[0012] 式中， 为第0个分析子带滤波器频谱核在变量x处的取值， 为第0个综合子带滤波器频谱核在变量x处的取值，s(x2)为分析提升滤波器，t(x2)为综合提升滤波器；

[0013] 步骤3、优化分析提升滤波器s(x2)，即求解以下优化问题：

[0014]    ③

[0015] 式中，λs是阻带截止频率，ch(x)是分析子带滤波器频谱核对应的频率向量，s是由分析提升滤波器系数构成的向量，Uh是分析子带滤波器频谱核对应的阻带能量矩阵，rh是分析子带滤波器频谱核对应的阻带能量向量；

[0016] 步骤4、固定分析提升滤波器s(x2)，优化综合提升滤波器t(x2)，即求解以下优化问题：

[0017]    ④

[0018] 式中，λs是阻带截止频率， 是第0个综合子带滤波器频谱核在变量x处的取值，t是由综合提升滤波器系数构成的向量，Ug是综合子带滤波器频谱核对应的阻带能量矩阵，rg是综合子带滤波器频谱核对应的阻带能量向量，cg(x)是综合子带滤波器频谱核对应的频率向量，ε是容忍因子；

[0019] 步骤5、将步骤3优化得到的分析提升滤波器s(x2)和步骤4优化得到的综合提升滤2
波器t(x)带入公式②中，得到第0个分析子带滤波器频谱核在变量x处的取值 和第0个综合子带滤波器频谱核在变量x处的取值

[0020] 步骤6、进行变量反替换λ＝x+1，得到第0个分析子带滤波器频谱核在特征根λ处的取值h0(λ)和第0个综合子带滤波器频谱核在特征根λ处的g0(λ)；

[0021] 步骤7、将步骤6所得的第0个分析子带滤波器频谱核在特征根λ处的取值h0(λ)和第0个综合子带滤波器频谱核在特征根λ处的g0(λ)带入公式⑤中，得到第1个分析子带滤波器频谱核在特征根λ处的取值h1(λ)和第1个综合子带滤波器频谱核在特征根λ处的g1(λ)，即[0022] h1(λ)＝g0(2-λ)，g1(λ)＝h0(2-λ)    ⑤

[0023] 与现有技术相比，本发明基于单级提升结构，能有效控制分析和综合频谱核的频率特性。其利用单级提升结构，将分析和综合频谱核表示为关于两个提升滤波器的函数，使其自动满足重构条件，从而将对分析和综合子带滤波器的优化转化为对两个提升滤波器的优化，即分析和综合滤波器的频率特性可描述为关于提升滤波器的函数，通过求解带约束的二次规划问题优化s(x),t(x)，以获得良好的频率特性。本发明方法在保证完全重构特性的前提下获得良好了的频率特性。

附图说明

[0024] 图1是图滤波器组的结构图。

[0025] 图2是子带滤波器长度较低时的优化结果图，其中(a)为现有方法设计所得的图滤波器，(b)为本发明方法设计所得的图滤波器。PRV表示重构特性。

[0026] 图3是子带滤波器长度较高时的优化结果图，其中(a)为现有算法设计所得的图滤波器，(b)为本发明算法设计所得的图滤波器。PRV表示重构特性。

[0027] 图4是本发明方法设计的图滤波器组对于密尼苏达交通网络的子带分解图，其中(a)为LL子带分解图，(b)为LH子带分解图，(c)为HL子带分解图，(d)为HH子带分解图。

具体实施方式

[0028] 一种基于提升结构的两通道双正交图滤波器组设计方法，其包含如下步骤：

[0029] 第一步：两通道图滤波器组的结构如图1所示。其分析和综合子带滤波器表示为关于分析和综合频谱核的函数：

[0030]

[0031] 式中，Hi为第i个分析子带滤波器，Gi为第i个综合子带滤波器，是hi(λ)为第i个分析子带滤波器频谱核在特征根λ处的取值，gi(λ)为第i个综合子带滤波器频谱核在特征根λ处的取值，Pλ为特征根λ对应的正交投影矩阵。λ为图G拉普拉斯矩阵的特征根，σ(G)图G拉普拉斯矩阵的特征根集合。i是子带滤波器的序号。

[0032] 如果满足h1(λ)＝g0(2-λ),g1(λ)＝h0(2-λ)，双正交条件为：

[0033]

[0034] 式中，h0(λ)为第0个分析子带滤波器频谱核在特征根λ处的取值，g0(λ)为第0个综合子带滤波器频谱核在特征根λ处的取值，h0(2-λ)为第0个分析子带滤波器频谱核在特征根2-λ处的取值，g0(2-λ)为第0个综合子带滤波器频谱核在特征根2-λ处的取值。进行变量替换x＝λ-1，并定义

[0035]

[0036] 式中， 为第0个分析子带滤波器频谱核在变量x处的取值， 为第0个综合子带滤波器频谱核在变量x处的取值。式(2)可写成：

[0037]

[0038] 式中， 为第0个分析子带滤波器频谱核在变量x处的取值， 为第0个综合子带滤波器频谱核在变量x处的取值。 为第0个分析子带滤波器频谱核在变量-x处的取值， 为第0个综合子带滤波器频谱核在变量-x处的取值。x是取值范围为[-1,1]的变量。

[0039] 第二步：引入单级提升结构，将 表示为关于提升滤波器的函数：

[0040]

[0041] 式中， 为第0个分析子带滤波器频谱核在变量x处的取值， 为第0个综合子带滤波器频谱核在变量x处的取值。s(x2)和t(x2)分别为分析和综合提升滤波器。其中，s(x),t(x)均为关于变量x的多项式：

[0042]

[0043] 式中，s(x)是分析提升滤波器，t(x)为综合提升滤波器，sn是提升滤波器s(x)的系数，Ls是提升滤波器s(x)的阶数，tn是提升滤波器t(x)的系数，Lt是提升滤波器t(x)的阶数。可验证，引入提升结构后，双正交条件(4)自动满足。

[0044] 第三步：优化s(x)，将问题归结为带约束的二次规划，目标函数为阻带能量，约束为通带平坦性

[0045]

[0046] 式中， 为第0个分析子带滤波器频谱核在变量x处的取值，sn是提升滤波器s(x)的系数，λs是阻带截止频率。将式(5)和(6)代入式(7)

[0047]

[0048] 式中，s是由提升滤波器s(x)系数构成的向量，Uh是分析子带滤波器频谱核对应的阻带能量矩阵，rh是分析子带滤波器频谱核对应的阻带能量向量，ch(x)是分析子带滤波器频谱核对应的频率向量。

[0049] 第四步：固定s(x)，优化t(x)，将问题归结为带约束的二次规划，目标函数为阻带能量，通带平坦约束 自动成立，考虑提升结构中，分析和综合频谱核的频率特性相互限制，故增加阻带波纹约束

[0050]

[0051] 将式(5)和(6)代入式(9)，可得

[0052]

[0053] 式中，ε是容忍因子，t是由提升滤波器t(x)系数构成的向量，Ug是综合子带滤波器频谱核对应的阻带能量矩阵，rg是综合子带滤波器频谱核对应的阻带能量向量，cg(x)是综合子带滤波器频谱核对应的频率向量。

[0054] 第五步：将第三、四步优化得到s(x)和t(x)代入式(5)，即得到 和[0055] 第六步：进行变量反替换λ＝x+1，得h0(λ)和g0(λ)。

[0056] 第七步：将h0(λ)和g0(λ)带入式(11)，得h1(λ)和g1(λ)

[0057] h1(λ)＝g0(2-λ)，g1(λ)＝h0(2-λ)    (11)

[0058] 上述s(x)和s(x2)只是对函数的变量进行了变换而已，本质上还是分析提升滤波器；同理，t(x)和t(x2)也只是对函数的变量进行了变换而已，本质上还是综合提升滤波器。

[0059] 为验证本方法的有效性，进行了仿真实验并与现有方法进行比较。仿真参数为：子带滤波器的通带和阻带截止频率设为：λp＝0.7,λs＝1.3。失真参数、子带滤波器长度和提升滤波器长度参数分两种情况，分别为：ε＝0.1,La＝10,Ls＝17,Ls＝4,Lt＝4,ε＝0.01,La＝18,Ls＝33,Ls＝8,Lt＝8。我们分别采用现有方法和本发明方法设计上述滤波器组。

[0060] 表1给出了两种方法的性能比较。设计所得的子带滤波器频率响应如图2和图3所示(两种长度情况)，图2是子带滤波器长度较低时的优化结果图，图3是子带滤波器长度较高时的优化结果图，其中(a)为现有方法设计所得的图滤波器，(b)为本发明方法设计所得的图滤波器。PRV表示重构特性。图4是本发明方法设计的图滤波器组对于密尼苏达交通网络的子带分解图，其中(a)为LL子带分解图，(b)为LH子带分解图，(c)为HL子带分解图，(d)为HH子带分解图。表2给出了本发明方法设计所得的提升滤波器的系数。

[0061] 表1

[0062]

[0063]

[0064] 表2

[0065]

[0066] 从仿真结果可以得出，本发明方法设计所得的子带滤波器在保证重构特性的前提下拥有良好的频率特性，优于现有的方法。