CN110210101B - 一种基于cem frm的动态非均匀窄过渡带滤波器组及设计方法 - Google Patents

一种基于cem frm的动态非均匀窄过渡带滤波器组及设计方法 Download PDF

Info

Publication number
CN110210101B
CN110210101B CN201910446571.0A CN201910446571A CN110210101B CN 110210101 B CN110210101 B CN 110210101B CN 201910446571 A CN201910446571 A CN 201910446571A CN 110210101 B CN110210101 B CN 110210101B
Authority
CN
China
Prior art keywords
filter
module
frm
cem
filter bank
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
CN201910446571.0A
Other languages
English (en)
Other versions
CN110210101A (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 CN201910446571.0A priority Critical patent/CN110210101B/zh
Publication of CN110210101A publication Critical patent/CN110210101A/zh
Application granted granted Critical
Publication of CN110210101B publication Critical patent/CN110210101B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design
    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00—Networks using digital techniques
    • H03H17/02—Frequency selective networks
    • H03H17/0211—Frequency selective networks using specific transformation algorithms, e.g. WALSH functions, Fermat transforms, Mersenne transforms, polynomial transforms, Hilbert transforms
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/38—Transceivers, i.e. devices in which transmitter and receiver form a structural unit and in which at least one part is used for functions of transmitting and receiving
    • H04B1/40—Circuits

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Algebra (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computing Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Signal Processing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
  • Noise Elimination (AREA)

Abstract

本发明涉及滤波器组技术领域，具体涉及一种基于CEM FRM的动态非均匀窄过渡带滤波器组及设计方法。本发明包括基于CEM FRM的分析滤波器组部分以及动态综合滤波器组部分。分析滤波器组部分包括延时模块、抽取模块、上支路延时模块、下支路原型滤波器模块、上支路屏蔽滤波器模块、下支路屏蔽滤波器模块、求和模块及IFFT模块。动态综合滤波器组部分包括：上支路屏蔽滤波器模块、下支路屏蔽滤波器模块、上支路延时模块、下支路原型滤波器模块、K倍插值模块、延时模块、求和模块及IFFT模块。设计方法包括基于CEMFRM的原型滤波器设计和屏蔽滤波器设计。本发明可以实现低复杂度的非均匀窄过渡带滤波器组，有效减少窄过渡带的非均匀滤波器组硬件资源消耗。

Description

一种基于CEM FRM的动态非均匀窄过渡带滤波器组及设计 方法
技术领域
本发明涉及滤波器组技术领域，具体涉及一种基于CEM FRM的动态非均匀窄过渡带滤波器组及设计方法。
背景技术
软件无线电被定义为是这样一种无线电，即其接收端的数字化是在天线后面的某一级，比如在宽带滤波、低噪声放大器和用来把射频信号下变频到中频的混频器及其中频放大器/滤波器等级连部件的后端进行的，对于发射机的数字化则正好相反。软件无线电刚开始是从军事通信领域提出来的，但经过了十几年的研究和推广，软件无线电已经应用于包括民用移动通信、雷达、电子战、测绘，甚至电视广播等无线电工程的各个领域。
由于多速率信号处理技术是软件无线电中的一个重要技术。在对射频信号直接进行带通采样之后，我们得到的数据流速率很高，导致后续的信号处理速度跟不上，尤其是在一些对实时性有要求的系统中。所以我们需要对A/D后的数据流进行降速处理，多速率信号处理技术为这种降速处理的实现提供了理论的依据。
滤波器组是多速率信号处理技术的核心。滤波器组是一组具有共同的输入或者共同的相加后的输出的滤波器。滤波器组用来实现对信号频率的各种分量进行分解，然后我们可以根据需要对各个子信号进行不同的处理和传输。根据滤波器组的结构，我们可以将其分为两类：一类是均匀滤波器组，另一类是非均匀滤波器组。现如今，随着软件无线电的广泛应用和多速率信号处理技术的发展，我们的研究重点已经从均匀滤波器组转变为非均匀滤波器组。
目前，设计非均匀滤波器组有两种方法：直接法和间接法。直接法是根据非均匀滤波器组中各个子滤波器的参数要求直接设计出所需要的子滤波器；间接法是利用分析滤波器组和综合滤波器组组合的方式来设计非均匀滤波器组。首先利用分析滤波器组设计出均匀滤波器组，然后利用综合滤波器组合并分析滤波器组中子滤波器从而实现非均匀滤波器组。
由于在利用滤波器组技术对信号进行处理的过程中，信号会产生误差失真，主要有相位、幅度、混叠这三种误差，为了去除这三种误差，使整个滤波器组系统做到精确或者近似精确重构，解决方案是设计特性、结构不同的滤波器组，所以这三种误差的消除也推动了滤波器组理论的不断发展。其中，采用间接法设计非均匀滤波器组时产生的混叠误差的主要原因是滤波器组存在过渡带。为了尽可能降低混叠误差，需要设计出足够窄的过渡带的滤波器组。但是，滤波器的过渡带越窄，滤波器组的复杂度就越高，工程实现就越困难，所以需要解决窄过渡带滤波器组的复杂度高的问题。
申请号为CN201610133236.1的专利，一种FRM窄过渡带滤波器组结构和文献《调制滤波器组及其设计方法研究》仅涉及基于FRM的均匀滤波器组，与本发明结构不同；申请号为CN201710244799.2的专利，一种5G系统中基于FRM技术的FBMC收发系统的设计方法，是首先利用FRM技术对原型滤波器进行设计，然后经过离散傅里叶变换调制成综合滤波器组和分析滤波器组，不是动态的非均匀滤波器组结构；文献《基于非均匀滤波器组的动态信道化滤波》仅仅是基于DFT滤波器组的设计方法研究，没有涉及FRM技术，与本发明设计方法不同。
发明内容
本发明的目的在于提供一种基于CEM FRM的动态非均匀窄过渡带滤波器组及设计方法，以降低窄过渡带滤波器组的结构设计复杂度，实现低复杂度的非均匀窄过渡带滤波器组。
第一方面，本发明实施例提供一种基于CEM FRM的动态非均匀窄过渡带滤波器组，包括：基于CEM FRM的分析滤波器组部分以及基于CEM FRM的动态综合滤波器组部分。
基于CEM FRM的分析滤波器组部分，包括：延时模块和D倍抽取模块、上支路延时模块
Figure GDA0002108378970000021
和下支路原型滤波器模块
Figure GDA0002108378970000022
上支路屏蔽滤波器模块
Figure GDA0002108378970000023
和下支路屏蔽滤波器模块
Figure GDA0002108378970000024
求和模块以及IFFT模块。
基于CEM FRM的动态综合滤波器组部分，包括：上支路屏蔽滤波器模块
Figure GDA0002108378970000025
和下支路屏蔽滤波器模块
Figure GDA0002108378970000026
上支路延时模块
Figure GDA0002108378970000027
和下支路原型滤波器模块
Figure GDA0002108378970000028
K倍插值模块和延时模块、求和模块以及IFFT模块。
第二方面，本发明实施例提供一种基于CEM FRM的动态非均匀窄过渡带滤波器组设计方法，包括：
步骤一：根据滤波器组的设计需求确定非均匀滤波器组的分析滤波器组的子滤波器个数M，抽取倍数D以及滤波器组中合成的基于CEM FRM的原型低通滤波器的通带截止频率ωp，阻带起始频率ωs，阻带衰减δa和插值倍数L；
步骤二：根据基于CEM FRM的原型低通滤波器的通带截止频率ωp，阻带起始频率ωs和阻带衰减δ设计合成基于CEM FRM原型低通滤波器中的原型半带滤波器，并对所述的原型半带滤波器进行复指数调制，得到调制后的原型半带滤波器的系数的虚部；
步骤三：根据基于CEM FRM的原型低通滤波器的通带截止频率ωp，阻带起始频率ωs和阻带衰减δ设计合成基于CEM FRM原型低通滤波器中的屏蔽滤波器，并对所述的屏蔽滤波器进行复指数调制，并对得到的调制后的屏蔽滤波器的系数的实部和虚部分别进行多相分解；
步骤四：分别对所述的调制后的原型半带滤波器的虚部以及调制后的屏蔽滤波器的系数的实部和虚部的多相形式进行抽取，根据基于CEM FRM的分析滤波器组结构的各个模块组成，得到基于CEM FRM的分析滤波器组；
步骤五：根据滤波器组的设计需求确定基于CEM FRM的分析滤波器组中需要合并的子滤波器个数K，对调制后的原型半带滤波器和屏蔽滤波器分别进行K倍抽取，得到抽取后的调制的原型半带滤波器的虚部，再对抽取后的调制的屏蔽滤波器的系数的实部和虚部分别进行多相分解，最后根据基于CEM FRM的动态综合滤波器组结构的各个模块组成，得到基于CEM FRM的动态综合滤波器组；
步骤六：将所述的基于CEM FRM的分析滤波器组和基于CEM FRM的动态综合滤波器组进行组合，得到基于CEM FRM的动态非均匀窄过渡带滤波器组；
本发明还包括这样一些结构特征：
所述步骤二，包括：
1)根据基于CEM FRM的滤波器的结构图，通过将通带截止频率ωp、阻带起始频率ωs和插值倍数L代入：
l＝ceil(ωpL/2π)
Figure GDA0002108378970000031
Figure GDA0002108378970000032
得到基于CEM FRM的滤波器中的原型半带滤波器Fa(z)的通带截止频率θ和阻带起始频率φ，上式中，l为正整数，ceil(ωpL/2π)表示大于或等于ωpL/2π的最小整数，由计算得到的通带截止频率θ、阻带起始频率φ和阻带衰减δa即可设计出原型半带滤波器Fa(z)；
2)令Na表示原型半带滤波器Fa(z)的阶数，对Fa(z)进行复指数调制，调制因子为ejn π/2，调制后得到的滤波器Fae(z)表达式为：
Figure GDA0002108378970000033
3)对滤波器Fae(z)进行L倍插值，得到插值后的F′ae(z)表达式为：
Figure GDA0002108378970000034
令F′ae(z)的单位采样响应为f′ae(n)，f′ae(n)为复数，则可表示为：
f′ae(n)＝f′ae，R(n)+f′ae，I(n)·j
那么，F′ae，I(z)的表达式为：
Figure GDA0002108378970000041
所述步骤三，包括：
1)屏蔽滤波器FMa(z)的通带截止频率ωpMa和阻带起始频率ωsMa分别为：
Figure GDA0002108378970000042
Figure GDA0002108378970000043
根据通带截止频率ωpMa，阻带起始频率ωsMa和阻带衰减δa即可设计出屏蔽滤波器FMa(z)；
2)令屏蔽滤波器FMa(z)的阶数为NMa，对FMa(z)进行复指数调制，调制因子为
Figure GDA0002108378970000044
调制后得到的滤波器FMae(z)表达式为：
Figure GDA0002108378970000045
当l为偶数时，屏蔽滤波器FMa(z)的调制因子
Figure GDA0002108378970000046
中的ω0的表达式为：
Figure GDA0002108378970000047
当l为奇数时，屏蔽滤波器FMa(z)的调制因子
Figure GDA0002108378970000048
中的ω0的表达式为：
Figure GDA0002108378970000049
令FMae(z)的单位采样响应为fMae(n)，fMae(n)为复数，则可以表示为：
fMae(n)＝fMae，R(n)+fMae，I(n)·j
那么，FMae，R(z)和FMae，I(z)的表达式分别为：
Figure GDA00021083789700000410
Figure GDA00021083789700000411
FMae，R(z)和FMae，I(z)的多相表达式分别为：
Figure GDA00021083789700000412
Figure GDA00021083789700000413
Figure GDA0002108378970000051
Figure GDA0002108378970000052
所述步骤四，包括：
1)设计延时和D倍抽取模块，根据输入信号的采样率和输出信号的采样率确定抽取模块中抽取倍数D，第n个子信道的延时模块为z-n+1；
2)设计上支路延时模块和下支路原型滤波器模块，第n个子信道的上支路延时模块中的延时为
Figure GDA0002108378970000053
下支路滤波器模块为2F′ae，I(z1/D)；
3)设计上支路屏蔽滤波器模块，下支路屏蔽滤波器模块及求和模块，第n个子信道的上支路屏蔽滤波器模块为FMa，R，n-1(zM/D)，下支路屏蔽滤波器模块为FMa，I，n-1(zM/D)，然后将下支路屏蔽滤波器模块的输出信号乘(-1)再与上支路屏蔽滤波器模块的输出信号相加；
4)设计IFFT模块，将第三步加和之后得到的输出信号作为M点IFFT模块的输入，得到的M个输出信号即为基于CEM FRM的分析滤波器组的M个输出信号；
所述步骤五，包括：
1)确定需要合并的基于CEM FRM的分析滤波器组的M个输出信号中的K个输出信号，其中，K≤M，K＝2n，n为正整数；
2)分别对滤波器FMae，I(z)进行M/K倍的抽取，抽取之后得到的滤波器为
Figure GDA0002108378970000054
则滤波器
Figure GDA0002108378970000055
和
Figure GDA0002108378970000056
的多相表达式分别为：
Figure GDA0002108378970000057
Figure GDA0002108378970000058
Figure GDA0002108378970000059
Figure GDA00021083789700000510
上式中，
Figure GDA00021083789700000511
和
Figure GDA00021083789700000512
分别是滤波器
Figure GDA00021083789700000513
和
Figure GDA00021083789700000514
的单位采样响应；
3)选择基于CEM FRM的分析滤波器组的M个输出信号中的K个输出信号作为基于CEM FRM的动态综合滤波器组的输入信号，其中，K≤M；
4)设计IFFT模块，对这K个输入信号之后进行K点IFFT变换；
5)设计上支路屏蔽滤波和下支路屏蔽滤波模块，将K点IFFT的K个输出信号中的每个子信号都分成两个支路，将每个子信号分别进行上支路屏蔽滤波和下支路屏蔽滤波，第k个子信道的上支路屏蔽滤波器模块为
Figure GDA0002108378970000061
下支路屏蔽滤波器模块为
Figure GDA0002108378970000062
6)设计上支路延时模块，下支路原型滤波器模块以及求和模块，第k个子信道的上支路和下支路屏蔽滤波器模块的输出信号将它们分别输入到上支路延时模块和下支路原型滤波器模块，其中，上支路延时模块的延时为
Figure GDA0002108378970000063
下支路的原型滤波器模块为2F′ae，I(z1/M)，下支路的原型滤波器模块的输出信号乘(-1)再与上支路延时模块的输出信号相加得到求和模块；
7)设计K倍插值模块和延时模块以及求和模块，将上述4)中得到的输出信号输入到插值模块中，插值倍数为K，最后输入到延时模块中，第k个子信道的延时模块的延时为z-k+1，最后将所有k个子信道的延时模块的输出进行加和，求得最后的输出信号y(n)；
本发明的有益效果在于：本发明利用基于CEM FRM的分析滤波器组对信号频谱进行均匀划分，然后利用基于CEM FRM的动态综合滤波器组对信号进行动态综合，可以实现低复杂度的非均匀窄过渡带滤波器组。
附图说明
图1为本发明的基于CEM FRM的动态非均匀窄过渡带滤波器组结构图；
图2为本发明的基于CEM FRM的分析滤波器组结构图；
图3为本发明的基于CEM FRM的动态综合滤波器组结构图；
图4为本发明的基于CEM FRM的原型滤波器幅频特性仿真图；
图5为本发明的基于CEM FRM的分析滤波器组幅频特性仿真图；
图6为本发明的基于CEM FRM的动态综合滤波器组幅频特性仿真图。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图对本发明做进一步描述：
本发明的技术方案是这样实现的：
一种基于CEM FRM的动态非均匀窄过渡带滤波器组结构，该结构包括基于CEM FRM的分析滤波器组结构和基于CEM FRM的动态综合滤波器组结构。
基于CEM FRM的分析滤波器组结构，该结构包含延时模块和D倍抽取模块，上支路延时模块
Figure GDA0002108378970000064
和下支路原型滤波器模块
Figure GDA0002108378970000065
上支路屏蔽滤波器模块
Figure GDA0002108378970000066
和下支路屏蔽滤波器模块
Figure GDA0002108378970000071
求和模块及IFFT模块组成。
基于CEM FRM的动态综合滤波器组结构，该结构包括IFFT模块，上支路屏蔽滤波器模块
Figure GDA0002108378970000072
和下支路屏蔽滤波器模块
Figure GDA0002108378970000073
上支路延时模块
Figure GDA0002108378970000074
和下支路原型滤波器模块
Figure GDA0002108378970000075
求和模块，K倍插值模块和延时模块以及求和模块组成。
基于CEM FRM的动态非均匀窄过渡带滤波器组的设计方法，包含如下步骤：
步骤1：确定基于CEM FRM的分析滤波器组的子滤波器个数M，抽取倍数D和滤波器组中合成的基于CEM FRM的原型低通滤波器的通带截止频率ωp，阻带起始频率ωs，阻带衰减δ和插值倍数L。
步骤2：根据基于CEM FRM的原型低通滤波器的通带截止频率ωp，阻带起始频率ωs和阻带衰减δ设计合成基于CEM FRM原型低通滤波器中的原型半带滤波器，并且对原型半带滤波器进行复指数调制，可以得到调制后的半带滤波器的系数的虚部。
步骤3：根据基于CEM FRM的原型低通滤波器的通带截止频率ωp，阻带起始频率ωs和阻带衰减δ设计合成基于CEM FRM原型低通滤波器中屏蔽滤波器，对屏蔽滤波器进行复指数调制，并且对调制后的屏蔽滤波器的系数的实部和虚部分别进行多相分解。
步骤4：分别对调制后的半带滤波器的虚部，调制后的屏蔽滤波器的系数的实部和虚部的多相形式进行抽取，并且根据图2中的基于CEM FRM的分析滤波器组结构的各个模块组成基于CEM FRM的分析滤波器组。
步骤5：确定基于CEM FRM的分析滤波器组中需要合并的子滤波器个数K，对调制后的基于CEM FRM原型低通滤波器中的原型半带滤波器和屏蔽滤波器分别进行K倍抽取。可以得到抽取后的调制的基于CEM FRM原型低通滤波器中的原型半带滤波器的虚部。然后再对抽取后的调制的屏蔽滤波器的系数的实部和虚部分别进行多相分解。最后根据图3中基于CEM FRM的动态综合滤波器组结构的各个模块组成基于CEM FRM的动态综合滤波器组。
步骤6：最后，将基于CEM FRM的分析滤波器组和基于CEM FRM的动态综合滤波器组进行组合，即可得到基于CEM FRM的动态非均匀窄过渡带滤波器组。
更为具体的说，本发明的一种基于CEM FRM的动态非均匀窄过渡带滤波器组设计方法，包括：
步骤1：根据滤波器组的设计需求确定非均匀滤波器组的分析滤波器组的子滤波器个数M，抽取倍数D以及滤波器组中合成的基于CEM FRM的原型低通滤波器的通带截止频率ωp，阻带起始频率ωs，阻带衰减δa和插值倍数L。
步骤2：
第一步：图3为基于CEM FRM的滤波器结构图，根据ωp，ωs和L由式(1)～(3)得到基于CEM FRM的滤波器中的原型半带滤波器Fa(z)的通带截止频率θ和阻带起始频率φ：
Figure GDA0002108378970000081
Figure GDA0002108378970000082
Figure GDA0002108378970000083
其中，l为正整数，ceil(ωpL/2π)表示大于或等于ωpL/2π的最小整数。由计算得到的通带截止频率θ，阻带起始频率φ和阻带衰减δa即可设计出原型半带滤波器Fa(z)。
第二步：令Na表示原型半带滤波器Fa(z)的阶数，对Fa(z)进行复指数调制，调制因子为ejnπ/2，调制后得到的滤波器Fae(z)表达式为：
Figure GDA0002108378970000084
第三步：对滤波器Fae(z)进行L倍插值，可以得到插值后的F′ae(z)表达式：
Figure GDA0002108378970000085
令F′ae(z)的单位采样响应为f′ae(n)。f′ae(n)为复数，那么它可以表示为：
f′ae(n)＝f′ae，R(n)+f′ae，I(n)·j#(6)
那么，F′ae，I(z)的表达式为：
Figure GDA0002108378970000086
步骤3：
第一步：屏蔽滤波器FMa(z)的通带截止频率ωpMa和阻带起始频率ωsMa分别为：
Figure GDA0002108378970000087
Figure GDA0002108378970000088
根据通带截止频率ωpMa，阻带起始频率ωsMa和阻带衰减δa即可设计出屏蔽滤波器FMa(z)。
第二步：令屏蔽滤波器FMa(z)的阶数为NMa，对FMa(z)进行复指数调制，调制因子为
Figure GDA0002108378970000089
调制后得到的滤波器FMae(z)表达式为：
Figure GDA0002108378970000091
当l为偶数时，屏蔽滤波器FMa(z)的调制因子
Figure GDA0002108378970000092
中的ω0的表达式为：
Figure GDA0002108378970000093
当l为奇数时，屏蔽滤波器FMa(z)的调制因子
Figure GDA0002108378970000094
中的ω0的表达式为：
Figure GDA0002108378970000095
令FMae(z)的单位采样响应为fMae(n)。fMae(n)为复数，那么它可以表示为：
fMae(n)＝fMae，R(n)+fMae，I(n)·j#(13)
那么，FMae，R(z)和FMae，I(z)的表达式分别为：
Figure GDA0002108378970000096
Figure GDA0002108378970000097
FMae，R(z)和FMae，I(z)的多相表达式分别为：
Figure GDA0002108378970000098
Figure GDA0002108378970000099
Figure GDA00021083789700000910
Figure GDA00021083789700000911
步骤4：
第一步：设计延时和D倍抽取模块。根据输入信号的采样率和输出信号的采样率确定抽取模块中抽取倍数D，第n个子信道的延时模块为z-n+1。
第二步：设计上支路延时模块和下支路原型滤波器模块。第n个子信道的上支路延时模块中的延时为
Figure GDA00021083789700000912
下支路滤波器模块为2F′ae，I(z1/D)。
第三步：设计上支路屏蔽滤波器模块，下支路屏蔽滤波器模块及求和模块。第n个子信道的上支路屏蔽滤波器模块为FMa，R，n-1(zM/D)，下支路屏蔽滤波器模块为FMa，I，n-1(zM/D)，然后将下支路屏蔽滤波器模块的输出信号乘-1再与上支路屏蔽滤波器模块的输出信号相加。
第四步：设计IFFT模块。将第三步加和之后得到的输出信号作为M点IFFT模块的输入，得到的M个输出信号即为基于CEM FRM的分析滤波器组的M个输出信号。
步骤5：
第一步：确定需要合并的基于CEM FRM的分析滤波器组的M个输出信号中的K个输出信号(K≤M，K＝2n，n为正整数)。
第二步：分别对滤波器FMae，I(z)进行M/K倍的抽取，抽取之后得到的滤波器为
Figure GDA0002108378970000101
那么，滤波器
Figure GDA0002108378970000102
和
Figure GDA0002108378970000103
的多相表达式分别为：
Figure GDA0002108378970000104
Figure GDA0002108378970000105
Figure GDA0002108378970000106
Figure GDA0002108378970000107
其中，
Figure GDA0002108378970000108
和
Figure GDA0002108378970000109
分别是滤波器
Figure GDA00021083789700001010
和
Figure GDA00021083789700001011
的单位采样响应。
第三步：选择基于CEM FRM的分析滤波器组的M个输出信号中的K个输出信号(K≤M)作为基于CEM FRM的动态综合滤波器组的输入信号。
第四步：设计IFFT模块。对这K个输入信号之后进行K点IFFT变换。
第五步：设计上支路屏蔽滤波和下支路屏蔽滤波模块。将K点IFFT的K个输出信号中的每个子信号都分成两个支路，将每个子信号分别进行上支路屏蔽滤波和下支路屏蔽滤波，第k个子信道的上支路屏蔽滤波器模块为
Figure GDA00021083789700001012
下支路屏蔽滤波器模块为
Figure GDA00021083789700001013
第六步：设计上支路延时模块，下支路原型滤波器模块以及求和模块。第k个子信道的上支路和下支路屏蔽滤波器模块的输出信号将它们分别输入到上支路延时模块和下支路原型滤波器模块，其中，上支路延时模块的延时为
Figure GDA00021083789700001014
下支路的原型滤波器模块为2F′ae，I(z1/M)，下支路的原型滤波器模块的输出信号乘-1再与上支路延时模块的输出信号相加得到求和模块。
第七步：设计K倍插值模块和延时模块以及求和模块。将第四步得到的输出信号输入到插值模块中，插值倍数为K。最后输入到延时模块中，第k个子信道的延时模块的延时为z-k+1，最后将所有k个子信道的延时模块的输出进行加和，即可求得最后的输出信号y(n)。
为了验证本发明的有效性，根据上述设计方法进行了仿真实验。仿真系统的抽样率设置为960MHz，首先设计非均匀滤波器组中的分析滤波器组部分，分析滤波器组的子滤波器个数为16，抽取倍数也为16(临界抽取)，基于CEM FRM的原型低通滤波器的归一化通带截止频率和阻带截止频率分别设置为0.0615和0.0635，阻带衰减设置为60dB，插值倍数L设置为32。首先，根据式(1)～(3)可以得到基于CEM FRM的原型低通滤波器的原型半带滤波器Fa(z)，它的通带截止频率和阻带起始频率分别设置为0.468和0.532，阻带衰减为60dB，滤波器阶数为102。然后，根据式(8)～(9)可以得到屏蔽滤波器FMa(z)的归一化通带截止频率和阻带起始频率分别为0.047875和0.1046875，滤波器阶数为384，调制因子f0为-0.015625(f0＝ω0/2π)。最终合成的FRM滤波器H(z)的幅频特性图如图4所示。然后根据步骤2～4，我们可以得到基于CEM FRM的分析滤波器组，其幅频特性图如图5所示。根据步骤5，我们分别在K为2和4时对动态综合滤波器组进行仿真，得到的基于CEM FRM的动态综合滤波器组幅频特性图如图6所示。
本发明中的一种基于复指数调制的频率响应屏蔽技术(CEM FRM)的动态非均匀窄过渡带滤波器组的结构，该结构包括基于CEM FRM的分析滤波器组结构和基于CEM FRM的动态综合滤波器组结构。基于CEM FRM的分析滤波器组结构，该结构包含延时和抽取模块，上支路延时模块，下支路原型滤波器模块，上支路屏蔽滤波器模块，下支路屏蔽滤波器模块，求和模块和IFFT模块。基于CEM FRM的动态综合滤波器组结构，该结构包括IFFT模块，上支路屏蔽滤波器模块，下支路屏蔽滤波器模块，上支路延时模块，下支路原型滤波器模块，求和模块以及延时和插值模块。
一种基于CEM FRM的动态非均匀窄过渡带滤波器组的结构，所述的基于CEM FRM的分析滤波器组结构设计方法如下：
步骤1：确定基于CEM FRM的分析滤波器组的子滤波器个数M，抽取倍数D和滤波器组中合成的基于CEM FRM的原型低通滤波器的通带截止频率ωp，阻带起始频率ωs，阻带衰减δ和插值倍数L。
步骤2：根据基于CEM FRM的原型低通滤波器的通带截止频率ωp，阻带起始频率ωs和阻带衰减δ设计合成基于CEM FRM原型低通滤波器中的原型半带滤波器，并且对原型半带滤波器进行复指数调制，可以得到调制后的半带滤波器的系数的虚部。
步骤3：根据基于CEM FRM的原型低通滤波器的通带截止频率ωp，阻带起始频率ωs和阻带衰减δ设计合成基于CEM FRM原型低通滤波器中屏蔽滤波器，对屏蔽滤波器进行复指数调制，并且对调制后的屏蔽滤波器的系数的实部和虚部分别进行多相分解。
步骤4：分别对调制后的半带滤波器的虚部，调制后的屏蔽滤波器的系数的实部和虚部的多相形式进行抽取，并且根据图2中的基于CEM FRM的分析滤波器组结构的各个模块组成基于CEM FRM的分析滤波器组。
一种基于CEM FRM的动态非均匀窄过渡带滤波器组的结构，所述的基于CEM FRM的动态综合滤波器组结构设计方法如下：
步骤1：确定基于CEM FRM的分析滤波器组的子滤波器个数M，抽取倍数D和滤波器组中合成的基于CEM FRM的原型低通滤波器的通带截止频率ωp，阻带起始频率ωs，阻带衰减δ和插值倍数L。
步骤2：根据基于CEM FRM的原型低通滤波器的通带截止频率ωp，阻带起始频率ωs和阻带衰减δ设计合成基于CEM FRM原型低通滤波器中的原型半带滤波器，并且对原型半带滤波器进行复指数调制，可以得到调制后的半带滤波器的系数的虚部。
步骤3：根据基于CEM FRM的原型低通滤波器的通带截止频率ωp，阻带起始频率ωs和阻带衰减δ设计合成基于CEM FRM原型低通滤波器中屏蔽滤波器，对屏蔽滤波器进行复指数调制，并且对调制后的屏蔽滤波器的系数的实部和虚部分别进行多相分解。
步骤4：确定基于CEM FRM的分析滤波器组中需要合并的子滤波器个数K，对调制后的基于CEM FRM原型低通滤波器中的原型半带滤波器和屏蔽滤波器分别进行K倍抽取。可以得到抽取后的调制的基于CEM FRM原型低通滤波器中的原型半带滤波器的虚部。然后再对抽取后的调制的屏蔽滤波器的系数的实部和虚部分别进行多相分解。最后根据基于CEMFRM的动态综合滤波器组结构中的各个模块组成基于CEM FRM的动态综合滤波器组。
步骤5：将基于CEM FRM的分析滤波器组和基于CEM FRM的动态综合滤波器组进行组合，即可得到基于CEM FRM的动态非均匀窄过渡带滤波器组。
综上，本发明涉及滤波器组技术领域，具体涉及一种基于CEM FRM的动态非均匀窄过渡带滤波器组及设计方法。本发明包括基于CEM FRM的分析滤波器组部分以及动态综合滤波器组部分。分析滤波器组部分包括延时模块、抽取模块、上支路延时模块、下支路原型滤波器模块、上支路屏蔽滤波器模块、下支路屏蔽滤波器模块、求和模块及IFFT模块。动态综合滤波器组部分包括：上支路屏蔽滤波器模块、下支路屏蔽滤波器模块、上支路延时模块、下支路原型滤波器模块、K倍插值模块、延时模块、求和模块及IFFT模块。设计方法包括基于CEM FRM的原型滤波器设计和屏蔽滤波器设计。本发明可以实现低复杂度的非均匀窄过渡带滤波器组，有效减少窄过渡带的非均匀滤波器组硬件资源消耗。

Claims (6)

1.一种基于CEM FRM的动态非均匀窄过渡带滤波器组，其特征在于，包括：基于CEM FRM的分析滤波器组部分以及基于CEM FRM的动态综合滤波器组部分；
基于CEM FRM的分析滤波器组部分，包括：延时模块和D倍抽取模块、上支路延时模块
Figure FDA0002073816470000011
和下支路原型滤波器模块
Figure FDA0002073816470000012
上支路屏蔽滤波器模块
Figure FDA0002073816470000013
和下支路屏蔽滤波器模块
Figure FDA0002073816470000014
求和模块以及IFFT模块；
基于CEM FRM的动态综合滤波器组部分，包括：上支路屏蔽滤波器模块
Figure FDA0002073816470000015
和下支路屏蔽滤波器模块
Figure FDA0002073816470000016
上支路延时模块
Figure FDA0002073816470000017
和下支路原型滤波器模块
Figure FDA0002073816470000018
K倍插值模块和延时模块、求和模块以及IFFT模块。
2.一种基于CEM FRM的动态非均匀窄过渡带滤波器组设计方法，包括权利要求1所述的基于CEM FRM的动态非均匀窄过渡带滤波器组，其特征在于，步骤如下：
步骤一：根据滤波器组的设计需求确定非均匀滤波器组的分析滤波器组的子滤波器个数M，抽取倍数D以及滤波器组中合成的基于CEM FRM的原型低通滤波器的通带截止频率ωp，阻带起始频率ωs，阻带衰减δa和插值倍数L；
步骤二：根据基于CEM FRM的原型低通滤波器的通带截止频率ωp，阻带起始频率ωs和阻带衰减δ设计合成基于CEM FRM原型低通滤波器中的原型半带滤波器，并对所述的原型半带滤波器进行复指数调制，得到调制后的原型半带滤波器的系数的虚部；
步骤三：根据基于CEM FRM的原型低通滤波器的通带截止频率ωp，阻带起始频率ωs和阻带衰减δ设计合成基于CEM FRM原型低通滤波器中的屏蔽滤波器，并对所述的屏蔽滤波器进行复指数调制，并对得到的调制后的屏蔽滤波器的系数的实部和虚部分别进行多相分解；
步骤四：分别对所述的调制后的原型半带滤波器的虚部以及调制后的屏蔽滤波器的系数的实部和虚部的多相形式进行抽取，根据基于CEM FRM的分析滤波器组结构的各个模块组成，得到基于CEM FRM的分析滤波器组；
步骤五：根据滤波器组的设计需求确定基于CEM FRM的分析滤波器组中需要合并的子滤波器个数K，对调制后的原型半带滤波器和屏蔽滤波器分别进行K倍抽取，得到抽取后的调制的原型半带滤波器的虚部，再对抽取后的调制的屏蔽滤波器的系数的实部和虚部分别进行多相分解，最后根据基于CEM FRM的动态综合滤波器组结构的各个模块组成，得到基于CEM FRM的动态综合滤波器组；
步骤六：将所述的基于CEM FRM的分析滤波器组和基于CEM FRM的动态综合滤波器组进行组合，得到基于CEM FRM的动态非均匀窄过渡带滤波器组。
3.根据权利要求2所述的一种基于CEM FRM的动态非均匀窄过渡带滤波器组设计方法，其特征在于，所述步骤二，包括：
1)根据基于CEM FRM的滤波器的结构图，通过将通带截止频率ωp、阻带起始频率ωs和插值倍数L代入：
l＝ceil(ωpL/2π)
Figure FDA0002073816470000021
Figure FDA0002073816470000022
得到基于CEM FRM的滤波器中的原型半带滤波器Fa(z)的通带截止频率θ和阻带起始频率φ，上式中，l为正整数，ceil(ωpL/2π)表示大于或等于ωpL/2π的最小整数，由计算得到的通带截止频率θ、阻带起始频率φ和阻带衰减δa即可设计出原型半带滤波器Fa(z)；
2)令Na表示原型半带滤波器Fa(z)的阶数，对Fa(z)进行复指数调制，调制因子为ejnπ/2，调制后得到的滤波器Fae(z)表达式为：
Figure FDA0002073816470000023
3)对滤波器Fae(z)进行L倍插值，得到插值后的F′ae(z)表达式为：
Figure FDA0002073816470000024
令F′ae(z)的单位采样响应为f′ae(n)，f′ae(n)为复数，则可表示为：
f′ae(n)＝f′ae，R(n)+f′ae，I(n)·j
那么，F′ae，I(z)的表达式为：
Figure FDA0002073816470000025
4.根据权利要求3所述的一种基于CEM FRM的动态非均匀窄过渡带滤波器组设计方法，其特征在于：所述步骤三，包括：
1)屏蔽滤波器FMa(z)的通带截止频率ωpMa和阻带起始频率ωsMa分别为：
Figure FDA0002073816470000026
Figure FDA0002073816470000027
根据通带截止频率ωpMa，阻带起始频率ωsMa和阻带衰减δa即可设计出屏蔽滤波器FMa(z)；
2)令屏蔽滤波器FMa(z)的阶数为NMa，对FMa(z)进行复指数调制，调制因子为
Figure FDA0002073816470000031
调制后得到的滤波器FMae(z)表达式为：
Figure FDA0002073816470000032
当l为偶数时，屏蔽滤波器FMa(z)的调制因子
Figure FDA0002073816470000033
中的ω0的表达式为：
Figure FDA0002073816470000034
当l为奇数时，屏蔽滤波器FMa(z)的调制因子
Figure FDA0002073816470000035
中的ω0的表达式为：
Figure FDA0002073816470000036
令FMae(z)的单位采样响应为fMae(n)，fMae(n)为复数，则可以表示为：
fMae(n)＝fMae，R(n)+fMae，I(n)·j
那么，FMae，R(z)和FMae，I(z)的表达式分别为：
Figure FDA0002073816470000037
Figure FDA0002073816470000038
FMae，R(z)和FMae，I(z)的多相表达式分别为：
Figure FDA0002073816470000039
Figure FDA00020738164700000310
Figure FDA00020738164700000311
Figure FDA00020738164700000312
5.根据权利要求4所述的一种基于CEM FRM的动态非均匀窄过渡带滤波器组设计方法，其特征在于：所述步骤四，包括：
1)设计延时和D倍抽取模块，根据输入信号的采样率和输出信号的采样率确定抽取模块中抽取倍数D，第n个子信道的延时模块为z-n+1；
2)设计上支路延时模块和下支路原型滤波器模块，第n个子信道的上支路延时模块中的延时为
Figure FDA0002073816470000041
下支路滤波器模块为2F′ae，I(z1/D)；
3)设计上支路屏蔽滤波器模块，下支路屏蔽滤波器模块及求和模块，第n个子信道的上支路屏蔽滤波器模块为FMa，R，n-1(zM/D)，下支路屏蔽滤波器模块为FMa，I，n-1(zM/D)，然后将下支路屏蔽滤波器模块的输出信号乘(-1)再与上支路屏蔽滤波器模块的输出信号相加；
4)设计IFFT模块，将第三步加和之后得到的输出信号作为M点IFFT模块的输入，得到的M个输出信号即为基于CEM FRM的分析滤波器组的M个输出信号。
6.根据权利要求5所述的一种基于CEM FRM的动态非均匀窄过渡带滤波器组设计方法，其特征在于：所述步骤五，包括：
1)确定需要合并的基于CEM FRM的分析滤波器组的M个输出信号中的K个输出信号，其中，K≤M，K＝2n，n为正整数；
2)分别对滤波器FMae，I(z)进行M/K倍的抽取，抽取之后得到的滤波器为
Figure FDA0002073816470000042
则滤波器
Figure FDA0002073816470000043
和
Figure FDA0002073816470000044
的多相表达式分别为：
Figure FDA0002073816470000045
Figure FDA0002073816470000046
Figure FDA0002073816470000047
Figure FDA0002073816470000048
上式中，
Figure FDA0002073816470000049
和
Figure FDA00020738164700000410
分别是滤波器
Figure FDA00020738164700000411
和
Figure FDA00020738164700000412
的单位采样响应；
3)选择基于CEM FRM的分析滤波器组的M个输出信号中的K个输出信号作为基于CEMFRM的动态综合滤波器组的输入信号，其中，K≤M；
4)设计IFFT模块，对这K个输入信号之后进行K点IFFT变换；
5)设计上支路屏蔽滤波和下支路屏蔽滤波模块，将K点IFFT的K个输出信号中的每个子信号都分成两个支路，将每个子信号分别进行上支路屏蔽滤波和下支路屏蔽滤波，第k个子信道的上支路屏蔽滤波器模块为
Figure FDA00020738164700000413
下支路屏蔽滤波器模块为
Figure FDA00020738164700000414
6)设计上支路延时模块，下支路原型滤波器模块以及求和模块，第k个子信道的上支路和下支路屏蔽滤波器模块的输出信号将它们分别输入到上支路延时模块和下支路原型滤波器模块，其中，上支路延时模块的延时为
Figure FDA0002073816470000051
下支路的原型滤波器模块为2F′ae，I(z1/M)，下支路的原型滤波器模块的输出信号乘(-1)再与上支路延时模块的输出信号相加得到求和模块；
7)设计K倍插值模块和延时模块以及求和模块，将上述4)中得到的输出信号输入到插值模块中，插值倍数为K，最后输入到延时模块中，第k个子信道的延时模块的延时为z-k+1，最后将所有k个子信道的延时模块的输出进行加和，求得最后的输出信号y(n)。
CN201910446571.0A 2019-05-27 2019-05-27 一种基于cem frm的动态非均匀窄过渡带滤波器组及设计方法 Active CN110210101B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910446571.0A CN110210101B (zh) 2019-05-27 2019-05-27 一种基于cem frm的动态非均匀窄过渡带滤波器组及设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910446571.0A CN110210101B (zh) 2019-05-27 2019-05-27 一种基于cem frm的动态非均匀窄过渡带滤波器组及设计方法

Publications (2)

Publication Number Publication Date
CN110210101A CN110210101A (zh) 2019-09-06
CN110210101B true CN110210101B (zh) 2022-08-02

Family

ID=67788883

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910446571.0A Active CN110210101B (zh) 2019-05-27 2019-05-27 一种基于cem frm的动态非均匀窄过渡带滤波器组及设计方法

Country Status (1)

Country Link
CN (1) CN110210101B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111010146B (zh) * 2019-12-13 2023-08-01 哈尔滨工程大学 一种基于快速滤波器组的信号重构结构及其设计方法
CN110957996B (zh) * 2019-12-13 2023-08-01 哈尔滨工程大学 一种基于abc算法的无乘法器frm滤波器组优化设计方法
CN113992190B (zh) * 2021-12-23 2022-04-01 南京先锋硕通无线技术有限公司 一种用于dvb-rcs2的双层滤波器组设计方法
CN115296647B (zh) * 2022-07-08 2023-04-14 中国人民解放军战略支援部队航天工程大学 基于frm技术的wola滤波器组及子带分割设计方法
CN115955379B (zh) * 2022-12-25 2024-06-28 哈尔滨工程大学 一种多尺度可配置的窄过渡带信道化器低复杂度实现方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102222221A (zh) * 2011-07-01 2011-10-19 哈尔滨工程大学 基于多光谱图像滤波输出统计分布规律的目标检测方法
CN202475373U (zh) * 2010-12-03 2012-10-03 约翰·梅扎林瓜联合有限公司 滤波器电路组件、同轴电缆接地块组件和连接器组件
CN104685331A (zh) * 2012-10-05 2015-06-03 法国大陆汽车公司 用于位于内燃机中的压电传感器的经改进的测量放大器电路
CN105811920A (zh) * 2016-03-09 2016-07-27 哈尔滨工程大学 一种frm窄过渡带滤波器组结构
CN107070836A (zh) * 2017-04-14 2017-08-18 蚌埠学院 一种5g系统中基于frm技术的fbmc收发系统的设计方法
CN108919202A (zh) * 2018-08-10 2018-11-30 杭州电子科技大学 一种基于认知机制的非均匀动态滤波器组的构建方法
CN109525256A (zh) * 2018-10-18 2019-03-26 哈尔滨工程大学 一种基于fpga的窄过渡带滤波器组的信道化发射结构

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7108192B2 (en) * 1999-09-17 2006-09-19 Silverbrook Research Pty Ltd Rotationally symmetric tags
KR100483023B1 (ko) * 2002-12-24 2005-04-15 한국전자통신연구원 광전송 시스템의 편광 모드 분산 보상 장치 및 그 보상 방법
US7447284B2 (en) * 2003-03-28 2008-11-04 Freescale Semiconductor, Inc. Method and apparatus for signal noise control
KR101547210B1 (ko) * 2013-12-05 2015-08-25 한국기초과학지원연구원 냉전자 소스원을 이용한 이온트랩 질량분석기

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN202475373U (zh) * 2010-12-03 2012-10-03 约翰·梅扎林瓜联合有限公司 滤波器电路组件、同轴电缆接地块组件和连接器组件
CN102222221A (zh) * 2011-07-01 2011-10-19 哈尔滨工程大学 基于多光谱图像滤波输出统计分布规律的目标检测方法
CN104685331A (zh) * 2012-10-05 2015-06-03 法国大陆汽车公司 用于位于内燃机中的压电传感器的经改进的测量放大器电路
CN105811920A (zh) * 2016-03-09 2016-07-27 哈尔滨工程大学 一种frm窄过渡带滤波器组结构
CN107070836A (zh) * 2017-04-14 2017-08-18 蚌埠学院 一种5g系统中基于frm技术的fbmc收发系统的设计方法
CN108919202A (zh) * 2018-08-10 2018-11-30 杭州电子科技大学 一种基于认知机制的非均匀动态滤波器组的构建方法
CN109525256A (zh) * 2018-10-18 2019-03-26 哈尔滨工程大学 一种基于fpga的窄过渡带滤波器组的信道化发射结构

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
A Non-Maximally Decimated Dynamic Reconfigurable Channelized Structure Based on Modulated Filter Bank;Wenxu Zhang等;《Int. J. Communications, Network and System Sciences》;20170814(第10期);第88-97页 *
A novel channelized FB architecture with narrow transition bandwidth based on CEM FRM;Tao Chen等;《Ann. Telecommun》;20150821;第71卷;第27-33页 *
一种基于CEM FRM技术的滤波器组设计新方法;张立志等;《信息工程大学学报》;20090615;第10卷(第2期);第177-180页 *

Also Published As

Publication number Publication date
CN110210101A (zh) 2019-09-06

Similar Documents

Publication Publication Date Title
CN110210101B (zh) 一种基于cem frm的动态非均匀窄过渡带滤波器组及设计方法
US7620673B2 (en) Complimentary discrete fourier transform processor
Chen et al. Non-maximally decimated analysis/synthesis filter banks: Applications in wideband digital filtering
CN105337587B (zh) 一种基于dft的非最大抽取系统综合滤波器组构造方法
CN109525256B (zh) 一种基于fpga的窄过渡带滤波器组的信道化发射结构
CN101881821B (zh) 一种分数阶傅里叶域信道化接收方法
CN107748354B (zh) 基于分析与综合的宽带数字波束形成装置
Sommen et al. On the relationship between uniform and recurrent nonuniform discrete-time sampling schemes
CN109962731B (zh) 近似完全重构星载低复杂度非均匀带宽数字信道化方法
Zhang et al. Design and FPGA implementation of a novel efficient FRM-based channelized receiver structure
CN115296647B (zh) 基于frm技术的wola滤波器组及子带分割设计方法
CN115955379B (zh) 一种多尺度可配置的窄过渡带信道化器低复杂度实现方法
Kuenzle et al. 3-D IIR filtering using decimated DFT-polyphase filter bank structures
Ogale et al. Cosine Modulated Non-Uniform Filter Banks.
US9306606B2 (en) Nonlinear filtering using polyphase filter banks
CN114966564A (zh) 一种宽带线性调频信号多通道并行去斜方法
Tang et al. Design and simulation of DRFM system based on digital channelized receiver and transmitter
Varghese et al. Design of computationally efficient ifir based filter structure for digital channelizer
Harris et al. Cascade of perfect reconstruction analysis and synthesis filter banks: The new architecture of next generation wideband receivers
Cecchi et al. A new approach to digital audio equalization
Vashkevich et al. Practical design of multi-channel oversampled warped cosine-modulated filter banks
Jiang et al. An improved wideband digital baseband conversion method
CN113507279B (zh) 一种高精度星间距离观测数据的降速率滤波方法
Lei et al. An Efficient Weighted Overlap-add Structure for Wideband Frequency Spectrum Analysis with Narrow Resolution Bandwidth
WO2014064408A2 (en) Method and apparatus for processing a signal

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