CN113203997B - 基于fpga的雷达超分辨测向方法、系统及应用 - Google Patents
基于fpga的雷达超分辨测向方法、系统及应用 Download PDFInfo
- Publication number
- CN113203997B CN113203997B CN202110406848.4A CN202110406848A CN113203997B CN 113203997 B CN113203997 B CN 113203997B CN 202110406848 A CN202110406848 A CN 202110406848A CN 113203997 B CN113203997 B CN 113203997B
- Authority
- CN
- China
- Prior art keywords
- matrix
- decomposition
- data
- flow pattern
- fpga
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 72
- 239000011159 matrix material Substances 0.000 claims abstract description 117
- 238000001228 spectrum Methods 0.000 claims abstract description 90
- 238000004364 calculation method Methods 0.000 claims abstract description 89
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 88
- 238000007667 floating Methods 0.000 claims abstract description 70
- 239000013598 vector Substances 0.000 claims abstract description 46
- 238000012545 processing Methods 0.000 claims abstract description 40
- 238000005259 measurement Methods 0.000 claims abstract description 26
- 230000003595 spectral effect Effects 0.000 claims description 52
- 239000013307 optical fiber Substances 0.000 claims description 7
- 238000011161 development Methods 0.000 claims description 2
- 230000007547 defect Effects 0.000 abstract description 8
- 238000010586 diagram Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000005855 radiation Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 238000002054 transplantation Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000010349 pulsation Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
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
-
- 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
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/14—Systems for determining direction or deviation from predetermined direction
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/04—Programme control other than numerical control, i.e. in sequence controllers or logic controllers
- G05B19/042—Programme control other than numerical control, i.e. in sequence controllers or logic controllers using digital processors
- G05B19/0423—Input/output
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Automation & Control Theory (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提出了一种基于FPGA的雷达超分辨测向方法、系统及应用,方法包括:获取雷达数据,写入硬浮点FPGA;计算雷达数据的协方差矩阵;将协方差矩阵转换为浮点模式,写入特征分解模块,进行基于QR分解的特征分解计算和基于施密特正交的QR分解计算,获取多个特征值和多个特征向量;根据预设的信源数门限,参考特征值进行信源数判定,结合特征向量构造噪声矩阵;根据预设的流型阵和噪声矩阵计算空间谱,基于预设的向量点乘器进行空间谱的计算和搜索,获取谱峰数据,进而获取雷达数据对应的空间坐标信息。本发明克服了现有技术中的信号处理平台处理能力不足及测量时间过长的缺陷,大大提升了DOA的计算效率,缩短了计算时间。
Description
技术领域
本发明涉及雷达实时DOA领域,特别涉及一种基于FPGA的雷达超分辨测向方法、系统及应用。
背景技术
MUSIC算法是一种基于矩阵特征空间分解的方法,被广泛应用在DOA(DOA,是电子、通信、雷达、声呐等研究领域的行业内用语,通过处理接收到的回波信号,获取目标的距离信息和方位信息)相关领域中。由于MUSIC算法包含非常大的计算量,现有系统的处理速度难以满足实际应用的需要,因此MUSIC算法只能普遍应用于一些实时性要求不太严格的领域,如对雷雨天气的研究、对慢速船舶的定位等,但在移动通信、电子侦察以及电子对抗等对实时性要求严格领域中受到极大限制。
由于整个MUSIC算法计算复杂度和灵活度都很大,信号处理平台(TIDSP和常规FPGA)在低功耗、轻重量限制及小型化要求下无法满足时间需求和精度需求,同时大多数的EW系统都要求浮点处理,所以国内大多采用DSP处理器运载MUSIC算法,处理速度慢,处理时间停留在ms量级,远远达不到实际应用的需要。
矩阵运算被广泛的运用于科学计算、雷达数字信号处理和图像处理等领域,随着矩阵运算量和复杂程度的不断增加,传统硬件平台逐渐无法满足系统实时计算的需求。FPGA具有运算速度快、灵活性高且可以并行计算的特点,这使其在矩阵运算中具有明显的优势。设计基于FPGA的矩阵求逆IP核及其实验平台,对提高矩阵求逆的效率和运算速度,具有重要的工程应用价值。然而,随着矩阵阶数增大,定点数据在计算过程中的截位误差累积,造成最终的数据精度大幅下降,一般的FPGA难以满足应用需求。
由此,需要有一种更好的方案来解决现有技术中信号处理平台(TIDSP和常规FPGA)在低功耗、轻重量限制及小型化要求下处理能力不足及测量时间过长的缺陷,从而提升DOA的计算效率,缩短计算时间。
发明内容
有鉴于此,本发明提出了一种基于硬浮点FPGA的雷达超分辨测向方法、系统及数据处理设备,具体方案如下:
一种基于硬浮点FPGA的雷达超分辨测向方法,包括:
获取雷达数据,通过多通道形式将所述雷达数据写入硬浮点FPGA,所述雷达数据为通道数×快拍数的定点矩阵;
通过所述硬浮点FPGA计算所述雷达数据的协方差矩阵;
将所述协方差矩阵由定点模式转换为浮点模式,写入预设的特征分解模块,所述特征分解模块对所述协方差矩阵进行基于QR分解的特征分解计算和基于施密特正交的QR分解计算,获取多个特征值和多个特征向量;
根据预设的信源数门限,参考所述特征值进行信源数判定,结合所述特征向量构造噪声矩阵;
获取预设的流型阵,根据所述流型阵和所述噪声矩阵计算空间谱,基于预设的向量点乘器进行所述空间谱的计算和搜索,获取谱峰数据;
根据所述谱峰数据获取所述雷达数据对应的空间坐标信息。
在一个具体实施例中,所述谱峰数据的获取过程包括,
初始谱峰获取:获取预设的流型阵,根据所述流型阵和所述噪声矩阵计算空间谱,基于预设的向量点乘器进行所述空间谱的计算和搜索,获取谱峰;
谱峰迭代:将所述谱峰作为粗测结果,基于所述粗测结果对所述流型阵进行汇聚和插值处理,获取更新后的流型阵,根据所述噪声矩阵和更新后的流型阵重新计算空间谱,对更新后的空间谱进行计算和搜索,获得更新后的谱峰,重复迭代,直至满足迭代停止条件;
谱峰数据获取:将满足迭代停止条件的谱峰作为所述谱峰数据。
在一个具体实施例中,所述空间谱的表达式为:
其中,Pmu表示空间谱,a表示流型阵,En表示噪声矩阵,h表示矩阵共轭转置。
在一个具体实施例中,所述硬浮点FPGA包括Intel Arria10系列,开发平台包括QUARTUS。
在一个具体实施例中,所述特征分解模块由多多个QR分解迭代构成,迭代的次数包括所述通道数的5倍;
特征分解计算时间的表达式为:
TQR=L×C×2k。
其中,TQR表示特征分解计算时间,L表示通道数,C表示迭代次数,k表示单次QR分解的时间。
在一个具体实施例中,所述协方差矩阵的表达式为:
A=XXT
其中,A表示协方差矩阵,X表示雷达数据;
对所述协方差矩阵进行QR分解,可得:
A=QR
其中,A表示协方差矩阵,Q表示正交矩阵,R表示上三角矩阵,RQ与A正交相似,具有相同的特征值。
在一个具体实施例中,所述流型阵包括128×128的矩阵;
每个谱峰获取的时间为128×128个周期。
一种基于硬浮点FPGA的雷达超分辨测向系统,包括如下,
数据获取单元:用于获取雷达数据,通过多通道形式将所述雷达数据写入硬浮点FPGA,所述雷达数据为通道数×快拍数的定点矩阵;
协方差计算单元:用于通过所述硬浮点FPGA中计算所述雷达数据的协方差矩阵;
特征分解单元:用于将所述协方差矩阵由定点模式转换为浮点模式,写入预设的特征分解模块,所述特征分解模块对所述协方差矩阵进行基于QR分解的特征分解计算和基于施密特正交的QR分解计算,获取多个特征值和多个特征向量;
噪声计算单元:用于根据预设的信源数门限,参考所述特征值进行信源数判定,结合所述特征向量构造噪声矩阵;
谱峰搜索计算单元:用于获取预设的流型阵,根据所述流型阵和所述噪声矩阵计算空间谱,基于预设的向量点乘器进行所述空间谱的计算和搜索,获取谱峰数据;
坐标获取单元:用于根据所述谱峰数据获取所述雷达数据对应的空间坐标信息。
在一个具体实施例中,所述谱峰搜索计算单元包括,
初始谱峰获取单元:用于获取预设的流型阵,根据所述流型阵和所述噪声矩阵计算空间谱,基于预设的向量点乘器进行所述空间谱的计算和搜索,获取谱峰;
谱峰迭代单元:用于将所述谱峰作为粗测结果,基于所述粗测结果对所述流型阵进行汇聚和插值处理,获取更新后的流型阵,根据所述噪声矩阵和更新后的流型阵重新计算空间谱,对更新后的空间谱进行计算和搜索,获得更新后的谱峰,重复迭代,直至满足迭代停止条件;
谱峰数据获取单元:用于将满足迭代停止条件的谱峰作为所述谱峰数据。
一种雷达数据处理设备,包括计算机、光纤和数据计算板卡,所述光纤分别连接所述计算机和所述数据计算板卡;
所述数据计算板卡上设置有上述所述的基于硬浮点FPGA的雷达超分辨测向系统。
有益效果:本发明提出了一种基于FPGA的雷达超分辨测向方法、系统及应用,克服了现有技术中的信号处理平台(DSP和常规FPGA)在低功耗、轻重量限制及小型化要求下处理能力不足及测量时间过长的缺陷。采用MUSIC算法,基于Intel硬浮点FPGA进行算法的工程化实现。其中以并行的QR迭代方法实现厄米特矩阵特征值分解,以并行的汇聚插值算法实现谱估计,同时基于特殊的硬浮点FPGA结构大大提升了DOA的计算效率,缩短了计算时间。
附图说明
图1为本发明实施例1基于硬浮点FPGA的雷达超分辨测向方法流程图;
图2为本发明实施例1硬浮点FPGA的结构图;
图3为本发明实施例1向量点乘模块图;
图4为本发明实施例1硬浮点IP的模块图;
图5为本发明实施例1谱峰示意图;
图6为本发明实施例2基于硬浮点FPGA的雷达超分辨测向系统结构图;
图7为本发明实施例3雷达超分辨测向数据处理设备结构图。
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本发明的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。
附图标记:1-数据获取单元;2-协方差计算单元;3-特征分解单元;4-噪声计算单元;5-谱峰搜索计算单元;6-坐标获取单元;12-计算机;14-光纤;16-数据计算板卡;51-初始谱峰获取单元;52-谱峰迭代单元;53-谱峰数据获取单元。
具体实施方式
在下文中,将更全面地描述本公开的各种实施例。本公开可具有各种实施例,并且可在其中做出调整和改变。然而,应理解:不存在将本公开的各种实施例限于在此公开的特定实施例的意图,而是应将本公开理解为涵盖落入本公开的各种实施例的精神和范围内的所有调整、等同物和/或可选方案。
本发明提供了一种基于硬浮点FPGA的雷达超分辨测向方法、系统及数据处理设备,克服了现有技术中的信号处理平台(DSP和常规FPGA)在低功耗、轻重量限制及小型化要求下处理能力不足及测量时间过长的缺陷。采用MUSIC算法,基于Intel硬浮点FPGA进行算法的工程化实现。其中以并行的QR迭代方法实现厄米特矩阵特征值分解,以并行的汇聚插值算法实现谱估计,同时基于特殊的硬浮点FPGA结构大大提升了DOA的计算效率,缩短了计算时间。
在雷达系统领域,尤其是在STAP、DBF、CFAR等雷达关键技术方面,目前使用FPGA方案的比较少,仅有少数研究所在从事这方面的研究。这主要是因为:1.国内FPGA在雷达方面的应用案例还比较少,缺乏借鉴经验;2.需要对FPGA的架构、工具以及算法有深入的了解,才能有效利用FPGA;3.算法由DSP、GPU、x86平台移植到FPGA平台存在较大的技术难度,平台之间的技术移植受限于平台的规则、技术的复杂度以及技术的适用性等方面,难以进行移植。不同厂家的FPGA芯片内部构造不同,尤其是高精度、高速率的FPGA,内部构造极为复杂,需要对该芯片程度极深的认识和理解。
本发明在某些特定领域,尤其是实时性要求严格的领域,具有很高的应用价值和推广价值。本发明可应用于弹载雷达导引头实时处理系统,使系统具备对高超音速辐射源进行连续、稳定的超分辨率测量的能力。
在本公开的各种实施例中使用的术语仅用于描述特定实施例的目的并且并非意在限制本公开的各种实施例。如在此所使用,单数形式意在也包括复数形式,除非上下文清楚地另有指示。除非另有限定,否则在这里使用的所有术语(包括技术术语和科学术语)具有与本公开的各种实施例所属领域普通技术人员通常理解的含义相同的含义。所述术语(诸如在一般使用的词典中限定的术语)将被解释为具有与在相关技术领域中的语境含义相同的含义并且将不被解释为具有理想化的含义或过于正式的含义,除非在本公开的各种实施例中被清楚地限定。
实施例1
本发明实施例1公开了一种基于硬浮点FPGA的雷达超分辨测向方法,具体步骤如说明书附图1所示,具体方案如下:
101、数据获取:获取雷达数据,通过多通道形式将雷达数据写入硬浮点FPGA,雷达数据为通道数×快拍数的定点矩阵;
102、协方差计算:通过硬浮点FPGA中计算雷达数据的协方差矩阵;
103、特征分解:将协方差矩阵由定点模式转换为浮点模式,写入预设的特征分解模块,特征分解模块对协方差矩阵进行基于QR分解的特征分解计算和基于施密特正交的QR分解计算,获取多个特征值和多个特征向量;
104、噪声构造:根据预设的信源数门限,参考特征值进行信源数判定,结合特征向量构造噪声矩阵;
105、谱峰数据获取:获取预设的流型阵,根据流型阵和噪声矩阵计算空间谱,基于预设的向量点乘器进行空间谱的计算和搜索,获取谱峰数据;
106、坐标获取:根据谱峰数据获取雷达数据对应的空间坐标信息。
随着矩阵阶数增大,定点数据在计算过程中的截位误差累积,造成最终结果精度大幅下降,难以满足应用需求,因此必须采用浮点处理。但常规FPGA进行浮点计算时,采用的软浮点(用逻辑构造乘加运算)结构,在运行过程中会降低矩阵求逆的效率。
针对这一弊端,本发明采用业界第一款硬浮点架构的FPGA(Intel Arria10)实现浮点复矩阵求逆运算,极大提升了求逆算法的效率
具体地,本实施例为了提升运算效率,设计选用搭载硬浮点DSP计算单元的IntelArria10系列FPGA实现,芯片内部结构如说明书附图2所示。该DSP单元采用芯片硬件定制浮点设计,遵循IEEE754标准,单个DSP即可实现一个浮点加法和一个浮点乘法,无需额外的逻辑资源。此外通过调整DSP内部流水寄存器级数,可以满足不同的latency和Fmax的需求,无需在timing上做过多努力,就能工作在极高的工作频率。在浮点计算密集的情况下,Arria10比起业界同等级的FPGA在资源、效率、运算速度方面均有明显优势。相对于传统FPGA,本实施例采用的雷达超分辨测向方法在求逆算法效率上能够提升2-3倍。
具体地,101、数据获取:获取雷达数据,通过多通道形式将雷达数据写入硬浮点FPGA,雷达数据为通道数×快拍数的定点矩阵。雷达数据通过高速串行总线发送,在写入硬浮点FPGA之前,还需要对雷达数据进行数据下变频处理。数字下变频,是指在超外差式接收机中经过混频后得到的中频信号比原始信号的频率低的一种混频方式。
具体地,102、协方差计算:通过硬浮点FPGA中计算雷达数据的协方差矩阵。计算单元dot(向量点乘器)的并行宽度通常为快拍数的8分之1。如果通道数为8,协方差矩阵计算时间为8×8×(512/8)=1024时钟周期。为了节省BRAM,需考虑:
(1)协方差矩阵计算为定点模式,计算完成后转为单精度浮点模式;
(2)计算所需的共轭对称矩阵用原始矩阵变换,无需存储。
具体的协方差矩阵计算包括:输入雷达数据X,快拍数为m,通道数为n,协方差矩阵地表达式为:
A=XXT
其中,雷达回波数据X为n×m的定点矩阵。
具体地,103、特征分解:将协方差矩阵由定点模式转换为浮点模式,写入预设的特征分解模块,特征分解模块对协方差矩阵进行基于QR分解的特征分解计算和基于施密特正交的QR分解计算,获取多个特征值和多个特征向量。
将协方差矩阵转为由定点模式转为浮点模式,写入特征分解模块的BRAM缓存。其中,特征分解模块由QRD(QR分解)迭代构成,迭代的次数可通过上位机设置,迭代次数通常为通道数的1.5倍。特征分解模块的启动信号由协方差矩阵计算完成信号经延时形成,即协方差矩阵计算完成后,延时预设时间即可启动特征分解模块。
特征分解包括基于QR分解的特征分解计算、基于施密特正交的QR分解计算和关键算子dot(向量点乘)的硬浮点实现。
基于QR分解的特征分解计算包括:对于非奇异方阵A,进行QR分解,可得A=QR,其中R为上三角矩阵,Q为正交矩阵,即QTQ=I,I为单位矩阵。
则可推导出:QTAQ=QTQRQ=RQ
也就是说,RQ与A是正交相似的,它们具有相同的特征值。因此可以采用如下的QR迭代格式:
Ak=QkRk,k∈N
Ak+1=Qk+1Rk+1,k∈N
取A1=A,经过第二个迭代式进行迭代,可以证明Ak收敛于矩阵A的特征值,即Uk=Q1Q2…Qk。
基于施密特正交的QR分解计算,本实施例通过Matlab实现,具体代码如下:
关键算子dot(向量点乘)的硬浮点实现涉及到N维方阵相乘,需要计算n×n次n维向量点乘。向量点乘的运算效率对整个方法的实现效率具有非常关键的作用,本实施例采用脉动架构实现向量点乘,可大大降低FPGA的资源占用,提升FPGA的运算效率。具体结构如说明书附图3所示。高效率的浮点算子(如浮点加,减,乘,向量点乘),都可以采用QUARTUS的硬浮点IP生成,如说明书附图4所示。
特征分解计算时间的表达式为:
TQR=L×C×2k。
其中,TQR表示特征分解计算时间,L表示通道数,C表示迭代次数,k表示单次QR分解的时间。如果通道数为8,特征分解计算时间约为8×1.5×2k(单次QRD时间)=24K时钟周期,完成特征分解计算后,输出8个特征值(D1,....D8)和8个特征向量(v1,....v8)。
具体地,104、噪声构造:根据预设的信源数门限,参考特征值进行信源数判定,结合特征向量构造噪声矩阵。105、谱峰数据获取:获取预设的流型阵,根据流型阵和噪声矩阵计算空间谱,基于预设的向量点乘器进行空间谱的计算和搜索,获取谱峰数据。
根据上位机设定的信源数门限,参考103获取的8个特征值结果,进行信源数判定,再选出相应的特征向量,形成噪声矩阵。然后读入DDR(双倍速率同步动态随机存储器)中的流型阵,计算空间谱。
整个空间谱的计算和搜索过程基于dot(向量点乘器)全部展开,每个谱峰点的计算/搜索时间为1个周期,整个谱峰生成,搜索时间为128×128=16384周期。
噪声空间计算具体包括:
将矩阵的特征值λ进行从小到大的排序:λ1≥λ2≥...≥λM>0,其中D个较大的特征值对应于信号,M-D个较小的特征值对应于噪声。用噪声向量构造一个噪声矩阵:
En=[vD+1,vD+2,...,vM]
空间谱的表达式为:
其中,Pmu表示空间谱,a表示流型阵,En表示噪声矩阵,h表示矩阵共轭转置。
空间谱生成也基于大量的向量点乘(dot),采用并行脉动架构可大大降低FPGA的资源占用。
谱峰搜索计算具体包括:采用并行九宫搜峰算法进行谱峰搜索。九宫搜峰算法如表1所示
表1九宫搜峰算法
P1 | P2 | P3 |
P4 | P5 | P6 |
P7 | P8 | P9 |
P5>(p1,p2,p3,p4,p6,p7,p8,p9)and p5>搜峰门限
其中,P5为所需谱峰,其坐标对应为辐射源的俯仰和方向角,谱峰搜索模块有大量并行浮点比较器构成。
特别地,在本实施例中,谱峰数据的获取需要经过迭代。一次完整的空间谱搜索和计算往往不能满足测量精度的需求。本实施例通过前一次搜索的峰值更新流型阵,进而实现谱峰的更新,大大提升了测量精度。具体流程如下:
初始谱峰获取:获取预设的流型阵,根据流型阵和噪声矩阵计算空间谱,基于预设的向量点乘器进行空间谱的计算和搜索,获取谱峰。此时的谱峰为第一次搜峰完成后的结果。
谱峰迭代:将谱峰作为粗测结果,基于粗测结果对流型阵进行汇聚和插值处理,获取更新后的流型阵,根据噪声矩阵和更新后的流型阵重新计算空间谱,对更新后的空间谱进行计算和搜索,获得更新后的谱峰,重复迭代,直至满足迭代停止条件。
谱峰数据获取:将满足迭代停止条件的谱峰作为谱峰数据。
例如,为进一步提升测量精度,在上一次搜峰(粗测)完成后,可以改变流型阵(基于粗测结果,对流型阵进行汇聚、插值),进行第二次测量(精测)。精测时间与粗测相同,测量精度的提升与插值的倍数相关。
106、坐标获取:根据谱峰数据获取雷达数据对应的空间坐标信息。峰值数据示意如说明书附图5所示。
本实施例采用的方法与现有技术中不同平台MUSIC性能对比如表2所示:(music规格,8通道数,512快拍数,128×128流型阵)
表2性能对比表
由上述实验结果可以看出,采用本实施例提供的方法,在处理时间和设备时延上远远领先现有技术。实际测试测试结果表明,本实施例提供的方法在满负荷工作时,功耗小于20W,远低于同样性能的DSP/GPU(提升3-5倍能效比)。本实施例提供的方法克服了现有的信号处理平台(DSP和常规FPGA)在低功耗、轻重量限制及小型化要求下处理能力不足及测量时间过长的缺陷。
本实施例提供了一种基于硬浮点FPGA的雷达超分辨测向方法,可应用于弹载实时处理系统中。该方法克服了现有技术中的信号处理平台(DSP和常规FPGA)在低功耗、轻重量限制及小型化要求下处理能力不足及测量时间过长的缺陷。采用MUSIC算法,基于Intel硬浮点FPGA进行算法的工程化实现。其中以并行的QR迭代方法实现厄米特矩阵特征值分解,以并行的汇聚插值算法实现谱估计,同时基于特殊的硬浮点FPGA结构大大提升了DOA的计算效率,缩短了计算时间。
实施例2
本发明实施例2公开了一种基于硬浮点FPGA的雷达超分辨测向系统,在实施例1的基础上,将实施例1的方法系统化,使其更具实际应用性。系统的结构图如说明书附图6所示,具体方案如下:
一种基于硬浮点FPGA的雷达超分辨测向系统,包括依次连接的数据获取单元1、协方差计算单元2、特征分解单元3、噪声计算单元4、谱峰搜索计算单元5和坐标获取单元6。
数据获取单元1:用于获取雷达数据,通过多通道形式将雷达数据写入硬浮点FPGA,雷达数据为通道数×快拍数的定点矩阵。数据在写入硬浮点FPGA之前,还包括对雷达数据进行数字下变频处理。
协方差计算单元2:用于通过硬浮点FPGA中计算雷达数据的协方差矩阵。协方差计算单元2计算协方差时为定点模式,计算完成后转为单精度浮点,计算协方差所需的共轭对称矩阵用原始矩阵变换,无需存储。
特征分解单元3:设置有特征分解模块,用于将协方差矩阵由定点模式转换为浮点模式,写入预设的特征分解模块,特征分解模块对协方差矩阵进行基于QR分解的特征分解计算和基于施密特正交的QR分解计算,获取多个特征值和多个特征向量。特征分解模块有QRD(QR分解)迭代构成,迭代的次数可有上位机设置,通常为通道数的1.5倍。特征分解模块的启动信号由协方差矩阵计算完成信号经延时形成。
在本实施例中,通道数为8,特征分解计算时间约为8×1.5×2k(单次QRD时间)=24K时钟周期,完成计算后,输出8个特征值(D1,....D8)和8个特征向量(v1,....v8)。
噪声计算单元4:用于根据预设的信源数门限,参考特征值进行信源数判定,结合特征向量构造噪声矩阵。根据上位机设定的信源数门限,参考8个特征值结果,进行信源数判定,再选出相应的特征向量,形成噪声矩阵。
谱峰搜索计算单元5:用于获取预设的流型阵,根据流型阵和噪声矩阵计算空间谱,基于预设的向量点乘器进行空间谱的计算和搜索,获取谱峰数据。读入DDR(双倍速率同步动态随机存储器)中的流型阵a(128*128矩阵),计算空间谱。整个空间谱的计算和搜索过程基于dot(向量点乘器)全部展开,每个谱峰点的计算/搜索时间为1个周期,整个谱峰生成及搜索时间为128×128=16384周期。
坐标获取单元6:用于根据谱峰数据获取雷达数据对应的空间坐标信息。
此外,谱峰搜索单元5包括初始谱峰获取单元51、谱峰迭代单元52和谱峰数据获取单元53。具体包括:
初始谱峰获取单元51:用于获取预设的流型阵,根据流型阵和噪声矩阵计算空间谱,基于预设的向量点乘器进行空间谱的计算和搜索,获取谱峰;
谱峰迭代单元52:用于将谱峰作为粗测结果,基于粗测结果对流型阵进行汇聚和插值处理,获取更新后的流型阵,根据噪声矩阵和更新后的流型阵重新计算空间谱,对更新后的空间谱进行计算和搜索,获得更新后的谱峰,重复迭代,直至满足迭代停止条件;
谱峰数据获取单元53:用于将满足迭代停止条件的谱峰作为谱峰数据。
本实施例提供了一种基于硬浮点FPGA的雷达超分辨测向系统,在实施例1的基础上,将实施例1的方法系统化,使其更具实际应用性。
实施例3
本发明实施例提供了一种雷达数据处理设备,结构如说明书附图7所示,具体方案如下:
一种雷达数据处理设备,包括计算机12、光纤14和数据计算板卡16,光纤14分别连接计算机12和数据计算板卡16。
数据计算板卡16上设置有实施例2所述的基于硬浮点FPGA的雷达超分辨测向系统。
本实施例提供了一种雷达数据处理设备,将实施例2的系统应用到具体场景。
本发明提供了一种基于FPGA的雷达超分辨测向方法、系统及用于,克服了现有技术中的信号处理平台(DSP和常规FPGA)在低功耗、轻重量限制及小型化要求下处理能力不足及测量时间过长的缺陷。采用MUSIC算法,基于Intel硬浮点FPGA进行算法的工程化实现。其中以并行的QR迭代方法实现厄米特矩阵特征值分解,以并行的汇聚插值算法实现谱估计,同时基于特殊的硬浮点FPGA结构大大提升了DOA的计算效率,缩短了计算时间。
本领域技术人员可以理解附图只是一个优选实施场景的示意图,附图中的模块或流程并不一定是实施本发明所必须的。本领域技术人员可以理解实施场景中的装置中的模块可以按照实施场景描述进行分布于实施场景的装置中,也可以进行相应变化位于不同于本实施场景的一个或多个装置中。上述实施场景的模块可以合并为一个模块,也可以进一步拆分成多个子模块。上述本发明序号仅仅为了描述,不代表实施场景的优劣。以上公开的仅为本发明的几个具体实施场景,但是,本发明并非局限于此,任何本领域的技术人员能思之的变化都应落入本发明的保护范围。
Claims (10)
1.一种基于硬浮点FPGA的雷达超分辨测向方法,其特征在于,包括,
获取雷达数据,通过多通道形式将所述雷达数据写入硬浮点FPGA,所述雷达数据为通道数×快拍数的定点矩阵;
通过所述硬浮点FPGA计算所述雷达数据的协方差矩阵;
将所述协方差矩阵由定点模式转换为浮点模式,写入预设的特征分解模块,所述特征分解模块对所述协方差矩阵进行基于QR分解的特征分解计算和基于施密特正交的QR分解计算,获取多个特征值和多个特征向量;
根据预设的信源数门限,参考所述特征值进行信源数判定,结合所述特征向量构造噪声矩阵;
获取预设的流型阵,根据所述流型阵和所述噪声矩阵计算空间谱,基于预设的向量点乘器进行所述空间谱的计算和搜索,获取谱峰数据;
根据所述谱峰数据获取所述雷达数据对应的空间坐标信息。
2.根据权利要求1所述的方法,其特征在于,所述谱峰数据的获取过程包括,
初始谱峰获取:获取预设的流型阵,根据所述流型阵和所述噪声矩阵计算空间谱,基于预设的向量点乘器进行所述空间谱的计算和搜索,获取谱峰;
谱峰迭代:将所述谱峰作为粗测结果,基于所述粗测结果对所述流型阵进行汇聚和插值处理,获取更新后的流型阵,根据所述噪声矩阵和更新后的流型阵重新计算空间谱,对更新后的空间谱进行计算和搜索,获得更新后的谱峰,重复迭代,直至满足迭代停止条件;
谱峰数据获取:将满足迭代停止条件的谱峰作为所述谱峰数据。
3.根据权利要求1所述的方法,其特征在于,所述空间谱的表达式为:
其中,Pmu表示空间谱,a表示流型阵,En表示噪声矩阵,h表示矩阵共轭转置。
4.根据权利要求1所述的方法,其特征在于,所述硬浮点FPGA包括Intel Arria10系列,开发平台包括QUARTUS。
5.根据权利要求1所述的方法,其特征在于,所述特征分解模块由多个QR分解迭代构成,迭代的次数为所述通道数的1.5倍;
特征分解计算时间的表达式为:
TQR=L×C×2k
其中,TQR表示特征分解计算时间,L表示通道数,C表示迭代次数,k表示单次QR分解的时间。
6.根据权利要求1所述的方法,其特征在于,所述协方差矩阵的表达式为:
A=XXT
其中,A表示协方差矩阵,X表示雷达数据;
对所述协方差矩阵进行QR分解,可得:
A=QR
其中,A表示协方差矩阵,Q表示正交矩阵,R表示上三角矩阵,RQ与A正交相似,具有相同的特征值。
7.根据权利要求2所述的方法,其特征在于,所述预设的流型阵包括128×128的矩阵;
每个谱峰获取的时间为128×128个周期。
8.一种基于硬浮点FPGA的雷达超分辨测向系统,其特征在于,包括如下,
数据获取单元:用于获取雷达数据,通过多通道形式将所述雷达数据写入硬浮点FPGA,所述雷达数据为通道数×快拍数的定点矩阵;
协方差计算单元:用于通过所述硬浮点FPGA中计算所述雷达数据的协方差矩阵;
特征分解单元:用于将所述协方差矩阵由定点模式转换为浮点模式,写入预设的特征分解模块,所述特征分解模块对所述协方差矩阵进行基于QR分解的特征分解计算和基于施密特正交的QR分解计算,获取多个特征值和多个特征向量;
噪声计算单元:用于根据预设的信源数门限,参考所述特征值进行信源数判定,结合所述特征向量构造噪声矩阵;
谱峰搜索计算单元:用于获取预设的流型阵,根据所述流型阵和所述噪声矩阵计算空间谱,基于预设的向量点乘器进行所述空间谱的计算和搜索,获取谱峰数据;
坐标获取单元:用于根据所述谱峰数据获取所述雷达数据对应的空间坐标信息。
9.根据权利要求8所述的系统,其特征在于,所述谱峰搜索计算单元包括,
初始谱峰获取单元:用于获取预设的流型阵,根据所述流型阵和所述噪声矩阵计算空间谱,基于预设的向量点乘器进行所述空间谱的计算和搜索,获取谱峰;
谱峰迭代单元:用于将所述谱峰作为粗测结果,基于所述粗测结果对所述流型阵进行汇聚和插值处理,获取更新后的流型阵,根据所述噪声矩阵和更新后的流型阵重新计算空间谱,对更新后的空间谱进行计算和搜索,获得更新后的谱峰,重复迭代,直至满足迭代停止条件;
谱峰数据获取单元:用于将满足迭代停止条件的谱峰作为所述谱峰数据。
10.一种雷达数据处理设备,其特征在于,包括计算机、光纤和数据计算板卡,所述光纤分别连接所述计算机和所述数据计算板卡;
所述数据计算板卡上设置有权利要求8或9所述的基于硬浮点FPGA的雷达超分辨测向系统。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110406848.4A CN113203997B (zh) | 2021-04-15 | 2021-04-15 | 基于fpga的雷达超分辨测向方法、系统及应用 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110406848.4A CN113203997B (zh) | 2021-04-15 | 2021-04-15 | 基于fpga的雷达超分辨测向方法、系统及应用 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113203997A CN113203997A (zh) | 2021-08-03 |
CN113203997B true CN113203997B (zh) | 2023-11-03 |
Family
ID=77027077
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110406848.4A Active CN113203997B (zh) | 2021-04-15 | 2021-04-15 | 基于fpga的雷达超分辨测向方法、系统及应用 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113203997B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115656961B (zh) * | 2022-12-26 | 2023-03-10 | 南京楚航科技有限公司 | 一种基于并行处理器的os-cfar处理方法及系统 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5990834A (en) * | 1997-08-29 | 1999-11-23 | Codar Ocean Sensors, Ltd. | Radar angle determination with music direction finding |
CN104133200A (zh) * | 2014-07-30 | 2014-11-05 | 西安电子科技大学 | 基于fpga的正交匹配追踪方法 |
CN104155648A (zh) * | 2014-08-26 | 2014-11-19 | 国家海洋局第一海洋研究所 | 基于阵列数据重排的高频地波雷达单次快拍music测向方法 |
CN106646345A (zh) * | 2016-09-09 | 2017-05-10 | 西北工业大学 | 基于数据域的测向误差补偿方法 |
CN107064901A (zh) * | 2017-04-27 | 2017-08-18 | 哈尔滨工业大学 | 一种舰载高频地波超视距雷达的目标方位估计方法 |
CN109633525A (zh) * | 2018-11-29 | 2019-04-16 | 上海无线电设备研究所 | 一种快速的均匀圆阵空间谱测向实现方法 |
CN111308430A (zh) * | 2020-02-25 | 2020-06-19 | 西安电子科技大学 | 基于fpga及多核dsp硬件架构的测向及抗干扰方法及系统 |
CN111896956A (zh) * | 2020-08-10 | 2020-11-06 | 西安电子科技大学 | 基于fpga和dsp的实时微波关联成像装置及方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6567034B1 (en) * | 2001-09-05 | 2003-05-20 | Lockheed Martin Corporation | Digital beamforming radar system and method with super-resolution multiple jammer location |
JP4679632B2 (ja) * | 2008-12-01 | 2011-04-27 | 三菱電機株式会社 | レーダ装置 |
JP5675285B2 (ja) * | 2010-11-10 | 2015-02-25 | 富士通テン株式会社 | レーダ装置 |
US20170054449A1 (en) * | 2015-08-19 | 2017-02-23 | Texas Instruments Incorporated | Method and System for Compression of Radar Signals |
JP6570675B2 (ja) * | 2018-02-15 | 2019-09-04 | 三菱電機株式会社 | レーダ信号処理装置 |
-
2021
- 2021-04-15 CN CN202110406848.4A patent/CN113203997B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5990834A (en) * | 1997-08-29 | 1999-11-23 | Codar Ocean Sensors, Ltd. | Radar angle determination with music direction finding |
CN104133200A (zh) * | 2014-07-30 | 2014-11-05 | 西安电子科技大学 | 基于fpga的正交匹配追踪方法 |
CN104155648A (zh) * | 2014-08-26 | 2014-11-19 | 国家海洋局第一海洋研究所 | 基于阵列数据重排的高频地波雷达单次快拍music测向方法 |
CN106646345A (zh) * | 2016-09-09 | 2017-05-10 | 西北工业大学 | 基于数据域的测向误差补偿方法 |
CN107064901A (zh) * | 2017-04-27 | 2017-08-18 | 哈尔滨工业大学 | 一种舰载高频地波超视距雷达的目标方位估计方法 |
CN109633525A (zh) * | 2018-11-29 | 2019-04-16 | 上海无线电设备研究所 | 一种快速的均匀圆阵空间谱测向实现方法 |
CN111308430A (zh) * | 2020-02-25 | 2020-06-19 | 西安电子科技大学 | 基于fpga及多核dsp硬件架构的测向及抗干扰方法及系统 |
CN111896956A (zh) * | 2020-08-10 | 2020-11-06 | 西安电子科技大学 | 基于fpga和dsp的实时微波关联成像装置及方法 |
Non-Patent Citations (1)
Title |
---|
基于CUDA架构的混合测向算法研究;金莉,宋万杰,施治国;雷达科学与技术;第16卷(第2期);第133-138页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113203997A (zh) | 2021-08-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110361691B (zh) | 基于非均匀阵列的相干信源doa估计fpga实现方法 | |
CN1211755C (zh) | 自适应滤波的系统和方法 | |
CN102680945B (zh) | 基于fpga的多普勒调频率估计方法 | |
CN107907855A (zh) | 一种互素阵列转化为均匀线阵的doa估计方法及装置 | |
CN110398732B (zh) | 低计算量自适应步长迭代搜索的目标方向检测方法 | |
CN104535067B (zh) | 一种基于分区搜索的脉冲信号到达时间快速计算方法 | |
CN114488064A (zh) | 一种距离-速度联合估计方法 | |
Li et al. | Hardware acceleration of MUSIC algorithm for sparse arrays and uniform linear arrays | |
CN113203997B (zh) | 基于fpga的雷达超分辨测向方法、系统及应用 | |
CN105608057A (zh) | 一种分时复用硬件资源的信号子空间分解的fpga实现模块及其fpga实现方法 | |
CN113391261A (zh) | 一种基于泰勒展开的增广互质阵列快速doa估计方法 | |
CN110471041B (zh) | 一种基于Vivado HLS的二维DOA估计方法 | |
CN111337896A (zh) | 一种实现动目标检测加速的方法 | |
CN113093111B (zh) | 基于压缩感知和遗传算法的均匀圆阵解调二维相干信号方法及系统 | |
CN114563760B (zh) | 一种基于sca阵型的二阶超波束形成方法、设备及介质 | |
CN113341371B (zh) | 一种基于l阵和二维esprit算法的doa估计方法 | |
Yang et al. | A efficient design of a real-time FFT architecture based on FPGA | |
CN109683128B (zh) | 冲击噪声环境下的单快拍测向方法 | |
CN113552533A (zh) | 一种空间谱估计方法、装置、电子设备及存储介质 | |
CN113238184A (zh) | 一种基于非圆信号的二维doa估计方法 | |
CN113221059B (zh) | 无需构造协方差矩阵的快速共轭梯度测向算法 | |
Jiang et al. | The DOA estimation for a mixture of uncorrelated and coherent sources via decomposing the coprime array | |
Zhao et al. | FPGA-Based CNN Quantization Method and Hardware Acceleration Module | |
CN115291160B (zh) | 二维doa估计方法、系统及计算机可读介质 | |
Ma et al. | Projection Matching Pursuit based DoA Estimation Scheme and its FPGA Implementation |
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 |