基于凸优化的非下采样图滤波器组的设计方法技术

本发明专利技术公开一种基于凸优化的非下采样图滤波器组的设计方法，利用直接构造法和优化方法进行非下采样图滤波器组的设计，同时全面考虑图滤波器组的频谱特性和完全重构条件。在图信号的去噪仿真实验中，与现有技术相比，本发明专利技术设计所得的非下采样图滤波器组的去噪性能更好。

Design method of filter bank for non down sampling graph based on convex optimization

The invention discloses a design method of a non sub sampled graph filter bank based on convex optimization. It uses the direct construction method and optimization method to design the non sub sampled graph filter banks, and considers the spectrum characteristics and the complete reconstruction conditions of the graph filter banks. Compared with the existing technology, the de-noising performance of the non - down - Sampling Filter Bank of the present invention is better than the existing technology in the de-noising simulation experiment of the graph signal.

【技术实现步骤摘要】
基于凸优化的非下采样图滤波器组的设计方法
本专利技术涉及图信号处理领域中的图滤波器组
，具体涉及一种基于凸优化的非下采样图滤波器组的设计方法。
技术介绍
在网络、计算机视觉和高维云数据等领域中，图提供了一个灵活的模型来表示数据。随着图信号处理的发展，越来越多的学者从事图信号处理领域的研究工作。图信号处理将传统信号处理中的诸多概念和理论拓展至图结构上，引申出了图傅里叶变换等重要概念。图信号既可表示小规模数据信号也可表示大规模数据信号，由于大规模图具有庞大的节点数和边数，图傅里叶变换是全局变换，它不适合于大规模图信号的处理。为此，许多学者提供了图小波和图滤波器组，其具备多尺度变换特性，适合于处理大规模图信号。近年来，图小波和图滤波器组已被广泛应用于图信号的多分辨分析，压缩，去噪和分类。Narang和Ortega最早提出完全重构两通道图小波滤波器组的设计方法，该方法在二分图中，针对上下采样运算引起的频谱混叠现象，设计出具有临界采样、混叠消除、完全重构特性的正交镜像图滤波器组，但该方法设计是针对于二分图或可以分解为二分图的图信号的。此后，Narang和Ortega构造出对于任意无向图的双正交的小波滤波器组，松弛正交条件，设计得到双正交图小波，其可以对任意图信号进行处理，但此双正交图小波的设计方法未考虑滤波器的频谱选择性。AkieSakiyama和YuichiTanaka提出图滤波器组的过采样拉普拉斯矩阵，在原始图的基础上附加节点和边构造出一个包含所有原始图信息的二分图，过采样对于图信号的处理有更大的设计自由，但通过对原始图附加节点和边，导致原始图结构的改变。随后他又提出M通道过采样图滤波器组的设计算法,M通道图滤波器组具有更多的子带划分，对于大规模的图信号处理有重大的意义。上述的图小波和图滤波器组的研究主要是基于临界采样和过采样进行滤波器组的设计，其结构中均含有图信号的下采样运算，这可能导致图滤波器组中对图信号采样时会要求特殊的图结构，如二分图、可二分图或循环图。而对于一般图的结构而言，基于图染色的采样模式并不精确，基于奇异值分解的采样模式是不适用于连通图的处理，而基于最大生成树的采样模式对一般图的采样运算也缺乏准确的定义。传统数字信号处理中，传统的非下采样小波滤波器组有着广泛的应用。目前，在图信号处理中难以准确定义一般图的图信号下采样运算且非下采样图滤波器的研究相对较少，因此对于非下采样图滤波器进一步的研究有待提出。
技术实现思路
本专利技术所要解决的是目前图滤波器组的设计方法难以准确定义一般图的图信号下采样运算的问题，提供一种基于凸优化的非下采样图滤波器组的设计方法。为解决上述问题，本专利技术是通过以下技术方案实现的：基于凸优化的非下采样图滤波器组的设计方法，具体包括如下步骤：步骤1，本专利技术首先采用样条图小波作为非下采样图滤波器组的分析滤波器组。根据样条图小波的定义，非下采样图滤波器的分析滤波器组可表示为；步骤2，根据步骤1中给定的分析滤波器组，利用非下采样图滤波器组的完全重构条件，直接构造出综合滤波器组：G0＝IN,G1＝IN步骤3，根据步骤1中给定的分析滤波器组，以综合滤波器组的阻带能量为目标函数，以完全重构条件为约束函数，约束优化问题是半正定规划问题，可以求解出综合滤波器。优化问题为式中，Es(g0)为综合低通滤波器的阻带能量；Es(g1)为综合高通滤波器的阻带能量；α,β为综合权值；h0,h1分别为分析低通和高通滤波器；g0,g1分别为综合低通和高通滤波器；cLh0(·),cLh1(·)为分析滤波器h0,h1的频率向量，cLg0(·),cLg1(·)为综合滤波器g0,g1的频率向量；Lh0,Lh1为分析滤波器h0,h1的长度，Lg0,Lg1为综合滤波器g0,g1的长度，N+1为频率离散点数；εr为重构误差容限；T代表转置；步骤4，根据前三步中直接构造法和优化设计法所得的两通道非下采样图滤波器组，对其低频分量采用级联的方式，从两通道非下采样图滤波器组扩展到多通道非下采样图滤波器组。上述步骤3中，采用半正定规划求解器SDPT3有效求解综合滤波器。上述步骤3中，权值α＝1与权值β＝0.1。上述步骤3中，λk的取值范围为[0,2]。与现有技术相比，本专利技术利用直接构造法和优化方法进行非下采样图滤波器组的设计，同时全面考虑图滤波器组的频谱特性和完全重构条件。在图信号的去噪仿真实验中，与现有技术相比，本专利技术设计所得的非下采样图滤波器组的去噪性能更好。附图说明图1为两通道非下采样图滤波器组的基本结构。图2为多通道非下采样图滤波器组的基本结构。图3为图滤波器的幅度响应。其中(a)为本专利技术中直接构造所得到图滤波器的幅度响应；(b)为本专利技术中优化设计所得到图滤波器的幅度响应。图4为本专利技术中所得到的非下采样图滤波器系统与已有方法得到的图滤波器系统进行去噪实验的仿真结果。其中(a)为输入信号；(b)为噪声信号(σ＝1/2)；(c)为基于现有方法1(临界采样双正交图滤波器组)设计的图滤波器组去噪后的输出信号；(d)为基于现有方法2(M通道过采样图滤波器组)设计的图滤波器组去噪后的输出信号；(e)为本专利技术直接构造方法设计所得的两通道非下采样图滤波器组去噪后的输出信号；(f)为本专利技术直接构造方法设计所得的三通道非下采样图滤波器组去噪后的输出信号；(g)为本专利技术优化设计方法所得的两通道非下采样图滤波器组去噪后的输出信号；(h)为本专利技术优化设计方法所得的三通道非下采样图滤波器组去噪后的输出信号。具体实施方式为了便于理解，下面以两通道和三通道非下采样图滤波器组为例，对本专利技术基于凸优化的非下采样图滤波器组的设计方法，进行进一步详细说明。根据图1中表示的两通道非下采样图滤波器组，其中两个线性滤波器H0,H1组成分析滤波器组，综合滤波器组由两个线性滤波器G0,G1组成。图信号f为图滤波器组的输入信号，输出信号为，f0,f1分别表示低频子带系数和高频子带系数。分析图滤波器组和综合图滤波器组可分别表示为σ(G)是由图G的拉普拉斯矩阵所有特征值构成的特征空间,Pλ,λ∈σ(G)是特征空间的投影矩阵，hi(λ),gi(λ)分别是分析子带滤波器和综合子带滤波器。两通道非下采样图滤波器组的输入输出关系可表示为T是图滤波器组的传递函数。两通道非下采样图滤波器组的完全重构条件可写为G0H0+G1H1＝IN(3)由式(1)、(3)得h0(λ)g0(λ)+h1(λ)g1(λ)＝1(4)多项式中λ∈[0,2]。第一步：采用样条图小波作为分析滤波器组。根据样条图小波的定义，任意图的样条图小波分析滤波器组可表示为HLP,HHP分别为任意图的样条图小波滤波器组的低通和高通滤波器。当n＝1时，定义非下采样图滤波器组的分析滤波器组为Λ是归一化拉普拉斯矩阵的正实数特征值λ∈[0,2]组成的对角矩阵，U是特征值λ对应的特征向量组成的矩阵，上标H表示共轭转置。根据公式(1)、(7)和(8)，可得非下采样图滤波器组的分析子带滤波器第二步：采用直接构造的方法，在图滤波器组的设计中，完全重构特性是首要考虑的。和传统滤波器组一样，重构误差是衡量图滤波器组整体性能的主要性能指标，当图滤波器组是无失真传输，重构误差为零时，则滤波器组是完全重构的，即设计的综合滤波器组能准确的重构出输入信本文档来自技高网...

【技术保护点】
基于凸优化的非下采样图滤波器组的设计方法，其特征是，包括步骤如下：步骤1，采用样条图小波作为非下采样图滤波器组的分析滤波器组，并根据样条图小波的定义，给定两通道下的分析滤波器组；步骤2，根据步骤1中给定的分析滤波器组，利用非下采样图滤波器组的完全重构条件，直接构造出两通道下的综合滤波器组；步骤3，基于给定的分析滤波器组，以综合滤波器组的阻带能量为目标函数，将完全重构条件为约束函数，构建优化问题来求解综合滤波器；其中优化问题为：

【技术特征摘要】
1.基于凸优化的非下采样图滤波器组的设计方法，其特征是，包括步骤如下：步骤1，采用样条图小波作为非下采样图滤波器组的分析滤波器组，并根据样条图小波的定义，给定两通道下的分析滤波器组；步骤2，根据步骤1中给定的分析滤波器组，利用非下采样图滤波器组的完全重构条件，直接构造出两通道下的综合滤波器组；步骤3，基于给定的分析滤波器组，以综合滤波器组的阻带能量为目标函数，将完全重构条件为约束函数，构建优化问题来求解综合滤波器；其中优化问题为：式中，Es(g0)为综合低...

【专利技术属性】
技术研发人员：蒋俊正，杨圣，欧阳缮，孙希延，纪元法，刘松辽，杨玉琳，曹想，赵海兵，杨杰，
申请(专利权)人：桂林电子科技大学，
类型：发明
国别省市：广西,45