CN108777569A - 基于多相滤波器的任意延时方法 - Google Patents
基于多相滤波器的任意延时方法 Download PDFInfo
- Publication number
- CN108777569A CN108777569A CN201810501912.5A CN201810501912A CN108777569A CN 108777569 A CN108777569 A CN 108777569A CN 201810501912 A CN201810501912 A CN 201810501912A CN 108777569 A CN108777569 A CN 108777569A
- Authority
- CN
- China
- Prior art keywords
- delay
- filter
- coefficient
- integer
- decimal
- 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.)
- Pending
Links
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H17/0009—Time-delay networks
Landscapes
- Filters That Use Time-Delay Elements (AREA)
Abstract
本发明涉及数字信号处理技术领域,特别涉及一种基于多相滤波器的任意延时方法。包括如下步骤:(A)将延时分解为整数延时和小数延时,获得整数延时d和小数延时Δ;(B)将原型FIR滤波器分解为I相滤波器,每相的滤波器系数个数为K;(C)根据多相滤波器的I和小数延时Δ获取对应的子滤波器或虚拟子滤波器的系数Px(k);(D)将整数延时d和小数延时对应的滤波器系数Px(k)配置到相应的位置;(E)将输入信号通过d级级联寄存器或者FIFO后输出获得整数倍延时d·Ts;(F)将经过整数倍延时的信号再通过系数为Px(k)的滤波器后输出即可。该延时方法精度高、性能好,可满足阵列信号延时等应用的要求;同时该方法运算量小,便于软件及FPGA实现。
Description
技术领域
本发明涉及数字信号处理技术领域,特别涉及一种基于多相滤波器的任意延时方法。
背景技术
在数字信号处理领域,对信号的延时处理是比较常用也是非常重要的,延时处理在相控阵天线的波束成型、延时相关、多通道相参校准、数字接收器的时间调整等方面均有广泛的应用。在这些应用中,大部分应用对延时的精度具有较高的要求,而传统的整数延时是采用数字方式进行延时,但现有的A/D采样速率是有限的,不可能无限高,因此延时精度受到限制。
参阅图1a、1b,为了提升延时精度,一般采用内插、滤波、延时、抽取的方式,通过增加信号的采样率来降低延时步进并保证输出采样率与输入采样率相一致,从而使得延时精度不再受限于原始采样时钟周期。其流程大致如下:
1、将延时分解为粗延时d·Ts和细延时Δ·Ts;
2、将信号进行d个采样点的粗延时;
3、进行细延时:
(1)将信号进行I倍内插提高采样率,缩小采样时间为
(2)计算当前采样率下细延时Δ·Ts对应的整数倍个T′s,即Δ·Ts≈d′·T′s;
(3)将内插后的信号进行整数倍d′采样点的延时;
4、对延时后的信号进行I倍抽取,获取原采样率下延时后的信号。
这种方法的优点是能够不受限于原始信号的采样率,其缺点也比较明显,即内插滤波的运算量大。
发明内容
本发明的目的在于提供一种基于多相滤波器的任意延时方法,可以实现高精度的任意小数延时。
为实现以上目的,本发明采用的技术方案为:一种基于多相滤波器的任意延时方法,包括如下步骤:(A)将延时分解为整数延时和小数延时:T=d·TS+Δ·TS,获得整数延时d和小数延时Δ;(B)将原型FIR滤波器分解为I相滤波器,每相的滤波器系数个数为K,经过多相分解之后的滤波器系数通过以下映射方式获得:Pi(k)=h(k·I+i);i=0,1,...,I-1;k=0,1,...,K-1,其中h(k·I+i)为FIR滤波器系数;(C)记x=Δ·I,若x∈{0,1,...,I-1},则直接根据步骤B中的映射方式获得滤波器系数Px(k);否则,通过牛顿法拟合滤波器系数,获取x值对应的虚拟子滤波器的系数Px(k);(D)将整数延时d和小数延时对应的滤波器系数Px(k)配置到相应的位置;(E)将输入信号通过d级级联寄存器或者FIFO后输出获得整数倍延时d·TS;(F)将经过整数倍延时的信号再通过系数为Px(k)的滤波器后输出即可。
与现有技术相比,本发明存在以下技术效果:本发明将延时分解为整数倍采样周期的延时和小数倍采样周期的延时,通过利用整数延时方法以及多相滤波器的分解、小数延时滤波器的牛顿法拟合,实现对输入信号进行任意延时,该延时方法精度高、性能好,可满足阵列信号延时等应用的要求;同时该方法运算量小,便于软件及FPGA实现。
附图说明
图1a、1b是现有技术中两种多相内插滤波器的结构示意图;
图2是本发明的信号任意延时实现框图;
图3是本发明整数延时的实现框图;
图4是本发明小数延时的实现框图,其中不需要拟合滤波器系数;
图5是本发明小数延时的实现框图,其中需要对滤波器系数进行拟合;
图6是将原型滤波器分解为I个多相结构的子滤波器示意图;
图7是在图6的基础上进行虚拟子滤波器系数拟合示意图。
具体实施方式
下面结合图2至图7,对本发明做进一步详细叙述。
参阅图2,一种基于多相滤波器的任意延时方法,包括如下步骤:(A)将延时分解为整数延时和小数延时:T=d·TS+Δ·TS,获得整数延时d和小数延时Δ,其中TS为采样周期;设计原型FIR滤波器h(n),此为固定滤波器系数,可先采用matlab等工具设计好。(B)参见图6,将原型FIR滤波器分解为I相滤波器,每相的滤波器系数个数为K,经过多相分解之后的滤波器系数通过以下映射方式获得:Pi(k)=h(k·I+i);i=0,1,...,I-1;k=0,1,...,K-1,;其中h(k·I+i)为FIR滤波器系数,这里的I和K都是正整数。(C)记x=Δ·I,若x∈{0,1,...,I-1},说明小数延时正好对应分解后的某一相子滤波器的输出,则直接根据步骤B中的映射方式获得滤波器系数Px(k),如图4所示;否则,通过牛顿法拟合滤波器系数,获取x值对应的虚拟子滤波器的系数Px(k),如图5所示,此时的小数延时正好对应拟合的虚拟子滤波器的输出。(D)将整数延时d和小数延时对应的滤波器系数Px(k)配置到相应的位置;(E)将输入信号通过d级级联寄存器或者FIFO后输出获得整数倍延时d·TS;(F)将经过整数倍延时的信号再通过系数为Px(k)的滤波器后输出即可。本发明将延时分解为整数倍采样周期的延时和小数倍采样周期的延时,通过利用整数延时方法以及多相滤波器的分解、小数延时滤波器的牛顿法拟合,实现对输入信号进行任意延时,该延时方法精度高、性能好,可满足阵列信号延时等应用的要求;同时该方法运算量小,便于软件及FPGA实现。
以上步骤中,步骤A-D是初始化的过程,可以在软件中完成,且只需要一次;对于信号的延时处理主要是通过整数延时以及一个滤波器完成的,即步骤E和步骤F。因此运算量比较小。
虚拟子滤波器的系数Px(k)可以有多种方式来拟合,本发明中优选地,所述的步骤C中,按如下方式获取x值对应的虚拟子滤波器Px(k)的系数:
(S1)计算滤波器系数的0~I-1阶差分:
零阶:P0(k)、P1(k)、...、PI-1(k);
一阶:ΔP0(k)=P1(k)-P0(k)、ΔP1(k)=P2(k)-P1(k)、...、
ΔPI-2(k)=PI-1(k)-PI-2(k);
二阶:Δ2P0(k)=ΔP1(k)-ΔP0(k)、Δ2P1(k)=ΔP2(k)-ΔP1(k)、...、Δ2PI-3(k)=ΔPI-2(k)-ΔPI-3(k);
……
I-1阶:ΔI-1P0(k)=ΔI-2P1(k)-ΔI-2P0(k);
(S2)根据下列公式计算j:
(S3)将步骤S1和S2中计算到的数据代入下式中计算得到Px(k):
式中k=0,1,...,K-1。通过上述步骤计算得到的虚拟子滤波器Px(k)的系数非常精确,这就保证了延时的高精度。
具体地,在设计原型滤波器时,一般原型FIR滤波器的系数个数M正好可以分解为两个整数相乘,这样就可以保证进行I相分解时,每相的滤波器个数相等。但也有可能存在不能分解的情况,即所述的步骤B中,若原型FIR滤波器的系数个数M不能够分解为两个整数I和K相乘,则需要对原型FIR滤波器系数补零来实现M′=I·K,其中M′为补零后的滤波器系数个数。
整数延时比较简单,可以通过采用移位寄存器或者FIFO的方式,直接对原始信号按采样点进行搬移即可。整数延时下,信号延时量为d·TS,其中d为正整数。在实际实现中,如果d比较小,可以采用d级移位寄存器进行实现,如图3所示。如果d过大,可以采用存储器代替寄存器进行长时间延时。即采用FIFO实现,其步骤如下:(E1)将输入信号缓存到FIFO中,FIFO的存储深度大于整数延时d,否则会溢出;(E2)等待整数延时d后,再将数据从FIFO中读取出来并输出即完成延时。
下面从原理上来详细阐述下以上流程。假设原型FIR滤波器的数字个数为M,将其分解为I相滤波器,每相的滤波器系数个数为K,则M=K·I。经过多相分解之后的滤波器系数可以通过以下映射方式获得:Pi(k)=h(k·I+i);i=0,1,...,I-1;k=0,1,...,K-1,其中h(k·I+i)为FIR滤波器系数。
因此,可以将该映射按自然顺序填充相继的列将一维序列h(n)(FIR滤波器系数)映射为I行K列的二维数组,每一行的K个值为一个子滤波器的系数,原FIR滤波器则分解为I个多相结构的子滤波器,如图6所示。
据传统的多相滤波器延时原理可知,输入的信号先提升采样率、延时后还需要再进行抽取,实际处理时没有必要对原始的信号采样率进行提升。因此,本发明在此基础上进行优化,参见多相内插滤波器第二种结构即图1b,由于该内插滤波器在开关转换之后输出采样率相比输入提升了I倍,同时由于输出是由开关输入I路原始采样率信号合成,因此开关转换之前的I路输出信号之间具有严格的延时关系:即具有不同延时,假设第0路信号的延时为0,则第一路信号延时为TS/I,第二路信号延时为2Ts/I,…,第I-1路信号延时为(I-1)Ts/I。
本发明中,当小数延时正好是{0、Ts/I、2Ts/I、...、(I-1)Ts/I}这些时间里的任何一个,通过开关直接固定连接在当前相应的位置即可实现对输入信号的固定延时,如图4所示。当Δ·Ts=i·TS/I时,将开关固定连接在第i路输出,如图4中实线所示,则得到对输入进行j·Ts/I延时的输出。该多相滤波器经过分解后,Pi(k)即为用于计算延时为i·Ts/I的子滤波器。由于小数延时是任意的,所以会存在小数延时不属于{0、Ts/I、2Ts/I、...、(I-1)Ts/I}这些时间里的任何一个,而是介于其中两个时间的中间,此时,通过虚拟出一个子滤波器,如图7所示,并通过拟合的方式,将该虚拟子滤波器Px(k)的系数计算出来,然后,将开关固定连接在该虚拟子滤波器上输出即可实现任意时间的小数延时。
Claims (4)
1.一种基于多相滤波器的任意延时方法,包括如下步骤:
(A)将延时分解为整数延时和小数延时:T=d·Ts+Δ·Ts,获得整数延时d和小数延时Δ;
(B)将原型FIR滤波器分解为I相滤波器,每相的滤波器系数个数为K,经过多相分解之后的滤波器系数通过以下映射方式获得:
Pi(k)=h(k·I+i);i=0,1,...,I-1;k=0,1,...,K-1
其中h(k·I+i)为FIR滤波器系数;
(C)记x=Δ·I,若x∈{0,1,...,I-1},则直接根据步骤B中的映射方式获得滤波器系数Px(k);否则,通过牛顿法拟合滤波器系数,获取x值对应的虚拟子滤波器的系数Px(k);
(D)将整数延时d和小数延时对应的滤波器系数Px(k)配置到相应的位置;
(E)将输入信号通过d级级联寄存器或者FIFO后输出获得整数倍延时d·Ts;
(F)将经过整数倍延时的信号再通过系数为Px(k)的滤波器后输出即可。
2.如权利要求1所述的基于多相滤波器的任意延时方法,其特征在于:所述的步骤C中,按如下方式获取x值对应的虚拟子滤波器Px(k)的系数:
(S1)计算滤波器系数的0~I-1阶差分:
零阶:P0(k)、P1(k)、...、PI-1(k);
一阶:ΔP0(k)=P1(k)-P0(k)、ΔP1(k)=P2(k)-P1(k)、...、ΔPI-2(k)=PI-1(k)-PI-2(k);
二阶:Δ2P0(k)=ΔP1(k)-ΔP0(k)、Δ2P1(k)=ΔP2(k)-ΔP1(k)、...、Δ2PI-3(k)=ΔPI-2(k)-ΔPI-3(k);
……
I-1阶:ΔI-1P0(k)=ΔI-2P1(k)-ΔI-2P0(k);
(S2)根据下列公式计算j;
(S3)将步骤S1和S2中计算到的数据代入下式中计算得到Px(k):
式中k=0,1,...,K-1。
3.如权利要求1所述的基于多相滤波器的任意延时方法,其特征在于:所述的步骤B中,若原型FIR滤波器的系数个数M不能够分解为两个整数I和K相乘,则需要对原型FIR滤波器系数补零来实现M′=I·K,其中M′为补零后的系数个数。
4.如权利要求1所述的基于多相滤波器的任意延时方法,其特征在于:所述的步骤E中,输入信号采用FIFO获得整数倍延时的步骤如下:
(E1)将输入信号缓存到FIFO中,FIFO的存储深度大于整数延时d;
(E2)等待整数延时d后,再将数据从FIFO中读取出来并输出即可。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810501912.5A CN108777569A (zh) | 2018-05-23 | 2018-05-23 | 基于多相滤波器的任意延时方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810501912.5A CN108777569A (zh) | 2018-05-23 | 2018-05-23 | 基于多相滤波器的任意延时方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108777569A true CN108777569A (zh) | 2018-11-09 |
Family
ID=64027472
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810501912.5A Pending CN108777569A (zh) | 2018-05-23 | 2018-05-23 | 基于多相滤波器的任意延时方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108777569A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110708039A (zh) * | 2019-09-17 | 2020-01-17 | 广西师范大学 | 一种farrow滤波器的系数拟合方法 |
CN111865270A (zh) * | 2020-07-07 | 2020-10-30 | 北京无线电计量测试研究所 | 一种信号延时方法和装置 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080298432A1 (en) * | 2007-06-04 | 2008-12-04 | Yugang Ma | Ranging system and method |
US20100249596A1 (en) * | 2009-03-24 | 2010-09-30 | Texas Instruments Incorporated | Receive beamformer for ultrasound having delay value sorting |
CN102291154A (zh) * | 2011-09-23 | 2011-12-21 | 电子科技大学 | 极坐标发射机 |
CN104407061A (zh) * | 2014-12-31 | 2015-03-11 | 南通友联数码技术开发有限公司 | 一种超声信号整数与小数精确延时系统及其方法 |
US20160277007A1 (en) * | 2015-03-20 | 2016-09-22 | Texas Instruments Incorporated | Systems and Methods of Variable Fractional Rate Digital Resampling |
CN106411379A (zh) * | 2016-09-29 | 2017-02-15 | 电子科技大学 | 一种降低硬件资源消耗的宽带波束形成设计方法 |
CN106549650A (zh) * | 2016-12-07 | 2017-03-29 | 西安电子科技大学 | 一种基于fpga的高阶fir滤波器的实现方法 |
CN107144821A (zh) * | 2017-04-07 | 2017-09-08 | 电子科技大学 | 宽带数字阵雷达中基于时延波束形成的高效接收通道 |
-
2018
- 2018-05-23 CN CN201810501912.5A patent/CN108777569A/zh active Pending
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080298432A1 (en) * | 2007-06-04 | 2008-12-04 | Yugang Ma | Ranging system and method |
US20100249596A1 (en) * | 2009-03-24 | 2010-09-30 | Texas Instruments Incorporated | Receive beamformer for ultrasound having delay value sorting |
CN102291154A (zh) * | 2011-09-23 | 2011-12-21 | 电子科技大学 | 极坐标发射机 |
CN104407061A (zh) * | 2014-12-31 | 2015-03-11 | 南通友联数码技术开发有限公司 | 一种超声信号整数与小数精确延时系统及其方法 |
US20160277007A1 (en) * | 2015-03-20 | 2016-09-22 | Texas Instruments Incorporated | Systems and Methods of Variable Fractional Rate Digital Resampling |
US9531343B2 (en) * | 2015-03-20 | 2016-12-27 | Texas Instruments Incorporated | Systems and methods of variable fractional rate digital resampling |
CN106411379A (zh) * | 2016-09-29 | 2017-02-15 | 电子科技大学 | 一种降低硬件资源消耗的宽带波束形成设计方法 |
CN106549650A (zh) * | 2016-12-07 | 2017-03-29 | 西安电子科技大学 | 一种基于fpga的高阶fir滤波器的实现方法 |
CN107144821A (zh) * | 2017-04-07 | 2017-09-08 | 电子科技大学 | 宽带数字阵雷达中基于时延波束形成的高效接收通道 |
Non-Patent Citations (6)
Title |
---|
(美)维纳•K•英格尔,约翰•G•普罗克斯著;刘树棠,陈志刚译: "《数字信号处理MATLAB版》", 31 July 2013, 西安交通大学出版社 * |
H. JOHANSSON AND F. HARRIS: "Polyphase Decomposition of Digital Fractional-Delay Filters", 《IEEE SIGNAL PROCESSING LETTERS》 * |
吴瑛: "《数字信号处理 第2版》", 31 December 2017, 西安电子科技大学出版社 * |
左军: "《数值计算方法与实验学习辅导》", 31 August 2015, 国防工业出版社 * |
王晨: "宽带数字阵列可重构接收通道结构研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
高 磊: "高精度可变延时控制技术研究", 《全球定位系统》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110708039A (zh) * | 2019-09-17 | 2020-01-17 | 广西师范大学 | 一种farrow滤波器的系数拟合方法 |
CN110708039B (zh) * | 2019-09-17 | 2023-08-08 | 广西师范大学 | 一种farrow滤波器的系数拟合方法 |
CN111865270A (zh) * | 2020-07-07 | 2020-10-30 | 北京无线电计量测试研究所 | 一种信号延时方法和装置 |
CN111865270B (zh) * | 2020-07-07 | 2023-04-25 | 北京无线电计量测试研究所 | 一种信号延时方法和装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104993827B (zh) | 模数转换器误差估计校正的装置及其方法 | |
CN108768343A (zh) | 基于多相滤波器的高精度延时方法 | |
JP3337424B2 (ja) | デジタル・フィルタ | |
CN100499371C (zh) | 一种可编程内插滤波器装置及其实现方法 | |
CN108777569A (zh) | 基于多相滤波器的任意延时方法 | |
CN105991137A (zh) | 可变分数速率数字重采样的系统和方法 | |
CN108427102A (zh) | 一种雷达基带杂波生成装置及方法 | |
CN101958697B (zh) | 一种多相滤波结构的实现方法及装置 | |
CN103117972A (zh) | 一种矢量信号分析方法和装置 | |
CN106506135A (zh) | 一种吉比特速率的全数字时域并行定时同步系统及方法 | |
CN106130507B (zh) | 高速并行位同步滤波器时变系数更新环路 | |
CN110365314A (zh) | 具有稀疏系数的可分离二维fir滤波器的设计方法 | |
CN105066990A (zh) | 一种适用于捷联惯性导航的高精度数字滤波器 | |
CN110166021A (zh) | 一种实现任意降采样率转换的数字信号处理方法 | |
CN101540749B (zh) | 可配置变换长度dft的前处理单元的实现方法及装置 | |
EP1274170A1 (fr) | Dispositif de conversion d'une séquence d'échantillons numériques | |
CN105634439A (zh) | 一种异步成形滤波器设计方法 | |
CN106385253B (zh) | 基于参数处理模块和锁相环级联的数字时间转换系统 | |
CN108648148A (zh) | 一种基于数字升采样再三次样条的数字图像任意点插值方法 | |
CN103595670B (zh) | 一种信号偏移校准方法和装置 | |
CN111131141B (zh) | 话音加、解密方法及装置 | |
CN107145457B (zh) | 基于片上ram提升多通道有效数据传输的装置及方法 | |
CN104579238B (zh) | 应用于全数字接收机的基于Schur算法的插值滤波器 | |
CN105007081B (zh) | 一种四相位数模转换方法及数模转换器 | |
CN102045078B (zh) | 一种基于fpga的软件接收机系统及实现方法 |
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 |