本发明属于外辐射源雷达的实时信号处理技术领域，涉及一种基于gpu加速的多通道fblms实现方法。
背景技术
在外辐射源雷达信号处理中，雷达回波信号中存在直达波和杂波等情况，需要对其进行直达波和杂波抑制处理。常见的杂波抑制方法有扩展关联法(extensivecancellationalgorithm,eca)、最小均方误差法(leastmeansquare,lms)、递归最小二乘法(recursiveleastsquares,rls)、归一化最小均方误差法(normalizedleastmeansquare,nlms)和频域分块最小均方误差法(frequencyblockleastmeansquare,fblms)等。eca和rls具有较好的直达波信号和多径信号抑制效果；但是eca算法涉及矩阵求逆，运算量太大，不适合实时系统；rls计算量为o(n2)，运算量也较大。nlms算法的杂波抑制效果略低于eca和rls算法，但计算量为o(n)，运算量较小；lms算法的计算量也为o(n)，但杂波抑制效果较差，不能满足要求；fblms算法由于采用了快速运算方法，大大降低了其计算复杂度，计算量较nlms和lms降低很多，改进的fblms算法在输入信号功率不断变化的情况下，可以达到更快的收敛速度，同时取得更小的稳态误差。所以，多通道下的fblms算法具有重要研究意义。
目前研究的fblms算法，主要分为基于cpu平台和gpu平台两种处理方式。基于cpu平台的实现方式，受限于cpu的实时处理能力，只能实现单通道的实时信号处理；基于gpu平台的实现方式，利用gpu的并行计算能力，若对数据进行分段而后进行并行处理，可以在一定程度上提高了运算速度，如李晓波、关欣等人提出的基于数据块交叉重组结构的直达波杂波抑制实现方法。但此方法没有针对fblms内部结构进行gpu算法优化，对数据进行分段会导致抑制比降低，同时也未提及多通道下的信号处理实现。多通道信号处理对于后期数字波束形成，目标方向估计具有重要意义。
技术实现要素：
本发明的目的，就是针对上述情况，提供一种基于gpu加速的多通道fblms实现方法，该方法针对多通道的回波信号，能以较快的速度进行数据处理，满足信号处理的实时性要求。
本发明的技术方案为：
第一步：信号分块。如图1所示，将待处理的参考信号x(n),n＝1,2,...,len和监测通道信号di(n),i＝1,2,...,m；n＝1,2,...,len进行分块。分块长度为l，分块后得参考信号x(k),k＝1,2,...,n以及监测通道信号di(k),i＝1,2,...,m；k＝1,2,...,n，其中len是数据长度，m为监测通道数，每个通道分块后共有n个数据块。
第二步：对分块后的参考信号x(k)进行扩展。如图2所示，每次数据块前面拼接上前一个数据块，使得每个数据块长度扩展为2l，需要说明的是，在第一块数据前补与块长度相同的零，标号为0；标号为0与1的两块数据分为第一组，长度为2l；接着标号为1与2的两块数据分为第二组，以此类推，得到n组长度为2l的数据。对分块数据进行重复扩展，保证在对消后拼接数据时不存在跳变点。
第三步：对扩展后的参考信号做fft变换，将时域信号换置频域。利用gpu并行处理的特性，对这n个数据块做同时做fft变换，将原本在fblms迭代中的参考信号频域转换拿到迭代过程之前实现。利用gpu并行的特性，多个数据块的fft变换可以一次完成，相比在迭代中完成n次fft所需时间要少，这是对fblms算法内部的优化。
第四步：对分块后的监测通道信号d(g,k),g＝1,2,...,m；k＝1,2,...n进行结构变形。如图2所示，将每个通道标号为g的块依次连接，记为分别记为g1,g2,...,gj,...gn分别表示监测通道i的第j块数据。以此类推，得到一组n块长度为ml的数据。
第五步：取出第三步中参考信号进行fft变换后的一组数据，记为ui(n),n＝1,2,...,2l；fblms滤波器在频率进行迭代，权向量长度为2l，每个通道都需要对应一个fblms滤波器，则m个通道的滤波器权值为ωi(n),n＝1,2,...,2lm，分为m段，每段权值的数据长度为2l。参考信号的ui与每个通道的ωi(n)分别相乘，得到yi(n)，n＝1,2,...,2lm。
第六步：按照第五步数据处理方式进行fblms算法运算，对yi进行ifft变换，得到yi(n),n＝1,2,...，2lm，yi是每一个通道经过滤波器滤波后的第i-1和第i块数据。实际使用中，我们需要对yi进行截取，仅需要每个通道数据中的后l个数据。经过截取，数据长度变为ml。
第七步：取出第四步中监测信号的第i块数据，记为di(n),n＝1,2,...,lm；yi(n)与di相减，得到ei(n),n＝1,2,...,lm，此时便是实现杂波对消后的所需要的输出信号。
第八步：对滤波器权向量进行更新迭代。记γ为遗忘因子，α为步进长度，p为ui的功率估计修正，p＝γ*p+(1-γ)*|ui|2，对ei进行fft变换，得到ei(n),n＝1,2,...,lm，对ei中的每一个通道的数据前插入长度为l的零数据块，插入后数据长度变为2ml，修正量
第九步：输出信号结构转换。各通道直达波抑制后的信号是以第四步输出的结构保存的，需要执行第四步的内过程，将其转换为各通道信号的形式。
本发明的有益效果为：
1.改进了算法的流程，对数据进行分块处理，降低了所需重复迭代次数；
2.针对fblms算法内部结构进行优化，在对参考数据进行分块和重叠后，对其进行fft变换，再进行权值迭代；在权系数迭代之前进行fft变换，减少了迭代过程中fft的运算量，提高了运算速度；
3.实现了对多通道回波信号的杂波抑制，利用gpu并行处理，满足实时性要求。在gpu显存充足的情况下，可以支持多通道的杂波抑制处理，并保持运算时间基本不变。
附图说明
图1为各通道数据分块；
图2为各通道数据块的结构变形；
图3为基于gpu的多通道fblms算法框图；
图4为输出信号结构转换示意图；
图5为使用直达波抑制之后的信号与参考信号无源相关后得到的距离多普勒图；
图6为距离维的相关峰；
图7为多普勒维的相关峰。
具体实施方式
下面通过实际数据的处理实例来验证本发明的优点。
以dtmb数字电视信号为外辐射源，dtmb信号带宽为7.56mhz，本次实验数据采样中心频率为688mhz，采样率为10mhz，接受机通道共有6个，1个作为参考通道，可作为监测通道的有5个。考虑外辐射源雷达目标信号微弱的特点，选择长累计，故实验迭代数据长度选择5mhz(即0.5s的数据)。
以下将做cpu串行处理方法，现有gpu处理方法以及本发明方法的性能比较。需要特别说明的是，由于现有gpu处理方法只适用于单通道，为更好比较单通道情况，特省略步骤第四步和第九步，简化为单通道模式以作比较。
结果如表1所示：
表1
图5是使用直达波抑制之后的信号与参考信号无源相关后得到的距离多普勒图，图中可清晰辨别出目标位置。图6和图7分别是距离维和多普勒维的相关峰，从中可清晰辨别目标信号相关峰。