CN101424738B - 基于fpga的实时数字脉冲压缩系统的处理方法 - Google Patents
基于fpga的实时数字脉冲压缩系统的处理方法 Download PDFInfo
- Publication number
- CN101424738B CN101424738B CN2008101976326A CN200810197632A CN101424738B CN 101424738 B CN101424738 B CN 101424738B CN 2008101976326 A CN2008101976326 A CN 2008101976326A CN 200810197632 A CN200810197632 A CN 200810197632A CN 101424738 B CN101424738 B CN 101424738B
- Authority
- CN
- China
- Prior art keywords
- ram
- fft
- pulse compression
- result
- dft
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Compression, Expansion, Code Conversion, And Decoders (AREA)
Abstract
本发明涉及一种基于FPGA的实时数字脉冲压缩系统的处理方法,采用的方法是:中频信号经过A/D采样后进行幅度判决,然后按以下步骤进行:FFT处理;匹配乘法处理:IFFT处理;输出脉压结果。本发明采用上述方法后,对4096点数据进行频域数字脉冲压缩,选取8个碟形算子、32个RAM(128*32bit)和48个ROM(934*32bit)进行运算,若系统处理速度为100MHz,完成整个压缩处理只需要49us,其中FFT用时23040ns,综合滤波用时2560ns,IFFT用时23040ns。可见,采用该优化设计方法可以大大提高处理速度,并且利用FPGA处理可以减小系统的功耗,在航天领域有较大的使用价值和较好的应用前景。
Description
技术领域
本发明属于雷达信号处理方法,涉及大时宽带宽积的线性调频数字脉冲压缩系统在频域内的实现方法,具体地讲是一种基于FPGA的实时数字脉冲压缩系统的处理方法。
背景技术
脉冲压缩技术因解决了雷达作用距离与分辨力之间的矛盾而成为现代雷达的一种重要体制。现今数字脉冲压缩系统大多采用频域处理方法,脉冲压缩处理的计算量大,对实时性要求较高,一般采取较多的DSP的硬件平台实现,硬件系统复杂,功耗大,不易广泛应用。利用FPGA处理大时宽带宽积的雷达脉冲信号时往往需要较长的处理时间,很难满足某些场合下的实时处理要求。
发明内容
本发明的目的是提供一种使大时宽带宽积的脉冲压缩技术易于工程化,通过一系列优化设计方法在有限的资源下实现实时压缩处理的基于FPGA的实时数字脉冲压缩系统的处理方法。
为了实现上述目的,本发明所采用的方法是:中频信号经过A/D采样后进行幅度判决,然后按以下步骤进行处理:
第一步骤:FFT处理
用高速时钟采样后得到有用信号,将其补0到4096点后放入32个RAM中进行后面的处理:先将4096点FFT输入数据生成4个1024点的序列,再分别对这4个序列进行1024点的FFT,而对于每一个1024点的FFT而言,又可以继续分成4个256点的FFT,这样继续一直分下去,通过6轮的计算最后就只剩4点的DFT了,直接计算即可;
第二步骤:匹配乘法(频谱搬移和综合滤波)处理
当FFT计算结束后,4096点的结果已经存放在32个RAM里面,接下来进行的就是频域的匹配乘法处理,匹配滤波器的系数已经存放在24个ROM里面,在计算过程中,直接将RAM里面的数据取出来与ROM取出来的系数直接相乘即可得到计算结果,并且存放于32个RAM里面;
第三步骤:IFFT处理
频域匹配乘法结束之后得到的32个RAM中的序列,输入序列按照基4抽取的规律分解成4组,整个4096点的DFT经过6轮的分解后会得到4点的DFT的计算,计算出第1轮的1024个4点的DFT后,将结果重新排列后进行第2轮的1024个4点的DFT计算,一直这样直到6轮计算结束,得到的结果就是最后脉冲压缩输出系列;
第四步骤:输出脉冲压缩结果
IFFT处理结束之后,脉冲压缩结果已经存放在32个RAM之中,需要将结果顺序从RAM中读出来送给输出,并且在输出结束之后将32个RAM内容清零已准备进行下一次的脉冲压缩处理。
本发明采用上述优化设计方法后,对4096点数据进行频域数字脉冲压缩,选取8个碟形算子、32个RAM(128*32bit)和48个ROM(934*32bit)进行运算,若系统处理速度为100MHz,完成整个压缩处理只需要49us,其中FFT用时23040ns,综合滤波用时2560ns,IFFT用时23040ns。可见,采用该优化设计方法可以大大提高处理速度,并且利用FPGA处理可以减小系统的功耗,在航天领域有较大的使用价值和较好的应用前景。
附图说明
图1为目前目标识别信号处理框图。
图2为本发明脉冲压缩整体流程图。
图3为本发明FFT处理流程图。
图4为本发明频域匹配乘法处理流程图。
图5为本发明IFFT处理流程图。
具体实施方式
下面结合附图对本发明作进一步的详细描述
脉冲压缩技术因解决了雷达作用距离与分辨力之间的矛盾,在精确制导领域具有广泛的应用前景。目前,目标自动识别常用步进跳频脉冲压缩体制、线性调频脉冲压缩体制和线性调频加步进跳频二次脉冲压缩体制,其信号处理过程如图1所示。其中,脉冲压缩主要用来提高信号检测性能,从而提高目标识的效果。
本发明基于FPGA处理平台,实现了实时数字脉冲压缩系统的优化设计。其处理流程如图2所示,中频信号经过A/D采样后进行幅度判决,然后进行、FFT、匹配乘法(频谱搬移和综合滤波)和IFFT等处理,输出脉冲压缩结果。
其中,FFT处理流程如图3所示。用高速时钟采样后得到有用信号,将其补0到4096点后放入32个RAM中进行后面的脉冲压缩处理。先将4096点FFT输入数据生成4个1024点的序列,再分别对这4个序列进行1024点的FFT。而对于每一个1024点的FFT而言,又可以继续分成4个256点的FFT,这样继续一直分下去,通过6轮的计算最后就只剩4点的DFT了,直接计算即可。
当FFT计算结束后,4096点的结果已经存放在32个RAM里面,接下来进行的就是频域的匹配乘法处理,匹配滤波器的系数已经存放在24个ROM里面,在计算过程中,直接将RAM里面的数据取出来与ROM取出来的系数直接相乘即可得到计算结果,并且存放于32个RAM里面。频域匹配乘法处理流程图如图4所示。
频域匹配乘法结束之后得到的32个RAM中的序列,输入序列按照基4抽取的规律分解成4组,整个4096点的DFT经过6轮的分解后会得到4点的DFT的计算。计算出第1轮的1024个4点的DFT后,将结果重新排列后进行第2轮的1024个4点的DFT计算,一直这样直到6轮计算结束,得到的结果就是最后脉冲压缩输出系列。IFFT处理的流程图如图5所示。
IFFT处理结束之后,脉冲压缩结果已经存放在32个RAM之中,需要将结果顺序从RAM中读出来送给输出,并且在输出结束之后将32个RAM内容清零已准备进行下一次的脉冲压缩处理。
整个脉冲压缩算法的具体实现上采用了一系列的优化设计方法。首先,优化脉冲压缩处理流程。系统始终都是一路在进行处理,只是在最后IFFT之后才将I,Q两路分离出来,这样,对于系统资源的需求大大减少,同时也可以减少总体的处理时间。
第二,采用流水线和并行处理技术。脉冲压缩算法在FPGA上实现时充分优化资源,压缩过程中采用并行的流水线处理,资源充分复用。
第三,优化时序控制方法。系统中有专门的RAM读地址控制和写地址控制模块,ROM读地址控制模块和碟形算子模块,由三级状态信号控制相应地址的产生,第一级状态信号有FFT状态、综合滤波状态、IFFT状态;第二级状态信号包括FFT和IFFT的轮次运算状态信号;第三级状态信号包括FFT和IFFT每轮运算状态下的N次运算状态信号。状态信号由RAM读地址控制模块产生,并由它传递给RAM写地址控制模块、ROM读地址控制模块和碟形算子模块,以控制模块的时序动作,这样便于时序控制和并行运算,加快处理速度。
第四,优化蝶形算子设计方法。在FFT、综合滤波和IFFT三种状态下都需要用到乘法器,系统通过状态信号来控制乘法器的输出数据,并且将计算结果进行相应的处理,就可以在不同的状态下使用同样的乘法器完成不同的功能,这样既节省了资源又加快了运行速度。另外,将4个乘法器输出的复数分成两组后分别相加或相减,再将结果相加,最后进行复数乘法,此方法相对于顺序相加所需时钟周期减少,而且时序稳定。
第五,采用滤波器合并技术。因为脉冲压缩处理需要进行低通滤波、匹配滤波和加权滤波,本系统在处理时把三个滤波器合并成一个综合滤波器,在频域滤波的时候可以有效减少处理时间,节省资源。
本发明已成功用于所研制的数字脉冲压缩样机,具用良好的压缩效果。对于时宽带宽积为500(时宽为10us、带宽为50MHz)、中频为150MHz的线性调频信号,在一片Xilinx公司的Virtex-4FPGA芯片上,采用该优化设计方法可以在111us内完成一个周期脉冲信号的压缩处理。压缩后的脉冲信号时宽为34ns,主旁瓣比优于-35dB,在节省处理时间的基础上又能达到所需要的技术指标和其他技术要求。
本说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
Claims (1)
1.一种基于FPGA的实时数字脉冲压缩系统的处理方法,采用的方法是:中频信号经过A/D采样后进行幅度判决,然后按以下步骤进行:
第一步骤:FFT处理
用高速时钟采样后得到有用信号,将其补0到4096点后放入32个RAM中进行后面的处理:先将4096点FFT输入数据生成4个1024点的序列,再分别对这4个序列进行1024点的FFT,而对于每一个1024点的FFT而言,又可以继续分成4个256点的FFT,这样继续一直分下去,通过6轮的计算最后就只剩4点的DFT了,直接计算即可;
第二步骤:匹配乘法处理
当FFT计算结束后,4096点的结果已经存放在32个RAM里面,接下来进行的就是频域的匹配乘法处理,匹配滤波器的系数已经存放在24个ROM里面,在计算过程中,直接将RAM里面的数据取出来与ROM取出来的系数直接相乘即可得到计算结果,并且存放于32个RAM里面;
第三步骤:IFFT处理
频域匹配乘法结束之后得到的32个RAM中的序列,输入序列按照基4抽取的规律分解成4组,整个4096点的DFT经过6轮的分解后会得到4点的DFT的计算,计算出第1轮的1024个4点的DFT后,将结果重新排列后进行第2轮的1024个4点的DFT计算,一直这样直到6轮计算结束,得到的结果就是最后脉冲压缩输出系列;
第四步骤:输出脉冲压缩结果
IFFT处理结束之后,脉冲压缩结果已经存放在32个RAM之中,需要将结果顺序从RAM中读出来送给输出,并且在输出结束之后将32个RAM内容清零已准备进行下一次的脉冲压缩处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008101976326A CN101424738B (zh) | 2008-11-13 | 2008-11-13 | 基于fpga的实时数字脉冲压缩系统的处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008101976326A CN101424738B (zh) | 2008-11-13 | 2008-11-13 | 基于fpga的实时数字脉冲压缩系统的处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101424738A CN101424738A (zh) | 2009-05-06 |
CN101424738B true CN101424738B (zh) | 2011-06-01 |
Family
ID=40615466
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008101976326A Expired - Fee Related CN101424738B (zh) | 2008-11-13 | 2008-11-13 | 基于fpga的实时数字脉冲压缩系统的处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101424738B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102353940B (zh) * | 2011-06-10 | 2013-06-12 | 西安电子科技大学 | 基于fpga的脉冲压缩优化方法 |
CN103926567B (zh) * | 2013-01-12 | 2017-03-22 | 安徽博微长安电子有限公司 | 高速实时脉冲压缩算法 |
CN103197287A (zh) * | 2013-03-25 | 2013-07-10 | 西安电子科技大学 | 高速实时频域脉冲压缩器及其处理方法 |
CN103529430B (zh) * | 2013-10-29 | 2015-11-04 | 北京无线电计量测试研究所 | 一种用于线性调频信号的匹配压缩装置及匹配压缩方法 |
CN103529429B (zh) * | 2013-10-29 | 2016-08-17 | 北京无线电计量测试研究所 | 用于线性调频信号的频域中的脉冲压缩方法 |
CN103808809B (zh) * | 2014-03-12 | 2016-02-17 | 哈尔滨工业大学 | 一种基于线性调频信号激励的超声回波信号实时处理装置及处理方法 |
CN106842144A (zh) * | 2015-12-03 | 2017-06-13 | 中国航空工业集团公司雷华电子技术研究所 | 并行多相结构数字脉压方法 |
CN107942298B (zh) * | 2017-11-03 | 2019-10-11 | 西安电子科技大学 | 一种用户可配置参数的低速目标检测器及方法 |
CN108132460B (zh) * | 2017-12-07 | 2021-07-16 | 西安电子工程研究所 | 基于频域通道均衡的脉冲压缩补偿算法 |
CN109614151B (zh) * | 2018-11-14 | 2023-02-28 | 上海无线电设备研究所 | 一种四核并行的大点数脉压数据处理方法 |
CN109490849A (zh) * | 2018-11-27 | 2019-03-19 | 重庆秦嵩科技有限公司 | 一种适用于大脉宽信号脉冲压缩的fpga实现方法 |
CN110058201B (zh) * | 2019-04-22 | 2022-11-22 | 西安电子工程研究所 | 一种基于fpga资源复用实现多波形多速率时域脉压的方法 |
CN110308427B (zh) * | 2019-06-30 | 2023-07-18 | 南京理工大学 | 基于fpga的lfm脉冲雷达频域脉冲压缩处理方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5991454A (en) * | 1997-10-06 | 1999-11-23 | Lockheed Martin Coporation | Data compression for TDOA/DD location system |
CN1982914A (zh) * | 2005-12-15 | 2007-06-20 | 四川川大智胜软件股份有限公司 | 基于译码和报文识别的雷达数据自适应无损压缩方法 |
CN101165510A (zh) * | 2006-10-18 | 2008-04-23 | 中国科学院电子学研究所 | 星载合成孔径雷达可变位数baq压缩系统及方法 |
-
2008
- 2008-11-13 CN CN2008101976326A patent/CN101424738B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5991454A (en) * | 1997-10-06 | 1999-11-23 | Lockheed Martin Coporation | Data compression for TDOA/DD location system |
CN1982914A (zh) * | 2005-12-15 | 2007-06-20 | 四川川大智胜软件股份有限公司 | 基于译码和报文识别的雷达数据自适应无损压缩方法 |
CN101165510A (zh) * | 2006-10-18 | 2008-04-23 | 中国科学院电子学研究所 | 星载合成孔径雷达可变位数baq压缩系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101424738A (zh) | 2009-05-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101424738B (zh) | 基于fpga的实时数字脉冲压缩系统的处理方法 | |
CN103973324B (zh) | 一种宽带数字接收机及其实时频谱处理方法 | |
CN102288941A (zh) | 一种基于fpga和dsp的中频lfm-pd雷达信号实时处理系统及处理方法 | |
CN100543492C (zh) | 多通道高频雷达接收机的数字信号处理方法 | |
CN107153206A (zh) | 一种基于fft的北斗卫星信号快速捕获方法 | |
CN103154932A (zh) | 用于分析信号、提供瞬时频率和短时傅里叶变换的方法以及用于分析信号的设备 | |
CN104467739A (zh) | 一种带宽、中心频点可调的数字滤波器及其实现方法 | |
CN101858938A (zh) | 基于自适应滤波原理的瞬时频率测量方法 | |
CN102882814A (zh) | 参数化模块化多通道数字下变频设计平台和方法 | |
CN111766597A (zh) | 调频连续波激光测距差拍信号测频方法及装置 | |
CN107544942A (zh) | 一种快速傅里叶变换的vlsi设计方法 | |
CN105807264A (zh) | 雷达脉冲重复频率检测与初始脉冲到达时间的估计方法 | |
CN103969508B (zh) | 一种实时高精密的电力谐波分析方法及装置 | |
CN208190613U (zh) | 一种基于fpga实现的分数阶积分器 | |
CN102111372B (zh) | 一种基于脉冲抵消方式的削峰方法 | |
CN112350691A (zh) | 一种基于多级滤波抽取和两级wola结构的地磁弱信号检测算法 | |
CN103744066A (zh) | 数字正交鉴相与匹配滤波的优化方法 | |
CN110208755A (zh) | 一种基于fpga的动态雷达回波数字下变频系统及方法 | |
CN204290909U (zh) | 一种带宽、中心频点可调的数字滤波器 | |
CN112526208B (zh) | 一种基于高耦合度迭代模型的电能质量测量系统及方法 | |
CN108919200A (zh) | 一种基于fpga的任意雷达波形在线设计方法 | |
CN103487806B (zh) | 一种基于时分复用的多普勒参数二次拟合方法 | |
CN102981156A (zh) | 一种超声成像后处理方法及装置 | |
Hu | Design and verification of FIR filter based on Matlab and DSP | |
CN106100588A (zh) | 一种基于fpga的可重构多通道数字下变频系统及方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110601 Termination date: 20141113 |
|
EXPY | Termination of patent right or utility model |