基于FPGA的箔条云雷达回波实时模拟方法技术

本发明专利技术公开了一种基于FPGA的动态箔条云雷达回波模拟方法，主要解决现有技术模拟箔条云雷达回波可移植性差，且只模拟单发箔条云回波的问题。其通过模拟运行在FPGA平台上，并对计算过程进行硬件优化来实现：1)计算箔条的初始姿态、每根箔条的初始直角坐标和每根箔条的雷达散射截面值；2)根据箔条高度计算得到大气密度；3)计算出箔条在一个采样时间点内的位移增量，进而计算出每根箔条的直角坐标；4)根据箔条位置的直角坐标与雷达坐标计算两者距离，根据该距离和箔条雷达散射截面值计算箔条云的雷达信号回波，完成一个时间点的计算；重复2)～4)得到不同时刻箔条云的回波。本发明专利技术可移植性好，可用于模拟多发箔条云的雷达回波。

【技术实现步骤摘要】
基于FPGA的箔条云雷达回波实时模拟方法
本专利技术属于雷达
，尤其涉及一种雷达对箔条云的回波实时模拟方法，可用于获得不同时刻箔条云中各个箔条的动态回波。
技术介绍
箔条云是一种使用最为广泛的无源干扰，由于它具有滞空时间长、作用面积大、回波强的特点，在防空、反舰等领域被广泛应用。箔条云回波计算主要分为三个部分，分别是箔条运动扩散计算、箔条雷达散射截面计算和箔条回波计算。其运动扩散与箔条释放平台的箔条抛洒速度、抛洒方向、箔条所在环境的风力、空气密度和粘度相关，是影响箔条的雷达回波特性的主要因素之一。箔条云内的箔条雷达散射截面统计值受箔条云内箔条总数、入射波频率因素影响。其回波特性与雷达的工作体制以及工作波形有关。在已有的技术中，大多数箔条云回波模拟都是基于PC平台进行的。中国专利CN108020819A，公开了一种弹载雷达跟踪到运动目标后的箔条云干扰回波模拟方法，该方法首先计算每个时刻的箔条云质心坐标，再求出每个距离门上的箔条数；然后得到单根箔条的天线增益、幅度、频率，求得单根箔条的干扰调制值；最后通过叠加所有箔条的回波值得到整个箔条云的回波。此专利实施是基于电脑平台，资源消耗大，系统结构复杂，此外只能对单发箔条云的回波进行仿真，不能对多发箔条的情况进行仿真。中国专利CN110850376A，公开了一种箔条干扰模拟方法及箔条干扰模拟器，该专利首先基于箔条电磁特性进行离线仿真，产生箔条干扰滤波器系数；之后在线产生具有箔条干扰特性的回波系数；最后产生箔条干扰模拟信号。此专利基于箔条的回波系数生成箔条回波，未考虑箔条的实际运动情况，不能准确反映箔条在空中扩散时雷达回波的变化情况；此外，基于电脑平台的仿真资源消耗大、移动性较差。
技术实现思路
本专利技术的目的在于针对上述现有技术的不足，提供一种基于FPGA的箔条云雷达回波实时模拟方法，以在准确反映箔条云在空中运动扩散的条件下，实现对多发箔条的雷达回波进行仿真，并减小硬件资源消耗。为实现上述目的，本专利技术包括如下步骤：(1)计算箔条的初始姿态，包括原点到箔条的连线与正z轴之间的天顶角θ、箔条在xOy平面的投影线与正x轴之间的方位角箔条转动角度ζ、水平扩散速度v，并将这些计算结果保存到RAM中；(2)计算每个箔条的初始直角坐标Lx0(i)、Ly0(i)、Lz0(i)，分别保存到RAM中，其中i为箔条的序号；(3)计算每个箔条的雷达散射截面RCS值，把结果保存到RAM中；(4)计算每条箔条所在高度z(i)处的大气密度ρ(i)；(5)计算每条箔条受力引起的三个方向上的位移增量dx1(i)、dy1(i)、dz1(i)；(6)计算箔条螺旋转动引起的水平位移增量dx2(i)、dy2(i)；(7)根据(5)得到的箔条位移增量和(6)得到的位移增量，计算本时间采样点的箔条直角坐标xt(i)、yt(i)、zt(i)；(8)根据雷达的空间坐标和每个箔条的空间坐标计算雷达到每个箔条的距离Rt(i)；(9)根据步骤(3)得到的每根箔条的雷达散射截面RCS值和步骤(8)得到的雷达到箔条的距离Rt(i)，计算箔条的窄带信号回波S(t,i)，并把结果存储到FPGA外部的DDR存储器；(10)循环步骤(4)到步骤(10)，计算不同时间采样点的箔条回波。本专利技术与现有技术相比，具有如下有优点：1)本专利技术由于基于FPGA平台进行模拟，因而资源消耗小，可移植性好。2)本专利技术由于在计算过程中采用多种硬件计算单元，效率相比于电脑更高。3)本专利技术从箔条的空间扩散模型出发，首先对箔条云在不同时刻的空间分布进行计算，然后根据空间分布计算每个时间采样点的箔条回波，能准确仿真空间中多发箔条云在不同时刻的回波。附图说明图1是本专利技术的实现流程图；图2是本专利技术中第1000个时间采样点时的箔条在空间中的分布；图3是本专利技术中在第100个时间采样点的箔条云回波信号；图4是图3箔条云回波幅度的直方图；图5是现有箔条云雷达回波仿真的箔条云回波幅度直方图。具体实施方式以下结合附图对本专利技术的实施例和效果做进一步详细描述。参照图1，本专利技术的具体实施步骤如下:步骤1.计算箔条的初始姿态、每个箔条的初始直角坐标和每个箔条的雷达散射截面值。这三个计算过程互不依赖、相互独立，在FPGA内采用并行方式同时计算，具体实现如下：1.1)计算箔条的初始姿态：在计算初始姿态参数时，需要生成均匀分布随机数和高斯分布随机数；本实例使用xorshiro128算法生成范围是(0,1)均匀分布随机数，模块名称为rand()，使用Box-Muller算法在FPGA内计算标准正态分布随机数z0～N(0,1)，并根据公式z＝z0σ+μ得到均值为μ、标准差为σ的高斯分布随机数z，模块名称为normrnd(μ,σ)；调用随机数模块，计算各初始姿态的如下参数：原点到第i个箔条的连线与正z轴之间的天顶角：箔条在xOy平面的投影线与正x轴之间的方位角：箔条转动角度：ζ(i)＝2π·rand(i)；水平扩散速度：v(i)＝normrnd(0,5)；把这些计算结果保存到RAM中；1.2)计算每个箔条的初始直角坐标Lx0(i)、Ly0(i)、Lz0(i)：首先，调用rand()在球坐标中计算箔条云中各个箔条的随机初始位置(r(i),s1(i),s2(i))：r(i)＝2rand(i)，s1(i)＝2πrand(i)，s2(i)＝2πrand(i)；然后，根据预先设定的抛洒箔条处的直角坐标(hx,hy,hz)和箔条的随机初始位置(r(i),s1(i),s2(i))，计算每个箔条的初始直角坐标：其中三角函数调用CORDICIP核计算；1.3)计算箔条云中各个箔条的雷达散射截面值：首先，计算单根箔条平均雷达散射截面值：σsingle_mean＝0.153λ2，其中λ为入射波长；其次，计算箔条云平均雷达散射截面值：σcloud_mean＝0.65Nσsingle_mean，其中N为单发箔条云中的箔条总数；然后，调用对数IP核计算箔条云中的各个箔条雷达散射截面值：σ(i)＝-σcloud_meanlog(1-rand(i))，并将计算结果存入RAM中，供后续步骤调用。步骤2.计算每条箔条所在高度zt(i)处的大气密度ρ(i)。2.1)计算每根箔条所在高度zt(i)对应的位势高度H(i)：其中，有效地球半径R′m＝8500km，zt(i)是每根箔条在时间采样点t时的高度；但t＝0时，由于箔条还未下落，所以z0(i)＝Lz0(i)；当t≥1时，由于箔条所在高度的大气密度会影响箔条的受力情况，进而改变箔条的瞬时速度，最终影响箔条在两个时间采样点中的位移增量，所以时间采样点t时的箔条位势高度采用t-1时的箔条高度zt-1(i)计算得出；...

【技术保护点】
1.一种基于FPGA的箔条云雷达回波实时模拟方法，其特征在于，包括如下：/n(1)计算箔条的初始姿态，包括原点到箔条的连线与正z轴之间的天顶角θ、箔条在xOy平面的投影线与正x轴之间的方位角

【技术特征摘要】
1.一种基于FPGA的箔条云雷达回波实时模拟方法，其特征在于，包括如下：
(1)计算箔条的初始姿态，包括原点到箔条的连线与正z轴之间的天顶角θ、箔条在xOy平面的投影线与正x轴之间的方位角箔条转动角度ζ、水平扩散速度v，并将这些计算结果保存到RAM中；
(2)计算每个箔条的初始直角坐标Lx0(i)、Ly0(i)、Lz0(i)，分别保存到RAM中，其中i为箔条的序号；
(3)计算每个箔条的雷达散射截面RCS值，把结果保存到RAM中；
(4)计算每条箔条所在高度z(i)处的大气密度ρ(i)；
(5)计算每条箔条受力引起的三个方向上的位移增量dx1(i)、dy1(i)、dz1(i)；
(6)计算箔条螺旋转动引起的水平位移增量dx2(i)、dy2(i)；
(7)根据(5)得到的箔条位移增量和(6)得到的位移增量，计算本时间采样点的箔条直角坐标xt(i)、yt(i)、zt(i)；
(8)根据雷达的空间坐标和每个箔条的空间坐标计算雷达到每个箔条的距离Rt(i)；
(9)根据步骤(3)得到的每根箔条的雷达散射截面RCS值和步骤(8)得到的雷达到箔条的距离Rt(i)，计算箔条的窄带信号回波S(t,i)，并把结果存储到FPGA外部的DDR存储器；
(10)循环步骤(4)到步骤(10)，计算不同时间采样点的箔条回波。


2.根据权利要求1所述的方法，其特征在于，所述(1)中计算箔条的初始状态，实现如下：
计算原点到第i个箔条的连线与正z轴之间的天顶角：

其中，rand()为生成取值(0,1)范围内的均匀分布随机数；
计算箔条在xOy平面的投影线与正x轴之间的方位角：
计算箔条转动角度：ζ(i)＝2π·rand(i)；
计算水平扩散速度：v(i)＝normrnd(0,5)，其中normrnd(0,5)为生成服从均值为0、标准差为5的正态分布的随机数。


3.根据权利要求1所述的方法，其特征在于，所述(2)中计算每个箔条的直角坐标Lx0(i)、Ly0(i)、Lz0(i)，实现如下：
(2a)在球坐标中计算箔条云中各个箔条的随机初始位置(r(i),s1(i),s2(i))：
r(i)＝2rand(i)，s1(i)＝2πrand(i)，s2(i)＝2πrand(i)；
(2b)根据预先设定的抛洒箔条处的直角坐标(hx,hy,hz)和箔条的随机初始位置(r(i),s1(i),s2(i))
计算每个箔条的初始直角坐标：





4.根据权利要求1所述的方法，其特征在于，所述(3)根据统计方法计算箔条云中各个箔条雷达散射截面值，按照如下步骤进行：
(3a)计算单根箔条的平均雷达散射截面值：
σsingle_mean＝0.1...

【专利技术属性】
技术研发人员：张民，张树豪，李金星，江旺强，魏鹏博，聂丁，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西;61