首页 > 专利 > 杭州电子科技大学 > 基于全通滤波器混合迭代技术的QMF组设计方法专利详情

基于全通滤波器混合迭代技术的QMF组设计方法   0    0

有效专利 查看PDF
专利申请流程有哪些步骤?
专利申请流程图
申请
申请号:指国家知识产权局受理一件专利申请时给予该专利申请的一个标示号码。唯一性原则。
申请日:提出专利申请之日。
2021-02-03
申请公布
申请公布指发明专利申请经初步审查合格后,自申请日(或优先权日)起18个月期满时的公布或根据申请人的请求提前进行的公布。
申请公布号:专利申请过程中,在尚未取得专利授权之前,国家专利局《专利公报》公开专利时的编号。
申请公布日:申请公开的日期,即在专利公报上予以公开的日期。
2021-07-09
授权
授权指对发明专利申请经实质审查没有发现驳回理由,授予发明专利权;或对实用新型或外观设计专利申请经初步审查没有发现驳回理由,授予实用新型专利权或外观设计专利权。
2021-10-08
预估到期
发明专利权的期限为二十年,实用新型专利权期限为十年,外观设计专利权期限为十五年,均自申请日起计算。专利届满后法律终止保护。
2041-02-03
基本信息
有效性 有效专利 专利类型 发明专利
申请号 CN202110148894.9 申请日 2021-02-03
公开/公告号 CN113014225B 公开/公告日 2021-10-08
授权日 2021-10-08 预估到期日 2041-02-03
申请年 2021年 公开/公告年 2021年
缴费截止日
分类号 H03H11/12H03H17/02 主分类号 H03H11/12
是否联合申请 独立申请 文献类型号 B
独权数量 1 从权数量 3
权利要求数量 4 非专利引证数量 1
引用专利数量 2 被引证专利数量 0
非专利引证 1、2008.02.07Hao Wang等.“Iterative technique foroptimizing two-channel: Quadraturemirror, bi-orthogonal and graph filterbank”《.Signal Processing》.2019,1-10. J.-H. Lee.“Design of two-channellinear-phase QMF banks based on real IIRall-pass filters”《.IEE Proc.-Vis. ImageSignal Process》.2003,第150卷(第5期),331-338. 李伟琪等“.外插脉冲响应 FIR 滤波器系数敏感度参数研究”《.通信技术》.2018,第51卷(第11期),2558-2564. 高翠翠等“.基于正交镜像滤波器组的子带脉冲压缩方法”《. 微型机与应用 》.2015,第34卷(第19期),27-29,42.;
引用专利 US2008033730A、US2008033730A 被引证专利
专利权维持 1 专利申请国编码 CN
专利事件 事务标签 公开、实质审查、授权
申请人信息
申请人 第一申请人
专利权人 杭州电子科技大学 当前专利权人 杭州电子科技大学
发明人 王浩、李祥振、张福洪、欧毅、严建阳、赵晨子、李伟琪、赵知劲 第一发明人 王浩
地址 浙江省杭州市下沙高教园区2号大街 邮编 310018
申请人数量 1 发明人数量 8
申请人所在省 浙江省 申请人所在市 浙江省杭州市
代理人信息
代理机构
专利代理机构是经省专利管理局审核,国家知识产权局批准设立,可以接受委托人的委托,在委托权限范围内以委托人的名义办理专利申请或其他专利事务的服务机构。
杭州君度专利代理事务所 代理人
专利代理师是代理他人进行专利申请和办理其他专利事务,取得一定资格的人。
杨舟涛
摘要
本发明公开了基于全通滤波器混合迭代技术的QMF组设计方法。该方法包括根据设计目标确定滤波器的个数、阶数、截止频率和每个滤波器的理想相位。然后计算滤波器的实际相位与实际相位差,得到整体失真函数的实际相位差。通过LSK处理和交替迭代的方法分别完成目标函数的分母与分子的线性化,再通过CVX工具包得到滤波器系数。最后判断结果是否满足设计要求,不满足时修正加权值后重新迭代。本方法通过数学处理,降低了设计过程中公式的复杂度与求解难度,结合了不同手段的非线性处理方法将滤波器组的整体相位误差与阻带幅度误差控制在要求范围内。
  • 摘要附图
    基于全通滤波器混合迭代技术的QMF组设计方法
  • 说明书附图:图1
    基于全通滤波器混合迭代技术的QMF组设计方法
  • 说明书附图:图2
    基于全通滤波器混合迭代技术的QMF组设计方法
  • 说明书附图:图3
    基于全通滤波器混合迭代技术的QMF组设计方法
  • 说明书附图:图4
    基于全通滤波器混合迭代技术的QMF组设计方法
  • 说明书附图:图5
    基于全通滤波器混合迭代技术的QMF组设计方法
法律状态
序号 法律状态公告日 法律状态 法律状态信息
1 2021-10-08 授权
2 2021-07-09 实质审查的生效 IPC(主分类): H03H 11/12 专利申请号: 202110148894.9 申请日: 2021.02.03
3 2021-06-22 公开
权利要求
权利要求书是申请文件最核心的部分,是申请人向国家申请保护他的发明创造及划定保护范围的文件。
1.基于全通滤波器混合迭代技术的QMF组设计方法,其特征在于;该方法具体包括以下步骤:
步骤一、确定基本参数
设置全频带上的频率点数为L,正交镜像滤波器组中的两个全通滤波器的阶数分别为N1、N2,其中N1=N2+1,设置三个初始加权值W1=1、W2=1、W3=W1+W2+1;令迭代初始系数k=0,全通滤波器的初始迭代系数ai(k)=0,i=1、2;全通滤波器组成的低通滤波器H0(z)的通带截止频率ωp、阻带截止频率ωs,组成的高通滤波器H1(z)的通带截止频率ωs、阻带截止频率ωp;
步骤二、确定实际相位差
s2.1、确定滤波器的理想相位
设置两个全通滤波器的理想相位分别为θd1(ω)=‑N1ω'+0.25ω'和θd2(ω)=‑N2ω'‑
0.25ω',其中ω'∈[0,π];
s2.2、求取滤波器的实际相位误差

全通滤波器的表达式为Ai(e ):

Ni分别表示两个全通滤波器的阶数N1、N2;ai(n)表示全通滤波器Ai(e )系数ai的第n个元素,n=0,2,…,Ni;
全通滤波器的实际相位为θi(ω):
全通滤波器的实际相位差为θei(ω):
θei(ω)=θi(ω)‑θdi(ω)
其中ω∈[0,ωp]
对公式(4)去三角化后:
其中ψi(ω)=0.5×(Niω+θdi(ω));
s2.3、计算整体失真传输函数的实际相位误差
滤波器组的整体失真传输函数表达式为 理想相位、实
际相位以及相位误差分别为:
θdT(ω)=θd1(2ω)+θd2(2ω)‑ω    (6)
θT(ω)=θ1(2ω)+θ2(2ω)‑ω    (7)
θeT(ω)=θT(ω)‑θdT(ω)    (8)
因此
对公式(9)去三角化后:
步骤三、求解全通滤波器系数ai(k)
s3.1、得到目标优化函数表式如式(11)所示
jω jω jω
其中a1、a2分别表示全通滤波器A1(e )、A2(e )的系数;a表示全通滤波器A1(e )与A2jω (k) jω
(e )的系数,a=[a1;a2];f1 表示第k次迭代时全通滤波器A1(e )的相位误差, 表示第jω jω
k迭代时全通滤波器A2(e )的相位误差, 表示第k次迭代时整体失真传输函数T(e )的jω
相位误差,M(a,ω)表示低通滤波器H0(e )的阻带幅度误差:
s3.2、对步骤3.1中的目标函数进行LSK处理,将其转换为线性问题
完成目标函数(11)中的分母部分的线性转化;
s3.3对目标函数分子部分由于两个全通滤波器系数存在互乘而导致的非线性问题采用交替迭代的方法进行处理:
所述交替迭代为:在第k次迭代过程中,将第k‑1次迭代得到的两个全通滤波器的系数代入公式(15)中,得到常数化后的分母;然后将常数化的分母的值与第k‑1次得到的全通滤波器A1的系数 的值代入公式(16)中,求解全通滤波器A2的系数a2;通过凸优化处理工具箱CVX得到一组滤波器系数,再将这组滤波器系数代入公式(15)中更新分母的值,并将更新后的分母与全通滤波器A2的系数代入公式(17)中,求解全通滤波器A1的系数,并通过凸优化处理工具箱CVX得到一组滤波器系数;
反复交替,直至迭代终止得到两个全通滤波器的系数a;
步骤四、目标判断
设置一个实数μ作为判断目标,当目标函数满足 时结束迭代,将步骤三
得到的滤波器系数作为最终得到的全通滤波器的系数;当目标函数不满足上述条件时,进(k) (k‑1)
入步骤五;其中E 是第k次迭代中目标函数的最大值,E 是第k‑1次迭代中目标函数的最大值;
步骤五、修正Wi
计算第k次迭代的相位误差 其中 表示第k次迭代的全通滤

波器Ai(e )的实际相位;然后计算相位误差 的包络
设置一个实数判断目标ε,如果不满足 使
然后令 返回步骤三;否则的话直接令
返回步骤三。

2.如权利要求1所述基于全通滤波器混合迭代技术的QMF组设计方法,其特征在于:所述正交镜像滤波器组由分析滤波器组和综合滤波器组构成;分析滤波器组由两个全通滤波jω jω
器A1(e )、A2(e )组成的低通滤波器H0(z)与高通滤波器H1(z)组成,综合滤波器组为分析滤波器组的复数共轭形式。

3.如权利要求1所述基于全通滤波器混合迭代技术的QMF组设计方法,其特征在于:全通滤波器的理想相位目标为保证IIR滤波器稳定。

4.如权利要求1或2所述基于全通滤波器混合迭代技术的QMF组设计方法,其特征在于:
滤波器组的整体失真传输函数的相位误差是分析滤波器组和综合滤波器组的综合误差,体现整个系统的相位误差情况。
说明书

技术领域

[0001] 本发明属于数字信号处理技术领域,具体涉及一种基于全通滤波器交替混合迭代技术的正交镜像滤波器组的设计方法。

背景技术

[0002] 数字信号处理系统经过近几十年的发展,已经在生活中的各种设备上代替了模拟系统,比如说移动通信系统、消费者电子产品和汽车电子产品,或者助听器。而滤波器作为数字系统里面不可或缺的基础性组件,其功能直接或间接性的影响到整个数字系统的性能。
[0003] 滤波器组的作用在最近二十年愈发显得重要,主要是因为数字技术的发展,单一通道对信号的处理已经不能满足人们的需求,滤波器组的出现可以帮助人们将一个信号划分为多个子带,对不同频率的子带作不同方式的处理,从而达到自己想要的结果,比如说,可以在频谱分析中作为跨多路复用器,或是处理子带的频域、时域信号。而处理子带信号需要通过一个分析‑合成滤波器组来获得重建的时域信号。
[0004] 近年来两通道QMFB(Quadrature Mirror Filter Bank,正交镜像滤波器组)被应用于越来越多的领域,例如语音和图像信号的子带编码,小波基设计等。由于QMFB的广泛应用,人们对于它的设计关注也越来越重视。
[0005] X.Zhang和H.Iwakura提出了一种利用数字全通滤波器(DAFs,Digital Allpass Filters)设计QMFB的方法。由全通滤波器和纯延迟部分组成的QMFB具有近似线性的相位响应,使用另一个附加的DAFs作为均衡器可以使相位畸变最小化。论文中采用一个DAFs和一个纯延迟共同设计一个QMFB,设计DAFs的方法是其于1992年提出的一个算法,可以简单描述为一种基于特征值问题设计等纹波相位响应的DAFs的方法,通过计算最大特征向量和迭代法得到了最优滤波系数。这种算法借用等波纹设计FIR的方法得到了全通滤波器,但是其计算极其复杂,并且该算法并没有考虑全通滤波器和纯延迟合成时的相位误差和得到的低通滤波器的阻带误差。
[0006] S.S.Lawson和A.Klouche‑Djedid提出了另外一种利用DAFs设计QMFB的方法,该方法具体表现为把全通滤波器的相位分解为分子和分母两部分,并只对分子部分进行优化,所采用的优化方法是根据正弦函数和余弦函数的正交性质,进行近似积分,从而得到全通滤波器的系数。

发明内容

[0007] 针对现有技术的不足,本发明提出了基于全通滤波器混合迭代技术的QMF组设计方法,基于全通滤波器混合迭代技术的QMF组设计方法,结合了两种不同的优化方法,对分子与分母部分都进行优化,降低了计算复杂程度与计算量,提高了滤波器性能。
[0008] 基于全通滤波器混合迭代技术的QMF组设计方法,具体包括以下步骤:
[0009] 步骤一、确定基本参数
[0010] 正交镜像滤波器组由分析滤波器组和综合滤波器组构成;分析滤波器组由低通滤波器H0(z)与高通滤波器H1(z)组成,低通滤波器H0(z)与高通滤波器H1(z)由两个全通滤波jω jω器A1(e )、A2(e )组成的,低通滤波器H0(z)的通带截止频率ωp、阻带截止频率ωs,高通滤波器H1(z)的通带截止频率ωs、阻带截止频率ωp;综合滤波器组为分析滤波器组的复数共轭形式。设置全频带上的频率点数为L,两个全通滤波器的阶数分别为N1、N2,其中N1=N2+1,三个初始加权值W1=1、W2=1、W3=W1+W2+1;令迭代初始系数k=0,全通滤波器的初始迭代系数ai(k)=0,i=1、2。
[0011] 步骤二、确定实际相位差
[0012] s2.1、确定滤波器的理想相位
[0013] 以IIR滤波器稳定为目标,设置两个全通滤波器的理想相位分别为θd1(ω)=‑N1ω'+0.25ω'和θd2(ω)=‑N2ω'‑0.25ω',其中ω'∈[0,π]。
[0014] s2.2、求取滤波器的实际相位误差
[0015] 全通滤波器的表达式为Ai(ejω):
[0016]
[0017] Ni分别表示两个全通滤波器的阶数N1、N2。ai(n)表示全通滤波器Ai(ejω)系数ai的第n个元素,n=0,2,…,Ni;
[0018] 全通滤波器的实际相位为θi(ω):
[0019]
[0020] 全通滤波器的实际相位差为θei(ω):
[0021]
[0022] 其中ω∈[0,ωp]。
[0023]
[0024] 对公式(4)去三角化后:
[0025]
[0026] 其中ψi(ω)=0.5×(Niω+θdi(ω))。
[0027] s2.3、计算整体失真传输函数的实际相位误差
[0028] 滤波器组的整体失真传输函数表达式为 整体失真传输函数的相位误差是分析滤波器组和综合滤波器组的综合误差,体现了整个系统的相位误差情况,其理想相位、实际相位以及相位误差分别为:
[0029] θdT(ω)=θd1(2ω)+θd2(2ω)‑ω                           (6)[0030] θT(ω)=θ1(2ω)+θ2(2ω)‑ω                             (7)[0031] θeT(ω)=θT(ω)‑θdT(ω)                              (8)[0032] 因此
[0033]
[0034] 对公式(9)去三角化后:
[0035]
[0036] 步骤三、求解全通滤波器系数ai(k)
[0037] s3.1、得到目标优化函数表式如式(11)所示
[0038]
[0039] 其中a1、a2分别表示全通滤波器A1(ejω)、A2(ejω)的系数;a表示全通滤波器A1(ejω)jω jω与A2(e )的系数,。 表示第k次迭代时全通滤波器A1(e )的相位误差, 表示第k迭代jω jω
时全通滤波器A2(e )的相位误差, 表示第k次迭代时整体失真传输函数T(e )的相位jω
误差,M(a,ω)表示低通滤波器H0(e )的阻带幅度误差:
[0040]
[0041]
[0042]
[0043] s3.2、对步骤3.1中的目标函数进行LSK处理,将其转换为线性问题[0044]
[0045] 完成目标函数(11)中的分母部分的线性转化。
[0046] s3.3对目标函数分子部分由于两个全通滤波器系数存在互乘而导致的非线性问题采用交替迭代的方法进行处理。所述交替迭代为:在第k次迭代过程中,将第k‑1次迭代得到的两个全通滤波器的系数代入公式(15)中,得到常数化后的分母。然后将常数化的分母的值与第k‑1次得到的全通滤波器A1的系数 的值代入公式(16)中,求解全通滤波器A2的系数a2;通过凸优化处理工具箱CVX得到一组滤波器系数,再将这组滤波器系数代入公式(15)中更新分母的值,并将更新后的分母与全通滤波器A2的系数代入公式(17)中,求解全通滤波器A1的系数,并通过凸优化处理工具箱CVX得到一组滤波器系数;
[0047]
[0048]
[0049]
[0050] 反复交替,直至迭代终止得到两个全通滤波器的系数a。
[0051] 步骤四、目标判断
[0052] 设置一个实数μ作为判断目标,当目标函数满足 时结束迭代,将步骤三得到的滤波器系数作为最终得到的全通滤波器的系数;当目标函数不满足上述条件(k) (k‑1)时,进入步骤五。其中E 是第k次迭代中目标函数的最大值,E 是第k‑1次迭代中目标函数的最大值。
[0053]
[0054] 步骤五、修正Wi
[0055] 计算第k次迭代的相位误差 其中 表示第k次迭代的jω
全通滤波器Ai(e )的实际相位;然后计算相位误差 的包络
[0056] 设置一个实数判断目标ε,如果不满足 使然后令 k=k+1,返回步骤三;否则的话直接令
k=k+1,返回步骤三。
[0057] 本发明具有以下有益效果:
[0058] 1、通过数学处理,降低了滤波器组设计过程中公式的复杂程度和求解的困难程度。
[0059] 2、在求解目标函数的过程中,结合了不同的非线性处理方法,保证了滤波器组的整体相位误差,并且将阻带幅度误差控制在一定范围内。

实施方案

[0065] 以下结合附图对本发明作进一步的解释说明;基于全通滤波器混合迭代技术的QMF组设计方法,具体包括以下步骤:
[0066] 步骤一、确定基本参数
[0067] 滤波器组中两个全通滤波器的阶数分别为N1=11,N2=10,全频带上的频率点数L=8N1+1,低通滤波器的通带截止频率ωp=0.4π、阻带截止频率ωs=0.6π;最小化重构误差,最大化阻带衰减。
[0068] 步骤二、确定实际相位差
[0069] 确定滤波器的理想相位,计算滤波器的实际相位误差,然后确定滤波器组失真传输函数的实际相位误差。
[0070] 步骤三、求解滤波器系数
[0071] 设置目标函数,分别通过LSK处理与交替迭代的方法将非线性的目标函数的分母与分子常数化,即将目标函数转换为线性问题,然后使用CVX工具包求解滤波器的系数。
[0072] 步骤四、目标判断
[0073] 设置一个实数μ作为判断目标,当第k次迭代中目标函数的最大值E(k)满足输出此时的滤波器系数作为最终设计得到的滤波器系数,否则进入步骤五进行加权值修正。
[0074] 步骤五、加权值修正
[0075] 计算第k次迭代的相位误差 其中 表示第k次迭代的jω
全通滤波器Ai(e )的实际相位;然后计算相位误差 的包络
[0076] 设置一个实数判断目标ε,如果不满足 使然后令 k=k+1,返回步骤三;否则的话直接令
k=k+1,返回步骤三。
[0077] 经过3次迭代后,得到最终滤波器组中两个滤波器的系数如下表1所示:
[0078] 表1
[0079]
[0080] 图1、图2分别为两个滤波器的相位误差图,图3~5分别为两通道QMFB幅频响应、重构误差频率响应和相位误差响应。
[0081] 计算该两通道QMFB的最大阻带衰减(PSR),通带最大群延迟误差(MVPGD),QMFB重构响应的最大相位误差(MVPR),最大群延迟误差(MVGR),以及重构响应误(MVFBR),与其他方法得到的两通道QMFB进行对比,对比结果如表2所示:
[0082] 表2本设计方法与其他设计方法的关键指标的比较
[0083]
[0084] 从表2的对比结果可以看出,本实施例所获得的两通道QMFB在上述各项关键指标上都要优于其他设计方法。
[0085] [1]Y.Jou,Z.Lin and F.Chen“,Design of all‑pass‑based quadrature mirror filter banks using a Lyapunov error criterion,”2017 International Conference on Intelligent Informatics and BiomedicalSciences(ICIIBMS),Okinawa,2017,pp.39‑42,
[0086] [2]Y.D.Jou,Y.S.Lin,and F.K.Chen“, An efficient least‑squares design of two‑channel quadrature mirror filters using IIR all‑pass filters,”International Journal of Science and Engineering (IJSE),vol.4,no.1,pp.221‑224,March 2014.
[0087] [3]J.‑H.Lee and Y.‑H.Yang,“Design of two‑channel linear phase QMF banks based on realIIR all‑pass filters,”IEE Proc.‑Vision,Image and Signal Processing,vol.150,no.5,pp.331‑338,Oct.2003.
[0088] [4]S.S.Lawson and A.Klouche‑Djedid“, Technique for design of two‑channel approximately linear phase QMF bank and its application to image compression,”IEEE Pro‑c.Vision,Image and Signal Processing,vol.148,no.2,pp.85‑92,Apr.2001.
[0089] [5]Y.Jou,Z.Lin and F.Chen“, Improved Design of Two‑Channel Quadrature Mirror Filter Using Explicit Toeplitz Expressions,”2018 IEEE 7th Global Conference on Consumer Electronics (GCCE),Nara,2018,pp.1‑2,doi:10.1109/GCCE.2018.8574708.

附图说明

[0060] 图1为实施例中全通滤波器A1的相位误差。
[0061] 图2为实施例中全通滤波器A2的相位误差。
[0062] 图3为实施例中两通道QMFB幅频响应图。
[0063] 图4为实施例中两通道QMFB重构误差频率响应图。
[0064] 图5为实施例中两通道QMFB重构相位误差。
版权所有:盲专网 ©2023 zlpt.xyz  蜀ICP备2023003576号