一种基于导弹运动参数的末制导前视雷达回波模拟方法
【技术领域】
[0001] 本发明属于雷达信号处理技术领域，尤其涉及一种末制导前视雷达回波模拟方 法。
【背景技术】
[0002] 在导弹飞行末端接近攻击目标时，导弹处于俯冲状态，实时获取导弹正前方攻击 目标位置信息和攻击目标附近环境，对导弹精确打击目标起着尤为重要的作用。传统的正 侧视单基SAR和斜前视DBS技术都不能对导弹正前方成像，然而前视雷达方位向回波信号 为目标散射特性和天线方向图的卷积的结果，可以通过对雷达方位向回波信号反解卷积， 恢复目标的方位角度信息，距离上则通过距离脉冲压缩技术来实现目标距离维的高分辨。 因此，利用末制导前视雷达获取的攻击目标方位角度信息和弹目距离，可以实现导弹俯冲 状态下对地面静止或移动目标的精确定位和精确打击。
[0003] 模拟末制导前视雷达回波信号，对评价末制导前视雷达成像算法、研究攻击目标 的散射特性效应、设计雷达参数和导弹运动参数具有重要的现实意义。由于导弹攻击目标 时，处于俯冲姿态，导弹与目标近似在一条直线上，弹载前视雷达的照射区域也位于导弹的 正前方，并且前视雷达的成像方法与传统的SAR和DBS技术成像方法不同，因而导致弹载前 视雷达的回波模型与常见的正侧视单基SAR、双基前视SAR、斜前视DBS等雷达的回波模型 也有着很大的不同。文献"Atime-domainrawsignalsimulatorforinterferometric SAR"（IEEEtransactionsongeoscienceandremotesensing, 2004,V. 42,NO. 9,pp 1811 - 1817)利用时域的方法来模拟生成星载干涉SAR的回波信号；文献"Arawsignal simulatorforbistaticSAR"（ChineseJournalofAeronautics, 2009,V.22,pp 434 - 443)采用了二维频域快速傅里叶变换来生成双基SAR的回波信号；文献"Motion platformforward-lookingreal-beamradarechomodeling"（Processdingsof 2011IEEECIEinternationalconferenceonradar,Vol. 2:1370 -1373)建立了运动平台 前视实波束雷达的回波模型。然而现有SAR回波生成方法或者已有的前视雷达回波模型， 不能直接应用于模拟末制导前视雷达的回波生成。
【发明内容】
[0004] 本发明是针对已有的回波信号生成方法难以直接应用于末制导前视雷达的回波 信号生成，提出一种基于导弹运动参数的末制导前视雷达回波模拟生成方法，利用导弹的 运动轨迹参数，模拟生成地面散射点的回波距离历史，通过先求解雷达照射地面区域中所 有散射点的雷达回波矩阵，再将所有散射点的雷达回波矩阵相加，获得整个雷达照射场景 的的回波矩阵，模拟了末制导前视雷达的回波生成。
[0005] 本发明的技术方案为：一种基于导弹运动参数的末制导前视雷达回波模拟方法， 具体包括以下步骤：
[0006] S1 :根据弹载雷达参数，将导弹运动参数的采样周期插值为脉冲重复周期，得到插 值后的导弹三维位置坐标向量；
[0007] S2:在弹载雷达照射的地面区域内设置地面散射点，根据散射点的位置信息和导 弹的运动参数，计算出散射点的斜距历史；
[0008] S3:由散射点的斜距历史、雷达脉冲时宽，计算出散射点回波距离向时域滤波器和 散射点的相位信息；
[0009] S4 :由散射点的方位角度信息和雷达天线方向图，构造出散射点回波方位向时域 滤波器；
[0010] S5:将散射点的散射系数乘以散射点相位矩阵、散射点回波距离向时域滤波器和 散射点回波方位向时域滤波器，得到该散射点的回波矩阵。
[0011] 进一步地，所述步骤S2包括以下分步骤：
[0012] S21 :设置弹载雷达照射场景的散射点矩阵，根据设置的散射点矩阵和弹载雷达参 数，求解出散射点矩阵对应的距离向时间矩阵和方位向时间矩阵；
[0013] S22 :从步骤S1中的导弹三维位置坐标向量中，截取与步骤S21中弹载雷达扫描 完弹载雷达照射的地面区域对应的导弹三维位置坐标向量，从散射点矩阵中取出一个散射 点，再利用散射点的位置信息和截取后的导弹三维坐标向量，求解出散射点的斜距历史矩 阵。
[0014] 进一步地，所述步骤S3包括以下分步骤：
[0015]S31:根据S22中的散射点的斜距历史矩阵和雷达发射的脉冲信号，得到雷达照射 到散射点后返回的延时矩阵和相位矩阵；
[0016] S32 :利用步骤S21中距离向时间矩阵和步骤S31中的延时矩阵，求解出散射点回 波距离维滤波器矩阵。
[0017] 更进一步地，所述距离维滤波器矩阵匕(队，Na)中任一元素值FJ队，NJ 的求 解方法如下：
[0018]
[0019] 其中，队为雷达回波距离向采样点，Na为雷达回波方位向采样点，Rt(N"Na)为距离 时间矩阵，Td(H)为延时矩阵，?；为雷达脉冲时宽，row为矩阵的行索引号，col为矩阵的 列索引号。
[0020] 进一步地，所述步骤S5包括以下分步骤：
[0021 ] S51 :将散射点的散射系数、步骤S31中散射点的相位矩阵、步骤S32中的散射点距 离维滤波器矩阵和步骤S4中的散射点方位向滤波器矩阵相乘，得到散射点最后的回波矩 阵；
[0022] S52 :重复步骤S22至步骤S51，得到雷达照射的地面区域中所有散射点的回波矩 阵，再将所有散射点的回波矩阵相加，得到雷达照射的地面区域场景的回波矩阵。
[0023] 本发明的有益效果：本发明提供一种基于导弹运动参数的末制导前视雷达回波模 拟方法，首先根据弹载雷达参数，将导弹运动参数的采样周期插值为脉冲重复周期；然后在 弹载雷达照射的地面区域内设置地面散射点，根据散射点的位置信息和导弹的运动参数， 得散射点的斜距历史；再由散射点的斜距历史、雷达脉冲时宽，得到散射点回波距离向时域 滤波器和散射点的相位信息，由散射点的方位角度信息和雷达天线方向图，构造散射点回 波方位向时域滤波器；最后根据散射点的散射系数、散射点相位信息、散射点回波距离向时 域滤波器和散射点回波方位向时域滤波器，得该散射点的回波信息，利用时域相干叠加方 法获得整个雷达照射场景的回波信息。本发明利用导弹的运动参数，模拟了弹载前视雷达 回波信号，为导弹运动参数设计和弹载雷达参数设计提供了依据。
【附图说明】
[0024] 图1是本发明末制导前视雷达成像几何模型。
[0025] 图2是本发明提供的方法的流程图。
[0026] 图3是本发明具体实施例中导弹运动轨迹三维坐标。
[0027] 图4是本发明具体实施例中9个散射点的地面分布。
[0028] 图5是本发明具体实施例中雷达天线方向图。
[0029] 图6是本发明具体实施例中9个散射点的二维时域回波。
[0030] 图7是本发明具体实施例中9个散射点经脉冲压缩处理后的结果。
[0031] 图8是本发明具体实施例中9个散射点经距离走动校正处理后的结果。
[0032] 图9是本发明具体实施例中9个散射点经反卷积处理后的结果。
【具体实施方式】
[0033] 下面结合附图和具体实施例对本发明所提出的方法做进一步的描述。
[0034] 本发明具体实例中采用的末制导前视雷达成像几何模型如图1所示，攻击目标为 三维坐标的坐标原点，x〇y平面位于水平地面，Z轴垂直于水平地面；具体实例的雷达参数 如表1所示。
[0035]
[0036] 如图2所示为本发明方案的流程图，具体包括以下步骤：
[0037] S1 :根据弹载雷达参数，将导弹运动参数的采样周期插值为脉冲重复周期，得到插 值后的导弹三维位置坐标向量；
[0038] 读取导弹运动参数，利用三次样条插值法，将导弹运动参数的采样周期变换为前 视雷达的脉冲重复周期prt，插值后的导弹三维位置坐标向量记为，导弹的三维 位置运动轨迹如图3所示。
[0039] S2:在弹载雷达照射的地面区域内设置地面散射点，根据散射点的位置信息和导 弹的运动参数，计算出散射点的斜距历史；
[0040] 所述步骤S2包括以下分步骤：
[0041] S21 :设置雷达照射场景的散射点矩阵，根据设置的散射点矩阵和弹载雷达参数， 求解出散射点矩阵对应的距离向时间矩阵和方位向时间矩阵；
[0042] 按照如图4所示，在xoy平面上设置9个散射点，9个散射点的散射系数都为1，记 为〇 ，地面场景方位向总点数为凡点，且Nx= 3,地面场景距离向总点数为Ny 点，且Ny= 3,9个散射点的方位角度和坐标在图5中已标出，方位向角度向量记为， 对应的方位时间向量为忍（/乂.；)，散射系数矩阵。（化乂'^^中第"列向量的方位角度大 小为中的第Xn个元素方位时间为的第Xn个元素x为 地面方位向坐标，y为地面距离坐标，地面场景对应的斜距范围为R_= 4500米，R_ = 5500米。通过散射系数矩阵σ(Nx，Ny)(x^Q)和雷达参数，求解出距离向时间矩阵Rt (队，Na) 和方位向时间矩阵At(H)，其中队为雷达回波距离向采样点，Na为雷达回波方位向采 样点，距离向时间矩阵Rt (Nr，NJ的每一列向量都为^ (X)，^ (%)的第m个元素为 ^(乂),"，方位向时间矩阵At(Nr，Na)的每一行向量都为1.(%)的第η个元素 为焉)《、4,(^?),、队和^求解表达式如下：
[0046] Nr=ceil((2X(Rnax-Rnin)/Vc+2XTr)/