CN105915193B - 一种用于多相滤波器的改进生成方法 - Google Patents

一种用于多相滤波器的改进生成方法 Download PDF

Info

Publication number
CN105915193B
CN105915193B CN201610446815.1A CN201610446815A CN105915193B CN 105915193 B CN105915193 B CN 105915193B CN 201610446815 A CN201610446815 A CN 201610446815A CN 105915193 B CN105915193 B CN 105915193B
Authority
CN
China
Prior art keywords
filter
subfilters
stopband
passband
subfilter
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.)
Expired - Fee Related
Application number
CN201610446815.1A
Other languages
English (en)
Other versions
CN105915193A (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201610446815.1A priority Critical patent/CN105915193B/zh
Publication of CN105915193A publication Critical patent/CN105915193A/zh
Application granted granted Critical
Publication of CN105915193B publication Critical patent/CN105915193B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H7/00—Multiple-port networks comprising only passive electrical elements as network components
    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H7/00—Multiple-port networks comprising only passive electrical elements as network components
    • H03H7/18—Networks for phase shifting
    • H03H7/21—Networks for phase shifting providing two or more phase shifted output signals, e.g. n-phase output

Landscapes

  • Filters That Use Time-Delay Elements (AREA)
  • Noise Elimination (AREA)

Abstract

本发明公开了一种用于多相滤波器的改进生成方法，属于无线电技术领域。本发明包括下列步骤：根据工作需求参数，基于FIR滤波器的窗函数法生成FIR滤波器；将FIR滤波器分成n个第一子滤波器，形成n路多相滤波器组，其中n≥2；再将FIR滤波器分为n个第二子滤波器，将调整通带后的第二子滤波器作为基带滤波器，两个相邻的第一子滤波器之间增加一个基带滤波器，所增加的基带滤波器的通带在所连接的两个第一子滤波器的阻带之间。本发明用于多相滤波器，可以有效消除滤波器组之间的谱间干扰。

Description

一种用于多相滤波器的改进生成方法
技术领域
本发明涉及软件无线电技术领域，具体来说，是涉及对多相滤波器进行改进。
背景技术
软件无线电系统的理想结构是低通采样结构，低通采样结构采样的射频信号的带宽最宽，能够并行处理的通道数最多。由于低通采样结构对所有波段内的信号采用的是同一采样频率，因此经数字变频后的基带信号都处于严重的过采样状态。过采样虽然有利于提高信号的信噪比，但其过高的输出数据率要求后续的信号处理要具有很高的处理速度。因此，如何保证在信号不失真的前提下有效的降低信号的数据率，从而降低后续同步、解调等信号处理的速度，成为数字信号处理的前端部分的重要任务。这种信号采样率变化的理论称为多速率信号处理理论。
多速率信号处理的目的就是在不使信号失真的前提下改变信号的采样率。减小信号采样率以减小数据冗余的过程称为信号的抽取；增大信号的采样率以增加数据冗余的过程称为信号的内插。所谓整数倍抽取是指把原始采样序列x(n)每隔(D-1)个数据取一个，以形成一个新序列xD(m)＝x(mD)。由于抽取后的频谱产生了严重的混叠，使得从抽取后的信号频域XD(ejω)中已经无法恢复出感兴趣的信号频谱分量。但是，如果首先采用一数字滤波器对原信号的频域X(ejω)进行滤波，使X(ejω)中只含有小于π/D的频率分量，再进行D倍抽取，则抽取后的频谱就不会发生混叠，经过抽取，数据流数率只有以前的1/D，大大降低了对后续处理(解调分析等)的速度要求。
滤波器是软件无线电系统处理中核心的部分，用于滤出杂散/噪声、信道选择等，基于数字滤波器的冲击响应的Z变换可以得到抽取处理中数字滤波器的多相结构。但是，在实际工程中，想要生成理想滤波器是不可能的，这就使得在多相滤波时，当信号落在两个滤波器之间时，两个滤波器可能都有信号输出，从而使得对信号的频段判定出现偏差，影响后续处理。
发明内容
本发明的目的在于解决现有技术所存在的问题，提供一种类似于“理想”的滤波器，从而可以解决在多相滤波器中可能出现的谱图干扰。
本发明的一种用于多相滤波器的改进生成方法，包括以下步骤：
步骤1：根据工作需求参数，基于FIR(Finite Impulse Response，有限长单位冲激响应滤波器，又称为非递归型滤波器)滤波器的窗函数法生成FIR滤波器；
步骤2：将FIR滤波器分成n(n≥2)个第一子滤波器，形成n路多相滤波器组，n的具体取值取决于工作需求；
步骤3：再将FIR滤波器分为n个第二子滤波器，将调整通带后的第二子滤波器作为基带滤波器，在两个相邻的第一子滤波器之间增加一个基带滤波器，所增加的基带滤波器的通带在所连接的两个第一子滤波器的阻带之间，且基带滤波器的通带WP与阻带WS满足2WP≥WS，即用对应的基带滤波器覆盖之前相邻通道的直接频段，通过逻辑判断重新划分滤波器的带宽。
因为在现有的方式生成的多相滤波器组中，因为滤波器不能做到理想，两个通道之间就会有重叠部分，而当信号出现在这个区域时，相邻两个通道都会判断信号出现，就会出现频谱干扰，而通过增加额外滤波器(本发明的基带滤波器)在两个第一子滤波器的相邻位置，就可以把这块区域的带宽重新划分，把其带宽变窄，在干扰区域的频带就落入了增加的基带滤波器中，原有通道间将不会出现谱间干扰，即将宽带带通滤波器通过逻辑划分成为较为窄带的相邻滤波器组。这个滤波器组就是“理想”的滤波器组。
进一步的，步骤1的具体方法为：
步骤101：在MATLAB中，输入工作需求参数，基于FIR滤波器的窗函数法生成FIR滤波器，得到FIR滤波器的工作参数；
步骤102：将FIR滤波器的工作参数保存至文本文档。
进一步的，步骤2的具体方法为：
基于FIR滤波器的工作参数，得到n个第一子滤波器的工作参数，并在FPGA的ROM核存储单元设置一个长度为n的数组，用于存储各第一子滤波器的工作参数。用FPGA中的ROM核存储保存有FIR滤波器的系数的文本，只要ROM的地址固定，每次输出的数据都将是不变的，从而简化多相滤波器的生成过程。
因此在生成两个相邻第一子滤波器之间的基带滤波器时，可以从FPGA的ROM核读取第一子滤波器的工作参数，作为指定相邻第一子滤波器之间的基带滤波器，并基于指定相邻第一子滤波器的阻带调整所述基带滤波器的通带，从而简化实现。例如将步骤1生成的FIR滤波器分为子滤波器0～子滤波器8(子滤波器0～子滤波器8的通带依次从低到高)，取子滤波器0作为子滤波器0与子滤波器1之间的基带滤波器(定义为基带滤波器0)，基于子滤波器0与子滤波器1的阻带调整基带滤波器0的通带，即其通带在子滤波器0与子滤波器1的阻带之间，同理，取子滤波器1作为子滤波器1与子滤波器2之间的基带滤波器并调整期通带，以此类推，可以得到其他相邻子滤波器之间的基带滤波器。
本发明的有益效果是，可以有效消除滤波器组之间的谱间干扰。
附图说明
图1是实施例的实现系统流程。
图2是多相滤波器组的实现示意图。
图3是在现有的多相滤波器组的通道上增加辅助通道的示意图。
图4是本发明实施例生成的“理想”多相滤波器组。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面结合实施方式和附图，对本发明作进一步地详细描述。
本实例提供了一种基于本发明的用于多相滤波器的改进生成方法而得到的多相滤波器，本实施例基于FPGA的数字电路，处理的信号是经过ADC采样的数字信号，采样率为1.25GHz，每个通道的带宽为156.25MHz。本实施例中只处理100MHz到500MHz的信号，即只需要使用多相滤波器中第一个到第三个通道。
如图1所示，是实施例的系统流程图，包括ADC数据采样以及对其采样数据进行补码转化，将处理完的数据送入多相滤波器通道，滤波之后进行8点并行FFT，接着对每一路通道进行1024点串行FFT，最后进行门限判定，确定信号是在哪个频段。
本实例中，用于多相滤波器的改进生成方法，包括以下步骤：
S1：基于FIR滤波器的窗函数法生成FIR滤波器；
基于工作需求参数(通带Wp＝π/8，阻带Ws＝3π/16，阻带纹波as＝65db)，使用MATLAB设计一个Kaiser窗低通FIR滤波器，就可以得到低通FIR滤波器的工作参数，并将该工作参数保持至文本文档。
由于只针对带宽为100MHz到500MHz的信号。在1.25GHz的采样率下，做8相多相滤波器，每个通道的通带带宽为156.25MHz，根据数字信号滤波器的特点，有：
通道(1)的频率范围为：78.125MHz-234.375MHz；
通道(2)的频率范围为：234.375MHz-390.625MHz；
通道(3)的频率范围为：390.625MHz-546.875MHz。
S2：将生成的低通FIR滤波器运用于n路多相滤波器组：参见图2，将步骤S1生成的低通FIR滤波器分成8个子滤波器(也称子带滤波器)：通带从低到高依次为子滤波器0～子滤波器7，也就是将128阶滤波器变成为8个16阶滤波器，在FPGA的ROM核存储单元设置一个长度为8的数组(图2中的Coe0～Coe7)，用于存储各子滤波器的工作参数，每一个Coe里面存储的就是子滤波器的工作参数，从而形成了8路多相滤波器组。即在处理时，将对应子滤波器的工作参数与采样数据进行卷积运算，得出的结果再进行8点FFT运算。
S3：在每两个相邻的子滤波器之间加一个基带滤波器，该基带滤波器的通带在连接的子滤波器的阻带之间，且基带滤波器的通带WP与阻带WS满足2WP≥WS：基于子滤波器生成指定相邻的子滤波器之间的基带滤波器，本实施例中，只需要使用8路多相滤波器中第一个到第三个通道，因此需要用到前面两个子带滤波器来增加额外的基带滤波器，即用基于第一个子滤波器(子滤波器0)生成的基带滤波器覆盖第一个通道和第二通道直接的频段，得到通道4，用基于第二个子滤波器生成的的基带滤波器覆盖第二个通道和第三通道直接的频段，得到通道5，覆盖之前的通道，如图3所示。从而得到本实施例所生成的“理想”多相滤波器组，如图4所示。
将新生成的8路多相滤波器组的频段范围求出，就可以很好的判断每个信号所处的频段了。8路多相滤波器组的频段范围为：
通道编号 频率 1024点FFT Index序号范围
通道1 [100M～195.3125M) [654～256)①
通道2 [273.4375M～351.5625M) [767～256)
通道3 [429.6875M～546.875M) [767～256)
通道4(辅助通道1) [195.3125M～273.4375M) [767～256)
通道5(辅助通道2) [351.5625M～429.6875M) [767～256)
①FFT输出的Index序号从0～1023周期出现。理论上通道1测量范围从117.1875M开始，但由于不使用其低频方向相邻滤波器，因此忽略通道0对100M～117.1875M的影响。
在数据经过1024点串行FFT计算之后，在图1所示的门限判别模块设立一个指数Index，对落在带内的信号的判别采用常规判别处理即可，如：
(1)对于通道1-3，检测谱线值大小是否高于预设门限以及Index是否在[767,256)范围内，如果在则认为该通道有信号出现，否者无。注意：就算该滤波器有高于门限的谱线，但是不在[767,256)范围内，同样认为没有信号。
(2)对通道4-5，按照步骤1来做，不同在于保留[256,767)范围内的信号。
上述判别算法中指数Index为串行FFT模块中FFT核中所输出的10bit数据。
以上所述，仅为本发明的具体实施方式，本说明书中所公开的任一特征，除非特别叙述，均可被其他等效或具有类似目的的替代特征加以替换；所公开的所有特征、或所有方法或过程中的步骤，除了互相排斥的特征和/或步骤以外，均可以任何方式组合。

Claims (4)

1.一种用于多相滤波器的改进生成方法，其特征在于，包括以下步骤：
步骤1：根据工作需求参数，使用MATLAB设计一个Kaiser窗低通FIR滤波器，得到低通FIR滤波器的工作参数，所述工作需求参数包括通带、阻带和阻带纹波；
步骤2：将FIR滤波器分成n个第一子滤波器，形成n路多相滤波器组，其中n≥2；
步骤3：再将FIR滤波器分为n个第二子滤波器，将调整通带后的第二子滤波器作为基带滤波器，在两个相邻的第一子滤波器之间增加一个基带滤波器，所增加的基带滤波器的通带在所连接的两个第一子滤波器的阻带之间，且基带滤波器的通带WP与阻带WS满足2WP≥WS。
2.如权利要求1所述的方法，其特征在于，将低通FIR滤波器的工作参数保存至文本文档。
3.如权利要求2所述的方法，其特征在于，步骤2中，基于生成的低通FIR滤波器的工作参数，得到n个第一子滤波器的工作参数，并在FPGA的ROM核存储单元设置一个长度为n的数组，用于存储各第一子滤波器的工作参数。
4.如权利要求3所述的方法，其特征在于，步骤3为，从FPGA的ROM核读取第一子滤波器的工作参数，作为指定相邻第一子滤波器之间的基带滤波器，并基于指定相邻第一子滤波器的阻带调整所述基带滤波器的通带：基带滤波器的通带在所连接的两个第一子滤波器的阻带之间，且基带滤波器的通带WP与阻带WS满足2WP≥WS。
CN201610446815.1A 2016-06-21 2016-06-21 一种用于多相滤波器的改进生成方法 Expired - Fee Related CN105915193B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610446815.1A CN105915193B (zh) 2016-06-21 2016-06-21 一种用于多相滤波器的改进生成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610446815.1A CN105915193B (zh) 2016-06-21 2016-06-21 一种用于多相滤波器的改进生成方法

Publications (2)

Publication Number Publication Date
CN105915193A CN105915193A (zh) 2016-08-31
CN105915193B true CN105915193B (zh) 2018-06-19

Family

ID=56758111

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610446815.1A Expired - Fee Related CN105915193B (zh) 2016-06-21 2016-06-21 一种用于多相滤波器的改进生成方法

Country Status (1)

Country Link
CN (1) CN105915193B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106549650B (zh) * 2016-12-07 2019-03-26 西安电子科技大学 一种基于fpga的高阶fir滤波器的实现方法
CN107241107B (zh) * 2017-06-28 2019-03-29 电子科技大学 一种数字信道化滤波器组实现方法
CN109561437B (zh) * 2017-09-26 2020-10-16 大唐移动通信设备有限公司 一种对信号进行处理的方法及装置
CN111414257A (zh) * 2020-03-30 2020-07-14 中国科学院新疆天文台 一种脉冲星信号多通道滤波方法、装置及存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8650236B1 (en) * 2009-08-04 2014-02-11 Altera Corporation High-rate interpolation or decimation filter in integrated circuit device
CN104301265A (zh) * 2014-04-11 2015-01-21 中国人民解放军国防科学技术大学 一种超宽带信号的数字正交解调实时处理方法
CN104506161A (zh) * 2014-10-11 2015-04-08 中国电子科技集团公司第十研究所 复系数希尔伯特带通滤波器分数倍采样率转换方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6870882B1 (en) * 1999-10-08 2005-03-22 At&T Corp. Finite-length equalization over multi-input multi-output channels
CN100576881C (zh) * 2007-12-26 2009-12-30 西安交通大学 一种基于多相滤波技术的自适应视频图像缩放引擎
CN102347768B (zh) * 2010-07-28 2014-03-12 中兴通讯股份有限公司 一种数字采样率的转换装置和方法
JP6296604B2 (ja) * 2014-04-17 2018-03-20 国立研究開発法人情報通信研究機構 Ｏｆｄｍ信号送信装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8650236B1 (en) * 2009-08-04 2014-02-11 Altera Corporation High-rate interpolation or decimation filter in integrated circuit device
CN104301265A (zh) * 2014-04-11 2015-01-21 中国人民解放军国防科学技术大学 一种超宽带信号的数字正交解调实时处理方法
CN104506161A (zh) * 2014-10-11 2015-04-08 中国电子科技集团公司第十研究所 复系数希尔伯特带通滤波器分数倍采样率转换方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于多相滤波器组的多参数雷达信号测量与FPGA实现;陈海峰;《中国优秀硕士学位论文全文数据库 信息科技辑》;20110815(第08期);I136-484 *

Also Published As

Publication number Publication date
CN105915193A (zh) 2016-08-31

Similar Documents

Publication Publication Date Title
CN105915193B (zh) 一种用于多相滤波器的改进生成方法
CN101977033B (zh) 用于井下仪器信号传输的数字滤波方法
CN101478525B (zh) 一种多载波分离的方法及多载波分离装置
CN107612865B (zh) 一种应用于低压电力线载波通信的信号降噪方法
CN104467739B (zh) 一种带宽、中心频点可调的数字滤波器及其实现方法
CN107411739A (zh) 基于双树复小波的脑电信号情绪识别特征提取方法
CN106980491A (zh) 一种a/d采样的改进的均值滤波算法
CN108880621B (zh) 一种电力线通信信号自适应滤波方法
CN109257068B (zh) 一种电力线通信信号自适应滤波方法
CN115296647B (zh) 基于frm技术的wola滤波器组及子带分割设计方法
CN115967404A (zh) 一种数字信道化方法及系统
CN109388882A (zh) 一种适用于专用集成电路设计的交错匹配滤波方法
CN110208755A (zh) 一种基于fpga的动态雷达回波数字下变频系统及方法
CN106872948B (zh) 基于两级数字下变频的雷达工作信道采集系统及其方法
Furtado et al. Optimized prototype filter based on the FRM approach for cosine-modulated filter banks
CN113341378A (zh) 基于频谱差分熵检测的自适应信道化接收方法
CN109039379B (zh) 电力线通信信号滤波方法
CN109474356B (zh) 宽带多通道信号能量检测系统及方法
CN103002197A (zh) 一种信号数据处理方法、装置及智能终端
CN114257285B (zh) 一种通感一体化基站感知信号滤波方法及系统
CN102104564B (zh) 一种通用超宽带接收信道化处理方法
CN109951174B (zh) 一种快速傅里叶变换滤波方法
CN108776259A (zh) 基于数字带宽限制技术的功率分析仪
CN101473535A (zh) 用于在数字域处理数字输入信号的方法及用于处理数字输入信号的数字滤波器电路
CN112688666B (zh) 一种仪器仪表二级均值滤波算法

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180619

Termination date: 20210621