在FPGA中SAR回波模拟中回波信号多普勒相位的两种实现方法
１ 引言
合成孔径雷达（synthetic aperture radar ，简称sar）作为一种高分辨微波成像雷达［１］，是地面信息获取的重要手段，它在国土测量、军事等领域发挥着重要作用。为确保所设计的sar系统能够满足各个用户的具体需求，原始数据仿真和成像处理以及图像指标的评估已经成为sar系统研制前的必需步骤。除此之外，在sar实时成像系统研制以及地面处理的过程中，我们同样需要大量的模拟回波数据［２］。因此，sar回波模拟技术能够为我们提供基本的研究手段和研究工具，它在sar的发展和应用中，具有及其重要的作用。
作为雷达与目标相对运动的直观反映，多普勒相位是sar回波模拟的原理依据，同样也是回波模拟精度的重要影响因素。在多普勒相位计算过程中存在数值开方运算，这使得其在fpga（现场可编程门阵列）实现中存在两种主要限制因素：其一，fpga硬件中不存在专门的开方器件；导致了开方运算在fpga实现中的复杂性。其二，为了保证原始数据的仿真精度及适应大斜距数据的仿真，导致sar回波模拟对fpga定点运算中的斜距动态范围要求大［３］。
因此，本文针对sar回波模拟中回波信号多普勒相位在fpga实现中存在的限制，采用泰勒级数展开和cordic算法两种常用方法对其进行实现与仿真。
２ sar回波模拟信号中的多普勒相位
sar的空间几何关系如图１，其中sar飞行路径的地面航迹方向称为方位方向，与其垂直的方向称为距离方向［４］。
图１ sar的空间几何关系
sar在运动过程中，以特定的脉冲重复周期
（pulse repetition time，简称prt）发射并接收脉冲串，雷达天线波束照射到地面以后，照射区域内的各个点目标（又被称作散射元）对入射波进行后向散射。发射信号经过目标和天线方向图的调制，成为携带目标信息和环境信息的sar回波。
sar的发射脉冲串一般为线性调频（chirp）信号，设定sar的发射脉冲串：
在本文中，要求斜距的近似误差最大为1/8波长，由表格１可见，目标点到sar平台的斜距量级为106，动态范围很大；由于fpga内部没有专用开方器件，使得多普勒相位运算过程复杂、速度低，在此使用泰勒级数展开与cordic算法。
３ 基于泰勒级数展开的多普勒相位计算方法
泰勒级数是幂级数的一种，如果有了某一函数的幂级数展开式，则我们就可用它进行近似计算，即在展开式有效的区间范围上，我们可按照已知精确度要求，利用这个幂级数展开式将该函数值近似地计算出来。
４ 基于cordic算法的多普勒相位计算方法
由于xilinx的cordic ip核对输入输出位宽的限制要求，对于斜距计算的通用方法来说，使其变得不可取。因此我们采用结合了ip核思想的cordic算法，cordic算法是一种数值型逼真计算方法，其fpga实现模块主要由各级寄存器、移位器、符号标志寄存器以及加法器组成，它在硬件电路实现上只用到了位移操作和加／减操作，这大大节约了fpga资源。
图２ 泰勒级数展开设计框图
假设sar平台与目标点的相对位置坐标为（ｘ，ｙ，z），经过两次cordic 算法运算即可得出目标点到sar平台的斜距值：
为了更方便地说明cordic算法，我们在此只针对sqrt1的求解过程进行详细的展开介绍。本文中，cordic ip核采用向量化模式，通过将任意方向的输入向量旋转一系列的预定角度最终达到与ｘ轴的对齐。故而此算法的最终结果即所有旋转角度的累加值，以及输入向量定标后的幅值（结果在ｘ分量中），同时ｙ分量的符号决定着下一次旋转的方向。在此模式中角度累加器初值为零，整个迭代运算结束后，结果为最终旋转角度［８］。于是以sqrt1为例，cordic算法中存在差分方程：
ｎ次旋转后最终达到与ｘ轴的对齐，此时
图３ cordic算法设计框图
由上图，输入信号x0，y0通过一系列移位操作与加减法操作最终得到xn，即sqrt1，同样可求出sqrt2即斜距值，最后斜距乘以波长系数４π／λ即可得到多普勒相位。
５ 仿真验证与性能分析
本文以xilinx virtex6sx315t为硬件平台，对两种多普勒相位计算方法进行了仿真验证与性能分析。
泰勒级数展开与cordic算法的斜距误差仿真结果如下图４：
图４ 斜距误差仿真结果图（左为泰勒级数展开法，右为cordic算法）
图５ 多普勒相位对比图（左为泰勒级数展开法，右为cordic算法）
图６ 脉冲压缩仿真图（左为泰勒级数展开法，右为cordic算法）
两种方法的性能分析包括资源消耗，fpga时序约束，以及斜距误差对脉压精度的影响。资源消耗如下：
由上表知，较泰勒级数展开法，cordic算法能够节省更多的资源，甚至在乘法器资源上是泰勒级数展开法的四分之一。
在fpga时序约束上，cordic算法最高能达到110mhz，泰勒级数展开法略高些：150mhz。
斜距误差对脉压精度的影响如下：
显然，在脉压精度上，两种方法均能满足精度要求，且泰勒级数展开法分辨率略高些。
６ 结论
本文针对sar回波模拟中回波信号多普勒相位在fpga实现中存在的主要限制因素：fpga硬件中不存在专门的开方器件；sar回波模拟对fpga定点运算中的斜距动态范围要求大。采用了两种常用的近似计算方法：泰勒级数展开以及cordic算法。并对其进行fpga仿真与验证，两种方法不仅很好的解决了斜距开方运算，仿真结果同时证明了其精度有效性，并在fpga资源消耗，时序约束，以及脉压精度上作出了具体评估与比较，在sar回波模拟中具有良好的应用前景。
如何读懂MOSFET的规格书datasheet
卫星机顶盒能不能装软件
电子市场增长稳定，下半年整体行业景气
如何让电动机启动控制回路_四种常见电动机控制回路详解
AI+医药研发的未来有多少想象空间
在FPGA中SAR回波模拟中回波信号多普勒相位的两种实现方法
思杰马克丁宣布大动作，代理商招募计划引爆行业热潮
双通道、快速、降压型控制器可实现0.4V至5.5V的动态电压调节
对于段码液晶屏的驱动IC我们该如何选择
MODBUS转PROFINET网关TS-180连接西门子PLC和工业称重仪表
2015年中国视频监控市场发展特点及未来展望
无人机图传在高速公路的应用场景概述
韩国4G网速变慢引发不满
用简单的退火方法来再生电极材料
功放机电路图系列一（六款模拟电路设计原理图详解）
行程开关工作原理是什么
TI德州仪器差分放大器的简单介绍
电力微气象监测预警系统
Apple Watch遭禁售 原因是侵犯专利
6a和5a快充的区别和优缺点
合成孔径雷达（synthetic aperture radar ，简称sar）作为一种高分辨微波成像雷达［１］，是地面信息获取的重要手段，它在国土测量、军事等领域发挥着重要作用。为确保所设计的sar系统能够满足各个用户的具体需求，原始数据仿真和成像处理以及图像指标的评估已经成为sar系统研制前的必需步骤。除此之外，在sar实时成像系统研制以及地面处理的过程中，我们同样需要大量的模拟回波数据［２］。因此，sar回波模拟技术能够为我们提供基本的研究手段和研究工具，它在sar的发展和应用中，具有及其重要的作用。
作为雷达与目标相对运动的直观反映，多普勒相位是sar回波模拟的原理依据，同样也是回波模拟精度的重要影响因素。在多普勒相位计算过程中存在数值开方运算，这使得其在fpga（现场可编程门阵列）实现中存在两种主要限制因素：其一，fpga硬件中不存在专门的开方器件；导致了开方运算在fpga实现中的复杂性。其二，为了保证原始数据的仿真精度及适应大斜距数据的仿真，导致sar回波模拟对fpga定点运算中的斜距动态范围要求大［３］。
因此，本文针对sar回波模拟中回波信号多普勒相位在fpga实现中存在的限制，采用泰勒级数展开和cordic算法两种常用方法对其进行实现与仿真。
２ sar回波模拟信号中的多普勒相位
sar的空间几何关系如图１，其中sar飞行路径的地面航迹方向称为方位方向，与其垂直的方向称为距离方向［４］。
图１ sar的空间几何关系
sar在运动过程中，以特定的脉冲重复周期
（pulse repetition time，简称prt）发射并接收脉冲串，雷达天线波束照射到地面以后，照射区域内的各个点目标（又被称作散射元）对入射波进行后向散射。发射信号经过目标和天线方向图的调制，成为携带目标信息和环境信息的sar回波。
sar的发射脉冲串一般为线性调频（chirp）信号，设定sar的发射脉冲串：
在本文中，要求斜距的近似误差最大为1/8波长，由表格１可见，目标点到sar平台的斜距量级为106，动态范围很大；由于fpga内部没有专用开方器件，使得多普勒相位运算过程复杂、速度低，在此使用泰勒级数展开与cordic算法。
３ 基于泰勒级数展开的多普勒相位计算方法
泰勒级数是幂级数的一种，如果有了某一函数的幂级数展开式，则我们就可用它进行近似计算，即在展开式有效的区间范围上，我们可按照已知精确度要求，利用这个幂级数展开式将该函数值近似地计算出来。
４ 基于cordic算法的多普勒相位计算方法
由于xilinx的cordic ip核对输入输出位宽的限制要求，对于斜距计算的通用方法来说，使其变得不可取。因此我们采用结合了ip核思想的cordic算法，cordic算法是一种数值型逼真计算方法，其fpga实现模块主要由各级寄存器、移位器、符号标志寄存器以及加法器组成，它在硬件电路实现上只用到了位移操作和加／减操作，这大大节约了fpga资源。
图２ 泰勒级数展开设计框图
假设sar平台与目标点的相对位置坐标为（ｘ，ｙ，z），经过两次cordic 算法运算即可得出目标点到sar平台的斜距值：
为了更方便地说明cordic算法，我们在此只针对sqrt1的求解过程进行详细的展开介绍。本文中，cordic ip核采用向量化模式，通过将任意方向的输入向量旋转一系列的预定角度最终达到与ｘ轴的对齐。故而此算法的最终结果即所有旋转角度的累加值，以及输入向量定标后的幅值（结果在ｘ分量中），同时ｙ分量的符号决定着下一次旋转的方向。在此模式中角度累加器初值为零，整个迭代运算结束后，结果为最终旋转角度［８］。于是以sqrt1为例，cordic算法中存在差分方程：
ｎ次旋转后最终达到与ｘ轴的对齐，此时
图３ cordic算法设计框图
由上图，输入信号x0，y0通过一系列移位操作与加减法操作最终得到xn，即sqrt1，同样可求出sqrt2即斜距值，最后斜距乘以波长系数４π／λ即可得到多普勒相位。
５ 仿真验证与性能分析
本文以xilinx virtex6sx315t为硬件平台，对两种多普勒相位计算方法进行了仿真验证与性能分析。
泰勒级数展开与cordic算法的斜距误差仿真结果如下图４：
图４ 斜距误差仿真结果图（左为泰勒级数展开法，右为cordic算法）
图５ 多普勒相位对比图（左为泰勒级数展开法，右为cordic算法）
图６ 脉冲压缩仿真图（左为泰勒级数展开法，右为cordic算法）
两种方法的性能分析包括资源消耗，fpga时序约束，以及斜距误差对脉压精度的影响。资源消耗如下：
由上表知，较泰勒级数展开法，cordic算法能够节省更多的资源，甚至在乘法器资源上是泰勒级数展开法的四分之一。
在fpga时序约束上，cordic算法最高能达到110mhz，泰勒级数展开法略高些：150mhz。
斜距误差对脉压精度的影响如下：
显然，在脉压精度上，两种方法均能满足精度要求，且泰勒级数展开法分辨率略高些。
６ 结论
本文针对sar回波模拟中回波信号多普勒相位在fpga实现中存在的主要限制因素：fpga硬件中不存在专门的开方器件；sar回波模拟对fpga定点运算中的斜距动态范围要求大。采用了两种常用的近似计算方法：泰勒级数展开以及cordic算法。并对其进行fpga仿真与验证，两种方法不仅很好的解决了斜距开方运算，仿真结果同时证明了其精度有效性，并在fpga资源消耗，时序约束，以及脉压精度上作出了具体评估与比较，在sar回波模拟中具有良好的应用前景。
如何读懂MOSFET的规格书datasheet
卫星机顶盒能不能装软件
电子市场增长稳定，下半年整体行业景气
如何让电动机启动控制回路_四种常见电动机控制回路详解
AI+医药研发的未来有多少想象空间
在FPGA中SAR回波模拟中回波信号多普勒相位的两种实现方法
思杰马克丁宣布大动作，代理商招募计划引爆行业热潮
双通道、快速、降压型控制器可实现0.4V至5.5V的动态电压调节
对于段码液晶屏的驱动IC我们该如何选择
MODBUS转PROFINET网关TS-180连接西门子PLC和工业称重仪表
2015年中国视频监控市场发展特点及未来展望
无人机图传在高速公路的应用场景概述
韩国4G网速变慢引发不满
用简单的退火方法来再生电极材料
功放机电路图系列一（六款模拟电路设计原理图详解）
行程开关工作原理是什么
TI德州仪器差分放大器的简单介绍
电力微气象监测预警系统
Apple Watch遭禁售 原因是侵犯专利
6a和5a快充的区别和优缺点