一种可变子带数字滤波器的制造方法

一种可变子带数字滤波器的制造方法
【技术领域】
[0001] 本发明涉及数字滤波器领域，特别是一种可变子带数字滤波器。
【背景技术】
[0002] 移动通信领域可用频带一般为受限资源，需要申请或者购买才能获得相应的运营 资格。在同一个制式内，频率资源也被分配给多个不同的运营商。无论中国还是国外，一家 运营商往往在一个制式上获得多个不连续频带的运营许可。例如GSM的890~893MHz(3MHz 带宽）和900~901MHz(1MHz带宽）二个子频带。不同国家有不同的无线频段的规划，同时 一个运营商不能获得整个频段的经营权。特别在国外，不同运营商的频段带宽不一样。此 外，一个制式里的多个不连续的频段的带宽不是固定的。所以为了适应现代无线移动网络 的运营特点，通信系统如基站、无线直放站、移动终端及手机需要具有多子带选频，带宽可 以实时设置的功能，以满足现场应用要求。模拟技术可以解决多网络融合和多子带，但不能 实现"带宽实时可设置"。
[0003] 利用数字信号处理技术灵活、精确的特点，可以有效地解决"带宽实时设置"的技 术难题。在数字信号处理技术中，数模转换器首先将连续模拟信号离散、量化和编码为数字 信号。然后，离散的数字信号经过数学算法的处理之后通过数模转换又还原为模拟信号。 数字信号处理的算法可以通过计算机或者数字信号处理器OSP)和专用集成电路（ASIC)， FPGA等以编程方式实现。然而，目前对实时可变子带滤波器的研发非常有限，还没有针对性 的完整的设计和实现方法。

【发明内容】

[0004] 本发明的主要目的在于提出一种高效的可变子带滤波器，有效实现可变子带滤波 器的同时，可以降低对硬件资源的占用。
[0005] 本发明采用如下技术方案：
[0006] -种可变子带数字滤波器，用于将[0, 2 31]频带分解为M个均匀子带进行选通滤 波，包括滤波输入端和滤波输出端，其特征在于：其结构包括直接型FIR滤波器或转置型 FIR滤波器。
[0007] 优选的，直接型FIR滤波器实现的结构如下，包括滤波输入端、滤波输出端、M-1个 单位延迟单元、M-1个合成累加器、M个可变系数乘法器和M个多相滤波器；该滤波输入端 一路连接第〇相多相滤波器输入端，另一路经M-1个单位延迟单元依次延迟构成M-1级延 迟支路；该M-1级延迟支路分别连接第一至第M-1相多相滤波器输入端；M个多相滤波器的 输出端分别连接到M个可变系数乘法器的一个输入端，而M个可变系数乘法器的另一输入 端分别连接M个可变系数Sp其中i= 0、1、2…M-1 ;该第0至第M-2个可变系数乘法器的 输出端分别连接第一至第M-1个合成累加器的一输入端，该第M-1个可变系数乘法器输出 端连接第M-1合成累加器另一输入端，第M-2合成累加器依次累加至第一合成累加器，该第 一合成累加器的输出端作为滤波输出端。
[0008] 优选的，转置型FIR滤波器实现的结构如下，包括M-1个单位延迟单元、M-1个合 成累加器、M个可变系数乘法器和M多相滤波器；该滤波输入端分成M支路，该M支路分别 连接第〇至第M-1个多相滤波器输入端；M个多相滤波器的输出端分别连接到M个可变系 数乘法器的一个输入端，而M个可变系数乘法器的另一输入端为M个可变系数Sp其中i= 0、1、2…M-1;该第0至第M-2个可变系数乘法器输出端分别连接第一至第M-1合成累加器 输入端，该第M-1个可变系数乘法器输出端经第M-1延迟单元连接第M-1合成累加器另一 输入端，该第M-2合成累加器至第一合成累加器的输出端分别经第M-2延迟单元至第一延 迟单元连接至上一级合成累加器输入端，该第一合成累加器输出端作为滤波输出端。
[0009] 优选的，所述M个多相滤波器为直接型FIR滤波器。
[0010] 优选的，所述多相滤波器包括子滤波输入端，子滤波输出端、N个常系数乘法器、 N-1个M延迟单元和N-1个加法器；该子滤波输入端一路连接第一个常系数乘法器输入端， 另一路经N-1个M延迟单元依次延迟构成N-1级延迟支路；该N-1级延迟支路分别连接第 二个至第N个常系数乘法器输入端，第一至第N个常系数乘法器的另一个输入端分别连接 滤波器系数b至hh，该第一至第N-1个常系数乘法器的输出端分别连接第一至第N-1合 加法器输入端，该第N个常系数乘法器的输出端连接第N-1加法器的另一输入端，第N-1加 法器依次累加至第一加法器，该第一加法器的输出端作为子滤波输出端。
[0011] 优选的，所述多相滤波器为转置型FIR滤波器。
[0012] 优选的，所述多相滤波器包括子滤波输入端、子滤波输出端、N个常系数乘法器、 N-1个M延迟单元和N-1个加法器；该子滤波输入端分成N支路，该N支路分别连接第一至 第N个常系数乘法器的一个输入端，该N个常系数乘法器的另一个输入端分别连接滤波器 系数&至hh，该第一至第N-1个常系数乘法器输出端分别连接第一至第N-1加法器输入 端，该第N个常系数乘法器输出端经第N-1个M延迟单元连接至第N-1加法器另一输入端， 该第N-2加法器至第一加法器的输出端分别经第N-2延迟单元至第一延迟单元连接至上一 级加法器输入端，该第一加法器输出端作为子滤波输出端。
[0013] 优选的，所述的合成累加器为采用进位保留和流水线结构以减少加法器的运算延 迟。
[0014] 优选的，所述的单位延迟单元、合成累加器、多相滤波器和可变系数乘法器均为采 用两通道对同相和正交信号同步进行滤波的复数运算单元。
[0015] 由上述对本发明的描述可知，与现有技术相比，本发明具有如下有益效果：
[0016] 一、本发明结构实现的可变子带滤波器可以任意组合不同的通道滤波器，形成2M 种通带结构的滤波器。
[0017] 二、本发明结构实现的可变子带滤波器不但可以实现全通功能，而且通带内的波 纹特性均与低通原型滤波器一致，信号经过滤波器后不会产生幅度失真
[0018] 三、本发明采用多相分解结构实现可变多子带滤波，可以有效减少所需硬件资源。
[0019] 四、本发明结构和一般可变数字滤波器的实现技术相比，可以减少所需改变参数 的个数，即一般可变滤波器所需改变的参数为数字滤波器系数的个数（一般为30以上），而 本发明的可变子带滤波器所需改变参数个数为子带数远远少于滤波器系数个数。
[0020] 五、本发明结构实现的可变子带滤波器使用的是FIR滤波器，具有线性相位特性， 即对带内信号的群延迟为常数。
【附图说明】
[0021] 图1是现有的可变子带数字滤波器的基本原理框图；
[0022] 图2是本发明可变子带数字滤波器的直接型结构（实施例一）；
[0023] 图3是可变子带数字滤波器的转置结构（实施例二）
[0024] 图4是直接型FIR多相滤波器的结构；
[0025] 图5是转置型FIR多相滤波器的结构；
[0026] 图6是第4带滤波器的冲击响应波形（N= 79);
[0027] 图7是第4带滤波器的幅频和相频响应（N= 79);
[0028] 图8是只选通第1个子带时，即，争时，可变子带滤波器的幅频响应；
[0029]图9是当选通第2和第3个子带时，即& =时，可变子带滤波器的幅频 响应。
【具体实施方式】
[0030] 以下通过【具体实施方式】对本发明作进一步的描述。
[0031] 实施例一
[0032] 一种可变子带数字滤波器，基于离散傅里叶变换（DFT)滤波器组实现的均匀子带 滤波器，用于将[0, 2 31 ]频带分解为M个均匀子带进行选通滤波。
[0033] 参照图2,其直接型FIR滤波器实现的结构包括滤波输入端x[n]、滤波输出端 y[n]、M-1个单位延迟单元z'M-1个合成累加器400、M个可变系数乘法器300和M个多相 滤波器匕(zM)…EM_i(zM)。该滤波输入端x[n] -路连接第0相多相滤波器输入端匕(zM)，另 一路经M-1个单位延迟单元依次延迟构成M-1级延迟支路；该M-1级延迟支路分别连接第 一至第M-1相多相滤波器Ei(zM)…Eh(zM)的输入端；M个多相滤波器EQ (zM)…Eh(zM)的输 出端分别连接到M个可变系数乘法器300的一个输入端，而M个可变系数乘法器300另一 输入端分别连接M个可变系数Si其中i=0、1、2…M-1。该第0至第M-2个可变系数乘法 器300的输出端分别连接第一至第M-1个合成累加器400的输入端，该第M-1个可变系数 乘法器300输出端连接第M-1合成累加器400另一输入端，第M-1合成累加器依次累加至 第一合成累加器，该第一合成累加器的输出端作为滤波输出端y[n]。
[0034] 单位延迟单元z'合成累加器400、多相滤波器EJ/)…EM_i(zM)和可变系数乘法 器的可变系数Si均为复数运算单元，即采用两通道对同相和正交信号同步执行。该单位延 迟单元夕由多位寄存器实现，其位宽等于输入信号位宽。合成累加器400为多输入加法 器，采用进位保留加法器将多个输入压缩为进位（Carry)与和值（Sum)两个位矢量，且仅在 接近输出端的第一个合成累加器相加得到输出结果。M个可变系数乘法器300另一端输入 的可变系数SpM-1彡i彡0,等于各个子带调制系数WikM之和，其中i为所连接多相滤波 器的序号，M-1多i多0，M为子带总数，k为子带序号，子带滤波器含有几个通带（通带为幅 度没有衰减的频率范围），Si就含有几项WikM，即：
[0035] 该实施例中的多相滤波器EJ/)为直接型FIR'滤汲器，i= 0、1…M-1。参照图4， 其包括子滤波输入端、子滤波输出端、N个常系数乘法器220、N-1个M延迟单元和N-l个加法器230。该子滤波输入端一路连接第一个常系数乘法器输入端，另一路经N-1个M延 迟单元z^M依次延迟构成N-1级延迟支路。该N-1级延迟支路分别连接第二个至第N个常 系数乘法器的一输入端，第一至第N个常系数乘法器的另一个输入端分别连接滤波器系数 hQ至hN_i，该第0至第N-2个常系数乘法器的输出端分别连接第一至第N-1合加法器输入 端，该第N-1个常系数乘法器的输出端连接第N-1加法器的另一输入端，第N-1加法器依次 累加至第一加法器，该第一加法器的输出端作为子滤波输出端。
[0036] 以下说明如何得到多相分解结构的可变子带滤波器各相滤波器 的系数。
[0037] 本发明可变子带滤波器原型低通滤波器为第M带（或者称为Nyquist)滤波器。该 原型低通滤波器为只有一个基本通带，即Si只含有一项k=0子带时的滤波器。因此，可变 子带滤波器的原型滤波器系数可由第M带（或者称为Nyquist)滤波器设计得到。第M带 (或者称为Nyquist)滤波器可以通过窗口法和优化方法设计得到。这里使用窗口法说明设 计过程。
[0038] 首先，根据可变子带滤波器阻带（通带）波纹（如60dB)和通带角频率和阻带角 频率的过渡带宽度（范围为〇至2Jr/M)，选择窗口形状。由于凯瑟窗可以通过改变参数得 到任意波纹指标，所以这里使用凯瑟窗。
[0039] 其次，通过对通带宽度为2 31 /M的低通滤波器的冲击响应波形进行抽样，抽样频 率为M。由通带角频率《为[-JT/M，jt/M]理想滤波器的傅里叶逆变换求出其冲击响应 函数为?
[0040]
[0041] 然后，第M带滤波器由以上冲击响应序列取关于时间序列号n= 0为中心的N= 2MJ-1点后（N为FIR滤波器长度），即MJ-1彡n彡-MJ+1，左移（延迟）（N-l)/2 =MJ-1 转换为因果FIR滤波器后，2MJ-2彡n彡0,各系数再乘以凯瑟窗序列Wk[n]后得到，其中下 标k表示凯瑟窗，即：
[0042]
[0043]
[0044]
[0045] 1。艿修止災塞迅凼数，n」以农不艿天卞m的希级数，此处m耿办坝即吋以满足精 度要求；a、0为控制滤波器性能的参数，0由a推出，通过a可以改变窗函数的特性。 a增大，滤波器幅度谱主瓣增宽，旁瓣幅度减小。例如，a=2.12时，滤波器幅度谱主瓣过 渡带宽为3Jr/N，通带波纹0. 27dB，阻带最小衰减30dB;a= 7. 865时，过渡带宽为10JT/ N，通带波纹2. 75e-4dB，阻带最小衰减80dB。附图6和7给出了第4带滤波器的冲击响应 波形和幅频特性与相频特性频谱图（a= 7. 9,J= 10)。
[0046] 最后将求出的第M带滤波器系数前面补一个0后长麼夺为2MT?可W分为M 组（或称M相）多相分解滤波器￡力）（ei[n])的系数，即：
，或
ejn]表示Ei的第n个系数，i为第i相。以上第4带滤波器经4相 分解后的系数如下表所示。
[0047]
[0048]
[0049] 以上设计了含有4个子带的可变子带滤波器，子带的选通可以通过改变Si实现。 例如，附图8呈现了当只选通第1个子带时，即$ =/￥{时，可变子带滤波器的幅频响应；附 图9呈现了当选通第2和第3个子带时，即& = /f2i+ e,'. 时，可变子带滤波器的幅频响 应。
[0050]实施二
[0051] 一种可变子带数字滤波器，用于将[0,231]频带分解为M个均匀子带进行选通滤 波。其与实施例一的区别如下：
[0052] 参照图3,其转置型FIR滤波器实现的结构包括M-1个单位延迟单元、M-1个合成 累加器400、M个可变系数乘法器300和M个多相滤波器& (/)…Rh(/);该滤波输入端分 成M支路，该M支路分别连接第0至第M个多相滤波器输入端&(/)…Rn(zM) ;M个多相滤 波器R。(/)…Rh(/)的输出端分别连接到M个可变系数乘法器300的一个输入端，而M个 可变系数乘法器300的另一输入端分别连接M个可变系数Sp其中i= 0、1、2…M-1。该第 0至第M-2个可变系数乘法器300输出端分别连接第一至第M-1合成累加器400的一输入 端，该第M-1个可变系数乘法器输出端经第M-1延迟单元连接第M-1合成累加器另一输入 端，该第M-2合成累加器至第一合成累加器的输出端分别经第M-2延迟单元至第一延迟单 元连接至上一级合成累加器输入端，该第一合成累加器输出端作为滤波输出端。
[0053] 该实施例中的单位延迟单元z'合成累加器400、多相滤波器RQ(zM)…l(zM)和 可变系数乘法器Si均为复数运算单元，即采用两通道对同相和正交信号同步执行。转置型 单位延迟单元寄存器位于多相滤波器输出端，因此位宽由多相滤波器输出结果的最大摆幅 和截断的精度设定。合成累加器为多输入加法器，采用进位保留加法器将多个输入压缩为 进位（Carry)与和值（Sum)两个位矢量，且仅在接近输出端的第一个合成累加器相加得到 输出结果。M个可变系数乘法器另一端输入的可变系数SpM-1彡i彡0,等于各个子带调 制系数WikM之和，其中i为所连接多相滤波器的序号，M-1多i多0，M为子带总数，k为子 带序号，子带滤波器含有几个通带（通带为幅度没有衰减的频率范围），Si就含有几项WikM， 即
[0054] 多相滤波器氏(/)为转置型FIR滤波器，其中i= 0、1、2…M-1。参照图5,它包括 子滤波输入端、子滤波输出端、N个常系数乘法器220、N-1个M延迟单元厂"和N-1个加法 器230。该子滤波输入端分成N支路，该N支路分别连接第一级至第N个常系数乘法器的 一个输入端，该N个常系数乘法器的另一个输入端分别连接滤波器系数心至h^，该第一至 第N-1个常系数乘法器输出端分别连接第一至第N-1加法器输入端，该第N个常系数乘法 器输出端经第N-1个M延迟单元连接至第N-1加法器另一输入端，该第N-2加法器至第一 加法器的输出端分别经第N-2延迟单元至第一延迟单元连接至上一级加法器输入端，该第 一加法器输出端作为子滤波输出端。
[0055] 上述仅为本发明的【具体实施方式】，但本发明的设计构思并不局限于此，凡利用此 构思对本发明进行非实质性的改动，均应属于侵犯本发明保护范围的行为。
【主权项】
1. 一种可变子带数字滤波器，用于将[〇, 2 31 ]频带分解为M个均匀子带进行选通滤波， 包括滤波输入端和滤波输出端，其特征在于：其结构包括直接型FIR滤波器或转置型FIR滤 波器。2. 如权利要求1所述的一种可变子带数字滤波器，其特征在于：直接型FIR滤波器实 现的结构如下，包括滤波输入端、滤波输出端、M-I个单位延迟单元、M-I个合成累加器、M个 可变系数乘法器和M个多相滤波器；该滤波输入端一路连接第0相多相滤波器输入端，另 一路经M-I个单位延迟单元依次延迟构成M-I级延迟支路；该M-I级延迟支路分别连接第 一至第M-I相多相滤波器输入端；M个多相滤波器的输出端分别连接到M个可变系数乘法 器的一个输入端，而M个可变系数乘法器的另一输入端分别连接M个可变系数S i,其中i = 0、1、2…M-I ;该第0至第M-2个可变系数乘法器的输出端分别连接第一至第M-I个合成累 加器的一输入端，该第M-I个可变系数乘法器输出端连接第M-I合成累加器另一输入端，第 M-2合成累加器依次累加至第一合成累加器，该第一合成累加器的输出端作为滤波输出端。3. 如权利要求1所述的一种可变子带数字滤波器，其特征在于：转置型FIR滤波器实 现的结构如下，包括M-I个单位延迟单元、M-I个合成累加器、M个可变系数乘法器和M多相 滤波器；该滤波输入端分成M支路，该M支路分别连接第0至第M-I个多相滤波器输入端； M个多相滤波器的输出端分别连接到M个可变系数乘法器的一个输入端，而M个可变系数乘 法器的另一输入端为M个可变系数S i,其中i = 0、1、2…M-I ;该第0至第M-2个可变系数 乘法器输出端分别连接第一至第M-I合成累加器输入端，该第M-I个可变系数乘法器输出 端经第M-I延迟单元连接第M-I合成累加器另一输入端，该第M-2合成累加器至第一合成 累加器的输出端分别经第M-2延迟单元至第一延迟单元连接至上一级合成累加器输入端， 该第一合成累加器输出端作为滤波输出端。4. 如权利要求2所述的一种可变子带数字滤波器，其特征在于：所述M个多相滤波器 为直接型FIR滤波器。5. 如权利要求4所述的一种可变子带数字滤波器，其特征在于：所述多相滤波器包括 子滤波输入端，子滤波输出端、N个常系数乘法器、N-I个M延迟单元和N-I个加法器；该子 滤波输入端一路连接第一个常系数乘法器输入端，另一路经N-I个M延迟单元依次延迟构 成N-I级延迟支路；该N-I级延迟支路分别连接第二个至第N个常系数乘法器输入端，第一 至第N个常系数乘法器的另一个输入端分别连接滤波器系数&至h N_i，该第一至第N-I个 常系数乘法器的输出端分别连接第一至第N-I合加法器输入端，该第N个常系数乘法器的 输出端连接第N-I加法器的另一输入端，第N-I加法器依次累加至第一加法器，该第一加法 器的输出端作为子滤波输出端。6. 如权利要求3所述的一种可变子带数字滤波器，其特征在于：所述多相滤波器为转 置型FIR滤波器。7. 如权利要求6所述的一种可变子带数字滤波器，其特征在于：所述多相滤波器包括 子滤波输入端、子滤波输出端、N个常系数乘法器、N-I个M延迟单元和N-I个加法器；该子 滤波输入端分成N支路，该N支路分别连接第一至第N个常系数乘法器的一个输入端，该N 个常系数乘法器的另一个输入端分别连接滤波器系数b至h N_i，该第一至第N-I个常系数 乘法器输出端分别连接第一至第N-I加法器输入端，该第N个常系数乘法器输出端经第N-I 个M延迟单元连接至第N-I加法器另一输入端，该第N-2加法器至第一加法器的输出端分 别经第N-2延迟单元至第一延迟单元连接至上一级加法器输入端，该第一加法器输出端作 为子滤波输出端。8. 如权利要求3所述的一种可变子带数字滤波器，其特征在于：所述的合成累加器为 采用进位保留和流水线结构以减少加法器的运算延迟。9. 如权利要求2或3所述的一种可变子带数字滤波器，其特征在于：所述的单位延迟 单元、合成累加器、多相滤波器和可变系数乘法器均为采用两通道对同相和正交信号同步 进行滤波的复数运算单元。
【专利摘要】一种可变子带数字滤波器，用于将[0,2π]频带分解为M个均匀子带进行选通滤波，其结构可以由直接型FIR滤波器或转置型FIR滤波器实现。本发明高效的可变子带滤波器，可以任意组合不同的通道滤波器，形成2M种通带结构的滤波器，有效实现可变子带滤波器的同时，可以降低对硬件资源的占用。
【IPC分类】H03H17/02
【公开号】CN104883157
【申请号】CN201510252334
【发明人】黄锐敏, 朱述伟, 李国刚, 凌朝东
【申请人】华侨大学
【公开日】2015年9月2日
【申请日】2015年5月18日