基于FPGA的多通道信道化方法与流程

阅读：评论：0

基于fpga的多通道信道化方法
技术领域
1.本发明涉及数字信道化领域，尤其涉及一种基于fpga的多通道信道化方法。

背景技术：

2.对于高采样率、大带宽的数据进行频谱计算时，由于dft(数字傅里叶变换)的点数太大，需要采用数字信道化，将带宽分割为多个子带进行分析，可以显著的降低dft点数。
3.目前数字信道化的结构如图1所示，为每个子信道配置一个dds(直接数字式频率合成器)ip核，对于多个通道的信道化，需要将该结构多次复用，具体的次数为通道数量，每个dds ip核对应不同频率来实现多通道信道化，会带来繁杂的数据缓存、顺序调整问题,并且也会消耗很多的fpga资源。

技术实现要素：

4.本发明要解决的技术问题就在于：针对现有技术存在的技术问题，本发明提供一种基于fpga的多通道信道化方法，能够有效的减小fpga的资源消耗并减少数据缓存、顺序调整问题。
5.为解决上述技术问题，本发明提出的技术方案为：
6.一种基于fpga的多通道信道化方法，包括以下步骤：
7.获取输入数据，将所述输入数据分为a路第一数据，每一路第一数据乘以预设的jm系数，然后经过滤波器，再乘以预设的ejm系数，最后进行并行快速傅里叶变换，得到a路的第二数据；
8.对于每路第二数据，根据预设的子信道第一中心频率将该路第二数据分为a路的第三数据，并根据预设的子信道第二中心频率将该路第二数据分为a路的第四数据，选取目标数量第三数据和对应的第四数据，每一路被选取第三数据和第四数据均乘以预设的jm系数，然后经过滤波器，再乘以预设的ejm系数，最后进行并行快速傅里叶变换，得到a路的第五数据。
9.可选的，a为4的倍数，进行并行快速傅里叶变换具体包括：采用按时间抽取的基-4类fft算法，对所有输入的第一数据，或者第三数据及第四数据分组后，乘以预设的旋转因子再进行蝶形运算，得到对应的第二数据或者第五数据；或者，采用按频率抽取的基-4类fft算法，对所有输入的第一数据，或者第三数据及第四数据分组后，进行蝶形运算再乘以预设的旋转因子，得到对应的第二数据或者第五数据。
10.进一步的，对于每路第二数据，根据预设的子信道第一中心频率将该路第二数据分为a路第三数据，并根据预设的子信道第二中心频率将该路第二数据分为a路第四数据之前，还包括：若该路序数在目标范围内，对该路的第二数据求共轭得到新的第二数据。
11.进一步的，子信道第一中心频率表达式如下：
12.ω
a,p
＝(p+1/4)2π/n
13.子信道第二中心频率表达式如下：
14.ω
b,p
＝(p+3/4)2π/n
15.上式中，p＝0,1,
…
,n-1，n为并行快速傅里叶变换的输入数据路数。
16.可选的，所述目标数量的第三数据和对应的第四数据为频谱在[-π/2,π/2]区间的第三数据和对应的第四数据。
[0017]
可选的，得到a路的第五数据后还包括：若当前通道的第五数据与上一通道的第五数据的存在频谱重叠的部分，丢弃当前通道的第五数据的频谱重叠部分。
[0018]
与现有技术相比，本发明的优点在于：
[0019]
本发明将多通道信道化按照所需通道数量拆分为两级信道化，每级信道化的流程相同，且信道化个数相同，相比目前的多通道信道化，可以显著降低资源消耗并减少数据缓存、顺序调整问题。
附图说明
[0020]
图1为传统的一路数字信道化流程图。
[0021]
图2为本发明实施例的多通道信道化流程图。
[0022]
图3为基4-fft算法的分解示意图。
[0023]
图4为基4-fft算法的蝶形计算示意图。
[0024]
图5为本发明实施例的并行快速傅里叶变换示意图。
[0025]
图6为本发明实施例中根据子信道第一中心频率和子信道第二中心频率划分的16信道数据示意图。
[0026]
图7为本发明实施例中多通道信道化的第二级详细流程图。
[0027]
图8为本发明实施例中第一级信道化输出的16个子信道中第k个子信道的等效频谱覆盖范围。
具体实施方式
[0028]
以下结合说明书附图和具体优选的实施例对本发明作进一步描述，但并不因此而限制本发明的保护范围。
[0029]
鉴于adc采样频率为4800msps，直接进行m＝256信道化，将会带来繁杂的数据缓存、顺序调整问题。为了避免这些问题，本实施例将m＝256拆为两级信道化，如图1所示，先进行第一级信道化：把1路数据分16路并行输入，每路乘以一个jm系数，然后过滤波器，再乘以一个ejm系数，最后并行输入dft模块进行并行快速傅里叶变换，得到的16路输出就是16路信道化的输出数据；
[0030]
然后进行第二级信道化：这16路数据每路1路分16路，其中1～8路的数据再次分为16路，每路乘以一个jm系数，然后过滤波器，再乘以一个ejm系数，最后并行输入dft模块进行并行快速傅里叶变换，通过原理推导，9～16路的数据求共轭后，再继续前述的步骤，求共轭为本领域常规技术手段，在此不再赘述，将一共得到16*16＝256路输出，就是256路信道化的数据，每一级信道化个数均为n＝16，即m＝n
×
n。
[0031]
图1中的的“n点fft”，n＝16，输入数据速率为300msps，工作时钟为300mhz。
[0032]
对于16路并行数据，如果直接进行dft运算，每次需要的计算量为n2次复数乘法和n(n-1)次复数加法。为加快dft的运算速度，考虑n为4的整数次幂情况，即n＝4m，m为正整
数。这种情况下采用基-4类fft(快速傅里叶变换)算法。基-4类fft算法又可按实现结构分为两子类：按时间抽取(dit,decimation in time)和按频率抽取(dif,decimation in frequency)，dit和dif，前者将输入按倒位序重新排列，输出几位自然顺序排列；后者的话，输入为自然顺序，输出为倒位序，即dit先乘以旋转因子后蝶形运算，dif先蝶形运算后乘以旋转因子。本实施例中采用dif基-4fft算法。dif基-4fft算法将一个n点dft运算被分解为4个n/4点dft运算，如果n/4可以被4整除，则可以进一步将每个n/4点dft运算分解为蝶形计算结构和四个n/16点dft，直到最后4点dft运算。如图4所示，每个蝶形计算模块由4个复数乘法器、8个复数加法器(减法可以用加法器变通实现)、4个复数移位寄存器组成。
[0033]
以n＝16为例，这个分解过程如图3和图5所示，16路并行被分为4组，每组包括4路并行数据，这4路并行数据案通过蝶形计算模块进行dft运算。最后一级分解，dft运算的频移因子简化为1、-1、j、-j等值因此不需要复数乘法器。从图中可以看出，当n＝16时，整个计算结构只需要9个复数乘法器、64个复数加法器。
[0034]
在基-4类fft算法中，蝶形计算模块的输出数据不需存储，直接流向下一级，整个算法不会产生数据堆积。采用这样的并行流水结构和同步时钟驱动，可以实现每个时钟节拍输出一次fft计算结果，便于后续进行第二级信道化和后续的接收信号频谱计算。
[0035]
为了避免信道化的时候有频率覆盖不完全，存在盲区，实现对复数信号的无盲区信道化接收，本实施例中针对第二级信道化输入的16路数据，即第一级信道化输出的16路数据，将每一路数据按照a、b两种混频设置子信道中心频率：
[0036]
a情况：ω
a,p
＝(p+1/4)2π/n,p＝0,1,
…
,n-1
[0037]
b情况：ω
b,p
＝(p+3/4)2π/n,p＝0,1,
…
,n-1
[0038]
其中n为16，即得到a情况的子信道第一中心频率的16路数据和b情况的子信道第二中心频率的16路数据，如图6所示，从图中可以看出，第一级信道化每个子信道的输出，其有用频谱集中在[-π/2,π/2]区间，因此从对应的子信道第一中心频率的16路数据和子信道第二中心频率的16路数据中，各选择[-π/2,π/2]区间的8个子信道输出，作为该第一级信道化子信道输出数据划分的16路并行数据。
[0039]
如图7所示，将前面被选取的[-π/2,π/2]区间的子信道数据每路乘以一个jm系数，然后过滤波器，再乘以一个ejm系数，最后并行输入dft模块进行前述的并行快速傅里叶变换，最终得到该第一级信道化子信道的16路并行数据，这些数据存在频谱重叠的部分，因此丢弃掉一部分频谱重叠的部分。该第一级信道化子信道经过第二级信道化之后的等效频谱覆盖范围如图8所示。图7中的n点fft，输入数据速率为18.75msps，工作时钟为300mhz，因此在fpga中实现时，可以复用图5所示并行快速傅里叶方法，以尽可能节省dsp48乘法器资源。
[0040]
根据上述技术构思，本实施例提出一种基于fpga的多通道信道化方法，如图1和图7所示，包括以下步骤：
[0041]
第一级信道化：获取输入数据，将所述输入数据分为a路第一数据，每一路第一数据乘以预设的jm系数，然后经过滤波器，再乘以预设的ejm系数，最后进行并行快速傅里叶变换，得到a路的第二数据；
[0042]
第二级信道化：对于每路第二数据，根据预设的子信道第一中心频率将该路第二数据分为a路的第三数据，并根据预设的子信道第二中心频率将该路第二数据分为a路的第
四数据，选取目标数量第三数据和对应的第四数据，每一路被选取第三数据和第四数据均乘以预设的jm系数，然后经过滤波器，再乘以预设的ejm系数，最后进行并行快速傅里叶变换，得到a路的第五数据。
[0043]
通过上述步骤，即可对于输入数据进行a*a个通道的信道化，频率划分比较细，直接就能到信号的具体位置。
[0044]
本实施例中，a为4的倍数，进行并行快速傅里叶变换具体包括：采用按时间抽取的基-4类fft算法，对所有输入的第一数据，或者第三数据及第四数据分组后，乘以预设的旋转因子再进行蝶形运算，得到对应的第二数据或者第五数据；或者，采用按频率抽取的基-4类fft算法，对所有输入的第一数据，或者第三数据及第四数据分组后，进行蝶形运算再乘以预设的旋转因子，得到对应的第二数据或者第五数据。
[0045]
本实施例中，对于每路第二数据，根据预设的子信道第一中心频率将该路第二数据分为a路第三数据，并根据预设的子信道第二中心频率将该路第二数据分为a路第四数据之前，还包括：若该路序数在目标范围内，对该路的第二数据求共轭得到新的第二数据，a＝16时，目标范围的序数为9～16。
[0046]
如前所述，子信道第一中心频率表达式如下：
[0047]
ω
a,p
＝(p+1/4)2π/n
[0048]
子信道第二中心频率表达式如下：
[0049]
ω
b,p
＝(p+3/4)2π/n
[0050]
上式中，p＝0,1,
…
,n-1，n为并行快速傅里叶变换的输入数据路数。
[0051]
本实施例中，所述目标数量的第三数据和对应的第四数据为频谱在[-π/2,π/2]区间的第三数据和对应的第四数据，a＝16时，目标数量的第三数据和对应的第四数据为频谱在[-π/2,π/2]区间的8个第三数据和8个对应的第四数据。
[0052]
本实施例中，得到a路的第五数据后还包括：若当前通道的第五数据与上一通道的第五数据的存在频谱重叠的部分，丢弃当前通道的第五数据的频谱重叠部分。
[0053]
综上所述，本实施例的方法将多通道的信道化拆分为两级信道化，减少了每一级信道化个数，且每一级信道化采用并行傅里叶变换的方法，具体采用按时间抽取的基-4类fft算法或者按频率抽取的基-4类fft算法，可以显著降低资源消耗并减少数据缓存、顺序调整问题。
[0054]
上述只是本发明的较佳实施例，并非对本发明作任何形式上的限制。虽然本发明已以较佳实施例揭露如上，然而并非用以限定本发明。因此，凡是未脱离本发明技术方案的内容，依据本发明技术实质对以上实施例所做的任何简单修改、等同变化及修饰，均应落在本发明技术方案保护的范围内。

技术特征：

1.一种基于fpga的多通道信道化方法，其特征在于，包括以下步骤：获取输入数据，将所述输入数据分为a路第一数据，每一路第一数据乘以预设的jm系数，然后经过滤波器，再乘以预设的ejm系数，最后进行并行快速傅里叶变换，得到a路的第二数据；对于每路第二数据，根据预设的子信道第一中心频率将该路第二数据分为a路的第三数据，并根据预设的子信道第二中心频率将该路第二数据分为a路的第四数据，选取目标数量第三数据和对应的第四数据，每一路被选取第三数据和第四数据均乘以预设的jm系数，然后经过滤波器，再乘以预设的ejm系数，最后进行并行快速傅里叶变换，得到a路的第五数据。2.根据权利要求1所述的基于fpga的多通道信道化方法，其特征在于，a为4的倍数，进行并行快速傅里叶变换具体包括：采用按时间抽取的基-4类fft算法，对所有输入的第一数据，或者第三数据及第四数据分组后，乘以预设的旋转因子再进行蝶形运算，得到对应的第二数据或者第五数据；或者，采用按频率抽取的基-4类fft算法，对所有输入的第一数据，或者第三数据及第四数据分组后，进行蝶形运算再乘以预设的旋转因子，得到对应的第二数据或者第五数据。3.根据权利要求1所述的基于fpga的多通道信道化方法，其特征在于，对于每路第二数据，根据预设的子信道第一中心频率将该路第二数据分为a路第三数据，并根据预设的子信道第二中心频率将该路第二数据分为a路第四数据之前，还包括：若该路序数在目标范围内，对该路的第二数据求共轭得到新的第二数据。4.根据权利要求1所述的基于fpga的多通道信道化方法，其特征在于，子信道第一中心频率表达式如下：ω
a,p
＝(p+1/4)2π/n子信道第二中心频率表达式如下：ω
b,p
＝(p+3/4)2π/n上式中，p＝0,1,
…
,n-1，n为并行快速傅里叶变换的输入数据路数。5.根据权利要求1所述的基于fpga的多通道信道化方法，其特征在于，所述目标数量的第三数据和对应的第四数据为频谱在[-π/2,π/2]区间的第三数据和对应的第四数据。6.根据权利要求1所述的基于fpga的多通道信道化方法，其特征在于，得到a路的第五数据后还包括：若当前通道的第五数据与上一通道的第五数据的存在频谱重叠的部分，丢弃当前通道的第五数据的频谱重叠部分。

技术总结

本发明公开了一种基于FPGA的多通道信道化方法，包括以下步骤：获取输入数据，将所述输入数据分为a路第一数据，每一路第一数据乘以预设的jm系数，然后经过滤波器，再乘以预设的ejm系数，最后进行并行快速傅里叶变换，得到a路的第二数据；对于每路第二数据，根据预设的子信道第一中心频率将该路第二数据分为a路的第三数据，并根据预设的子信道第二中心频率将该路第二数据分为a路的第四数据，选取目标数量第三数据和对应的第四数据，每一路被选取第三数据和第四数据均乘以预设的jm系数，然后经过滤波器，再乘以预设的ejm系数，最后进行并行快速傅里叶变换，得到a路的第五数据。本发明能够有效的减小FPGA的资源消耗并减少数据缓存、顺序调整问题。顺序调整问题。顺序调整问题。