CN108919227A - 一种基于gpu加速的多通道fblms实现方法 - Google Patents
一种基于gpu加速的多通道fblms实现方法 Download PDFInfo
- Publication number
- CN108919227A CN108919227A CN201810940572.6A CN201810940572A CN108919227A CN 108919227 A CN108919227 A CN 108919227A CN 201810940572 A CN201810940572 A CN 201810940572A CN 108919227 A CN108919227 A CN 108919227A
- Authority
- CN
- China
- Prior art keywords
- channel
- data
- signal
- block
- length
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明属于外辐射源雷达的实时信号处理技术领域,涉及一种基于GPU加速的多通道FBLMS实现方法。本发明针对FBLMS算法内部结构进行优化,在对参考数据进行分块和重叠后,对其进行FFT变换,再进行权值迭代;在权系数迭代之前进行FFT变换,减少了迭代过程中FFT的运算量,提高了运算速度。实现了对多通道回波信号的杂波抑制,利用GPU并行处理,满足实时性要求。在GPU显存充足的情况下,可以支持多通道的杂波抑制处理,并保持运算时间基本不变。
Description
技术领域
本发明属于外辐射源雷达的实时信号处理技术领域,涉及一种基于GPU加速的多通道FBLMS实现方法。
背景技术
在外辐射源雷达信号处理中,雷达回波信号中存在直达波和杂波等情况,需要对其进行直达波和杂波抑制处理。常见的杂波抑制方法有扩展关联法(ExtensiveCancellation Algorithm,ECA)、最小均方误差法(Least Mean Square,LMS)、递归最小二乘法(Recursive Least Squares,RLS)、归一化最小均方误差法(Normalized Least MeanSquare,NLMS)和频域分块最小均方误差法(Frequency Block Least Mean Square,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,修正量滤波器权向量的更新方程为ω(n+1)=ω(n)+△ω。转到第五步,迭代N次。
第九步:输出信号结构转换。各通道直达波抑制后的信号是以第四步输出的结构保存的,需要执行第四步的内过程,将其转换为各通道信号的形式。
本发明的有益效果为:
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分别是距离维和多普勒维的相关峰,从中可清晰辨别目标信号相关峰。
Claims (1)
1.一种基于GPU加速的多通道FBLMS实现方法,其特征在于,包括以下步骤:
S1、信号分块:
根据滤波器长度L对个通道数据分块,即将待处理的参考信号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个数据块;
S2、对分块后的参考信号x(k)进行扩展:
每个数据块前面拼接上前一个数据块,使得每个数据块长度扩展为2L,即第一块数据前补与块长度相同的零,标号为0,第二块数据前拼接上第一块数据,标号为1,以此类推,然后将标号为0与1的两块数据分为第一组,长度为2L;接着标号为1与2的两块数据分为第二组,以此类推,得到N组长度为2L的数据;
S3、对扩展后的参考信号做FFT变换,将时域信号变换为频域信号:
利用GPU并行处理的特性,对步骤S2得到的N个数据块同时做FFT变换;
S4、对分块后的监测通道信号d(g,k),g=1,2,...,M;k=1,2,...N进行结构变形:
将每个通道标号为g的块依次连接,分别记为g1,g2,...,gj,...gN,分别表示监测通道i的第j块数据;以此类推,得到一组长度为ML的数据;
S5、取出步骤S3中参考信号进行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;
S6、按照步骤S5数据处理方式进行FBLMS算法运算,对Yi进行IFFT变换,得到yi(n),n=1,2,...,2LM,yi是每一个通道经过滤波器滤波后的第i-1和第i块数据;
S7、取出步骤S4中监测信号的第i块数据,记为Di(n),n=1,2,...,LM;yi(n)与Di相减,得到ei(n),n=1,2,...,LM,此时便是实现杂波对消后的所需要的输出信号;
S8、对滤波器权向量进行更新迭代:
记γ为遗忘因子,α为步进长度,P为Ui的功率估计修正,P=γ*P+(1-γ)*|Ui|2,对ei进行FFT变换,得到Ei(n),n=1,2,...,LM,对Ei中的每一个通道的数据前插入长度为L的零数据块,插入后数据长度变为2ML,修正量滤波器权向量的更新方程为ω(n+1)=ω(n)+△ω,转到步骤S5,迭代N次;
S9、输出信号结构转换:各通道直达波抑制后的信号是以步骤S4输出的结构保存的,执行步骤S4的内过程,将其转换为各通道信号的形式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810940572.6A CN108919227B (zh) | 2018-08-17 | 2018-08-17 | 一种基于gpu加速的多通道fblms实现方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810940572.6A CN108919227B (zh) | 2018-08-17 | 2018-08-17 | 一种基于gpu加速的多通道fblms实现方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108919227A true CN108919227A (zh) | 2018-11-30 |
CN108919227B CN108919227B (zh) | 2021-12-31 |
Family
ID=64405339
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810940572.6A Active CN108919227B (zh) | 2018-08-17 | 2018-08-17 | 一种基于gpu加速的多通道fblms实现方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108919227B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110208752A (zh) * | 2019-06-27 | 2019-09-06 | 电子科技大学 | 一种基于gpu的雷达mti/mtd实现方法 |
CN110208751A (zh) * | 2019-06-21 | 2019-09-06 | 电子科技大学 | 一种基于gpu加速的雷达信号旁瓣对消实现方法 |
CN111131099A (zh) * | 2019-12-05 | 2020-05-08 | 北京航空航天大学杭州创新研究院 | 一种面向同频全双工的低复杂度fblms自干扰抵消实现方法 |
WO2021196165A1 (zh) * | 2020-04-03 | 2021-10-07 | 华为技术有限公司 | 频率分析方法、装置及雷达 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6798754B1 (en) * | 1997-11-13 | 2004-09-28 | National University Of Singapore | Acoustic echo cancellation equipped with howling suppressor and double-talk detector |
CN101509975A (zh) * | 2009-02-11 | 2009-08-19 | 北京航空航天大学 | 一种基于双子孔径不同重叠度的运动目标检测方法 |
US20130063299A1 (en) * | 2010-02-16 | 2013-03-14 | Cavitid Inc. | Systems, Methods and Apparatuses for Remote Device Detection |
JP2013113722A (ja) * | 2011-11-29 | 2013-06-10 | Mitsubishi Electric Corp | レーダ装置 |
CN103197300A (zh) * | 2013-03-26 | 2013-07-10 | 中国科学院电子学研究所 | 一种基于gpu的外辐射源雷达直达波杂波对消实时处理方法 |
CN106802409A (zh) * | 2017-01-19 | 2017-06-06 | 电子信息系统复杂电磁环境效应国家重点实验室 | 基于多cpu处理的外辐射源雷达实时信号处理方法 |
CN106936407A (zh) * | 2017-01-12 | 2017-07-07 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 频域块最小均方自适应滤波方法 |
CN107071197A (zh) * | 2017-05-16 | 2017-08-18 | 中山大学花都产业科技研究院 | 一种基于全相位多延迟分块频域的回音消除方法及系统 |
CN107171984A (zh) * | 2017-05-09 | 2017-09-15 | 东南大学 | 一种异步多载波系统频域信道估计方法 |
CN107657243A (zh) * | 2017-10-11 | 2018-02-02 | 电子科技大学 | 基于遗传算法优化的神经网络雷达一维距离像目标识别方法 |
US20180131403A1 (en) * | 2015-10-20 | 2018-05-10 | The Aerospace Corporation | Circuits and methods for reducing an interference signal that spectrally overlaps a desired signal |
-
2018
- 2018-08-17 CN CN201810940572.6A patent/CN108919227B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6798754B1 (en) * | 1997-11-13 | 2004-09-28 | National University Of Singapore | Acoustic echo cancellation equipped with howling suppressor and double-talk detector |
CN101509975A (zh) * | 2009-02-11 | 2009-08-19 | 北京航空航天大学 | 一种基于双子孔径不同重叠度的运动目标检测方法 |
US20130063299A1 (en) * | 2010-02-16 | 2013-03-14 | Cavitid Inc. | Systems, Methods and Apparatuses for Remote Device Detection |
JP2013113722A (ja) * | 2011-11-29 | 2013-06-10 | Mitsubishi Electric Corp | レーダ装置 |
CN103197300A (zh) * | 2013-03-26 | 2013-07-10 | 中国科学院电子学研究所 | 一种基于gpu的外辐射源雷达直达波杂波对消实时处理方法 |
US20180131403A1 (en) * | 2015-10-20 | 2018-05-10 | The Aerospace Corporation | Circuits and methods for reducing an interference signal that spectrally overlaps a desired signal |
CN106936407A (zh) * | 2017-01-12 | 2017-07-07 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 频域块最小均方自适应滤波方法 |
CN106802409A (zh) * | 2017-01-19 | 2017-06-06 | 电子信息系统复杂电磁环境效应国家重点实验室 | 基于多cpu处理的外辐射源雷达实时信号处理方法 |
CN107171984A (zh) * | 2017-05-09 | 2017-09-15 | 东南大学 | 一种异步多载波系统频域信道估计方法 |
CN107071197A (zh) * | 2017-05-16 | 2017-08-18 | 中山大学花都产业科技研究院 | 一种基于全相位多延迟分块频域的回音消除方法及系统 |
CN107657243A (zh) * | 2017-10-11 | 2018-02-02 | 电子科技大学 | 基于遗传算法优化的神经网络雷达一维距离像目标识别方法 |
Non-Patent Citations (4)
Title |
---|
JORGE LORENTE: "《The frequency partitioned block modified filtered-x NLMS with orthogonal correction factors for multichannel Active Noise Control》", 《DIGITAL SIGNAL PROCESSING》 * |
李晓波; 关欣; 仲利华: "《基于GPU的外辐射源雷达信号处理实时实现方法》", 《系统工程与电子技术》 * |
杨鹏程; 吕晓德; 柴致海: "《基于频域分块RDS-LMS算法的机载外辐射源雷达杂波对消 》", 《电子与信息学报》 * |
辛召强; 沈晓峰: "《利用CUDA快速实现IMM目标跟踪》", 《雷达科学与技术》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110208751A (zh) * | 2019-06-21 | 2019-09-06 | 电子科技大学 | 一种基于gpu加速的雷达信号旁瓣对消实现方法 |
CN110208752A (zh) * | 2019-06-27 | 2019-09-06 | 电子科技大学 | 一种基于gpu的雷达mti/mtd实现方法 |
CN111131099A (zh) * | 2019-12-05 | 2020-05-08 | 北京航空航天大学杭州创新研究院 | 一种面向同频全双工的低复杂度fblms自干扰抵消实现方法 |
CN111131099B (zh) * | 2019-12-05 | 2022-05-31 | 北京航空航天大学杭州创新研究院 | 一种面向同频全双工的低复杂度fblms自干扰抵消实现方法 |
WO2021196165A1 (zh) * | 2020-04-03 | 2021-10-07 | 华为技术有限公司 | 频率分析方法、装置及雷达 |
Also Published As
Publication number | Publication date |
---|---|
CN108919227B (zh) | 2021-12-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108919227A (zh) | 一种基于gpu加速的多通道fblms实现方法 | |
Zennaro | P-stability properties of Runge-Kutta methods for delay differential equations | |
CN109977914B (zh) | 基于vmd的自适应降噪方法 | |
CN104345306B (zh) | 基于Khatri‑Rao子空间的目标波达角估计方法 | |
CN103197300B (zh) | 一种基于gpu的外辐射源雷达直达波杂波对消实时处理方法 | |
CN113702909A (zh) | 一种基于声音信号到达时间差的声源定位解析解计算方法及装置 | |
CN107105111A (zh) | 一种组合步长成比例仿射投影回声消除方法 | |
CN104698446A (zh) | 基于动目标检测的雷达回波滤波方法 | |
EP1425853B1 (en) | Method and apparatus for generating a set of filter coefficients for a time updated adaptive filter | |
CN106802409A (zh) | 基于多cpu处理的外辐射源雷达实时信号处理方法 | |
Valagiannopoulos | An overview of the Watson transformation presented through a simple example | |
CN108919239B (zh) | 一种基于双组合自适应滤波的超声直达波对消方法 | |
Tian et al. | Long-time coherent integration and motion parameters estimation of radar moving target with unknown entry/departure time based on SAF-WLVT | |
CN112748898B (zh) | 一种复数向量运算装置及运算方法 | |
CN111965616A (zh) | 基于cpu和gpu的时分mimo雷达信号处理方法 | |
Jiang et al. | An efficient ADBF algorithm based on Keystone transform for wideband array system | |
CN109787676B (zh) | 一种高动态下的零陷展宽方法 | |
CN113805139A (zh) | 基于聚焦变换的宽带信号稀疏表示波达方向估计方法 | |
CN112698263A (zh) | 一种基于正交传播算子的单基地互质mimo阵列doa估计算法 | |
CN112731283A (zh) | 基于多级维纳滤波器的高亚音速飞行目标声学测向方法 | |
CN111796251A (zh) | 利用gpu对时分mimo雷达信号进行处理的方法 | |
KR20080006509A (ko) | 적응 원형 마이크로폰 어레이를 이용한 다자간 회의시스템의 개발 | |
JP2004156920A (ja) | 自己相関行列算出方法、tdl信号処理装置、アダプティブ・アレイ・アンテナ及びサイドローブ・キャンセラ | |
CN112433135B (zh) | 一种基于马顿核的整体最小二乘的局部放电信号定位方法 | |
Jianshu et al. | New recursive algorithm for matrix inversion |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |