首页 > 专利 > 杭州电子科技大学 > 图滤波器的滤波器组频率选择性误差交替优化设计方法专利详情

图滤波器的滤波器组频率选择性误差交替优化设计方法   0    0

有效专利 查看PDF
专利申请流程有哪些步骤?
专利申请流程图
申请
申请号:指国家知识产权局受理一件专利申请时给予该专利申请的一个标示号码。唯一性原则。
申请日:提出专利申请之日。
2021-08-18
申请公布
申请公布指发明专利申请经初步审查合格后,自申请日(或优先权日)起18个月期满时的公布或根据申请人的请求提前进行的公布。
申请公布号:专利申请过程中,在尚未取得专利授权之前,国家专利局《专利公报》公开专利时的编号。
申请公布日:申请公开的日期,即在专利公报上予以公开的日期。
2021-11-26
授权
授权指对发明专利申请经实质审查没有发现驳回理由,授予发明专利权;或对实用新型或外观设计专利申请经初步审查没有发现驳回理由,授予实用新型专利权或外观设计专利权。
2022-08-23
预估到期
发明专利权的期限为二十年,实用新型专利权期限为十年,外观设计专利权期限为十五年,均自申请日起计算。专利届满后法律终止保护。
2041-08-18
基本信息
有效性 有效专利 专利类型 发明专利
申请号 CN202110949942.4 申请日 2021-08-18
公开/公告号 CN113630104B 公开/公告日 2022-08-23
授权日 2022-08-23 预估到期日 2041-08-18
申请年 2021年 公开/公告年 2022年
缴费截止日
分类号 H03H17/00G06F30/367 主分类号 H03H17/00
是否联合申请 独立申请 文献类型号 B
独权数量 1 从权数量 2
权利要求数量 3 非专利引证数量 1
引用专利数量 5 被引证专利数量 0
非专利引证 1、2004.03.04杨圣《.非时变和两维可分时变非下采样图滤波器组的设计》《.CNKI》.2019,D. B. H. Tay, Y. Tanaka and A. Sakiyama《.Near Orthogonal OversampledGraph Filter Banks》《. IEEE SignalProcessing Letters》.2016,;
引用专利 US2010283639A、AU2020103329A、US2015317808A、WO2008151380A、US2004042557A 被引证专利
专利权维持 1 专利申请国编码 CN
专利事件 事务标签 公开、实质审查、授权
申请人信息
申请人 第一申请人
专利权人 杭州电子科技大学 当前专利权人 杭州电子科技大学
发明人 王浩、欧毅、严建阳 第一发明人 王浩
地址 浙江省杭州市下沙高教园区2号大街 邮编 310018
申请人数量 1 发明人数量 3
申请人所在省 浙江省 申请人所在市 浙江省杭州市
代理人信息
代理机构
专利代理机构是经省专利管理局审核,国家知识产权局批准设立,可以接受委托人的委托,在委托权限范围内以委托人的名义办理专利申请或其他专利事务的服务机构。
杭州君度专利代理事务所 代理人
专利代理师是代理他人进行专利申请和办理其他专利事务,取得一定资格的人。
杨舟涛
摘要
本发明公开了图滤波器的滤波器组频率选择性误差交替优化设计方法。该方法通过提升滤波器系数作为媒介,在保证滤波器组完美重构性能不被破坏的情况下,使用交替优化设计方法将滤波器组的频率选择性误差指标进行优化,以得到一组既保证重构误差值小,又将频率选择性误差减到尽可能最小的系数。相较于以往通过牺牲频率选择性误差值来获取较好重构误差值的滤波器组设计方法,本发明在已经具备极佳重构误差值的基础上,优化图滤波器组的频率选择性误差值,同时保证重构误差值小的指标能够非常好的保留,在两个指标中找到了一个极佳的平衡点,使得优化出的图滤波器组不仅在滤波性能上同时也在图像美观程度上有了进一步的改观。
  • 摘要附图
    图滤波器的滤波器组频率选择性误差交替优化设计方法
  • 说明书附图:图1
    图滤波器的滤波器组频率选择性误差交替优化设计方法
  • 说明书附图:图2
    图滤波器的滤波器组频率选择性误差交替优化设计方法
  • 说明书附图:图3
    图滤波器的滤波器组频率选择性误差交替优化设计方法
  • 说明书附图:图4
    图滤波器的滤波器组频率选择性误差交替优化设计方法
  • 说明书附图:图5
    图滤波器的滤波器组频率选择性误差交替优化设计方法
  • 说明书附图:图6
    图滤波器的滤波器组频率选择性误差交替优化设计方法
法律状态
序号 法律状态公告日 法律状态 法律状态信息
1 2022-08-23 授权
2 2021-11-26 实质审查的生效 IPC(主分类): H03H 17/00 专利申请号: 202110949942.4 申请日: 2021.08.18
3 2021-11-09 公开
权利要求
权利要求书是申请文件最核心的部分,是申请人向国家申请保护他的发明创造及划定保护范围的文件。
1.图滤波器的滤波器组频率选择性误差交替优化设计方法,其特征在于:该方法包括以下步骤:
步骤一、根据设计要求,分别确定图滤波器组中低通滤波器h0(x)的通带截止频率ωp,阻带截止频率ωs,提升滤波器的系数sn、tn的最大长度Ls、Lt;
步骤二、由提升滤波器的系数sn、tn求出低通滤波器h0(x)与低通滤波器g0(x)的表达式,进而求出此时滤波器组的频率选择性误差以及重构误差;
2‑
1.低通滤波器h0(x)与低通滤波器g0(x)的表达式
低通滤波器h0(x)与低通滤波器g0(x)的表达式分别为:
其中s(x)和t(x)分别为:
2‑
2.确定近似正交图滤波器组的完美重构条件
对于近似正交图滤波器组结构,低通滤波器h0(x)、g0(x)与高通滤波器h1(x)、g1(x)存在着严格的等式关系:
h1(x)=g0(2‑x),g1(x)=h0(2‑x)
对于图滤波器组结构,完美重构条件为:
2
h0(x)g0(x)+h1(x)g1(x)=c
令 并将低通滤波器与高通滤波器之间的等式关系代入到完美重构条件中,得到
对于近似正交图滤波器组的完美重构条件公式:
h0(x)g0(x)+h0(2‑x)g0(2‑x)=2
2‑
3.求近似正交图滤波器组的实际重构误差
定义实际条件与完美重构条件间 的差值E(x)为实际重构误差:
E(x)=|h0(x)g0(x)+h0(2‑x)g0(2‑x)‑2|
其中最大重构误差Emax为:
Emax=||E(x)||∞
2‑
4.求近似正交图滤波器组的实际频率选择性误差
在给定低通滤波器h0(x)的通带截止频率ωp,阻带截止频率ωs后,低通滤波器h0(x)的通带误差ε_h0_pass和阻带误差ε_h0_stop,g0(x)的通带误差ε_g0_pass和阻带误差ε_g0_stop分别为:
T
ε_g0_stop=||[g0(xs=ws) g0(xs+1) … g0(2)]||∞
步骤三、设计交替优化算法
3‑
1.设计优化目标函数为:
其中s、t分别为提升滤波器初始系数
3‑
2.固定s系数,优化t系数
先固定初始系数 对系数 进行优化,优化目标函数为:
其中 为经过一次优化以后新的t系数,令新的系数替代原来的系数;
3‑
3.固定t系数,优化s系数
再固定优化好的系数 对系数 进行优化,优化目标函数为:
其中 为经过一次优化以后新的s系数,令新的系数替代原来的系数;
步骤四、将最终输出的最大重构误差值以及最大频率选择性误差值作为最终优化结果,同时输出交替优化结束时得到的提升滤波器系数s、t,优化结束。

2.如权利要求1 所述图滤波器的滤波器组频率选择性误差交替优化设计方法,其特征在于:图滤波器组的理想重构条件为:在x∈[0,2]的频域上满足:当x=0时, 当x=1 时,h0(x)=1。

3.如权利要求1所述图滤波器的滤波器组频率选择性误差交替优化设计方法,其特征在于:目标函数是一个最小最大化问题,同时也是一个线性问题,使用CVX函数进行优化求解。
说明书

技术领域

[0001] 本发明属于图信号处理技术领域,具体涉及一种改进的基于图滤波器的两通道近似正交图滤波器组频率选择性误差交替优化设计方法。

背景技术

[0002] 滤波器组被广泛应用于通信、语音编码、音频编码和图像信号处理,而且M通道的均匀滤波器组的理论和设计已经得到广泛地研究。在这些应用中,经常要将信号分解成一组子带信号,各占据原信号频带的一个子频带,通常互不交叠或适当交叠,有时有需要将这些子带信号重新合成为所需要的信号,而要完成这些任务就需要数字滤波器组。前者是要设计分析滤波器组,它是单输入—多输出;后者是要设计综合滤波器组,它是多输入—单输出。不同于传统滤波器,传统滤波器组在设计时,往往频域性能要求较高,而图滤波器组在设计时对频域性能要求较低。同时在应用方面,图滤波器组往往应用于图信号的过滤。本发明研究的就是近似正交图滤波器组(Bi‑orthogonal Graph Wavelet Filter Banks)。
[0003] 文献(Jun‑Zheng Jiang and Fang Zhou and Peng‑Lang Shui.Lifting‑based design of two‑channel biorthogonal graph filter bank[J].IET Signal Processing,2016,10(6):670‑675.)提出一种设计两通道近似正交图滤波器组的新结构,在重构误差和频率选择性两个图滤波器组的指标之中找到了一个平衡点,使得设计出的图滤波器组重构误差很小的同时频率选择性也有一定的保证。但基于该结构设计出的图滤波器组仍有进一步优化的空间,使得两个指标的平衡点更为均衡,同时在不大程度影响重构误差的基础上使得画出的图滤波器组图像更为美观直接。

发明内容

[0004] 本发明针对现有技术的不足,提出了一种基于图滤波器的两通道近似正交图滤波器组频率选择性误差交替优化设计方法。该方法具体包含以下步骤:
[0005] 步骤一、根据设计要求,分别确定图滤波器组中低通滤波器h0(x)的通带截止频率ωp,阻带截止频率ωs,提升滤波器的系数sn、tn的最大长度Ls、Lt。
[0006] 步骤二、由提升滤波器的系数sn、tn求出低通滤波器h0(x)与低通滤波器g0(x)的表达式,进而求出此时滤波器组的频率选择性误差以及重构误差。
[0007] 2‑1.由提升滤波器的系数sn、tn求出低通滤波器h0(x)与低通滤波器g0(x)的表达式低通滤波器h0(x)与低通滤波器g0(x)的表达式分别为:
[0008]
[0009]
[0010] 其中s(x)和t(x)分别为:
[0011]
[0012]
[0013] 2‑2.确定近似正交图滤波器组的完美重构条件
[0014] 图滤波器组的理想重构条件为:在x∈[0,2]的频域上满足:当x=0时,当x=1时,h0(x)=1。同时,对于近似正交图滤波器组结构,低通滤波器h0(x)、g0(x)与高通滤波器h1(x)、g1(x)存在着严格的等式关系:
[0015] h1(x)=g0(2‑x),g1(x)=h0(2‑x)
[0016] 对于图滤波器组结构,完美重构条件为:
[0017] h0(x)g0(x)+h1(x)g1(x)=c2
[0018] 令 并将低通滤波器与高通滤波器之间的等式关系代入到上式,得到对于近似正交图滤波器组的完美重构条件公式:
[0019] h0(x)g0(x)+h0(2‑x)g0(2‑x)=2
[0020] 2‑3.求近似正交图滤波器组的实际重构误差
[0021] 实际设计过程中往往并不能达到完美重构条件,因此两者之间存在一个差值。定义两者之间的差值E(x)为实际重构误差:
[0022] E(x)=|h0(x)g0(x)+h0(2‑x)g0(2‑x)‑2|
[0023] 在以重构误差数值为核心比较图滤波器组性能优劣的情况下,只需要取重构误差中最大的值进行比较,即最大重构误差Emax:
[0024] Emax=||E(x)||∞
[0025] 2‑4.求近似正交图滤波器组的实际频率选择性误差
[0026] 由于设计出的滤波器通带与阻带应平坦才能达到很好的滤波效果,因此在考虑滤波器组指标时频率选择性也是一个重要衡量指标。在给定低通滤波器h0(x)通带截止频率ωp,阻带截止频率ωs后,低通滤波器h0(x)的通带误差ε_h0_pass和阻带误差ε_h0_stop,g0(x)的通带误差ε_g0_pass和阻带误差ε_g0_stop分别为:
[0027]
[0028] ε_h0_stop=‖[h0(xs=ws) h0(xs+1) … h0(2)]T‖∞
[0029]
[0030] ε_g0_stop=‖[g0(xs=ws) g0(xs+1) … g0(2)]T‖∞
[0031] 步骤三、设计交替优化算法
[0032] 3‑1.设计优化目标函数为:
[0033]
[0034] 其中s、t分别为提升滤波器初始系数 该目标函数是一个最小最大化问题,同时也是一个线性问题,可以使用CVX函数进行优化求解。
[0035] 3‑2.固定s系数,优化t系数
[0036] 先固定初始系数 对系数 进行优化,优化目标函数可以写为:
[0037]
[0038] 其中 为经过一次优化以后新的t系数,使用CVX函数进行优化后,令新的系数替代原来的系数。
[0039] 3‑3.固定t系数,优化s系数
[0040] 再固定优化好的系数 对系数 进行优化,优化目标函数可以写为:
[0041]
[0042] 其中 为经过一次优化以后新的s系数,使用CVX函数进行优化后,令新的系数替代原来的系数。
[0043] 步骤四、将最终输出的最大重构误差值以及最大频率选择性误差值作为最终优化结果,同时输出交替优化结束时得到的提升滤波器系数s、t,优化结束。
[0044] 本发明具有以下有益效果:
[0045] 相较于以往通过牺牲频率选择性误差值来获取较好重构误差值的滤波器组设计方法,本发明在已经具备极佳重构误差值的基础上,优化图滤波器组的频率选择性误差值,同时保证重构误差值小的指标能够非常好的保留,在两个指标中找到了一个极佳的平衡点,使得优化出的图滤波器组不仅在滤波性能上同时也在图像美观程度上有了进一步的改观。

实施方案

[0052] 以下结合附图对本发明作进一步说明。
[0053] 实施例1:
[0054] 设置图滤波器组中低通滤波器h0(x)的通带截止频率ωp=0.6,阻带截止频率ωs=1.4,提升滤波器的系数sn、tn的最大长度Ls=5、Lt=4,进行计算机模拟仿真。
[0055] 经过1次迭代优化后的提升滤波器的系数如表1所示,优化前后的低通滤波器频谱图像如图1、2所示,优化前后的最大重构误差值与频率选择性误差值如表2所示。
[0056]
[0057] 表1
[0058]
[0059] 表2
[0060] 实施例2:
[0061] 设置图滤波器组中低通滤波器h0(x)的通带截止频率ωp=0.6,阻带截止频率ωs=1.4,提升滤波器的系数sn、tn的最大长度Ls=7、Lt=5,进行计算机模拟仿真。
[0062] 经过1次迭代优化后的提升滤波器的系数如表3所示,优化前后的低通滤波器频谱图像如图3、4所示。优化前后的最大重构误差值与频率选择性误差值如表4所示。
[0063]
[0064] 表3
[0065]
[0066] 表4
[0067] 实施例3:
[0068] 设置图滤波器组中低通滤波器h0(x)的通带截止频率ωp=0.6,阻带截止频率ωs=1.4,提升滤波器的系数sn、tn的最大长度Ls=9、Lt=8,进行计算机模拟仿真。
[0069] 经过1次迭代优化后的提升滤波器的系数如表5所示,优化前后的低通滤波器频谱图像如图5、6所示。优化前后的最大重构误差值与频率选择性误差值如表6所示。为了贴合该滤波器组设计结构,在画图时采用 的换元方法,则此时
[0070]
[0071] 表5
[0072]
[0073] 表6
[0074] 从表2、4、6的对比结果可以看出,本方法得到的滤波器组的频率选择性误差值在优化后得到了显著的改善,同时仅使用微量的重构误差值为代价,没有破坏图滤波器组的完美重构条件。本方法与对比文件1(CN 113096052 A)中的结果对比如表7所示。
[0075]
[0076] 表7
[0077] 从表7可以看出,本方法仅使用微小重构误差为代价,在频率选择性误差上得到了极其良好的优化结果,不仅没有破坏图滤波器组的完美重构条件,也让滤波器在性能和美观程度上都得到了极大的改良。

附图说明

[0046] 图1为实施例1中迭代优化前后低通滤波器h0(x)谱图的比较
[0047] 图2为实施例1中迭代优化前后低通滤波器g0(x)谱图的比较
[0048] 图3为实施例2中迭代优化前后低通滤波器h0(x)谱图的比较
[0049] 图4为实施例2中迭代优化前后低通滤波器g0(x)谱图的比较
[0050] 图5为实施例3中迭代优化前后低通滤波器h0(x)谱图的比较
[0051] 图6为实施例3中迭代优化前后低通滤波器g0(x)谱图的比较
版权所有:盲专网 ©2023 zlpt.xyz  蜀ICP备2023003576号