CN114626006B - 雷达成像中cs算法补偿因子实时生成的fpga实现方法 - Google Patents

雷达成像中cs算法补偿因子实时生成的fpga实现方法 Download PDF

Info

Publication number
CN114626006B
CN114626006B CN202210274459.5A CN202210274459A CN114626006B CN 114626006 B CN114626006 B CN 114626006B CN 202210274459 A CN202210274459 A CN 202210274459A CN 114626006 B CN114626006 B CN 114626006B
Authority
CN
China
Prior art keywords
factor
real
algorithm
multiplication
fft
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
Application number
CN202210274459.5A
Other languages
English (en)
Other versions
CN114626006A (zh
Inventor
闵锐
李晋
徐浩典
曹雨欣
余雷
皮亦鸣
杨晓波
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN202210274459.5A priority Critical patent/CN114626006B/zh
Publication of CN114626006A publication Critical patent/CN114626006A/zh
Application granted granted Critical
Publication of CN114626006B publication Critical patent/CN114626006B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • G06F17/141—Discrete Fourier transforms
    • G06F17/142—Fast Fourier transforms, e.g. using a Cooley-Tukey type algorithm
    • 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
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9021—SAR image post-processing techniques
    • 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/40—Means for monitoring or calibrating
    • G01S7/4052—Means for monitoring or calibrating by simulation of echoes
    • 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
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F9/00—Arrangements for program control, e.g. control units
    • G06F9/06—Arrangements for program control, e.g. control units using stored programs, i.e. using an internal store of processing equipment to receive or retain programs
    • G06F9/30—Arrangements for executing machine instructions, e.g. instruction decode
    • G06F9/30003—Arrangements for executing specific machine instructions
    • G06F9/30007—Arrangements for executing specific machine instructions to perform operations on data operands
    • G06F9/3001—Arithmetic instructions
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F9/00—Arrangements for program control, e.g. control units
    • G06F9/06—Arrangements for program control, e.g. control units using stored programs, i.e. using an internal store of processing equipment to receive or retain programs
    • G06F9/30—Arrangements for executing machine instructions, e.g. instruction decode
    • G06F9/30003—Arrangements for executing specific machine instructions
    • G06F9/30007—Arrangements for executing specific machine instructions to perform operations on data operands
    • G06F9/30025—Format conversion instructions, e.g. Floating-Point to Integer, decimal conversion
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F9/00—Arrangements for program control, e.g. control units
    • G06F9/06—Arrangements for program control, e.g. control units using stored programs, i.e. using an internal store of processing equipment to receive or retain programs
    • G06F9/30—Arrangements for executing machine instructions, e.g. instruction decode
    • G06F9/30003—Arrangements for executing specific machine instructions
    • G06F9/30007—Arrangements for executing specific machine instructions to perform operations on data operands
    • G06F9/30036—Instructions to perform operations on packed data, e.g. vector, tile or matrix operations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F9/00—Arrangements for program control, e.g. control units
    • G06F9/06—Arrangements for program control, e.g. control units using stored programs, i.e. using an internal store of processing equipment to receive or retain programs
    • G06F9/30—Arrangements for executing machine instructions, e.g. instruction decode
    • G06F9/38—Concurrent instruction execution, e.g. pipeline or look ahead
    • G06F9/3867—Concurrent instruction execution, e.g. pipeline or look ahead using instruction pipelines
    • G06F9/3869—Implementation aspects, e.g. pipeline latches; pipeline synchronisation and clocking
    • 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)
  • Theoretical Computer Science (AREA)
  • Software Systems (AREA)
  • Mathematical Physics (AREA)
  • Remote Sensing (AREA)
  • General Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Discrete Mathematics (AREA)
  • Electromagnetism (AREA)
  • Computing Systems (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明属于雷达成像信号和FPGA技术领域，具体涉及一种雷达成像中CS算法补偿因子实时生成的FPGA实现方法。本发明主要是提出了基于FPGA实现的CS算法补偿因子实时生成模块，该模块由21级流水线结构组成，每一层流水结构主要利用Floating Point IP核和Cordic IP核对输入参数进行快速计算，实现CS算法补偿因子的实时生成，解决了在高速时钟下，传统DSP方案计算时间长，无法满足***成像实时性要求的问题。

Description

雷达成像中CS算法补偿因子实时生成的FPGA实现方法
技术领域
本发明属于雷达成像信号和FPGA技术领域，具体涉及一种雷达成像中CS算法补偿因子实时生成的FPGA实现方法。
背景技术
合成孔径雷达(Synthetic Aperture Radar,SAR)具有全天时，全天候的工作特点以及高分辨的成像精度，其已经在多云雾地区遥感测图、军事侦察、国民经济建设等方面发挥着巨大作用。近些年来，随着硬件制造水平的快速发展，SAR实时成像***设计收到越来越多的研究。而SAR成像算法中CS算法利用相位相乘替代了矩阵差值，避免了复杂的运算，还能较好地保持图像的相位信息，具有良好的成像效果。CS算法涉及到大量相位因子的计算，由于***需要满足实时性的要求，因此CS算法相位因子计算效率会影响到CS算法成像处理的快慢。
在早期基于CS算法的SAR成像处理***中，多采用DSP作为处理器进行浮点数运算，但随着要求的不断提高，DSP逐渐不能满足实时性的要求。如今FPGA能够快速完成加减乘除、三角函数等成像处理中常用的基本运算，因此基于FPGA的CS算法乘法因子实时生成方案具有极高的数据处理速率。
中国科学院计算技术研究所简方军、褚超、韩承德发表的论文“Chirp Scaling成像算法中Cs和Bro数值计算与硬件设计”提出了一种新的易于硬件实现的多项式逼近计算方法，并给出了该方法的计算精度及误差分析。该文提出的方法降低了计算复杂度，实测成像质量与原计算相当，但该***对方位样本点数16384、距离样本点数16384的SAR图像进行处理，需要花费8秒完成三个因子的全部计算，无法满足计算的实时性要求。
乘法因子生成模块接收景中心距离R_nc、雷达有效速度Vr、发射脉冲时宽Tr、雷达工作频率f0、距离向采样点数Nr、方位向采样点数Na、波束斜视角theta、距离向过采样系数alpha、方位向过采样系数belta、带宽Bw、天线方位向长度La，用于如下三个相位因子的计算：
Figure BDA0003555282720000011
其中s_sc用于第一次相位相乘实现Chirp Scaling操作，Km为变换到距离多普勒域的距离调频率，D_fn_ref_Vr为参考频率处的徙动因子，D_fn_Vr_mtx为形成矩阵，tr_mtx为距离时间轴矩阵，R_ref为参考目标斜距，c为光速。
Figure BDA0003555282720000021
其中H1用于第二次相位相乘实现距离压缩、SRC和一致RCMC操作，fr_mtx为距离频率轴矩阵。
Figure BDA0003555282720000022
Figure BDA0003555282720000023
其中Haz*H2用于第三次相位相乘实现方位压缩与相位校正操作，Haz为方位向匹配滤波器，H2为附加相位矫正项，R0_RCMC为随距离变化的斜距。
发明内容
针对上述问题，为了满足计算的实时性要求，本发明提供了一种雷达成像中CS算法补偿因子实时生成的FPGA实现方法。
本发明的技术方案为：
雷达成像中CS算法补偿因子实时生成的FPGA实现方法，其特征在于，包括以下步骤：
S1、获取雷达回波数据；
S2、采用CS算法对雷达回波数据进行处理，获得雷达成像图像，CS算法对雷达回波数据的处理过程包括四次FFT和三次相位相乘，第一次FFT为方位向FFT，将数据变换到距离多普勒域，采用Chirp Scaling相位因子进行第一次相位相乘，使所有目标的徙动轨迹一致化；进行第二次FFT，第二次FFT为距离向FFT，将数据变换到二维频域，采用RCMC和距离向压缩因子进行第二次相位相乘，完成距离压缩、SRC和一致RCMC；进行第三次FFT，第三次FFT为距离向IFFT，将数据变回距离多普勒域，采用相位校正和方位压缩因子进行第三次相位相乘；进行第四次FFT，第四次FFT为方位向IFFT，获得输出图像；
所述的Chirp Scaling相位因子、RCMC和距离向压缩因子、相位校正和方位压缩因子为CS算法乘法因子，CS算法乘法因子为通过CS算法补偿因子实时生成模块实时生成，CS算法补偿因子实时生成模块为基于FPGA的方式实现，具体为：采用21级流水线的结构，每一层流水结构利用Floating Point IP核和Cordic IP核对输入的雷达信号参数进行计算，从而获得包括Chirp Scaling相位因子、RCMC和距离向压缩因子、相位校正和方位压缩因子在内的CS算法乘法因子；CS算法补偿因子实时生成模块的实时生成控制方式为，通过控制信号端口接收FFT模块输出的有效信号标志位，在输出端口设置对应的因子输出有效信号标志位，从而在每一次相位相乘时根据上一次FFT的具体操作获得下一次相位相乘需要的对应CS算法乘法因子。
具体的，CS算法乘法因子包括两种数据输出顺序，分别为：
a.以方位向顺序输出：
a1,1、a1,2、a1,3、a1,4、
a2,1、a2,2、a2,3、a2,4
……
a8191,1、a8191,2、a8191,3、a8191,4
a8192,1、a8192,2、a8192,3、a8192,4
四路因子，直到当前四列方位向因子输出完毕，再输出下四列方位向因子：
a1,5、a1,6、a1,7、a1,8
a2,5、a2,6、a2,7、a2,8
……
a8191,5、a8191,6、a8191,7、a8191,8
a8192,5、a8192,6、a8192,7、a8192,8
b.以距离向顺序输出：
a1,1、a2,1、a3,1、a4,1
a1,2、a2,2、a3,2、a4,2
……
a1，8191、a2，8191、a3,8191、a4,8191
a1，8192、a2，8192,、a3，8192,、a4，8192
四路因子，直到当前四列距离向因子输出完毕，再输出下四列距离向因子：
a1,1、a2,1、a3,1、a4,1
a1,2、a2,2、a3,2、a4,2
……
a1，8191、a2，8191、a3,8191、a4,8191
a1，8192、a2，8192、a3，8192、a4，8192
直至完成一帧数据对应的因子输出为止。
上述方案中，CS算法乘法因子实时生成模块参数输入位宽均为32bits，输出因子位宽为32bit，模块中间部分均以64bits位宽进行数据处理，该模块时钟为300MHz。乘法因子生成模块结构为21级流水线结构，每一级流水线有固定的首次延时，每一级流水线均例化了大量Floating Point IP核与Cordic IP核用于数学运算，Floating Point IP核与Cordic IP核均采用pipeline流水线结构，整个***首次延时为21级流水线每一级首次延时的总和。其中Floating Point IP核例化为Add、Subtract、Multiply、Divide、Square-root、Fixed-to-Float、Float-to-Fixed、Float-to-Float等模式以实现数据的加减乘除、开方、平方、定点数浮点数相互转化等功能，Cordic IP核例化为Sin and Cos模式以实现求正余弦功能。
本发明的有益效果为，本发明实现CS算法补偿因子的实时生成，解决了在高速时钟下，传统DSP方案计算时间长，无法满足***成像实时性要求的问题。
附图说明
图1为本发明的方法流程图；
图2为乘法因子生成模块示意图；
图3为输出因子选择方式示意图；
图4为四路因子生成模块示意图；
图5为乘法因子方位向输出顺序示意图；
图6为乘法因子距离向输出顺序示意图。
具体实施方式
下面结合附图对本发明进行详细的描述；
本发明提出一种雷达成像信号处理中基于CS算法补偿因子实时生成的FPGA实现方式，该乘法因子生成模块接收雷达参数，并在FFT模块等控制信号控制下，对参数进行计算，生成乘法因子并输给复乘模块进行相位相乘。
具体流程如图1所述，CS算法主要步骤包括四次FFT和三次相位相乘，三次相位相乘所需因子由乘法因子实时生成模块产生，具体为：
通过方位向FFT将数据变换到距离多普勒域后，乘法因子实时生成模块输出ChirpScaling相位因子，通过相位相乘实现Chirp Scaling操作，使所有目标的徙动轨迹一致化。这是第一步相位相乘。
通过距离向FFT将数据变换到二维频域后；乘法因子实时生成模块输出第二个相位因子相位相乘，完成对数据的距离压缩、SRC和一致RCMC。这是第二步相位相乘。
通过距离向IFFT将数据变回距离多普勒域后；乘法因子实时生成模块输出第三个相位因子相位相乘进行相位相乘，实现方位压缩与相位校正，这是最后一次相位相乘。
乘法因子生成模块框架如图2所示，整个乘法因子生成模块结构为21级流水线结构，每一级流水线有固定的首次延时，每一级流水线均例化了大量Floating Point IP核与Cordic IP核用于数学运算，Floating Point IP核与Cordic IP核均采用pipeline流水线结构，整个***首次延时为21级流水线每一级首次延时的总和。其中Floating Point IP核例化为Add、Subtract、Multiply、Divide、Square-root、Fixed-to-Float、Float-to-Fixed、Float-to-Float等模式以实现数据的加减乘除、开方、平方、定点数浮点数相互转化等功能，Cordic IP核例化为Sin and Cos模式以实现求正余弦功能。
乘法因子生成模块接收景中心距离R_nc、雷达有效速度Vr、发射脉冲时宽Tr、雷达工作频率f0、距离向采样点数Nr、方位向采样点数Na、波束斜视角theta、距离向过采样系数alpha、方位向过采样系数belta、带宽Bw、天线方位向长度La，用于三个相位因子的计算。
三次相位相乘所需的相位因子均由CS算法乘法因子实时生成模块产生，CS算法乘法因子实时生成模块输入端口包括参数输入端口、控制信号输入端口，参数输入端口用于接收雷达参数，控制信号端口用于接收FFT模块输出有效信号标志位，用于控制乘法因子实时生成模块，当控制信号输入端口有效时，将参数输入端口使能接收参数，第一级流水结构进行计算，并且输出使能信号，作为下一级流水结构输入使能信号。最后一级流水线结构输出端口即CS算法乘法因子实时生成模块输出端口，包括因子输出端口以及因子输出有效信号标志位，每次输出四路因子及其有效信号用于后续模块进行处理。
乘法因子实时生成模块参数输入位宽均为32bits，输出因子位宽为32bit，模块中间部分均以64bits位宽进行数据处理。整个***在参数输入、因子输出时数据为单精度浮点数形式，在模块中间处理流程部分，基本为双精度浮点数形式。
整个乘法因子实时生成模块采用分时复用方式，在三次相位相乘时刻对应输出相应因子，其具体分时复用方式如图3所示：
在数据进行第一次方位向FFT后，输出数据有效时，乘法因子实时生成模块输出Chirp Scaling相位因子进行第一次相位相乘；在数据进行第一次距离向FFT后，输出数据有效时，乘法因子实时生成模块输出距离压缩和距离徙动矫正因子进行第二次相位相乘；在数据进行第一次距离向IFFT后，输出数据有效时，乘法因子实时生成模块输出方位压缩和残余相位补偿因子进行第三次相位相乘。
在CS算法硬件实现方案中，数据流均以四路并行数据进行传输处理，因此对应的乘法因子同样以四路因子流形式进行输出，相邻两路之间存在一个固定的差值，利用该差值生成思路因子的模块架构如图4所示：该子模块接受第一路因子和固定差值，将差值与第一路因子相加得下一路因子，同时将第一路因子打拍至同步。按照这个规律重复四次，即可得到四路并行的因子。
乘法因子实时生成模块会按照FFT处理后的数据输出顺序生成并输出对应的乘法因子，其具体顺序主要有以下两种方式：
乘法因子以方位向输出顺序如图5所示，首先输出a1,1、a1,2、a1,3、a1,4；a2,1、a2,2、a2,3、a2,4……a8191,1、a8191,2、a8191,3、a8191,4；a8192,1、a8192,2、a8192,3、a8192,4四路因子，直到当前四列方位向因子输出完毕，再输出下四列方位向因子a1,5、a1,6、a1,7、a1,8；a2,5、a2,6、a2,7、a2,8……a8191,5、a8191,6、a8191,7、a8191,8；a8192,5、a8192,6、a8192,7、a8192,8直至完成一帧数据对应的因子输出为止。
乘法因子以距离向输出顺序如图6所示，首先输出a1,1、a2,1、a3,1、a4,1；a1,2、a2,2、a3,2、a4,2……a1，8191、a2，8191、a3,8191、a4,8191；a1，8192、a2，8192,、a3，8192,、a4，8192,四路因子，直到当前四列距离向因子输出完毕，再输出下四列距离向因子a1,1、a2,1、a3,1、a4,1；a1,2、a2,2、a3,2、a4,2……a1，8191、a2，8191、a3,8191、a4,8191；a1，8192、a2，8192、a3，8192、a4，8192直至完成一帧数据对应的因子输出为止。

Claims (2)

1.雷达成像中CS算法补偿因子实时生成的FPGA实现方法，其特征在于，包括以下步骤：
S1、获取雷达回波数据；
S2、采用CS算法对雷达回波数据进行处理，获得雷达成像图像，CS算法对雷达回波数据的处理过程包括四次FFT和三次相位相乘，第一次FFT为方位向FFT，将数据变换到距离多普勒域，采用Chirp Scaling相位因子进行第一次相位相乘，使所有目标的徙动轨迹一致化；进行第二次FFT，第二次FFT为距离向FFT，将数据变换到二维频域，采用RCMC和距离向压缩因子进行第二次相位相乘，完成距离压缩、SRC和一致RCMC；进行第三次FFT，第三次FFT为距离向IFFT，将数据变回距离多普勒域，采用相位校正和方位压缩因子进行第三次相位相乘；进行第四次FFT，第四次FFT为方位向IFFT，获得输出图像；
所述的Chirp Scaling相位因子、RCMC和距离向压缩因子、相位校正和方位压缩因子为CS算法乘法因子，CS算法乘法因子为通过CS算法补偿因子实时生成模块实时生成，CS算法补偿因子实时生成模块为基于FPGA的方式实现，具体为：采用21级流水线的结构，每一层流水结构利用Floating Point IP核和Cordic IP核对输入的雷达信号参数进行计算，从而获得包括Chirp Scaling相位因子、RCMC和距离向压缩因子、相位校正和方位压缩因子在内的CS算法乘法因子；CS算法补偿因子实时生成模块的实时生成控制方式为，通过控制信号端口接收FFT模块输出的有效信号标志位，在输出端口设置对应的因子输出有效信号标志位，从而在每一次相位相乘时根据上一次FFT的具体操作获得下一次相位相乘需要的对应CS算法乘法因子。
2.根据权利要求1所述的雷达成像中CS算法补偿因子实时生成的FPGA实现方法，其特征在于，CS算法乘法因子包括两种数据输出顺序，分别为：
a.以方位向顺序输出：
a1,1、a1,2、a1,3、a1,4、
a2,1、a2,2、a2,3、a2,4
……
a8191,1、a8191,2、a8191,3、a8191,4
a8192,1、a8192,2、a8192,3、a8192,4
四路因子，直到当前四列方位向因子输出完毕，再输出下四列方位向因子：
a1,5、a1,6、a1,7、a1,8
a2,5、a2,6、a2,7、a2,8
……
a8191,5、a8191,6、a8191,7、a8191,8
a8192,5、a8192,6、a8192,7、a8192,8
b.以距离向顺序输出：
a1,1、a2,1、a3,1、a4,1
a1,2、a2,2、a3,2、a4,2
……
a1，8191、a2，8191、a3,8191、a4,8191
a1，8192、a2，8192、a3，8192、a4，8192
四路因子，直到当前四列距离向因子输出完毕，再输出下四列距离向因子：
a5,1、a6,1、a7,1、a8,1
a5,2、a6,2、a7,2、a8,2
……
a5，8191、a6，8191、a7,8191、a8,8191
a5，8192、a6，8192、a7，8192、a8，8192
直至完成一帧数据对应的因子输出为止。
CN202210274459.5A 2022-03-21 2022-03-21 雷达成像中cs算法补偿因子实时生成的fpga实现方法 Active CN114626006B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210274459.5A CN114626006B (zh) 2022-03-21 2022-03-21 雷达成像中cs算法补偿因子实时生成的fpga实现方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210274459.5A CN114626006B (zh) 2022-03-21 2022-03-21 雷达成像中cs算法补偿因子实时生成的fpga实现方法

Publications (2)

Publication Number Publication Date
CN114626006A CN114626006A (zh) 2022-06-14
CN114626006B true CN114626006B (zh) 2023-03-14

Family

ID=81904900

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210274459.5A Active CN114626006B (zh) 2022-03-21 2022-03-21 雷达成像中cs算法补偿因子实时生成的fpga实现方法

Country Status (1)

Country Link
CN (1) CN114626006B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116185335B (zh) * 2023-04-28 2023-06-23 中国电子科技集团公司第五十四研究所 一种BP成像算法的System Generator实现方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101430380A (zh) * 2008-12-19 2009-05-13 北京航空航天大学 基于非均匀采样的大斜视角机载sar聚束模式成像方法
CN103197317A (zh) * 2013-03-26 2013-07-10 西安电子科技大学 基于fpga的sar成像方法
CN103954963A (zh) * 2014-05-05 2014-07-30 北京理工大学 一种基于cs算法的步进频sar成像方法
CN104215962A (zh) * 2014-08-04 2014-12-17 北京理工大学 基于单片FPGA的Chirp Scaling成像方法
CN107918124A (zh) * 2017-10-26 2018-04-17 西安电子科技大学 带有方位空变校正的机载大斜视高分辨sar成像方法
CN110658520A (zh) * 2019-08-19 2020-01-07 中国科学院电子学研究所 一种合成孔径雷达成像***及方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10042046B2 (en) * 2015-07-07 2018-08-07 Mitsubishi Electric Research Laboratories, Inc. System and method for radar imaging using distributed arrays and compressive sensing

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101430380A (zh) * 2008-12-19 2009-05-13 北京航空航天大学 基于非均匀采样的大斜视角机载sar聚束模式成像方法
CN103197317A (zh) * 2013-03-26 2013-07-10 西安电子科技大学 基于fpga的sar成像方法
CN103954963A (zh) * 2014-05-05 2014-07-30 北京理工大学 一种基于cs算法的步进频sar成像方法
CN104215962A (zh) * 2014-08-04 2014-12-17 北京理工大学 基于单片FPGA的Chirp Scaling成像方法
CN107918124A (zh) * 2017-10-26 2018-04-17 西安电子科技大学 带有方位空变校正的机载大斜视高分辨sar成像方法
CN110658520A (zh) * 2019-08-19 2020-01-07 中国科学院电子学研究所 一种合成孔径雷达成像***及方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Jerome L.V.M.Stanislaus等.Low-complexity FPGA implementation of compressive sensing reconstruction.2013,第671-675页. *
徐涛 等.一种基于OSA与CS原理的SAR成像算法研究.2013,第49卷(第07期),第126-131,170页. *

Also Published As

Publication number Publication date
CN114626006A (zh) 2022-06-14

Similar Documents

Publication Publication Date Title
CN103197317B (zh) 基于fpga的sar成像方法
CN101369018B (zh) 一种星机联合双基地合成孔径雷达频域成像方法
CN110596646B (zh) 一种基于mimo体制的提高雷达角度分辨率的布局及方法
CN107229048B (zh) 一种高分宽幅sar动目标速度估计与成像方法
CN108387877B (zh) 一种多输入多输出雷达的运动目标相位修正方法
US4292634A (en) Real-time multiple-look synthetic aperture radar processor for spacecraft applications
CN108120980A (zh) 一种星载sar多模式成像信号处理算法的fpga的实现方法
CN106896360A (zh) 一种sar信号处理算法的fpga实现方法
CN114626006B (zh) 雷达成像中cs算法补偿因子实时生成的fpga实现方法
CN109633622B (zh) 基于时间反转二阶Keystone变换的雷达机动目标跨单元积累方法
CN114545411B (zh) 一种基于工程实现的极坐标格式多模高分辨sar成像方法
CN109444901A (zh) 一种异构环境下多子阵sas子孔径成像方法
Zheng et al. Fast parameter estimation algorithm for cubic phase signal based on quantifying effects of Doppler frequency shift
Zhang et al. Influence of the stop-and-hop assumption on synthetic aperture sonar imagery
CN110618411A (zh) 一种机载雷达杂波实时信号产生方法
Yu et al. Ground moving target motion parameter estimation using Radon modified Lv's distribution
CN112327301A (zh) 基于子孔径grft的低信噪比下参数化平动补偿快速方法
Yang et al. GF-3 data real-time processing method based on multi-satellite distributed data processing system
Zhou et al. Unambiguous reconstruction for multichannel nonuniform sampling SAR signal based on image fusion
CN104215962A (zh) 基于单片FPGA的Chirp Scaling成像方法
Liu et al. Parameter estimation of moving targets in the SAR system with a low PRF sampling rate
CN115308703A (zh) 一种高频mimo雷达辅助校准源及校准方法
CN113376632B (zh) 一种基于预处理和改进pfa的大斜视机载sar成像方法
Wang et al. Study on motion compensation method for W-band UAV MISAR real-time imaging
CN111220957B (zh) 一种基于fpga的距离徙动补偿***

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