CN111010146B - 一种基于快速滤波器组的信号重构结构及其设计方法 - Google Patents

一种基于快速滤波器组的信号重构结构及其设计方法 Download PDF

Info

Publication number
CN111010146B
CN111010146B CN201911280893.9A CN201911280893A CN111010146B CN 111010146 B CN111010146 B CN 111010146B CN 201911280893 A CN201911280893 A CN 201911280893A CN 111010146 B CN111010146 B CN 111010146B
Authority
CN
China
Prior art keywords
filter
stage
sub
channel
passband
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.)
Active
Application number
CN201911280893.9A
Other languages
English (en)
Other versions
CN111010146A (zh
Inventor
张文旭
崔鑫磊
张春光
姚雨双
马丹
秦涛
代雪飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN201911280893.9A priority Critical patent/CN111010146B/zh
Publication of CN111010146A publication Critical patent/CN111010146A/zh
Application granted granted Critical
Publication of CN111010146B publication Critical patent/CN111010146B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • H03H17/0211Frequency selective networks using specific transformation algorithms, e.g. WALSH functions, Fermat transforms, Mersenne transforms, polynomial transforms, Hilbert transforms

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Compression, Expansion, Code Conversion, And Decoders (AREA)
  • Surface Acoustic Wave Elements And Circuit Networks Thereof (AREA)

Abstract

本发明公开了基于快速滤波器组的信号重构结构及其设计方法,属于滤波器组技术领域。包括基于FFB的分析滤波器组(AFFB)结构和基于FFB的综合滤波器组(SFFB)结构。AFFB和SFFB结构都由L级滤波器组构成,AFFB将输入信号频谱均匀划分为N个子信道(N=2L),根据检测处理单元,通过SFFB对相邻子信道的信号进行重构还原成原来的信号。对各级的原型滤波器进行设计,每个滤波器可由对各级原型滤波器和互补滤波器的传递函数进行因子替换得到。FFB技术实现了低复杂度和良好频响特性的统一,适用于具有窄过渡带特性的滤波器组设计,使其具有良好的相邻子通道合并特性。本发明给出的结构与其他设计方法相比,具有更低的实现复杂度,节省了大量的硬件资源。

Description

一种基于快速滤波器组的信号重构结构及其设计方法
技术领域
本发明属于滤波器组技术领域,具体涉及基于快速滤波器组的信号重构结构及其设计方法。
背景技术
信号重构技术是近年来发展十分迅速的一种信号处理技术,它主要研究如何将处理过的信号还原成原来的信号。信号重构技术已广泛应用于光学通讯、图像处理、语音处理、地球物理信号处理、电子显微学、天文学、古地磁学、x射线结晶学等学科领域,并越来越受到人们的重视。
多速率信号处理是信号处理理论中的重要组成部分,凭借自身的许多优点受到了广大学者们的关注,并且应用于多个领域之中,而多速率信号处理理论的核心是由多个频率特性相关的滤波器构成的滤波器组。滤波器组是一组具有共同的输入或者共同的相加后的输出的滤波器。一般情况下,滤波器组用来实现对信号频率的各种分量进行分解,然后我们可以根据需要对各个子信号进行不同的处理或传输,在另一端再用一组滤波器将处理后的子信号进行综合。用来对信号进行分解的滤波器组我们称之为分析滤波器组,用来对信号进行综合的滤波器组我们称之为综合滤波器组。随着多速率信号处理技术的发展,滤波器组在数字通信、图像处理和信号处理等领域都得到了广泛的应用。
对于滤波器组的重构结构,比较常见的有基于树型结构、基于Farrow结构以及基于分析和综合滤波器组几种形式。在基于分析和综合滤波器组的形式中,目前主要运用的滤波器有基于可完全重构或近似完全重构的余弦调制滤波器组、正弦调制滤波器组、复指数调制滤波器组、调制DFT滤波器组和DFT滤波器组。由于FFB是一种非常有效的低复杂度滤波器组设计方法,本文利用AFFB和SFFB方法实现信号的重构。
在FFB研究方面,文献《快速滤波器组技术及在实时频谱分析中的应用研究》、《快速滤波器组在频率估计中的应用与研究》和文献《Low-complexity reconfigurable fastfilter bank for multi-standard wireless receivers》都是对基于FFB的分析滤波器组部分进行研究,与本发明设计结构不同;在滤波器组重构部分,文献《基于FRM的低复杂度信号重构系统研究》、专利《一种基于时域稀疏性的FRI信号重构方法》和专利《一种多通道声信号重构方法》都不是采用快速滤波器组实现信号重构,与本发明所用方法不同。
发明内容
本发明的目的在于提供实现低复杂度滤波器组信号重构的基于快速滤波器组的信号重构结构及其设计方法。
本发明的目的通过以下技术方案来实现:
基于快速滤波器组的信号重构结构及其设计方法,包括AFFB和SFFB结构及其设计方法,步骤如下:
针对AFFB结构及其设计方法,其设计步骤为:
步骤1:确定子信道数N,得到滤波器级数L。对于N通道AFFB,其对应的总级数L=log2 N,各级编号采用升序次序,即其各级编号k分别为0,…,L-2,L-1。全局通带边界频率为ωp,阻带边界频率为ωs,通带纹波δp,阻带衰减δs
步骤2:令φn(n=0,1,...,N-1)为通道n的通带边界频率,则有:
φ0=ωp (1)
AFFB的通道频率响应由第0级的子滤波器的频率响应决定,令θk(k=0,1,...,L-1)为第k级原型滤波器的通带边界频率,θk,m(k=0,1,...,L-1,m=0,1,...,2k-1)为第k级的第m个子滤波器的通带边界频率,则θ0,0与φ0满足:
θ0,0=φ0 (2)
步骤3:假设每个通道的过渡带宽为ωt,有ωt=ωsp。计算得出每个通道的带宽B以及通道0的通带边界频率φ0:
Q0=2L-1-k|k=0=2L-1 (3)
Q0=2L-1-k|k=0=2L-1 (4)
步骤4:对0级滤波器进行插值,插值系数为Q0=2L-1-k|k=0=2L-1,由升采样的频域特性可知,子滤波器插值后的频谱相较于原型滤波器压缩了2L-1倍,因此可得第0级原型滤波器通带边界θ0与子滤波器的通带边界θ0,0之间的关系式:
步骤5:结合(2)、(4)、(5)可得θ0关于通道数N以及通道过渡带宽ωt的表达式:
步骤6:由于滤波器的复杂度随着过渡带的变宽而降低,因此为了降低FFB的总体复杂度,设计各级原型滤波器主传递函数和互补滤波器的传递函数和/>(k=0,1,...,L-1)的过渡带应尽量地宽,各级原型滤波器通带边界频率θk的宽松边界:
结合式(5)、(6)、(7)得到各级原型滤波器的通带边界频率θk(k=1,...,L-1)的表达式:
步骤7:对于给定的全局通带纹波δp和阻带衰减δs,由于每个子滤波器对最终得到的通道频率响应的纹波幅度的影响十分复杂,因此很难确定各级的最优通带纹波。所以假设各级的子滤波器均为等纹波滤波器,通带纹波δ定义为:
δ=min(δp/Ls) (10)
步骤8:Hk,m(z)表示第k(k=0,1,...,L-1)级的第m(m=0,1,...,2k-1)个子滤波器,和/>分别表示Hk,m(z)上下支路的主传递函数和互补传递函数。为了更清楚的表示,用mn,k代替m表示n通道的第k级的第m个子滤波器的编号,则/>为mn,k的二进制编码的倒置编码对应的十进制编码所对应的值,则/>和/>可通过将原型滤波器对/>和/>的z因子替换为/>得到,其中/>
步骤9:AFFB结构中通道n(n=0,1,...,N-1)的频率响应是由各级的子滤波器级联而成,满足
其中αAN为归一化因子,为各级子滤波器的频率响应。
针对SFFB结构及其设计方法,其设计步骤为:
步骤1:确定子信道数N,得到滤波器级数L。对于N通道SFFB,其对应的总级数L=log2 N,各级编号采用降序次序,即其各级编号k分别为L-1,L-2,…,0。与AFFB采用相同的全局设计参数,全局通带边界频率为ωp,阻带边界频率为ωs,通带纹波δp,阻带衰减δs
步骤2:令φn(n=0,1,...,N-1)为通道n的通带边界频率,则有:
φ0=ωp (12)
SFFB的通道频率响应由第0级的子滤波器的频率响应决定,令θk(k=0,1,...,L-1)为第k级原型滤波器的通带边界频率,θk,m(k=0,1,...,L-1,m=0,1,...,2k-1)为第k级的第m个子滤波器的通带边界频率,则θ0,0与φ0满足:
θ0,0=φ0 (13)
步骤3:假设每个通道的过渡带宽为ωt,有ωt=ωsp。计算得出每个通道的带宽B以及通道0的通带边界频率φ0:
Q0=2L-1-k|k=0=2L-1 (14)
Q0=2L-1-k|k=0=2L-1 (15)
步骤4:对0级滤波器进行插值,插值系数为Q0=2L-1-k|k=0=2L-1,由升采样的频域特性可知,子滤波器插值后的频谱相较于原型滤波器压缩了2L-1倍,因此可得第0级原型滤波器通带边界θ0与子滤波器的通带边界θ0,0之间的关系式:
步骤5:结合(13)、(15)、(16)可得θ0关于通道数N以及通道过渡带宽ωt的表达式:
步骤6:由于滤波器的复杂度随着过渡带的变宽而降低,因此为了降低FFB的总体复杂度,设计各级原型滤波器主传递函数和互补滤波器的传递函数和/>(k=0,1,...,L-1)的过渡带应尽量地宽,各级原型滤波器通带边界频率θk的宽松边界:
结合式(16)、(17)、(18)得到各级原型滤波器的通带边界频率θk(k=1,...,L-1)的表达式:
步骤7:对于给定的全局通带纹波δp和阻带衰减δs,由于每个子滤波器对最终得到的通道频率响应的纹波幅度的影响十分复杂,因此很难确定各级的最优通带纹波。所以假设各级的子滤波器均为等纹波滤波器,通带纹波δ定义为:
δ=min(δp/Ls) (21)
步骤8:SFFB中各子滤波器可视作两输入单输出子系统,位于第k(k=0,1,...,L-1)级的第m(m=0,1,...,2k-1)个子滤波器传递函数为相应的,用/>分别表示/>上下支路的主传递函数和互补传递函数,则/>和/>满足式:
仍然表示mn,k的二进制编码的倒置编码对应的十进制编码所对应的值,Qk=2L-1-k仍然表示第k级的插值系数,则第k级的第mn,k个(mn,k=0,1,...,2k-1)子滤波器对的频率响应满足式(23)、(24):
步骤9:SFFB结构中通道的频率响应由各级的子滤波器级联而成,满足
其中为归一化因子,/>为SFFB各级子滤波器的传递函数。
本发明的有益效果在于:
本发明提供一种基于快速滤波器组的信号重构结构及其设计方法,利用AFFB和SFFB方法实现信号重构。AFFB将输入信号频谱均匀划分,SFFB再将划分后的子信道的信号重构还原成原来的信号。本发明给出的结构在相同指标的情况下,与其他设计方法相比,具有更低的实现复杂度,节省了大量的硬件资源。
附图说明
图1为本发明中基于FFB的信号重构整体结构图;
图2为本发明中基于FFB的分析滤波器组结构图;
图3为本发明中基于FFB的综合滤波器组结构图;
图4为本发明中基于FFB的滤波器组的频率响应图;
图5为本发明的输入信号幅频特性仿真图;
图6为本发明的经过AFFB结构后子信道的信号幅频特性仿真图;
图7为本发明的经过SFFB结构后重构信号幅频特性仿真图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步说明:
基于快速滤波器组的信号重构结构及其设计方法,其特征在于该结构采用基于FFB的分析滤波器组(AFFB,Analysis Fast Filter Bank)结构和基于FFB的综合滤波器组(SFFB,Synthesis Fast Filter Bank)结构组合的方式实现信号的重构,基于快速滤波器组的信号重构结构及其设计方法,包括AFFB和SFFB结构及其设计方法;
AFFB和SFFB结构都由L级滤波器组构成,AFFB将输入信号频谱均匀划分为N个子信道,根据检测处理单元,通过SFFB对相邻子信道的信号进行重构还原成原来的信号,对各级的原型滤波器进行设计,每个滤波器可由对各级原型滤波器和互补滤波器的传递函数进行因子替换得到。
图1为基于FFB的信号重构整体结构图,左半部分为AFFB结构,更详细的信道划分如图2所示,右半部分为SFFB结构,更详细的信道划分如图3所示。其中间环节部分为信号的检测单元。
图2为基于FFB的分析滤波器组结构图,设计方法如下:
步骤1:确定子信道数N,得到滤波器级数L。对于N通道AFFB,其对应的总级数L=log2 N,各级编号采用升序次序,即其各级编号k分别为0,…,L-2,L-1。全局通带边界频率为ωp,阻带边界频率为ωs,通带纹波δp,阻带衰减δs
步骤2:令φn(n=0,1,...,N-1)为通道n的通带边界频率,则有:
φ0=ωp (1)
AFFB的通道频率响应由第0级的子滤波器的频率响应决定,令θk(k=0,1,...,L-1)为第k级原型滤波器的通带边界频率,θk,m(k=0,1,...,L-1,m=0,1,...,2k-1)为第k级的第m个子滤波器的通带边界频率,则θ0,0与φ0满足:
θ0,0=φ0 (2)
步骤3:假设每个通道的过渡带宽为ωt,有ωt=ωsp。计算得出每个通道的带宽B以及通道0的通带边界频率φ0:
Q0=2L-1-k|k=0=2L-1 (3)
Q0=2L-1-k|k=0=2L-1 (4)
步骤4:对0级滤波器进行插值,插值系数为Q0=2L-1-k|k=0=2L-1,由升采样的频域特性可知,子滤波器插值后的频谱相较于原型滤波器压缩了2L-1倍,因此可得第0级原型滤波器通带边界θ0与子滤波器的通带边界θ0,0之间的关系式:
步骤5:结合(2)、(4)、(5)可得θ0关于通道数N以及通道过渡带宽ωt的表达式:
步骤6:由于滤波器的复杂度随着过渡带的变宽而降低,因此为了降低FFB的总体复杂度,设计各级原型滤波器主传递函数和互补滤波器的传递函数和/>(k=0,1,...,L-1)的过渡带应尽量地宽,各级原型滤波器通带边界频率θk的宽松边界:
结合式(5)、(6)、(7)得到各级原型滤波器的通带边界频率θk(k=1,...,L-1)的表达式:
步骤7:对于给定的全局通带纹波δp和阻带衰减δs,由于每个子滤波器对最终得到的通道频率响应的纹波幅度的影响十分复杂,因此很难确定各级的最优通带纹波。所以假设各级的子滤波器均为等纹波滤波器,通带纹波δ定义为:
δ=min(δp/Ls) (10)
步骤8:Hk,m(z)表示第k(k=0,1,...,L-1)级的第m(m=0,1,...,2k-1)个子滤波器,和/>分别表示Hk,m(z)上下支路的主传递函数和互补传递函数。为了更清楚的表示,用mn,k代替m表示n通道的第k级的第m个子滤波器的编号,则/>为mn,k的二进制编码的倒置编码对应的十进制编码所对应的值,则/>和/>可通过将原型滤波器对/>和/>的z因子替换为/>得到,其中/>
步骤9:AFFB结构中通道n(n=0,1,...,N-1)的频率响应是由各级的子滤波器级联而成,满足
其中αAN为归一化因子,为各级子滤波器的频率响应。
图3为基于FFB的综合滤波器组结构图,设计方法如下:
步骤1:确定子信道数N,得到滤波器级数L。对于N通道SFFB,其对应的总级数L=log2 N,各级编号采用降序次序,即其各级编号k分别为L-1,L-2,…,0。与AFFB采用相同的全局设计参数,全局通带边界频率为ωp,阻带边界频率为ωs,通带纹波δp,阻带衰减δs
步骤2:令φn(n=0,1,...,N-1)为通道n的通带边界频率,则有:
φ0=ωp (12)
SFFB的通道频率响应由第0级的子滤波器的频率响应决定,令θk(k=0,1,...,L-1)为第k级原型滤波器的通带边界频率,θk,m(k=0,1,...,L-1,m=0,1,...,2k-1)为第k级的第m个子滤波器的通带边界频率,则θ0,0与φ0满足:
θ0,0=φ0 (13)
步骤3:假设每个通道的过渡带宽为ωt,有ωt=ωsp。计算得出每个通道的带宽B以及通道0的通带边界频率φ0:
Q0=2L-1-k|k=0=2L-1 (14)
Q0=2L-1-k|k=0=2L-1 (15)
步骤4:对0级滤波器进行插值,插值系数为Q0=2L-1-k|k=0=2L-1,由升采样的频域特性可知,子滤波器插值后的频谱相较于原型滤波器压缩了2L-1倍,因此可得第0级原型滤波器通带边界θ0与子滤波器的通带边界θ00之间的关系式:
步骤5:结合(13)、(15)、(16)可得θ0关于通道数N以及通道过渡带宽ωt的表达式:
步骤6:由于滤波器的复杂度随着过渡带的变宽而降低,因此为了降低FFB的总体复杂度,设计各级原型滤波器主传递函数和互补滤波器的传递函数和/>(k=0,1,...,L-1)的过渡带应尽量地宽,各级原型滤波器通带边界频率θk的宽松边界:
结合式(16)、(17)、(18)得到各级原型滤波器的通带边界频率θk(k=1,...,L-1)的表达式:
步骤7:对于给定的全局通带纹波δp和阻带衰减δs,由于每个子滤波器对最终得到的通道频率响应的纹波幅度的影响十分复杂,因此很难确定各级的最优通带纹波。所以假设各级的子滤波器均为等纹波滤波器,通带纹波δ定义为:
δ=min(δp/Ls) (21)
步骤8:SFFB中各子滤波器可视作两输入单输出子系统,位于第k(k=0,1,...,L-1)级的第m(m=0,1,...,2k-1)个子滤波器传递函数为相应的,用/>分别表示/>上下支路的主传递函数和互补传递函数,则/>和/>满足式:
仍然表示mn,k的二进制编码的倒置编码对应的十进制编码所对应的值,Qk=2L-1-k仍然表示第k级的插值系数,则第k级的第mn,k个(mn,k=0,1,...,2k-1)子滤波器对的频率响应满足式(23)、(24):
步骤9:SFFB结构中通道的频率响应由各级的子滤波器级联而成,满足
其中为归一化因子,/>为SFFB各级子滤波器的传递函数。
为了验证本发明的有效性,根据上述设计方法进行了仿真实验。
设计的全局指标如下:(1)通道数N=16,FFB级数L=4,(2)全局通带边界频率ωp=0.0612π,ωs=0.0638π,过渡带宽度Δb=0.0026π,(3)通带纹波δp小于0.04dB,阻带衰减δs大于80dB。由全局指标确定滤波器的各级原型滤波器的通带频率θk分别为0.4896π、0.2552π、0.1276π和0.0638π。系统采样频率设定为960KHz,输入信号为2个实信号同时输入,分别为110KHz-180KHz的线性调频信号和300KHz的正弦信号。图4为基于FFB的滤波器组的频率响应图,图5为输入信号幅频特性仿真图;图6为信号经过AFFB结构后子信道的幅频特性仿真图;图7为信号经过SFFB结构后重构信号幅频特性仿真图。
表1为在滤波器过渡带宽相等的情况下(Δb=0.0026π),16通道基于FFB的信号重构结构与基于低通滤波器组(LPFB)、多相滤波器组(PPFB)、频率响应屏蔽技术的滤波器组(FRM)几种常用方法的重构结构的复杂度以及延时的对比。
表1复杂度及延时对比
通过表1的重构滤波器组复杂度数据表数据对比可知,在滤波器过渡带宽相等的情况下(Δb=0.0026π),本发明所推导出的基于FFB的信号重构结构相比于基于多相结构的信号重构结构、基于FRM的信号重构结构分别节省了66.79%和62.16%的复杂度,延时比多相结构增加了3.85倍,但是比基于FRM的重构结构的延时减少了56.95%。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (2)

1.一种基于快速滤波器组的信号重构结构的设计方法,其特征在于:采用基于FFB的分析滤波器组—AFFB结构和基于FFB的综合滤波器组—SFFB结构组合的方式实现信号的重构;
AFFB结构和SFFB结构都由L级滤波器组构成,AFFB结构用于将输入信号频谱均匀划分为N个子信道,N=2L;SFFB结构用于对相邻子信道的信号进行重构,还原成原来的信号;
所述基于FFB的分析滤波器组—AFFB结构中每级的滤波器设计方法如下:
步骤1.1:确定子信道数N和滤波器总级数L;
对于N通道AFFB结构,其对应的滤波器总级数L=log2N,各级编号采用升序次序,即其各级编号k分别为0,...,L-2,L-1,全局通带边界频率为ωp,阻带边界频率为ωs,通带纹波δp,阻带衰减δs
步骤1.2:令φn为通道n的通带边界频率,则有:
φ0=ωp (1)
AFFB结构的通道频率响应由第0级的子滤波器的频率响应决定,令θk为第k级原型滤波器的通带边界频率,θk,m为第k级的第m个子滤波器的通带边界频率,则θ0,0与φ0满足:
θ0,0=φ0 (2)
其中,n=0,1,...,N-1;m=0,1,...,2k-1;
步骤1.3:假设每个通道的过渡带宽为ωt,有ωt=ωsp,计算得出每个通道的带宽B以及通道0的通带边界频率φ0
步骤1.4:对0级滤波器进行插值,插值系数Q0为:
Q0=2L-1-k|k=0=2L-1 (4)
由升采样的频域特性可知,子滤波器插值后的频谱相较于原型滤波器压缩了2L-1倍,因此可得第0级原型滤波器通带边界θ0与子滤波器的通带边界θ0,0之间的关系式:
步骤1.5:结合(2)、(4)、(5)可得θ0关于通道数N以及通道过渡带宽ωt的表达式:
步骤1.6:由于滤波器的复杂度随着过渡带的变宽而降低,因此为了降低FFB的总体复杂度,设计各级原型滤波器主传递函数和互补滤波器的传递函数和/>的过渡带应尽量地宽,各级原型滤波器通带边界频率θk的宽松边界:
结合式(5)、(6)、(7)得到各级原型滤波器的通带边界频率θk的表达式:
步骤1.7:对于给定的全局通带纹波δp和阻带衰减δs,由于每个子滤波器对最终得到的通道频率响应的纹波幅度的影响十分复杂,因此很难确定各级的最优通带纹波,所以假设各级的子滤波器均为等纹波滤波器,通带纹波δ定义为:
步骤1.8:Hk,m(z)表示第k级的第m个子滤波器,和/>分别表示Hk,m(z)上下支路的主传递函数和互补传递函数,为了更清楚的表示,用mn,k代替m表示n通道的第k级的第m个子滤波器的编号,则/>为mn,k的二进制编码的倒置编码对应的十进制编码所对应的值,则/>和/>可通过将原型滤波器对/>和/>的z因子替换为得到,其中/>
步骤1.9:AFFB结构中通道n的频率响应是由各级的子滤波器级联而成,满足:
其中αAN为归一化因子,为各级子滤波器的频率响应。
2.根据权利要求1所述的一种基于快速滤波器组的信号重构结构的设计方法,其特征在于:所述的基于FFB的综合滤波器组—SFFB结构中每级的滤波器设计方法如下:
步骤2.1:确定子信道数N和滤波器总级数L;
对于N通道SFFB结构,其对应的滤波器总级数L=log2N,各级编号采用降序次序,即其各级编号k分别为L-1,L-2,…,0,与AFFB结构采用相同的全局设计参数,全局通带边界频率为ωp,阻带边界频率为ωs,通带纹波δp,阻带衰减δs
步骤2.2:令φn为通道n的通带边界频率,则有:
φ0=ωp (12)
SFFB结构的通道频率响应由第0级的子滤波器的频率响应决定,令θk为第k级原型滤波器的通带边界频率,θk,m为第k级的第m个子滤波器的通带边界频率,则θ0,0与φ0满足:
θ0,0=φ0 (13)
步骤2.3:假设每个通道的过渡带宽为ωt,有ωt=ωsp,计算得出每个通道的带宽B以及通道0的通带边界频率φ0
Q0=2L-1-k|k=0=2L-1 (15)
步骤2.4:对0级滤波器进行插值,插值系数Q0为:
Q0=2L-1-k|k=0=2L-1 (15)
由升采样的频域特性可知,子滤波器插值后的频谱相较于原型滤波器压缩了2L-1倍,因此可得第0级原型滤波器通带边界θ0与子滤波器的通带边界θ0,0之间的关系式:
步骤2.5:结合(13)、(15)、(16)可得θ0关于通道数N以及通道过渡带宽ωt的表达式:
步骤2.6:由于滤波器的复杂度随着过渡带的变宽而降低,因此为了降低FFB的总体复杂度,设计各级原型滤波器主传递函数和互补滤波器的传递函数和/>的过渡带应尽量地宽,各级原型滤波器通带边界频率θk的宽松边界:
结合式(16)、(17)、(18)得到各级原型滤波器的通带边界频率θk的表达式:
步骤2.7:对于给定的全局通带纹波δp和阻带衰减δs,由于每个子滤波器对最终得到的通道频率响应的纹波幅度的影响十分复杂,因此很难确定各级的最优通带纹波,所以假设各级的子滤波器均为等纹波滤波器,通带纹波δ定义为:
步骤2.8:SFFB中各子滤波器可视作两输入单输出子系统,位于第k级的第m个子滤波器传递函数为相应的,用/>和/>分别表示/>上下支路的主传递函数和互补传递函数,则/>和/>满足式:
仍然表示mn,k的二进制编码的倒置编码对应的十进制编码所对应的值,Qk=2L -1-k仍然表示第k级的插值系数,则第k级的第mn,k个子滤波器对的频率响应满足式(23)、(24):
步骤2.9:SFFB结构中通道的频率响应由各级的子滤波器级联而成,满足
其中为归一化因子,/>为SFFB各级子滤波器的传递函数。
CN201911280893.9A 2019-12-13 2019-12-13 一种基于快速滤波器组的信号重构结构及其设计方法 Active CN111010146B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911280893.9A CN111010146B (zh) 2019-12-13 2019-12-13 一种基于快速滤波器组的信号重构结构及其设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911280893.9A CN111010146B (zh) 2019-12-13 2019-12-13 一种基于快速滤波器组的信号重构结构及其设计方法

Publications (2)

Publication Number Publication Date
CN111010146A CN111010146A (zh) 2020-04-14
CN111010146B true CN111010146B (zh) 2023-08-01

Family

ID=70114673

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911280893.9A Active CN111010146B (zh) 2019-12-13 2019-12-13 一种基于快速滤波器组的信号重构结构及其设计方法

Country Status (1)

Country Link
CN (1) CN111010146B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104579243A (zh) * 2015-02-11 2015-04-29 桂林电子科技大学 一种dft调制滤波器组的快速设计方法
CN105281707A (zh) * 2015-09-09 2016-01-27 哈尔滨工程大学 一种动态可重构滤波器组低复杂度的实现方法
CN106953649A (zh) * 2017-03-03 2017-07-14 东南大学 一种基于外插脉冲响应法实现的快速滤波器组的非均匀数字信道化方法
CN107294512A (zh) * 2017-05-25 2017-10-24 哈尔滨工程大学 一种基于树型结构的非均匀滤波器组滤波方法
CN110210101A (zh) * 2019-05-27 2019-09-06 哈尔滨工程大学 一种基于cem frm的动态非均匀窄过渡带滤波器组及设计方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103997314B (zh) * 2014-06-05 2015-02-25 山东大学 一种改进的二级frm滤波器设计方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104579243A (zh) * 2015-02-11 2015-04-29 桂林电子科技大学 一种dft调制滤波器组的快速设计方法
CN105281707A (zh) * 2015-09-09 2016-01-27 哈尔滨工程大学 一种动态可重构滤波器组低复杂度的实现方法
CN106953649A (zh) * 2017-03-03 2017-07-14 东南大学 一种基于外插脉冲响应法实现的快速滤波器组的非均匀数字信道化方法
CN107294512A (zh) * 2017-05-25 2017-10-24 哈尔滨工程大学 一种基于树型结构的非均匀滤波器组滤波方法
CN110210101A (zh) * 2019-05-27 2019-09-06 哈尔滨工程大学 一种基于cem frm的动态非均匀窄过渡带滤波器组及设计方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
Haridas N等.Low-complexity technique to get arbitrary variation in the bandwidth of a digital FIR filter.《IET Signal Processing》.2016,第372-377页. *
Zhang Wenxu等.esign and FPGA Implementation of a Novel Efficient FRM-based Channelized Receiver Structure.《IEEE ACCESS》.2019,第114778-114787页. *
Zhang Wenxu等.Low-complexity channelizer based on FRM for passive radar multi-channel wideband receiver.《Circuits, Systems, and Signal Processing》.2019,第420-438页. *
刘小蒙 ; 邵高平 ; 汪洋 ; 岳强 ; .一种改进的动态信道化滤波方法.信息工程大学学报.2016,(第04期),第448-453页. *
唐鹏飞 ; 林钱强 ; 袁斌 ; 陈曾平 ; .一种新的动态信道化接收机设计方法.国防科技大学学报.2013,(第03期),第164-169页. *
张文旭 ; 陈亚静 ; 陈涛 ; 窦峥 ; .基于FRM的偶型排列信道化改进结构.北京邮电大学学报.2016,(第03期),第27-33页. *
陈涛 ; 岳玮 ; 刘颜琼 ; 司锡才 ; .宽带数字信道化接收机部分信道重构技术.哈尔滨工程大学学报.2011,(第12期),第1610-1616页. *

Also Published As

Publication number Publication date
CN111010146A (zh) 2020-04-14

Similar Documents

Publication Publication Date Title
Evangelista et al. Discrete frequency warped wavelets: Theory and applications
CN102318189B (zh) 低延迟调制滤波器组
CN105281707B (zh) 一种动态可重构滤波器组的实现方法
CN104579243B (zh) 一种dft调制滤波器组的快速设计方法
Kumar et al. Design of nearly perfect reconstructed non-uniform filter bank by constrained equiripple FIR technique
Sudharman et al. Design of power efficient variable bandwidth non-maximally decimated FRM filters for wideband channelizer
Evangelista et al. Frequency-warped filter banks and wavelet transforms: A discrete-time approach via Laguerre expansion
CN105279350A (zh) 近似完全重构非均匀余弦调制滤波器组的设计方法
CN107294512B (zh) 一种基于树型结构的非均匀滤波器组滤波方法
Bregovic et al. A general-purpose optimization approach for designing two-channel FIR filterbanks
Chan et al. Theory and design of a class of cosine-modulated non-uniform filter banks
JPH10285031A (ja) 帯域合成フィルタバンク及びフィルタリング方法並びに帯域分割フィルタバンク及びフィルタリング方法並びに符号化装置並びに復号化装置
CN104486711A (zh) 用于数字助听器的低复杂度可调滤波器组及其工作方法
CN111010146B (zh) 一种基于快速滤波器组的信号重构结构及其设计方法
CN110957996A (zh) 一种基于abc算法的无乘法器frm滤波器组优化设计方法
Sudharman et al. Design of reconfigurable FRM channelizer using resource shared non-maximally decimated masking filters
CN115296647B (zh) 基于frm技术的wola滤波器组及子带分割设计方法
Bindiya et al. Design of totally multiplier-less sharp transition width tree structured filter banks for non-uniform discrete multitone system
CN106160702A (zh) 近似完全重构单原型dft调制滤波器组的设计方法
CN106712836B (zh) 基于近似精确重构滤波器组的信道化器设计方法及装置
EP1421579B1 (en) Audio coding with non-uniform filter bank
CN107241082B (zh) 基于凸优化松弛的dft调制滤波器组的设计方法
Maneesha et al. A Channel Combiner approach for the design of Near Perfect Reconstruction non uniform Filter Banks
Fettweis et al. Realization of filter banks by means of wave digital filters
CN114331926B (zh) 基于换元思想的两通道图滤波器组系数设计优化方法

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