[go: up one dir, main page]

CN105634439A - 一种异步成形滤波器设计方法 - Google Patents

一种异步成形滤波器设计方法 Download PDF

Info

Publication number
CN105634439A
CN105634439A CN201510971094.1A CN201510971094A CN105634439A CN 105634439 A CN105634439 A CN 105634439A CN 201510971094 A CN201510971094 A CN 201510971094A CN 105634439 A CN105634439 A CN 105634439A
Authority
CN
China
Prior art keywords
filter
shaping
shaping filter
representing
coefficient
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.)
Granted
Application number
CN201510971094.1A
Other languages
English (en)
Other versions
CN105634439B (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.)
Nanjing Panda Handa Technology Co Ltd
Original Assignee
Panda Electronics Group Co Ltd
Nanjing Panda Handa Technology Co Ltd
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 Panda Electronics Group Co Ltd, Nanjing Panda Handa Technology Co Ltd filed Critical Panda Electronics Group Co Ltd
Priority to CN201510971094.1A priority Critical patent/CN105634439B/zh
Publication of CN105634439A publication Critical patent/CN105634439A/zh
Application granted granted Critical
Publication of CN105634439B publication Critical patent/CN105634439B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H2017/0072Theoretical filter design

Landscapes

  • Complex Calculations (AREA)
  • Synchronisation In Digital Transmission Systems (AREA)

Abstract

本发明公开了一种异步成形滤波器设计方法,包括查找表生成模块、NCO控制模块和成形滤波器模块,其中,查找表生成模块用于产生成形滤波器系数并制作查找表,NCO控制模块用于产生输入信号速率及样点相位信息,成形滤波器模块用于实现从输入信号速率到基带处理时钟的成形输出。本发明对基带处理时钟和所处理信号速率不受限于整数倍关系,具备较强的通用性;且方法实现简单,时序控制简单,占用资源较小,易于工程实现和模块化处理,便于移植。

Description

一种异步成形滤波器设计方法
技术领域
本发明涉及一种通用异步成形滤波器设计方法,属于无线通信物理层技术领域。
背景技术
在数字传输系统中,系统的传输带宽的有限的,而信号的频谱宽度是无限的,在有限的带宽中传输带宽无限的信号,会造成信号时域波形的失真,引起码间串扰,为了有效利用频谱、消除码间串扰,需要对发送信号进行成形滤波。
传统的成形滤波器设计通常采用整数倍成形方式,即滤波输出采样速率为输入信号速率的整数倍,这样设计的优点是易于硬件实现,缺点是要求基带处理时钟必须是所处理信号速率的整数倍,而通常情况下,所给时钟源与信号速率都不具备整数倍关系,需要通过外接数字频率合成器设备来达到分数倍变频的目的,增加了硬件开销。由此可见,设计一种异步成形滤波器,使得对于任意速率的基带信号,经过成形滤波器后输出的采样速率都为基带处理时钟,是非常有必要的。
中国专利CN204258746和CN1042189提出了一种分数倍内插成型滤波器实现方法,该滤波器包括数据同步模块、分段滤波模块、重采样模块、数控振荡器控制模块。数控振荡器控制模块根据基带调制速率输入产生慢溢出标志送给数据同步模块产生慢同步数据输出,慢同步数据输出经过分段滤波模块产生滤波器组输出,滤波器组输出和快溢出标志和分数时延滤波器时延经过重采样模块产生成型滤波数据。该发明通过在整数倍成型滤波器后接重采样模块实现分数倍成型滤波输出,需要的乘法器和加法器数量较多,消耗资源量大,且时序控制困难,时延较大。
发明内容
本发明目的是,针对常规成形滤波器设计只能实现整数倍成形,应用范围受限,增加硬件开销等问题,提出了一种异步成形滤波器设计方法,可实现任意输入信号速率到基带处理时钟的成形滤波输出,不受限于整数倍关系,可在任意基带处理时钟,任意信号速率下使用,具备较强的通用性,当系统时钟源与信号速率不具备整数倍关系时,节省了硬件开销;方法实现简单,工程实现复杂度低,对于可编程逻辑器件等开发平台,只需要一个NCO控制模块和成形滤波器即可实现异步成形输出,无需通过重采样模块实现异步时序,时序控制简单,占用资源较小;接口简单,易于模块化实现,便于移植。
本发明为解决上述技术问题提出的技术方案是:
一种异步成形滤波器设计方法,包括查找表生成模块、NCO控制模块和成形滤波器模块,其中,查找表生成模块用于产生成形滤波器系数并制作查找表,NCO控制模块用于产生输入信号速率及样点相位信息,成形滤波器模块用于实现从输入信号速率到基带处理时钟的成形输出,包括以下步骤:
步骤1,根据给定的滤波器参数通过查找表生成模块生成成形滤波器的成形系数,同时查找表生成模块根据成形滤波器的成形系数制作查找表;
步骤2,根据输入信号速率通过NCO控制模块产生输入信号速率时钟及样点相位信息;
步骤3,根据输入数据,通过成形滤波器模块实现从输入信号速率到基带处理时钟的成形输出。
所述步骤1中根据给定的滤波器参数通过查找表生成模块生成成形滤波器,同时查找表生成模块根据成形滤波器的成形系数制作查找表的方法,如下:
步骤11,根据基带处理时钟频率fm,输入信号速率fs,得到成形滤波器的成形倍数R:
R=2K,
其中,表示向上取整,K表示成形倍数R的位宽,fm表示基带处理时钟频率,fs表示输入信号速率。根据调制器系统给定滚降系数α和滤波器相关长度NT,以及成形倍数R即可制作成形滤波器;
步骤12,根据给定的量化精度对步骤11得到的成形滤波器进行成形系数量化,将浮点型成形系数量化成合适位宽的整数;
步骤13,对步骤12中的成形滤波器的成形系数进行截短,得到截短后的滤波器系数F:
步骤14,根据滤波器相关长度NT和步骤11得到的成形倍数R对步骤13中得到的截短后的滤波器系数F进行分组,生成成形滤波器系数矩阵:
f j i = F ( j - 1 ) R + i , i = 1 , 2 , ... , R j = 1 , 2 , ... , N T ;
其中,表示系数矩阵中的第i行第j列元素,R表示成形倍数,NT表示滤波器相关长度。
步骤15,根据步骤14得到的成形滤波器系数矩阵获取每列最大量化位宽:
其中,Mj表示第j列最大量化位宽,符号表示向上取整,|·|表示取绝对值,表示i=1,2,…,R个元素中的最大值,表示系数矩阵中的第i行第j列元素;
步骤16,按照步骤15得到的每列最大量化位宽,将步骤14得到的成形滤波器系数矩阵中的系数表示为二进制数,若则将其转换为补码形式;将各行对应列合并,生成二进制数组:
其中,Bi表示数组中的第i个元素,表示将转换为位宽为Mj的二进制数,&表示将比特位并置,元素Bi的位宽M为Mj表示第j列最大量化位宽,R表示数组长度;
步骤17,将步骤16生成的二进制数组按照所用开发系统规定的查找表格式生成二进制查找表:
所述步骤2中根据输入信号速率通过NCO控制模块产生输入信号速率时钟及样点相位信息的方法如下:
步骤21,根据基带处理时钟频率fm、输入信号速率fs生成频率控制字
其中,表示频率控制字,fm表示基带处理时钟频率、fs表示输入信号速率,N为NCO控制模块的累加精度;
步骤22,根据步骤21得到的频率控制字通过频率积分实现瞬时相位输出:
其中,表示第n个基带处理时钟nTm时的瞬时相位,mod表示取余函数,表示频率控制字,N为NCO控制模块的累加精度;
步骤23,根据步骤22得到的瞬时相位通过NCO控制模块溢出标志产生输入信号速率时钟的时钟使能信号:
时,NCO控制模块中的瞬时值溢出,提取溢出标志,将该溢出标志作为输入信号速率时钟的时钟使能信号;
步骤24,通过NCO控制模块高位获取需要的成形样点相位信息:
根据步骤11得到的成形倍数R的位宽K,得到第n个基带处理时钟时的成形样点相位Pn
其中,Pn表示第n个基带处理时钟时的成形样点相位,表示第n个基带处理时钟nTm时的瞬时相位,表示提取位宽为N的二进制数的第a至b位。
所述步骤3中根据输入数据,通过成形滤波器模块实现从输入信号速率到基带处理时钟的成形输出的方法如下:
步骤31,设置深度为2NT的数据寄存器,对输入数据进行缓存:
数据寄存器采用FIFO结构,利用NCO控制模块给出的输入信号速率时钟的时钟使能信号,在时钟使能信号有效时读取输入数据并更新寄存器,该数据寄存器的深度为2NT
步骤32,根据NCO控制模块给出的成形样点相位信息和生成的查找表,提取2NT个成形滤波器系数:
使用步骤1中生成的二进制查找表和步骤24中NCO控制模块给出的成形样点相位信息通过以下方法获取成形滤波器的全部系数,将第n个基带处理时钟时的读取地址表示为:
{ A 1 = P n A 2 = n o t P n ;
式中,A1表示查找表的读取地址1,用以提取前半段的成形滤波器系数,A2表示查找表的读取地址2,用以提取后半段的成形滤波器系数。Pn表示第n个基带处理时钟时的成形样点相位,not表示按位取反;
则可得到成形滤波器的第k个系数Ck
C k = { ( T A 1 ) M ( M - &Sigma; m = 0 k - 1 M m - 1 ) &RightArrow; ( M - &Sigma; m = 0 k - 1 M m ) 1 &le; k &le; N T ( T A 2 ) M ( &Sigma; m = 2 N T - k + 1 N T + 1 M m - 1 ) &RightArrow; ( &Sigma; m = 2 N T - k + 2 N T + 1 M m ) N T < k &le; 2 N T ;
式中,Ck表示成形滤波器的第k个系数,A1表示查找表的读取地址1,A2表示查找表的读取地址2,TA表示查找表中地址A对应的元素,M表示元素Bi的位宽,NT表示滤波器相关长度,M0等于0,表示提取位宽为M的二进制数的第a至b位;
步骤33,将缓存的数据与成形滤器系数对应位置相乘;
对于单比特输入数据,直接根据输入比特将成形滤波器系数取正或取负即可;对于多进制输入数据,需要2NT个乘法器;
步骤34,根据步骤32得到的成形滤波器的第k个系数Ck,可得到第n+nd个基带处理时钟时成形输出数据:
y ( n + n d ) = &Sigma; k = 1 2 N T C k &CenterDot; D k ;
其中,y(n+nd)表示第n+nd个基带处理时钟时的成形输出数据,nd表示由乘法器和加法器造成的时延,Ck表示成形滤波器的第k个系数,Dk表示在第n个基带处理时钟时数据寄存器中的第k个数据;所需加法器个数为2NT-1。
优选的:所述步骤11中成形滤波器为平方根升余弦滤波器,该平方根升余弦滤波器阶数为2RNT+1。
优选的:所述步骤13中对成形滤波器的成形系数的右半部分系数进行截去,进而得到截短后的成形系数F:
F = F 1 , F 2 , ... , F RN T ;
其中,F为截短后的滤波器系数向量,为相应的滤波器系数值。
优选的:所述步骤21中的NCO控制模块的累加精度N大于步骤11中的成形倍数R的位宽K。
优选的:所述步骤32中成形滤波器模块读取地址A1和A2进行延时处理,其中,A1为查找表读取地址1,A2为查找表读取地址2。
本发明的一种异步成形滤波器设计方法,相比现有技术,具有以下有益效果:
1,对基带处理时钟和所处理信号速率关系不做要求,不受限于整数倍关系。节省了硬件开销,无需增加额外的时钟模块来产生分数倍的时钟。
2,优化了查找表的设计,合理分配资源,节约了存储量。满足多通路多速率的信号系统处理要求,可实现多路速率不同的扩频通信信号以同一采样速率成形并复接。
3,对于时变速率的基带信号可采用统一的成形滤波结构,降低了软硬件设计的复杂度。
4,本发明只需一个异步成形滤波器即可,无需通过重采样模块实现异步时序,方法实现简单,所需模块较少,工程实现复杂度低。
5,本发明所需的乘法器、加法器数量明显减小,占用资源较少,无需复杂的时序控制逻辑,利于工程实现。
6,本发明接口更加简单明了,易于模块化实现,便于移植。
附图说明
图1为本发明通用异步成形滤波器实现框图。
图2为本发明查找表生成模块具体实现框图。
图3为本发明NCO控制模块具体实现框图。
图4为本发明成形滤波器模块具体实现框图。
具体实施方式
附图非限制性地公开了本发明一个优选实施例的结构示意图,以下将结合附图详细地说明本发明的技术方案。
实施例
本实施例的一种异步成形滤波器设计方法,如图1所示,包括查找表生成模块、NCO控制模块和成形滤波器模块,其中,查找表生成模块用于产生成形滤波器系数并制作查找表,NCO控制模块用于产生输入信号速率及样点相位信息,成形滤波器模块用于实现从输入信号速率到基带处理时钟的成形输出,具体包括以下步骤。
步骤1,根据给定的滤波器参数通过查找表生成模块生成成形滤波器的成形系数,同时查找表生成模块根据成形滤波器的成形系数制作查找表。如图2所示,其具体实施过程如下:
步骤11,针对给定的滤波器参数,生成成形滤波器
成形滤波器通常采用平方根升余弦滤波器,该平方根升余弦滤波器阶数为2RNT+1。根据调制器系统给定滚降系数α以及滤波器相关长度NT,基带处理时钟频率fm,输入信号速率fs,得到成形滤波器的成形倍数R:在设计滤波器时,滚降系数α一般由调制器系统给定,成型倍数R一般选为2的幂次方,在本发明中,R的选择要求大于基带处理时钟频率与输入信号速率之比,则R满足
R=2K,
其中,表示向上取整,K表示成形倍数R的位宽,fm表示基带处理时钟频率,fs表示输入信号速率。滤波器相关长度NT可任意选取,NT越大,输出波形越接近理论平方根升余弦波形,但所需的乘法器和加法器也越多,设计时可根据实际应用需求和硬件成本进行折中选择。根据以上假设,得到的滤波器阶数为2RNT+1。
步骤12,根据给定的量化精度对步骤11得到的成形滤波器进行成形系数量化
根据系统精度及位宽要求,将浮点型成形滤波器系数量化成合适(对应)的位宽。
步骤13,对步骤12中的成形滤波器的成形系数进行截短,得到截短后的滤波器系数F
由滤波器的性质可知,系数是左右对称的,为节约存储量,截去右半部分的系数。进而得到截短后的成形系数F:
F = F 1 , F 2 , ... , F RN T ;
其中,F为截短后的滤波器系数向量;为相应的滤波器系数值。
步骤14,生成成形滤波器系数矩阵:
按照所设的成形倍数和相关长度,对滤波器系数进行分组,生成滤波器系数矩阵,根据滤波器相关长度NT和步骤11得到的成形倍数R对步骤13中得到的截短后的滤波器系数F进行分组,生成成形滤波器系数矩阵:
F j i = F ( j - 1 ) R + i ; i = 1 , 2 , ... , R j = 1 , 2 , ... , N T ;
其中,表示系数矩阵中的第i行第j列元素,R表示成形倍数,NT表示滤波器相关长度。
步骤15,获取每列最大量化位宽
根据步骤14得到的成形滤波器系数矩阵获取每列最大量化位宽:
其中,Mj表示第j列最大量化位宽,符号表示向上取整,|·|表示取绝对值,表示i=1,2,…,R个元素中的最大值,表示系数矩阵中的第i行第j列元素。
步骤16,生成二进制数组
按照步骤15得到的每列最大量化位宽,将步骤14得到的成形滤波器系数矩阵中的系数表示为二进制数,若则将其转换为补码形式;将各行对应列合并,生成二进制数组:
其中,Bi表示数组中的第i个元素,表示将转换为位宽为Mj的二进制数,&表示将比特位并置,元素Bi的位宽M为Mj表示第j列最大量化位宽,R表示数组长度。
步骤17,生成查找表
按照所用开发系统规定的查找表格式,生成二进制查找表,因此将步骤16生成的二进制数组按照所用开发系统规定的查找表格式生成二进制查找表:
步骤2,根据输入信号速率通过NCO控制模块产生输入信号速率时钟及样点相位信息;如图3所示,其具体实施过程如下:
步骤21,根据基带处理时钟频率fm、输入信号速率fs生成频率控制字
其中,表示频率控制字,fm表示基带处理时钟频率、fs表示输入信号速率,N为NCO控制模块的累加精度;N为NCO的累加精度,可任意选取,N越大,得到的时钟越精确,需要的存储量也越大,设计时可根据实际应用需求进行选择。本发明中,N的选择要求大于步骤11中的K,其余不做限定。
步骤22,通过频率积分实现瞬时相位输出:
根据步骤21得到的频率控制字通过频率积分实现瞬时相位输出:
其中,表示第n个基带处理时钟nTm时的瞬时相位,mod表示取余函数,表示频率控制字,N为NCO控制模块的累加精度;
步骤23,通过NCO溢出标志产生输入信号速率时钟:
根据步骤22得到的瞬时相位通过NCO控制模块溢出标志产生输入信号速率时钟的时钟使能信号:
时,NCO控制模块中的瞬时值溢出,提取溢出标志,将该溢出标志作为输入信号速率时钟的时钟使能信号;
步骤24,通过NCO控制模块高位获取需要的成形样点相位信息:
根据步骤11得到的成形倍数R的位宽K,得到第n个基带处理时钟时的成形样点相位Pn
其中,Pn表示第n个基带处理时钟时的成形样点相位,表示第n个基带处理时钟nTm时的瞬时相位,表示提取位宽为N的二进制数的第a至b位。
步骤3,根据输入数据,通过成形滤波器模块实现从输入信号速率到基带处理时钟的成形输出,如图4所示,其具体实施过程如下:
步骤31,设置深度为2NT的数据寄存器,对输入数据进行缓存:
数据寄存器采用FIFO结构,利用NCO控制模块给出的输入信号速率时钟的时钟使能信号,在时钟使能信号有效时读取输入数据并更新寄存器,该数据寄存器的深度为2NT
步骤32,根据NCO控制模块给出的成形样点相位信息和生成的查找表,提取2NT个成形滤波器系数:
在查找表设计过程中,截去了滤波器的右半部分,使用时可通过NCO给出的成形样点相位信息获取成形滤波器的全部系数,使用步骤1中生成的二进制查找表通过步骤24中NCO控制模块给出的成形样点相位信息获取成形滤波器的全部系数,将第n个基带处理时钟时的读取地址表示为:
A 1 = P n A 2 = n o t P n ;
式中,A1表示查找表的读取地址1,用以提取前半段的成形滤波器系数,A2表示查找表的读取地址2,用以提取后半段的成形滤波器系数。Pn表示第n个基带处理时钟时的成形样点相位,not表示按位取反。
成形滤波器的第k个系数Ck
C k = { ( T A 1 ) M ( M - &Sigma; m = 0 k - 1 M m - 1 ) &RightArrow; ( M - &Sigma; m = 0 k - 1 M m ) 1 &le; k &le; N T ( T A 2 ) M ( &Sigma; m = 2 N T - k + 1 N T + 1 M m - 1 ) &RightArrow; ( &Sigma; m = 2 N T - k + 2 N T + 1 M m ) N T < k &le; 2 N T ;
式中,Ck表示成形滤波器的第k个系数,A1表示查找表的读取地址1,A2表示查找表的读取地址2,TA表示查找表中地址A对应的元素,M表示元素Bi的位宽,NT表示滤波器相关长度,M0等于0,表示提取位宽为M的二进制数的第a至b位。
需要注意的是,通过NCO产生数据速率时钟和通过数据速率时钟更新数据寄存器都存在一定的延时,为保证滤波器严格与数据对齐,需将读取地址A1和A2进行延时处理。
步骤33,将缓存的数据与成形滤器系数对应位置相乘;
对于单比特输入数据,直接根据输入比特将成形滤波器系数取正或取负即可;对于多进制输入数据,需要2NT个乘法器;
步骤34,根据步骤32得到的成形滤波器的第k个系数Ck,可得到第n+nd个基带处理时钟时成形输出数据:
y ( n + n d ) = &Sigma; i = 1 2 N T C k &CenterDot; D k ;
其中,y(n+nd)表示第n+nd个基带处理时钟时的成形输出数据,nd表示由乘法器和加法器造成的时延,Ck表示成形滤波器的第k个系数,Dk表示在第n个基带处理时钟时数据寄存器中的第k个数据;所需加法器个数为2NT-1。
由此,通过2NT个乘法器和2NT-1个加法器实现了从任意输入信号速率到基带处理时钟的成形滤波输出,不受限于整数倍关系,具备较强的通用性;且方法实现简单,时序控制简单,占用资源较小,易于工程实现和模块化处理,便于移植。
上面结合附图所描述的本发明优选具体实施例仅用于说明本发明的实施方式,而不是作为对前述发明目的和所附权利要求内容和范围的限制,凡是依据本发明的技术实质对以上实施例所做的任何简单修改、等同变化与修饰,均仍属本发明技术和权利保护范畴。

Claims (8)

1.一种异步成形滤波器设计方法,包括查找表生成模块、NCO控制模块和成形滤波器模块,其中,查找表生成模块用于产生成形滤波器系数并制作查找表,NCO控制模块用于产生输入信号速率及样点相位信息,成形滤波器模块用于实现从输入信号速率到基带处理时钟的成形输出,其特征在于,包括以下步骤:
步骤1,根据给定的滤波器参数通过查找表生成模块生成成形滤波器的成形系数,同时查找表生成模块根据成形滤波器的成形系数制作查找表;
步骤2,根据输入信号速率通过NCO控制模块产生输入信号速率时钟及样点相位信息;
步骤3,根据输入数据,通过成形滤波器模块实现从输入信号速率到基带处理时钟的成形输出。
2.根据权利要求1所述异步成形滤波器设计方法,其特征在于:所述步骤1中根据给定的滤波器参数通过查找表生成模块生成成形滤波器,同时查找表生成模块根据成形滤波器的成形系数制作查找表的方法,如下:
步骤11,根据基带处理时钟频率fm,输入信号速率fs,得到成形滤波器的成形倍数R:
其中,表示向上取整,K表示成形倍数R的位宽,fm表示基带处理时钟频率,fs表示输入信号速率;根据调制器系统给定滚降系数α和滤波器相关长度NT,以及成形倍数R即可制作成形滤波器;
步骤12,根据给定的量化精度对步骤11得到的成形滤波器进行成形系数量化,将浮点型成形系数量化成对应位宽的整数;
步骤13,对步骤12中的成形滤波器的成形系数进行截短,得到截短后的滤波器系数F:
步骤14,根据滤波器相关长度NT和步骤11得到的成形倍数R对步骤13中得到的截短后的滤波器系数F进行分组,生成成形滤波器系数矩阵:
F j i = F ( j - 1 ) R + i , i = 1 , 2 , ... , R j = 1 , 2 , ... , N T ;
其中,表示系数矩阵中的第i行第j列元素,R表示成形倍数,NT表示滤波器相关长度;
步骤15,根据步骤14得到的成形滤波器系数矩阵获取每列最大量化位宽:
其中,Mj表示第j列最大量化位宽,符号表示向上取整,|·|表示取绝对值,表示i=1,2,…,R个元素中的最大值,表示系数矩阵中的第i行第j列元素;
步骤16,按照步骤15得到的每列最大量化位宽,将步骤14得到的成形滤波器系数矩阵中的系数表示为二进制数,若则将其转换为补码形式;将各行对应列合并,生成二进制数组:
其中,Bi表示数组中的第i个元素,表示将转换为位宽为Mj的二进制数,&表示将比特位并置,元素Bi的位宽M为Mj表示第j列最大量化位宽,R表示数组长度;
步骤17,将步骤16生成的二进制数组按照所用开发系统规定的查找表格式生成二进制查找表:
3.根据权利要求2所述异步成形滤波器设计方法,其特征在于:所述步骤2中根据输入信号速率通过NCO控制模块产生输入信号速率时钟及样点相位信息的方法如下:
步骤21,根据基带处理时钟频率fm、输入信号速率fs生成频率控制字
其中,表示频率控制字,fm表示基带处理时钟频率、fs表示输入信号速率,N为NCO控制模块的累加精度;
步骤22,根据步骤21得到的频率控制字通过频率积分实现瞬时相位输出:
其中,表示第n个基带处理时钟nTm时的瞬时相位,mod表示取余函数,表示频率控制字,N为NCO控制模块的累加精度;
步骤23,根据步骤22得到的瞬时相位通过NCO控制模块溢出标志产生输入信号速率时钟的时钟使能信号:
时,NCO控制模块中的瞬时值溢出,提取溢出标志,将该溢出标志作为输入信号速率时钟的时钟使能信号;
步骤24,通过NCO控制模块高位获取需要的成形样点相位信息:
根据步骤11得到的成形倍数R的位宽K,得到第n个基带处理时钟时的成形样点相位Pn
其中,Pn表示第n个基带处理时钟时的成形样点相位,表示第n个基带处理时钟nTm时的瞬时相位,表示提取位宽为N的二进制数的第a至b位,N为NCO控制模块的累加精度。
4.根据权利要求3所述异步成形滤波器设计方法,其特征在于:所述步骤3中根据输入数据,通过成形滤波器模块实现从输入信号速率到基带处理时钟的成形输出的方法如下:
步骤31,设置深度为2NT的数据寄存器,对输入数据进行缓存:
数据寄存器采用FIFO结构,利用NCO控制模块给出的输入信号速率时钟的时钟使能信号,在时钟使能信号有效时读取输入数据并更新寄存器,该数据寄存器的深度为2NT
步骤32,根据NCO控制模块给出的成形样点相位信息和生成的查找表,提取2NT个成形滤波器系数:
使用步骤1中生成的二进制查找表和步骤24中NCO控制模块给出的成形样点相位信息通过以下方法获取成形滤波器的全部系数,将第n个基带处理时钟时的读取地址表示为:
A 1 = P n A 2 = notP n ;
式中,A1表示查找表的读取地址1,用以提取前半段的成形滤波器系数,A2表示查找表的读取地址2,用以提取后半段的成形滤波器系数;Pn表示第n个基带处理时钟时的成形样点相位,not表示按位取反;
则可得到成形滤波器的第k个系数Ck
C k = ( T A 1 ) M ( M - &Sigma; m = 0 k - 1 M m - 1 ) &RightArrow; ( M - &Sigma; m = 0 k M m ) 1 &le; k &le; N T ( T A 2 ) M ( &Sigma; m = 2 N T - k + 1 N T + 1 M m - 1 ) &RightArrow; ( &Sigma; m = 2 N T - k + 2 N T + 1 M m ) N T &le; k &le; 2 N T ;
式中,Ck表示成形滤波器的第k个系数,A1表示查找表的读取地址1,A2表示查找表的读取地址2,TA表示查找表中地址A对应的元素,M表示元素Bi的位宽,NT表示滤波器相关长度,M0等于0,表示提取位宽为M的二进制数的第a至b位;
步骤33,将缓存的数据与成形滤器系数对应位置相乘;
对于单比特输入数据,直接根据输入比特将成形滤波器系数取正或取负即可;对于多进制输入数据,需要2NT个乘法器;
步骤34,根据步骤32得到的成形滤波器的第k个系数Ck,可得到第n+nd个基带处理时钟时成形输出数据:
y ( n + n d ) = &Sigma; k = 1 2 N T C k &CenterDot; D k ;
其中,y(n+nd)表示第n+nd个基带处理时钟时的成形输出数据,nd表示由乘法器和加法器造成的时延,Ck表示成形滤波器的第k个系数,Dk表示在第n个基带处理时钟时数据寄存器中的第k个数据;所需加法器个数为2NT-1。
5.根据权利要求4所述异步成形滤波器设计方法,其特征在于:所述步骤11中成形滤波器为平方根升余弦滤波器,该平方根升余弦滤波器阶数为2RNT+1。
6.根据权利要求5所述异步成形滤波器设计方法,其特征在于:所述步骤13中对成形滤波器的成形系数的右半部分系数进行截去,进而得到截短后的成形系数F:
F = F 1 , F 2 , ... , F RN T ;
其中,F为截短后的滤波器系数向量,F1,F2,…,为相应的滤波器系数值。
7.根据权利要求6所述异步成形滤波器设计方法,其特征在于:所述步骤21中的NCO控制模块的累加精度N大于步骤11中的成形倍数R的位宽K。
8.根据权利要求7所述异步成形滤波器设计方法,其特征在于:所述步骤32中成形滤波器模块读取地址A1和A2进行延时处理,其中,A1为查找表读取地址1,A2为查找表读取地址2。
CN201510971094.1A 2015-12-22 2015-12-22 一种异步成形滤波器设计方法 Active CN105634439B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510971094.1A CN105634439B (zh) 2015-12-22 2015-12-22 一种异步成形滤波器设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510971094.1A CN105634439B (zh) 2015-12-22 2015-12-22 一种异步成形滤波器设计方法

Publications (2)

Publication Number Publication Date
CN105634439A true CN105634439A (zh) 2016-06-01
CN105634439B CN105634439B (zh) 2018-08-14

Family

ID=56049059

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510971094.1A Active CN105634439B (zh) 2015-12-22 2015-12-22 一种异步成形滤波器设计方法

Country Status (1)

Country Link
CN (1) CN105634439B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107704656A (zh) * 2017-09-06 2018-02-16 南京理工大学 基于fpga的参数化成型滤波器的设计与实现方法
CN110264412A (zh) * 2019-05-16 2019-09-20 北京奇艺世纪科技有限公司 图像处理方法、装置、终端设备以及存储介质
CN113708741A (zh) * 2021-05-28 2021-11-26 翱捷科技股份有限公司 一种系数可配置数字滤波器的设计方法及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1349335A2 (en) * 2002-03-26 2003-10-01 Kabushiki Kaisha Kenwood Selection of interpolation method for multicarrier reception
CN201663654U (zh) * 2009-11-20 2010-12-01 深圳市同洲电子股份有限公司 一种非整数倍插值装置及信号调制系统

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1349335A2 (en) * 2002-03-26 2003-10-01 Kabushiki Kaisha Kenwood Selection of interpolation method for multicarrier reception
CN201663654U (zh) * 2009-11-20 2010-12-01 深圳市同洲电子股份有限公司 一种非整数倍插值装置及信号调制系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
齐巍 等: "高动态测量系统中的自适应成形滤波器设计", 《西安交通大学学报》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107704656A (zh) * 2017-09-06 2018-02-16 南京理工大学 基于fpga的参数化成型滤波器的设计与实现方法
CN110264412A (zh) * 2019-05-16 2019-09-20 北京奇艺世纪科技有限公司 图像处理方法、装置、终端设备以及存储介质
CN113708741A (zh) * 2021-05-28 2021-11-26 翱捷科技股份有限公司 一种系数可配置数字滤波器的设计方法及装置
CN113708741B (zh) * 2021-05-28 2023-12-01 翱捷科技股份有限公司 一种系数可配置数字滤波器的设计方法及装置

Also Published As

Publication number Publication date
CN105634439B (zh) 2018-08-14

Similar Documents

Publication Publication Date Title
CN105634439B (zh) 一种异步成形滤波器设计方法
US8131790B2 (en) Decimation filter
CN107704656A (zh) 基于fpga的参数化成型滤波器的设计与实现方法
CN107196881A (zh) 一种高动态脉冲成形信号模拟方法和装置
JP2003516692A (ja) 正準符号桁数字係数表現を用いる再構成可能な有限応答フィルタ
EP0132988B1 (en) Digital synthesis technique for pulses having predetermined time and frequency domain characteristics
KR100340048B1 (ko) 승산기를 사용하지 않는 유한 임펄스 응답 필터 장치
CN102983839A (zh) 一种基于fpga实现gmsk信号发生器的方法
CN110247642B (zh) 一种fir滤波方法及滤波器
CN102983838A (zh) 一种基于fpga实现高斯滤波器数字逻辑电路的方法
CN102510325A (zh) 一种数字分路系统
CN113556101B (zh) Iir滤波器及其数据处理方法
CN114124034A (zh) 一种基于fpga的高速并行插值滤波器设计
US8242829B1 (en) Multichannel interpolator
EP1355421B1 (en) Digital matched filter
CN103457573A (zh) 一种数据率自适应的高斯滤波器
CN116111982A (zh) 一种基于fpga的iir滤波方法、存储介质及装置
CN118971901B (zh) 一种用于信号智能识别的可配置信道化实现方法及系统
CN116132232B (zh) 一种用于高速数字通信的多路并行上采样方法
JP5557339B2 (ja) デシメーションフィルタ及びデシメーション処理方法
KR100312413B1 (ko) 룩-업테이블방식을이용한디지털파형정형필터링장치
CN110635780A (zh) 一种基于fpga的变速率基带脉冲成形滤波器实现方法及滤波器
CN113872566B (zh) 带宽连续可调的调制滤波装置和方法
CN101894095B (zh) 快速哈达玛变换装置及方法
KR20030010143A (ko) 고속 클럭 활용을 위한 에프아이알 필터의 곱셈 장치 및그 방법

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
TR01 Transfer of patent right

Effective date of registration: 20240918

Address after: No. 3, Lianhe Village, Qinhuai District, Nanjing, Jiangsu 210007

Patentee after: NANJING PANDA HANDA TECHNOLOGY Co.,Ltd.

Country or region after: China

Address before: 210002 No. 301 East Zhongshan Road, Jiangsu, Nanjing

Patentee before: PANDA ELECTRONICS GROUP Co.,Ltd.

Country or region before: China

Patentee before: NANJING PANDA HANDA TECHNOLOGY Co.,Ltd.

TR01 Transfer of patent right