CN104539261B - 一种任意采样率转换的内插滤波处理方法 - Google Patents
一种任意采样率转换的内插滤波处理方法 Download PDFInfo
- Publication number
- CN104539261B CN104539261B CN201410745446.7A CN201410745446A CN104539261B CN 104539261 B CN104539261 B CN 104539261B CN 201410745446 A CN201410745446 A CN 201410745446A CN 104539261 B CN104539261 B CN 104539261B
- Authority
- CN
- China
- Prior art keywords
- fir filter
- signal
- filter
- coefficient
- clock
- 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
- Complex Calculations (AREA)
Abstract
本发明涉及一种任意采样率转换的内插滤波处理方法,该方法步骤包括:(1)根据输入信号的带宽要求和滤波精度系数,产生FIR滤波器抽头系数表;(2)用NCO产生1倍数据时钟信号,并得到该时钟信号的时钟相位Phase;(3)采用1倍数据时钟信号进行信号接收;(4)根据时钟相位Phase在FIR滤波器系数表内提取相应的滤波系数;(5)将提取得到的滤波系数应用于FIR滤波器,对接收信号进行滤波处理,即可将输入信号的采样率转换到***工作频率fs,实现采样率转换的内插滤波处理;该方法中采用的插值滤波器由一个横向FIR滤波器和一个装订有滤波器系数的ROM表组成,实现简单,耗费硬件资源少,实用性强。
Description
技术领域
本发明涉及数字滤波技术领域,特别涉及一种任意采样率转换的内插滤波处理方法。
背景技术
在现有数字内插滤波器设计中,将信号采样率变到***工作时钟,为了适应数据速率可变的要求,一种方法是根据传输数据的速率,采用重新配置时钟芯片来改变DA工作时钟,这种方法受硬件平台限制;另外一种方法是针对特定几个速率采用多级滤波器进行整倍数插值滤波,这种方法滤波器级数多、结构复杂、耗费资源巨大,并且数据速率不能连续可变,实用性不强。
发明内容
本发明的目的在于克服现有技术的不足,提供了一种任意采样率转换的内插滤波处理方法,该方法中采用的插值滤波器由一个横向FIR滤波器和一个装订有滤波器系数的ROM表组成,在每个***工作时钟fs周期内,根据量化后的时钟相位phase查找内插滤波器系数表,实时更新FIR滤波器系数并进行滤波运算,得到采样率为fs的插值后的信号,该内插滤波方法实现简单,耗费硬件资源少,实用性强。
本发明的上述目的主要是通过如下技术方案予以实现的:
一种任意采样率转换的内插滤波处理方法,包括以下步骤:
1、一种任意采样率转换的内插滤波处理方法,其特征在于包括以下步骤:
(1)、根据输入信号的带宽Bw和设定的信号分辨率参数K,生成K组FIR滤波器系数,并将所述K组滤波器系数保存在FIR滤波器抽头系数表内;其中,所述信号分辨率参数K=2L,L为正整数;
(2)、根据输入信号的码率Rb和***时钟频率fs,采用NCO产生1倍数据时钟信号,所述数据时钟信号的频率Rs=Rb,所述数据时钟信号的时钟相位Phase等于所述NCO累加值的高L位截取值,即所述时钟相位Phase的取值范围为0~K-1;其中,所述NCO的位宽为M且M≥L;
(3)、按照步骤(2)产生的1倍数据时钟信号对输入信号进行接收,得到接收信号;
(4)、根据信号接收时钟的时钟相位Phase的取值,在步骤(1)产生的FIR滤波器抽头系数表中提取FIR滤波器的抽头系数,其中,如果所述时钟相位Phase=k-1,则在所述FIR滤波器抽头系数表中提取得到第k组低通FIR滤波器系数h′k(n),其中,k=1,2,…,K,n=1,2,…,N,N为所述FIR滤波器的抽头个数;
(5)、将步骤(4)提取得到的FIR滤波器系数应用于FIR滤波器中,对步骤(3)中得到的接收信号进行内插滤波处理,得到内插滤波后信号。
上述的任意采样率转换的内插滤波处理方法,在步骤(1)中,根据输入信号的带宽Bw和设定的信号分辨率参数K,生成K组FIR滤波器系数,并将所述K组滤波器系数保存在FIR滤波器抽头系数表内,具体实现过程如下:
(1a)、根据信号带宽Bw,生成一组满足所述信号带宽内滤波要求的FIR滤波器的抽头系数h(n),其中,n=1,2,…,N,N为所述FIR滤波器的抽头个数;
(1b)、根据设定的信号分辨率参数K得到K个比例函数,其中,第k个所述比例函数为k=1,2,…,K,n=1,2,…,N;
(1c)、对步骤(1a)得到的FIR滤波器抽头系数h(n)进行N点离散傅里叶变换,得到FIR滤波器响应函数H(jw)=DFT(h(n)),其中,DFT()代表离散傅里叶变换;所述响应函数H(jw)的N个值分别为h1,h2,…,hn,…,hN;
(1d)、根据步骤(1b)得到的K个比例函数和步骤(1c)得到的响应函数H(jw),得到K个新的响应函数,其中,第k个所述新的响应函数为H′k(jw),H′k(jw)的N个值分别为h1′,h2′,…,hn′,…,hN′,其中,hn′=hn×fk(n),k=1,2,…,K,n=1,2,…,N;
(1e)、对步骤(1d)得到的K个新的响应函数进行逆傅里叶变换得K组滤波器系数,其中,第k组所述滤波器系数为h′k(n)=IDFT(H′k(jw)),IDFT代表逆离散傅里叶变换,k=1,2,…,K;
(1f)、将步骤(1e)中得到K组滤波器系数依次存入数据表中,得到FIR滤波器抽头系数表。
上述的任意采样率转换的内插滤波处理方法,FIR滤波器抽头系数表装订在ROM表中。
本发明与现有技术相比的优点在于:
(1)、本发明由一个横向FIR滤波器和一个装订有滤波器系数的ROM构成分数倍插值滤波器,通过该分数倍插值滤波器可以实现不同速率信号的采样率由信号频率到***工作频率的的采样率转换,实现结构简单,且可以满足输入信号的连续可变要求,实用性强。
(2)本发明的分数倍插值滤波器由一个横向FIR滤波器和一个装订有滤波器系数的ROM表组成,在每个***工作时钟fs周期内,根据量化后的时钟相位phase查找内部滤波器系数表,实时更新FIR滤波器系数并进行滤波运算,得到采样率为fs的插值后的信号。该滤波器实现方法简单,所用硬件资源少,可以有效克服现有技术受硬件平台限制、本发明方法中的插值滤波器结构简单、运算量小;
(3)本发明中使用的装订滤波器系数的ROM表大小可根据输入信号分辨率的要求灵活调整,即能满足性能要求,又能节省资源,灵活性强。
附图说明
图1为本发明的内插滤波处理方法的处理流程框图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步详细的描述:
如图1所示的本发明处理流程框图所示,本发明的任意采样率转换的内插滤波处理方法,包括以下步骤:
(1)、根据输入信号的带宽Bw和设定的信号分辨率参数K,生成K组FIR滤波器系数,并将所述K组滤波器系数保存在FIR滤波器抽头系数表内;其中,所述信号分辨率参数K=2L,L为正整数;具体实现步骤如下:
(1a)、根据信号带宽Bw,生成一组满足所述信号带宽内滤波要求的FIR滤波器的抽头系数h(n),其中,n=1,2,…,N,N为所述FIR滤波器的抽头个数;
(1b)、根据设定的精度系数K得到K个比例函数,其中,第k个所述比例函数为k=1,2,…,K,n=1,2,…,N;
(1c)、对步骤(1a)得到的FIR滤波器抽头系数h(n)进行N点离散傅里叶变换,得到FIR滤波器响应函数H(jw)=DFT(h(n)),其中,DFT()代表离散傅里叶变换;所述响应函数H(jw)的N个值分别为h1,h2,…,hn,…,hN;
(1d)、根据步骤(1b)得到的K个比例函数和步骤(1c)得到的响应函数H(jw),得到K个新的响应函数,其中,第k个所述新的响应函数为H′k(jw),H′k(jw)的N个值分别为h1′,h2′,…,hn′,…,hN′,其中,hn′=hn×fk(n),k=1,2,…,K,n=1,2,…,N;
(1e)、对步骤(1d)得到的K个新的响应函数进行逆傅里叶变换得K组滤波器系数,其中,第k组所述滤波器系数为h′k(n)=IDFT(H′k(jw)),IDFT代表逆离散傅里叶变换,k=1,2,…,K;
(1f)、将步骤(1e)中得到K组滤波器系数依次存入数据表中,得到FIR滤波器抽头系数表,该FIR滤波器抽头系数表可以装订在ROM表,便于硬件实现,可以与横向FIR滤波器联合使用,实现可变采样率的插值滤波处理。
在本实施例中,FIR滤波器设定为32阶的横向滤波器,即具有32个延迟单元,共N=33个抽头系数;且设定的信号分辨率参数K=256,即正整数L=8;
(2)、根据输入信号的码率Rb和***时钟频率fs,采用NCO产生1倍数据时钟信号,所述数据时钟信号的频率Rs=Rb,所述数据时钟信号的时钟相位Phase等于所述NCO累加值的高L位截取值,即所述时钟相位Phase的取值范围为0~K-1;其中,所述NCO的位宽为M且M≥L;该NCO的工作频率为fs;
本实施例中,NCO的位宽为M=32,在设定的输入信号分辨率参数K=256时,时钟相位phase为NCO累加值的高8位截取值,也就是说时钟相位的取值为0~255;在本发明中,就是采用该时钟相位得到FIR滤波器系数的提取序号,实现插值滤波的系数提取。
(3)、按照步骤(2)产生的1倍数据时钟信号对输入信号进行接收,得到接收信号;即该接收信号的采样率为Rs,本发明的目的就是将该信号的采样率统一转换为***采样率,并且保证信号不失真;
(4)、根据信号接收时钟的相位Phase的取值,在步骤(1)产生的FIR滤波器抽头系数表中提取FIR滤波器的抽头系数,其中,如果所述时钟相位Phase=k-1,则在所述FIR滤波器抽头系数表中提取得到第k组低通FIR滤波器系数h′k(n),其中,k=1,2,…,K,n=1,2,…,N,N为所述FIR滤波器的抽头个数;
(5)、将步骤(4)提取得到的FIR滤波器系数应用于FIR滤波器中,对步骤(3)中得到的接收信号进行内插滤波处理,得到内插滤波后信号,该信号的采样率为fs,即经过以上处理将采样率为Rs的信号内插滤波为采样率为fs的信号,从而实现采样率的转换。
在本实施例中,本发明的插值滤波器由一个32阶横向FIR滤波器和一个装订有滤波器系数的ROM表组成。在每个***工作时钟fs周期内,根据量化后的时钟相位phase,在ROM表中查找滤波器系数,实时更新FIR滤波器系数并进行滤波运算,得到采样率为fs的插值后的信号。如果对每个滤波系数按16位进行量化,时钟相位phase按8位进行量化,则装订FIR滤波器系数的ROM表的存储量为33*16*28=135168bit。
本发明采用可变系数的分数倍插值滤波器,实现不同速率信号插值滤波,整个功能可采用FPGA实现。
以上所述,仅为本发明最佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。
本发明说明书中未作详细描述的内容属于本领域专业技术人员公知技术。
Claims (2)
1.一种任意采样率转换的内插滤波处理方法,其特征在于包括以下步骤:
(1)、根据输入信号的带宽Bw和设定的信号分辨率参数K,生成K组FIR滤波器系数,并将所述K组滤波器系数保存在FIR滤波器抽头系数表内;其中,所述信号分辨率参数K=2L,L为正整数;具体实现过程如下:
(1a)、根据信号带宽Bw,生成一组满足所述信号带宽内滤波要求的FIR滤波器的抽头系数h(n),其中,n=1,2,…,N,N为所述FIR滤波器的抽头个数;
(1b)、根据设定的信号分辨率参数K得到K个比例函数,其中,第k个所述比例函数为
(1c)、对步骤(1a)得到的FIR滤波器抽头系数h(n)进行N点离散傅里叶变换,得到FIR滤波器响应函数H(jw)=DFT(h(n)),其中,DFT()代表离散傅里叶变换;所述响应函数H(jw)的N个值分别为h1,h2,…,hn,…,hN;
(1d)、根据步骤(1b)得到的K个比例函数和步骤(1c)得到的响应函数H(jw),得到K个新的响应函数,其中,第k个所述新的响应函数为H′k(jw),Hk′(jw)的N个值分别为h′1,h′2,…,hn′,…,hN′,其中,hn′=hn×fk(n),k=1,2,…,K,n=1,2,…,N;
(1e)、对步骤(1d)得到的K个新的响应函数进行逆傅里叶变换得K组滤波器系数,其中,第k组所述滤波器系数为h′k(n)=IDFT(H′k(jw)),IDFT代表逆离散傅里叶变换,k=1,2,…,K;
(1f)、将步骤(1e)中得到K组滤波器系数依次存入数据表中,得到FIR滤波器抽头系数表;
(2)、根据输入信号的码率Rb和***时钟频率fs,采用NCO产生1倍数据时钟信号,所述数据时钟信号的频率Rs=Rb,所述数据时钟信号的时钟相位Phase等于NCO累加值的高L位截取值,即所述时钟相位Phase的取值范围为0~K-1;其中,所述NCO的位宽为M且M≥L;
(3)、按照步骤(2)产生的1倍数据时钟信号对输入信号进行接收,得到接收信号;
(4)、根据信号接收时钟的时钟相位Phase的取值,在步骤(1)产生的FIR滤波器抽头系数表中提取FIR滤波器的抽头系数,其中,如果所述时钟相位Phase=k-1,则在所述FIR滤波器抽头系数表中提取得到第k组低通FIR滤波器系数h′k(n),其中,k=1,2,…,K,n=1,2,…,N,N为所述FIR滤波器的抽头个数;
(5)、将步骤(4)提取得到的FIR滤波器系数应用于FIR滤波器中,对步骤(3)中得到的接收信号进行内插滤波处理,得到内插滤波后信号。
2.根据权利要求1所述的一种任意采样率转换的内插滤波处理方法,其特征在于:FIR滤波器抽头系数表装订在ROM表中。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410745446.7A CN104539261B (zh) | 2014-12-08 | 2014-12-08 | 一种任意采样率转换的内插滤波处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410745446.7A CN104539261B (zh) | 2014-12-08 | 2014-12-08 | 一种任意采样率转换的内插滤波处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104539261A CN104539261A (zh) | 2015-04-22 |
CN104539261B true CN104539261B (zh) | 2017-08-29 |
Family
ID=52854748
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410745446.7A Active CN104539261B (zh) | 2014-12-08 | 2014-12-08 | 一种任意采样率转换的内插滤波处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104539261B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105306008A (zh) * | 2015-11-13 | 2016-02-03 | 北京遥测技术研究所 | 一种任意采样率转换的内插滤波处理方法 |
CN108871385B (zh) * | 2017-05-12 | 2021-09-07 | 西门子公司 | 编码器、电机、编码器数据处理方法及存储介质 |
CN108366320A (zh) * | 2018-01-08 | 2018-08-03 | 联创汽车电子有限公司 | 车载前馈主动降噪*** |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1192609A (zh) * | 1997-01-10 | 1998-09-09 | 哈里公司 | 重新抽样电路与相关方法 |
CN1226888C (zh) * | 2001-09-07 | 2005-11-09 | 三星Thales株式会社 | 使用查找表的块内插滤波器结构 |
CN103650445A (zh) * | 2011-03-09 | 2014-03-19 | 联邦科学技术研究组织 | 用于通信***的任意采样率转换 |
-
2014
- 2014-12-08 CN CN201410745446.7A patent/CN104539261B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1192609A (zh) * | 1997-01-10 | 1998-09-09 | 哈里公司 | 重新抽样电路与相关方法 |
CN1226888C (zh) * | 2001-09-07 | 2005-11-09 | 三星Thales株式会社 | 使用查找表的块内插滤波器结构 |
CN103650445A (zh) * | 2011-03-09 | 2014-03-19 | 联邦科学技术研究组织 | 用于通信***的任意采样率转换 |
Non-Patent Citations (3)
Title |
---|
"一种用于符号同步的奈奎斯特滤波器优化设计";叶双应 等;《电子器件》;20060630;第29卷(第2期);第493-496页 * |
"基于FPGA的EDSL***QAM解调器的研究与实现";杨霖;《中国优秀硕士学位论文全文数据库 信息科技辑》;20091115(第11期);第四章第4.2-4.3节 * |
"基于Gardner算法的位定时同步研究";朱雪阳;《中国优秀硕士学位论文全文数据库 信息科技辑》;20100915(第09期);第3章 * |
Also Published As
Publication number | Publication date |
---|---|
CN104539261A (zh) | 2015-04-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104539262B (zh) | 一种连续可变速率的数字成型滤波处理方法 | |
CN104539261B (zh) | 一种任意采样率转换的内插滤波处理方法 | |
CN101567701B (zh) | 一种高效多路数字下变频器*** | |
CN104320088B (zh) | 一种数字下变频电路 | |
CN104467739B (zh) | 一种带宽、中心频点可调的数字滤波器及其实现方法 | |
CN101257482A (zh) | 数字基带可变速率转换调制***的实现方法和实现装置 | |
CN101719769A (zh) | 一种任意抽取滤波装置 | |
CN105281707B (zh) | 一种动态可重构滤波器组的实现方法 | |
CN101478525A (zh) | 一种多载波分离的方法及多载波分离装置 | |
CN105337587B (zh) | 一种基于dft的非最大抽取***综合滤波器组构造方法 | |
CN103117972A (zh) | 一种矢量信号分析方法和装置 | |
CN104393854A (zh) | 基于fpga的时分复用级联积分梳状抽取滤波器及其实现方法 | |
CN106059530A (zh) | 一种频率响应与系数量化位数弱相关的半带滤波器结构 | |
CN110365351B (zh) | 一种基于滤波器复用架构实现信道化接收的方法及装置 | |
CN104883157B (zh) | 一种可变子带数字滤波器 | |
CN106849904A (zh) | 数字滤波设备 | |
CN102694563A (zh) | 一种用于声音信号的数字解调下变频***及方法 | |
CN104301085B (zh) | 一种资源复用的完全重构子带分析器 | |
CN102891662B (zh) | 一种通用的速率下变换、上变换装置及方法 | |
CN103956987A (zh) | 一种数字滤波器设计方法 | |
CN106972833B (zh) | 一种可任意倍数重采样的数字上变频器 | |
CN106130581A (zh) | 一种多相滤波宽带数字信道化接收机改进*** | |
CN105306008A (zh) | 一种任意采样率转换的内插滤波处理方法 | |
CN103001586A (zh) | 一种宽带双通道数字下变频器 | |
CN209267558U (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 |