PCM数据滚降滤波处理方法技术

本发明专利技术所述的PCM数据滚降滤波处理方法，根据冲击响应函数的时间取值采取时域叠加计算以进行滚降滤波处理，从而实现降低对FPGA数据处理平台通用化设计的技术要求、通过少量的加法运算降低占用计算资源和缩短数据处理时间的设计目的。包括以下实施步骤：步骤1、冲击响应函数h(t)时域离散化；步骤2、基带数据与自然冲击函数序列转化；步骤3、时域波形拟合；步骤4、时域叠加计算，根据当前时刻之前或之后的若干个数据符号的等效时刻，对滚降滤波冲击响应函数进行取样，并将所有数据符号的影响进行叠加计算。叠加计算。叠加计算。

【技术实现步骤摘要】
PCM数据滚降滤波处理方法


[0001]本专利技术涉及一种基于FPGA硬件架构实现的PCM数据滚降滤波处理方法，属于无线通信领域。

技术介绍

[0002]在无线通信系统调制发送基带数据信号之前，为有效地利用频谱和防止码间干扰问题，普遍地对基带数据进行滚降滤波处理。
[0003]特别对于采用FPGA作为数据处理硬件平台，基于该硬件架构的简便性特点，现有滚降滤波采取整数倍插值的滤波方式，如在先公开的专利CN201610596082.X，其处理过程包括：(1)采用基带数据速率整数倍N的高速时钟对基带数据进行高倍抽样，对抽样后的数据进行补零，即第1点为原始基带数据，第2至N点数据填充零；(2)对高倍抽样后的数据通过数字滤波器进行滚降滤波处理。
[0004]如果高倍抽取时钟与数据处理平台主工作时钟不满足倍数关系，进一步地对步骤(2)滚降滤波后的数据进行插值滤波处理，如另一在先公开专利CN201410469317.X，其处理过程是将滚降滤波后数据进行Farrow滤波插值计算，将高倍抽样时钟域数据转换到数据处理平台工作时钟域。
[0005]如上述方案的现有技术均基于滚降滤波器采用整数倍插值的滤波方法，存在的缺陷与不足之处在于：
[0006]1、基带数据速率与数据处理时钟速率存在整数倍关系，在实践中基带数据速率随应用场景不同变化范围很大，导致整数倍的滚降滤波器数据处理时钟频率也随之变化，难以实现平台硬件的通用化。以常用的16倍高倍抽取为例，当基带数据速率达10Mbps时，数据处理时钟需要达到160MHz；当基带数据速率达3.069Mbps时，数据处理时钟需要达到49.104MHz。不同的数据处理时钟频率需求对FPGA数据处理平台的通用化设计形成了较大困扰。
[0007]2、滚降滤波器需要进行大量乘法和加法运算，占用较多的计算资源。滚降滤波过程通过时域卷积计算过程完成，对于M阶滚降滤波器，需要M次乘法运算和M﹡(M
‑
1)/2次加法运算。
[0008]3、如果高倍抽取时钟与数据处理平台主工作时钟不满足倍数关系，采用Farrow插值滤波处理，需要实时更新时间间隔信息并计算插值系数，实现过程繁杂、占用时间较长；以4点Farrow插值计算为例，如下式所示的插值系数计算，式中u为时间间隔，f(n)为插值系数。
[0009]f(
‑
2)＝
‑
1/6*u3+1/2*u2‑
1/3*u
[0010]f(
‑
1)＝1/2*u3‑
u2‑
1/2*u+1
[0011]f(0)＝
‑
1/2*u3+1/2*u2+u
[0012]f(1)＝1/6*u3‑
1/6*u
[0013]如上式所示，Farrow插值滤波涉及到立方、平方、数值乘/加等多类计算，计算过程
中还必须对时间间隔u进行实时计时更新。
[0014]有鉴于此，特提出本专利申请。

技术实现思路

[0015]本专利技术所述的PCM数据滚降滤波处理方法，在于解决上述现有技术存在的问题而根据冲击响应函数的时间取值采取时域叠加计算以进行滚降滤波处理，从而实现降低对FPGA数据处理平台通用化设计的技术要求、通过少量的加法运算降低占用计算资源和缩短数据处理时间的设计目的。
[0016]为实现上述设计目的，所述的PCM数据滚降滤波处理方法，是根据冲击响应函数h(t)的时间取值进行时域叠加计算后得到。
[0017]具体地，包括以下实施步骤：
[0018]步骤1、冲击响应函数h(t)时域离散化
[0019]对冲击响应函数的有效时间周期按M段(M＞＝2,M为整数)进行离散化；
[0020]步骤2、基带数据与自然冲击函数序列转化
[0021]基带数据为基带脉冲数据序列，将基带数据进行差分计算并转化为自然冲击函数序列；
[0022]步骤3、时域波形拟合
[0023]将当前采样时刻之前数据符号对应的冲击响应函数波形、当前采样时刻之后的时间符号对应的冲击函数响应波形在当前时刻进行拟合；包括，
[0024](a)当前时刻按数据符号周期定位
[0025](b)定位后的时刻进行周期延拓
[0026](c)周期延拓后的时刻离散化
[0027]在当前时刻之前或之后数据符号的等效时刻，根据上述步骤1确定的时间量化精度进行离散化处理；
[0028]离散后的等效时刻t
ni
＝
‑
M/2，
‑
M/2+1,...,0,1,...,M/2；
[0029]在等效时刻的拟合数据值为h
d
(t
ni
)；
[0030]步骤4、时域叠加计算
[0031]根据当前时刻之前或之后的若干个数据符号的等效时刻，对滚降滤波冲击响应函数进行取样，并将所有数据符号的影响进行叠加计算。
[0032]进一步地，在步骤2中，转化后的自然冲击函数序列表示为，
[0033][0034]其中，i＝0，1，...，分别对应i*t
s
时刻的采样数据，t
s
为抽样周期。
[0035]进一步地，在步骤4中，时域叠加计算的处理结果为，
[0036][0037]其中，为离散后的等效时刻，N0为周期延拓的符号位数，T
b
为符号周期，t
i
'为当前时刻数据符号周期定位后的时间，d
n
为步骤2中转化的数据符号序列。
[0038]综上内容，本申请所述的PCM数据滚降滤波处理方法具有如下优点：
[0039]1、与现有技术相比，本申请对FPGA数据处理平台的通用化设计要求较低，究其原因是进行滚降滤波处理时无需要求基带数据速率与数据处理时钟速率存在整数倍关系，申请方案中对原始输入的PCM数据直接进行时域抽样，抽样速率不小于数据码率的2倍即可，从而能够适用于所有FPGA数据处理平台。
[0040]2、本申请采取采样时域拟合方式，无需进行时域卷积运算，即滚降滤波过程采用前/后数据符号冲击响应函数波形的叠加。由于省略了乘法计算、且加法计算次数仅与前后关联数据符号的数线性相关，因此既明显地减少了大量计算资源的占用，同时简化了计算过程、有效地缩短处理周期。
附图说明
[0041]以下附图是本申请具体实施方式的举例说明。
[0042]图1是滚降滤波器冲击响应函数的时域连续波形示意图；
[0043]图2是滚降滤波器冲击响应函数的时域离散波形示意图；
[0044]图3是待处理基带PCM数据图；
[0045]图4是恢复出的冲击函数序列波形图；
[0046]图5是时域叠加计算示意图。
具体实施方式
[0047]实施例1，所述PCM数据滚降滤波处理方法采取的设计思路如下：
[0048]基带数据脉冲本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种PCM数据滚降滤波处理方法，其特征在于：根据冲击响应函数h(t)的时间取值进行时域叠加计算后得到；包括以下实施步骤，步骤1、冲击响应函数h(t)时域离散化对冲击响应函数的有效时间周期按M段(M＞＝2,M为整数)进行离散化；步骤2、基带数据与自然冲击函数序列转化基带数据为基带脉冲数据序列，将基带数据进行差分计算并转化为自然冲击函数序列；步骤3、时域波形拟合将当前采样时刻之前数据符号对应的冲击响应函数波形、当前采样时刻之后的时间符号对应的冲击函数响应波形在当前时刻进行拟合；包括，(a)当前时刻按数据符号周期定位(b)定位后的时刻进行周期延拓(c)周期延拓后的时刻离散化在当前时刻之前或之后数据符号的等效时刻，根据上述步骤1确定的时间量化精度进行离散化处理；离散后的等效时刻t
ni
＝
‑
M/2，
‑
M/2+1,...,0,1,...,...

【专利技术属性】
技术研发人员：童亚钦，邵飞，康峰源，张伟，李汶芯，
申请(专利权)人：山东航天电子技术研究所，
类型：发明
国别省市：