CN105281707B - 一种动态可重构滤波器组的实现方法 - Google Patents
一种动态可重构滤波器组的实现方法 Download PDFInfo
- Publication number
- CN105281707B CN105281707B CN201510570781.2A CN201510570781A CN105281707B CN 105281707 B CN105281707 B CN 105281707B CN 201510570781 A CN201510570781 A CN 201510570781A CN 105281707 B CN105281707 B CN 105281707B
- Authority
- CN
- China
- Prior art keywords
- filter
- sub
- band
- filter group
- bands
- 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
Links
Landscapes
- Compression, Expansion, Code Conversion, And Decoders (AREA)
- Complex Calculations (AREA)
- Filters That Use Time-Delay Elements (AREA)
Abstract
本发明提供的是一种动态可重构滤波器组的低复杂度实现方法。该方法依照以下操作步骤进行:设计预划分子带带宽的、阶数为划分子带数目整数倍的有限长单位脉冲响应FIR低通原型滤波器;通过将设计生成的低通原型滤波器系数进行抽取构造各分支子带滤波器;将信号进行一系列的单位延迟、抽取、子带滤波、离散傅里叶运算等分析部分的处理分离得到子带数目的基带信号;根据经分析滤波器组后得到的低频子带信号确定需重构的子带数目,对低通原型滤波器系数进行动态抽取、变换,构造综合滤波器组的各子带滤波器;得到基带的重构信号。本发明的设计方法有效降低了动态可重构滤波器组的复杂度,且其动态实现显著提高了***的灵活性。
Description
技术领域
本发明涉及的是多速率数字滤波器组领域,并提供了一种动态可重构滤波器组低复杂度实现方法,该方法可应用于软件无线电、数字通信***、语音信号处理、图像处理以及电子侦察***等相关领域中。
背景技术
当前滤波器组在多速率数字信号处理和数字通信***中占据有非常重要的地位。滤波器组广泛应用于语音信号处理、视频信号处理、图像信号处理、数字通信、信道化接收等多种领域。
目前,随着滤波器组划分子带数目的增加对有限的硬件资源提出了较高的实现难度,这均要求滤波器组***能够具有更高的灵活性和较低的复杂度。在相关文献中均有提出降低滤波器组设计复杂度的方法,如文献《FIR数字滤波器组的低复杂度实现方法及装置》中提出通过构造滤波器组的内核滤波器来降低复杂度,文献《分析滤波器组部分和分析滤波方法》中提出了优化改进滤波器组的方法等。本发明在较低复杂度实现滤波器组的基础上还实现了其对信号的动态可重构功能,通过动态重构的方式提高了重构滤波器组的灵活性,大大降低了实现***的复杂度。
发明内容
本发明的目的是提供一种具有较高的灵活性和较低的复杂度的方法设计实现滤波器组的动态可重构,通过检测分析滤波子带信号动态确定重构子带滤波器系数,实现***的灵活性设计。
本发明提供了一种动态重构滤波器组的低复杂度实现方法,主要操作步骤是:
101.设计预划分子带带宽的、阶数为划分子带数目整数倍的有限长单位脉冲响应FIR低通原型滤波器;
102.通过将设计生成的低通原型滤波器系数进行抽取,构造各分支子带滤波器;
103.将信号进行一系列的单位延迟、抽取、子带滤波、离散傅里叶运算等分析部分的处理分离得到子带数目的基带信号;
104.根据经分析滤波器组后得到的低频子带信号确定需重构的子带数目,对低通原型滤波器系数进行动态抽取、变换,构造综合滤波器组的各子带滤波器;
105.对经分析部分分离出的子带数目的基带信号进行一系列的离散傅里叶运算、子带滤波、复因子调制、内插、单位延迟等综合部分的处理得到基带的重构信号;
其中,步骤101中具体操作如下:
所述低通原型滤波器h(n)的通带带宽为2π/M,级数基于以下公式:
N=βM,
在此M为划分子带数目,n=0,1,…,N-1,β为整数,
其中所述低通原型滤波器的优化是通过最小化整个可重构滤波器组实现时出现的混叠误差、幅度失真误差所构成的目标函数Od(m)来实现的,
Od(m)=mOα+(1-m)Oβ
其中m为加权,Oα为混叠误差,Oβ为幅度失真误差。
其中,步骤102中具体操作如下:
构造子带滤波器,各子带滤波器系数是对原型滤波器的抽取得到的,其中子带系数是基于以下公式:
hl(n)=h(Mn+l),
在此M为划分子带数目以及抽取倍数,l为子带序号索引,n为滤波器系数索引,而hl(n)为属于子带序号索引l的滤波器的系数,其中l=0,1,…M-1,n=0,1,…N/M-1。
其中,步骤103中具体操作如下:
分离实信号得到基带信号的分析滤波器组部分,包括可实现延迟和抽取功能的移位寄存器装置(301)、子带滤波器组(302)、离散傅里叶DFT变换装置(303),其中在离散傅里叶DFT变换装置的输出端,获得经分析滤波器组后的子带低频信号。
其中,步骤104中具体操作如下:
所述的实信号在经分析滤波器组后得到的低频子带信号,可经过对其检测分析实现动态重构,包括检测判定装置(201)及动态控制移位寄存器(202)。
其中,步骤105中具体操作如下:
并行FFT变换装置(401),复调制因子模块(402),子带滤波器组(403)、插值器(404)、移位寄存及加法器装置(405),其中经过移位寄存及加法器装置(405)后,获得经过整个动态重构滤波器组***后的低频输入信号。
本发明设计的动态重构滤波器组低复杂度方法具有如下的优点:与传统滤波器组结构相比,此结构在很大程度上降低了复杂度,减小了运算量,提高了运算效率,节约利用了有限的硬件资源,而且具有较好的灵活性,能够根据情况实现动态重构。
附图说明
图1本发明动态可重构滤波器组低复杂度实现方法流程图;
图2本发明实现滤波器组动态重构的方法框图;
图3本发明分析滤波部分的实现框图;
图4本发明综合滤波部分的实现框图。
具体实施方式
下面结合附图举例对本发明做更详细地描述:
图1是本发明动态可重构滤波器组低复杂度实现方法流程图,所述方法包括:
101.设计预划分子带带宽的、阶数为划分子带数目整数倍的有限长单位脉冲响应FIR低通原型滤波器;
详细的实施说明如下:
设计低通原型滤波器h(n)的通带带宽为2π/M,级数基于以下公式:
N=βM,
在此M为划分子带数目,n=0,1,…,N-1,β为整数,
按照图3,经过分析滤波后,DFT变换装置(303)的输出子带低频信号由以下公式给出:
其中,WM k=ej2πk/M为输入信号的调制因子,Ak(z)为k索引下的子带输出,H(z)为原型滤波器h(n)对应的z变换,k为子带输出号索引,k=0,1,…M-1
按照图4,经过综合滤波后,移位寄存器及加法器(405)的输出信号由以下公式给出:
其中,H(z)为原型滤波器h(n)对应的z变换,F(z)为原型滤波器f(n)对应的z变换,k1为综合子带的起始序号索引。
其中所述低通原型滤波器的优化是通过优化减小原型滤波器阻带衰减不理想引入的混叠误差和幅度失真误差Oβ构成的目标函数Od(m)来实现的,其中Od(m)为:
Od(m)=mOα+(1-m)Oβ
其中m为加权,Oα为混叠误差,Oβ为幅度失真误差。
102.通过将设计生成的低通原型滤波器系数进行抽取,构造各分支子带滤波器;
详细的实施说明如下:
构造子带滤波器,各子带滤波器系数是对原型滤波器的抽取得到的,其中子带系数是基于以下公式:
hl(n)=h(Mn+l),
其中M为划分子带数目以及抽取倍数,l为子带序号索引,n为滤波器系数索引,hl(n)为属于子带序号索引l的滤波器的系数,l=0,1,…M-1,n=0,1,…N/M-1。
103.将信号进行一系列的单位延迟、抽取、子带滤波、离散傅里叶运算等分析部分的处理分离得到子带数目的基带信号;
详细的实施说明如下:
分离实信号得到基带信号的分析滤波器组部分,包括可实现延迟和抽取功能的移位寄存器、子带滤波器组、离散傅里叶变换装置,其中在离散傅里叶变换装置的输出端,获得经分析滤波器组后的子带低频信号。
按照图3,经过分析滤波后,DFT变换装置(303)的输出子带低频信号由以下公式给出:
其中,WM k=ej2πk/M为输入信号的调制因子,Ak(z)为k索引下的子带输出,H(z)为原型滤波器h(n)对应的z变换,k为子带输出号索引,k=0,1,…M-1
104.根据经分析滤波器组后得到的低频子带信号确定需重构的子带数目,对低通原型滤波器系数进行动态抽取、变换,构造综合滤波器组的各子带滤波器;
详细的实施说明如下:
按照图4,经过综合滤波后,移位寄存器及加法器(405)的输出信号由以下公式给出:
其中,H(z)为原型滤波器h(n)对应的z变换,F(z)为原型滤波器f(n)对应的z变换,k1为综合子带的起始序号索引
将上式转化为矩阵形式关系如下:
设计的原型滤波器通过优化使其满足条件T=[1 0 … 0],
移位寄存器及加法器(405)的输出信号转化为:
即整个***完成的仅是对输入信号X(z)的尺度变换和频谱搬移,实现了重构。
按照图2,检测判定装置(201)给出需要综合的子带数目N,当N≠2的整次幂时,取大于N的最接近的2的整次幂数,以用快速变换FFT提高***运算的效率。利用N动态控制制移位寄存器(202)装置即可得到动态综合滤波各子带滤波器系数。
105.对经分析部分分离出的子带数目的基带信号进行一系列的离散傅里叶运算、子带滤波、复因子调制、内插、单位延迟等综合部分的处理得到基带的重构信号;
详细的实施说明如下:
并行FFT变换装置(401),复调制因子模块(402),子带滤波器组(403)、插值器(404)、移位寄存及加法器装置(405),其中经过移位寄存及加法器装置(405)后,获得经过整个动态重构滤波器组后的低频输入信号。
按照图4,经过综合滤波后,移位寄存器及加法器(405)的输出信号由以下公式给出:
即整个***完成的仅是对输入信号X(z)的尺度变换和频谱搬移,实现了重构。
Claims (3)
1.一种动态可重构滤波器组的实现方法,其特征在于包括以下步骤:
101:预划分子带带宽的有限长单位脉冲响应FIR低通原型滤波器,阶数为划分子带数目整数倍;
102:将生成的低通原型滤波器系数进行抽取构造各分支子带滤波器;
103:将信号进行单位延迟、抽取、子带滤波、离散傅里叶运算得到子带数目的基带信号;
104:根据分析滤波器组后得到的低频子带信号确定需重构的子带数目,对低通原型滤波器系数进行动态抽取、变换,构造综合滤波器组的各子带滤波器;
105:对分析部分分离出的子带数目的基带信号进行离散傅里叶运算、子带滤波、复因子调制、内插、单位延迟处理得到基带的重构信号。
2.根据权利要求1所述的动态可重构滤波器组的实现方法,其特征在于,分析滤波器组部分,其中所述子带数目为M,对有限长单位脉冲响应FIR低通原型滤波器,其特征在于,包括:
所述低通原型滤波器h(n)的通带带宽为2π/M,阶数基于以下公式:
N=βM,
其中n=0,1,…,N-1和β为整数,
所述低通原型滤波器的优化是通过最小化整个可重构滤波器组实现时出现的混叠误差、幅度失真误差所构成的目标函数Od(m)来实现的,
Od(m)=mOα+(1-m)Oβ,
其中m为加权,Oα为混叠误差,Oβ为幅度失真误差。
3.根据权利要求1所述的动态可重构滤波器组的实现方法,其特征在于,所述的构造各分支子带滤波器是M个子带滤波器通过抽取器对低通原型滤波器的系数进行M倍抽取,利用M组抽取得到系数构造,其中子带系数是基于以下公式:
hl(n)=h(Mn+l),
其中l为子带序号索引,n为滤波器系数索引,hl(n)为属于子带序号索引l的滤波器的系数,其中l=0,1,…M-1,n=0,1,…N/M-1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510570781.2A CN105281707B (zh) | 2015-09-09 | 2015-09-09 | 一种动态可重构滤波器组的实现方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510570781.2A CN105281707B (zh) | 2015-09-09 | 2015-09-09 | 一种动态可重构滤波器组的实现方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105281707A CN105281707A (zh) | 2016-01-27 |
CN105281707B true CN105281707B (zh) | 2018-12-25 |
Family
ID=55150155
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510570781.2A Active CN105281707B (zh) | 2015-09-09 | 2015-09-09 | 一种动态可重构滤波器组的实现方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105281707B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107294512B (zh) * | 2017-05-25 | 2020-07-28 | 哈尔滨工程大学 | 一种基于树型结构的非均匀滤波器组滤波方法 |
CN107241107B (zh) * | 2017-06-28 | 2019-03-29 | 电子科技大学 | 一种数字信道化滤波器组实现方法 |
CN108900461B (zh) * | 2018-06-28 | 2020-10-13 | 东南大学 | 一种基于大规模mimo的无线通信***宽带信号设计方法 |
CN110149118A (zh) * | 2019-04-28 | 2019-08-20 | 北京航空航天大学 | 一种基于非均匀滤波器组的动态信道化方法 |
CN111010146B (zh) * | 2019-12-13 | 2023-08-01 | 哈尔滨工程大学 | 一种基于快速滤波器组的信号重构结构及其设计方法 |
CN111277245B (zh) * | 2020-03-24 | 2023-03-10 | 西安电子科技大学 | 一种滤波正交频分复用***的低阶子带滤波器设计方法 |
CN112039497A (zh) * | 2020-08-21 | 2020-12-04 | 安徽蓝讯电子科技有限公司 | 一种滤波器组多载波的子带滤波器 |
CN113992190B (zh) * | 2021-12-23 | 2022-04-01 | 南京先锋硕通无线技术有限公司 | 一种用于dvb-rcs2的双层滤波器组设计方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101303764A (zh) * | 2008-05-16 | 2008-11-12 | 西安电子科技大学 | 基于非下采样轮廓波的多传感器图像自适应融合方法 |
CN101635564A (zh) * | 2001-04-02 | 2010-01-27 | 编码技术股份公司 | 利用复指数调制滤波器组的混叠减小 |
US20100228555A1 (en) * | 2005-09-16 | 2010-09-09 | Per Ekstrand | Partially Complex Modulated Filter Bank |
-
2015
- 2015-09-09 CN CN201510570781.2A patent/CN105281707B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101635564A (zh) * | 2001-04-02 | 2010-01-27 | 编码技术股份公司 | 利用复指数调制滤波器组的混叠减小 |
US20100228555A1 (en) * | 2005-09-16 | 2010-09-09 | Per Ekstrand | Partially Complex Modulated Filter Bank |
CN101303764A (zh) * | 2008-05-16 | 2008-11-12 | 西安电子科技大学 | 基于非下采样轮廓波的多传感器图像自适应融合方法 |
Non-Patent Citations (1)
Title |
---|
DFT调制滤波器组的设计算法研究;蒋俊正;《中国博士学位论文全文数据库》;20111215(第12期);第30-31、86-87页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105281707A (zh) | 2016-01-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105281707B (zh) | 一种动态可重构滤波器组的实现方法 | |
CN102318189B (zh) | 低延迟调制滤波器组 | |
CN107241107B (zh) | 一种数字信道化滤波器组实现方法 | |
CN101257482B (zh) | 数字基带可变速率转换调制***的实现方法和实现装置 | |
CN105337587B (zh) | 一种基于dft的非最大抽取***综合滤波器组构造方法 | |
CN107749762A (zh) | 一种多相滤波数字信道化实现方法 | |
US7620673B2 (en) | Complimentary discrete fourier transform processor | |
RU2012142732A (ru) | Устройство и способ обработки входного звукового сигнала с помощью каскадированного банка фильтров | |
CN101567701B (zh) | 一种高效多路数字下变频器*** | |
CN102355273B (zh) | 数字信道化方法及电路 | |
CN105279350B (zh) | 近似完全重构非均匀余弦调制滤波器组的设计方法 | |
Ambede et al. | Design and realization of variable digital filters for software-defined radio channelizers using an improved coefficient decimation method | |
Sakthivel et al. | Low complexity reconfigurable channelizers using non-uniform filter banks | |
Sudharman et al. | Design of power efficient variable bandwidth non-maximally decimated FRM filters for wideband channelizer | |
CN110210101A (zh) | 一种基于cem frm的动态非均匀窄过渡带滤波器组及设计方法 | |
CN107294512B (zh) | 一种基于树型结构的非均匀滤波器组滤波方法 | |
CN105811920B (zh) | 一种frm窄过渡带滤波器组结构 | |
CN104883157B (zh) | 一种可变子带数字滤波器 | |
CN110957996B (zh) | 一种基于abc算法的无乘法器frm滤波器组优化设计方法 | |
CN102891662A (zh) | 一种通用的速率下变换、上变换装置及方法 | |
CN115296647B (zh) | 基于frm技术的wola滤波器组及子带分割设计方法 | |
CN106712836B (zh) | 基于近似精确重构滤波器组的信道化器设计方法及装置 | |
Deep et al. | Spectral parameter approximation based tunable digital filters on Zynq SoC | |
CN109962731B (zh) | 近似完全重构星载低复杂度非均匀带宽数字信道化方法 | |
Amir et al. | Reconfigurable filter bank structures for low complexity digital channelizer using fractional interpolation and MFIR filters with cosine modulation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |