CN114331926A - 基于换元思想的两通道图滤波器组系数设计优化方法 - Google Patents
基于换元思想的两通道图滤波器组系数设计优化方法 Download PDFInfo
- Publication number
- CN114331926A CN114331926A CN202111638212.9A CN202111638212A CN114331926A CN 114331926 A CN114331926 A CN 114331926A CN 202111638212 A CN202111638212 A CN 202111638212A CN 114331926 A CN114331926 A CN 114331926A
- Authority
- CN
- China
- Prior art keywords
- filter bank
- filter
- omega
- low
- pass
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Complex Calculations (AREA)
Abstract
本发明公开了基于换元思想的两通道图滤波器组系数优化设计方法。通过两个不同值域的变量作为媒介,在保证滤波器组完美重构性能不被破坏的情况下,使用换元方法将高阶数滤波器组的频率选择性误差指标进行优化,以得到在高阶数时双正交滤波器组重构误差与频率选择性误差都小的系数组。相较于现有的双正交图滤波器组设计方法只能设计低阶数的不足,本发明在高阶数滤波器组上实现创新,提出一种高阶低阶均可使用的设计方法,使得在高阶数时图滤波器组依然能够得到较好的设计效果。
Description
技术领域
本发明属于图信号处理技术领域,具体涉及一种改进的基于换元思想的图滤波器组系数设计优化方法。
背景技术
滤波器组被广泛应用于通信、语音编码、音频编码和图像信号处理,M通道的均匀滤波器组的理论和设计已经得到广泛地研究。在这些应用中,经常要将信号分解成一组子带信号,每个子带信号各占据原信号频带的一个子频带,这些子带信号之间的频带通常互不交叠或适当交叠,有时有需要将这些子带信号重新组合成为一个信号,而要完成上述过程就需要数字滤波器组。前者是要设计分析滤波器组,它是单输入—多输出;后者是要设计综合滤波器组,它是多输入—单输出。不同于传统滤波器在设计时,往往频域性能要求较高,图滤波器组在设计时对频域性能要求较低。在应用方面,图滤波器组往往应用于图信号的过滤。
Jun-Zheng Jiang,Fang Zhou,Peng-Lang Shui提出了一种设计两通道近似正交图滤波器组的新结构,在重构误差和频率选择性两个图滤波器组的指标之中找到了一个平衡点,使得设计出的图滤波器组重构误差很小的同时频率选择性也有一定的保证。但该方法仅在低阶数图滤波器设计时有效,在高阶数图滤波器设计时频谱图像效果十分差。因此需要针对高阶数的图滤波器提出一种新的设计、优化方法。
发明内容
本发明针对现有技术的不足,提出了一种基于换元思想的两通道图滤波器组系数优化设计方法,
基于换元思想的两通道图滤波器组系数设计优化方法,具体包含以下几个步骤:
步骤一、根据设计要求,分别确定图滤波器组中低通滤波器h0(λ)的通带截止频率λp,阻带截止频率λs,滤波器最大长度Lh与Lg。
步骤二、确定换元前变量λ和换元后变量ω之间的关系。
2-1.λ为图滤波器中常用的自变量,定义域为λ∈[0,2],在λ域下,图滤波器常用的表达式可以写为:
其中a(i)为低通滤波器h0(λ)的系数;而ω为传统滤波器中常用的自变量,定义域为ω∈[0,π],自变量λ与自变量ω的关系为:
x=λ-1,x∈[-1,1]
x=cosθ,θ∈[π,0]
ω=π-θ,ω∈[0,π]
联立上述三式,可以得到:
λ=cos(π-ω)+1,ω∈[0,π]
将该公式展开,可得:
λ=1-cosω,ω∈[0,π]
即可以构建出在ω域下的图滤波器的表达式为:
其中b(i)为低通滤波器h0(ω)的系数,同时由该关系式可以根据λ域下的通带截止频率λp和阻带截止频率λs推导出ω域下的通带截止频率ωp和阻带截止频率ωs值。
2-2.通过传统滤波器中Kaiser窗生成低通滤波器的方法,生成一个通带截止频率为ωp,阻带截止频率为ωs,滤波器最大长度为Lh的低通滤波器h0(ω),同时令低通滤波器g0(ω)的初始值与h0(ω)相等。
2-3.确定在ω域下的完美重构条件,图滤波器组的理想重构条件在为x∈[0,2]的频域上应当满足:当λ=0时,当λ=1时,h0(λ)=1。同时,对于近似正交图滤波器组结构,低通滤波器h0(λ),g0(λ)与高通滤波器h1(λ),g1(λ)存在着严格的等式关系:
h1(λ)=g0(2-λ),g1(λ)=h0(2-λ)
对于任意图滤波器组结构,完美重构条件可以用如下公式表示:
h0(λ)g0(λ)+h1(λ)g1(λ)=c2
h0(λ)g0(λ)+h0(2-λ)g0(2-λ)=2
因此在ω域下,近似正交图滤波器组的完美重构条件公式可以改写为:
h0(ω)g0(ω)+h0(π-ω)g0(π-ω)=2
而实际设计过程中往往并不能达到完美重构条件,因此两者之间存在一个差值。定义实际重构误差E(ω)为:
E(ω)=|h0(ω)g0(ω)+h0(π-ω)g0(π-ω)-2|
最大重构误差定义为:
Emax=|E(ω)|inf
2-4.由于设计出的滤波器通带与阻带应平坦才能达到较好的滤波效果,因此在考虑滤波器组指标时频率选择性也是一个重要衡量指标。确定在ω域下的频率选择性误差关系式,低通滤波器h0(λ)的通带误差ε_h0_pass和阻带误差ε_h0_stop,g0(λ)的通带误差ε_g0_pass和阻带误差ε_g0_stop可为:
ε_h0_stop=||[h0(λs) h0(λs+1) … h0(2)]T||∞
ε_g0_stop=||[g0(λs) g0(λs+1) … g0(2)]T||∞
在ω域下为:
ε_h0_stop=||[h0(ωs) h0(ωs+1) … h0(π)]T||∞
ε_g0_stop=||[g0(ωs) g0(ωs+1) … g0(π)]T||∞
2-5.将实际重构误差E(ω)分别对h0(ω)和g0(ω)的系数hs与gs进行求导,得到如下两个梯度表达式:
步骤三、设计信赖域迭代优化算法
E(ω)0=|h0(hs;ω)g0(gs;ω)+h0(hs;π-ω)g0(gs;π-ω)-2|
则第1次迭代时的优化目标函数为:
同时,对于修改变量与应满足::由于该优化目标函数是一个最小最大化问题,同时也是一个线性问题,可以使用CVX函数进行优化求解。求解得到第1次迭代时修改变量与的具体值,再对h0(ω)和g0(ω)的系数hs与gs进行重新赋值操作:
优化目标函数为:
步骤四、设定迭代次数K,重复K次步骤三的算法优化求解过程,完成两通道图滤波器组系数设计优化。
本发明具有的优点是:
现有的双正交图滤波器组设计方法在低阶数滤波器时有很好的效果,但是来到高阶数时便会破坏滤波器本身的频谱特性,造成频率选择性误差非常大。本申请能够设计高阶数的双正交图滤波器组,在双正交图滤波器组的完美重构条件不被破坏的同时,频率选择性指标也能兼顾。
附图说明
图1为实施例1中通过本方法优化设计后得到的滤波器h0(ω)和h1(ω)图像。
图2为实施例1中通过现有技术优化设计后得到的滤波器h0(λ)和h1(λ)图像。
图3为实施例2中通过本方法优化设计后得到的滤波器h0(ω)和h1(ω)图像。
图4为实施例2中通过现有技术优化设计后得到的滤波器h0(λ)和h1(λ)图像。
具体实施方式
以下结合附图对本发明作进一步说明。
实施例1
本实施例中,设置低通滤波器h0(ω)的通带截止频率λp=0.6,阻带截止频率λs=1.4;最大长度Lh=50,Lg=50;迭代次数K=500。
通过本方法优化后,得到的重构误差值和频率选择性误差指标如表1所示:
表1
得到的低通滤波器h0(ω)和g0(ω)如图1所示,其中h1(ω)=g0(π-ω)。
通过现有技术[1]优化后得到的低通滤波器h0(ω)和g0(ω)如图2所示,与本方法得到重构误差值和频率选择性误差比较结果如表2所示:
表2
实施例2
本实施例中,设置低通滤波器h0(ω)的通带截止频率λp=0.6,阻带截止频率λs=1.4;最大长度Lh=50,Lg=60;迭代次数K=500。
通过本方法优化后,得到的重构误差值和频率选择性误差指标如表3所示:
表3
得到的低通滤波器h0(ω)和g0(ω)如图3所示。
通过现有技术优化后得到的低通滤波器h0(ω)和g0(ω)如图4所示,与本方法得到重构误差值和频率选择性误差比较结果如表4所示:
表4
从比较结果可以看出,本方法在设计高阶数滤波器的效果明显好于现有,能够解决现有的双正交图滤波器组设计方法不能设计高阶数滤波器的问题,做到保证图滤波器组的完美重构条件满足的情况下,同时兼顾频率选择性,整体提升双正交图滤波器组的性能和频谱图的美观程度。
[1]《Lifting-based design of two-channel biorthogonal graph filterbank》,Jun-Zheng Jiang, Fang Zhou,Peng-Lang Shui。
Claims (4)
1.基于换元思想的两通道图滤波器组系数设计优化方法,对由低通滤波器h0(x)、g0(x)与高通滤波器h1(x)、g1(x)组成的近似正交图滤波器组进行优化,在λ域下,低通滤波器h0(λ)的通带截止频率为λp,阻带截止频率为λs;其特征在于:具体包括以下步骤:
步骤一、对于λ∈[0,2],在λ域下图滤波器h(λ)的表达式为:
其中a(i)为低通滤波器h0(λ)的系数;
定义ω∈[0,π],将与图滤波器h(λ)相关的自变量λ写为:
即:
λ=1-cosω,ω∈[0,π]
因此在ω域下图滤波器h(ω)的表达式为:
其中b(i)为低通滤波器h0(ω)的系数;
步骤二、根据λ域下低通滤波器h0(λ)的通带截止频率λp和阻带截止频率λs,使用步骤一中图滤波器h(ω)的表达式计算得到ω域下低通滤波器h0(ω)的通带截止频率ωp和阻带截止频率ωs;低通滤波器h0(ω)的通带误差ε_h0_pass和阻带误差ε_h0_stop,g0(ω)的通带误差ε_g0_pass和阻带误差ε_g0_stop分别为:
ε_h0_stop=||[h0(ωs) h0(ωs+1)…h0(π)]T||∞
ε_g0_stop=||[g0(ωs) g0(ωs+1)…g0(π)]T||∞
ω域下的近似正交图滤波器组的完美重构条件为:
h0(ω)g0(ω)+h0(π-ω)g0(π-ω)=2
实际重构误差E(ω)为:
E(ω)=|h0(ω)g0(ω)+h0(π-ω)g0(π-ω)-2|
分别对两个低通滤波器的系数hs、gs进行求导,得到两个梯度表达式为:
步骤三、设计信赖域迭代优化算法,对两个低通滤波器的系数hs、gs进行迭代优化;
s3.2、当E(ω)i<E(ω)i-1时,修改两个低通滤波器的系数为h′s、g′s,通过s3.1进行第i+1次迭代:
步骤四、重复步骤三K次,完成两通道图滤波器组系数设计优化。
2.如权利要求1所述基于换元思想的两通道图滤波器组系数设计优化方法,其特征在于:近似正交图滤波器组结构的低通滤波器h0(λ),g0(λ)与高通滤波器h1(λ),g1(λ)的等式关系为:
h1(λ)=g0(2-λ),g1(λ)=h0(2-λ)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111638212.9A CN114331926B (zh) | 2021-12-29 | 2021-12-29 | 基于换元思想的两通道图滤波器组系数设计优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111638212.9A CN114331926B (zh) | 2021-12-29 | 2021-12-29 | 基于换元思想的两通道图滤波器组系数设计优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114331926A true CN114331926A (zh) | 2022-04-12 |
CN114331926B CN114331926B (zh) | 2022-06-10 |
Family
ID=81016834
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111638212.9A Active CN114331926B (zh) | 2021-12-29 | 2021-12-29 | 基于换元思想的两通道图滤波器组系数设计优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114331926B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5644648A (en) * | 1991-12-23 | 1997-07-01 | Lucent Technologies Inc. | Method and apparatus for connected and degraded text recognition |
CN1794300A (zh) * | 2005-12-20 | 2006-06-28 | 上海电力学院 | 基于最优小波滤波器组的多传感器图像融合方法 |
CN104506164A (zh) * | 2014-12-29 | 2015-04-08 | 桂林电子科技大学 | 基于两步法的图滤波器组优化设计方法 |
CN105243241A (zh) * | 2015-11-12 | 2016-01-13 | 桂林电子科技大学 | 基于提升结构的两通道双正交图滤波器组设计方法 |
CN107239623A (zh) * | 2017-06-08 | 2017-10-10 | 桂林电子科技大学 | 基于凸优化的m通道过采样图滤波器组的优化设计方法 |
CN113096052A (zh) * | 2021-03-02 | 2021-07-09 | 杭州电子科技大学 | 基于交替迭代的两通道近似正交图滤波器组系数优化方法 |
-
2021
- 2021-12-29 CN CN202111638212.9A patent/CN114331926B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5644648A (en) * | 1991-12-23 | 1997-07-01 | Lucent Technologies Inc. | Method and apparatus for connected and degraded text recognition |
CN1794300A (zh) * | 2005-12-20 | 2006-06-28 | 上海电力学院 | 基于最优小波滤波器组的多传感器图像融合方法 |
CN104506164A (zh) * | 2014-12-29 | 2015-04-08 | 桂林电子科技大学 | 基于两步法的图滤波器组优化设计方法 |
CN105243241A (zh) * | 2015-11-12 | 2016-01-13 | 桂林电子科技大学 | 基于提升结构的两通道双正交图滤波器组设计方法 |
CN107239623A (zh) * | 2017-06-08 | 2017-10-10 | 桂林电子科技大学 | 基于凸优化的m通道过采样图滤波器组的优化设计方法 |
CN113096052A (zh) * | 2021-03-02 | 2021-07-09 | 杭州电子科技大学 | 基于交替迭代的两通道近似正交图滤波器组系数优化方法 |
Non-Patent Citations (1)
Title |
---|
CHIEN-CHENG TSENG,ET AL: ""Designs of Bipartite Graph Filter Bank Using Graph"", 《2018 INTERNATIONAL SYMPOSIUM ON COMPUTER, CONSUMER AND CONTROL (IS3C)》 * |
Also Published As
Publication number | Publication date |
---|---|
CN114331926B (zh) | 2022-06-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR100565821B1 (ko) | 복소 지수 변조 필터 뱅크를 사용한 엘리어싱 저감 | |
Evangelista et al. | Discrete frequency warped wavelets: Theory and applications | |
KR101812003B1 (ko) | 고 주파수 재구성을 위한 복소 지수 변조 필터 뱅크 | |
CN107256537B (zh) | 一种设计两通道正交图滤波器组的设计方法 | |
CN107294512B (zh) | 一种基于树型结构的非均匀滤波器组滤波方法 | |
CN105243241B (zh) | 基于提升结构的两通道双正交图滤波器组设计方法 | |
CN103384901B (zh) | 改进变换域中的滤波 | |
CN107239623B (zh) | 基于凸优化的m通道过采样图滤波器组的优化设计方法 | |
CN113096052A (zh) | 基于交替迭代的两通道近似正交图滤波器组系数优化方法 | |
CN114331926B (zh) | 基于换元思想的两通道图滤波器组系数设计优化方法 | |
Furtado et al. | Optimized prototype filter based on the FRM approach for cosine-modulated filter banks | |
CN105719255A (zh) | 基于拉氏结构的可变带宽线性相位滤波器方法 | |
CN113630104B (zh) | 图滤波器的滤波器组频率选择性误差交替优化设计方法 | |
CN111010144B (zh) | 改进的两通道iir的qmfb设计方法 | |
CN111289106A (zh) | 一种基于数字滤波的光谱降噪方法 | |
Matei | Analytic design of uniform circular filter banks | |
Cecchi et al. | A new approach to digital audio equalization | |
Kumar et al. | An efficient closed-form design method for nearly perfect reconstruction of non-uniform filter bank | |
CN113014225B (zh) | 基于全通滤波器混合迭代技术的qmf组设计方法 | |
Virgulti et al. | IIR filter approximation of an innovative digital audio equalizer | |
Sumanth et al. | Performance Analysis of Non-uniform Filter Bank with Equi ripple Band pass Filter | |
Anurag et al. | An efficient iterative method for nearly perfect reconstruction non-uniform filter bank | |
CN115221823A (zh) | 一种两通道正交镜像滤波器组的系数稀疏优化设计方法 | |
CN111010146A (zh) | 一种基于快速滤波器组的信号重构结构及其设计方法 | |
Dutta Roy et al. | Design of maximally flat and monotonic FIR filters using the Bernstein polynomial |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |