调频连续波激光引信精确定距幅度平方相干算法_刘海明.pdf

第 卷 第期 年月探 测 与 控 制 学 报 收稿日期：作者简介：刘海明（），男，甘肃武威人，硕士研究生。调频连续波激光引信精确定距幅度平方相干算法刘海明，贺伟（西安邮电大学通信与信息工程学院，陕西 西安 ）摘要：在云雾环境中，粒子会对激光进行散射和吸收，使引信回波信号受到很强的干扰，目标信号难以提取和识别，这给激光引信在设定距离处准确输出起爆信号造成了一定的困难。针对这一问题基于蒙特卡洛仿真方法建立了云雾环境下调频连续波激光引信的多次散射模型。该模型模拟了光子的多径传播和多次散射；同时，利用幅度平方相干（）算法实现调频连续波（）激光引信的精确定距，并和快速傅里叶变换（）算法处理结果进行对比。仿真结果表明，算法在云雾环境中有良好的抗干扰能力，识别目标的正确率要高于 算法。该算法为激光引信抗云雾干扰提供了一种新方法。关键词：激光引信；多次散射；云雾干扰；算法中图分类号：文献标志码：文章编号：（），（，）：，（），（）（）：；引言由于激光具有高亮度、单色性、发散角小、相干性好等特性，可以建立具有良好时间、空间分辨率和高测距精度的探测系统。激光引信具有波束窄、旁瓣小、抗干扰能力强、精准起爆等优点。现阶段脉冲激光引信以其抗电磁干扰能力强、结构简单、成本低等优点被广泛应用。然而，脉冲激光引信容易受到云雾和烟尘等悬浮粒子物的干扰，难以准确识别目标信号。与脉冲激光相比，调频连续波（，）激光引信是依赖于发射信号和接收信号之间的差频信号，通过频域信号处理方法获取目标距离信息，受回波信号的时域特征影响较小，所以它在相同干扰条件下目标识别能力强于脉冲激光引信。因此，研究在云雾环境中 激光引信对提升引信抗干扰能力具有十分重要的意义。国内外学者针对云雾对激光引信干扰的问题已开展了大量的研究。文献 研究了气溶胶粒径分布对不同波长激光回波信号的影响。文献 研究了不同雾霾浓度下散射介质的消光系数对激光后向散射的影响。文献 基于单气溶胶粒子后向散射模型，对云雾环境下的目标回波信号展宽现象进行了研究。文献 通过仿真相干激光引信在多种云雾干扰场景下的回波信号，分析了回波信号的时域和频域特性。文献 提出了 激光引信的原型。文献 研究了激光引信的烟雾后向散射干扰信号特性并提出了基于频谱幅值归一化阈值的检测方法。文献 研究了激光引信的云雾回波的频域特性，并用样机进行了初步试验研究。文献 基于 算法设计了四象限激光引信，有效地提升了引信的抗干扰能力，然而其复杂和昂贵的缺点限制了它的广泛应用。本文提出一种新的基于固定距离的幅度平方相干（，）算法，在不增加额外探测器的情况下，依然具有良好的探测精度。云雾中激光散射的原理大气环境中，在相对湿度和风速等自然条件下，空气中的水蒸气会凝结为水滴悬浮在空中，使得能见度下降。通常云雾的粒径在 之间，粒径分布（）由对数正态分布描述，定义为（）（），（）式（）中，是粒子半径，为粒径分布的标准偏差，为云雾粒子的平均半径。激光在云雾中传播发生光散射，可以看作是激光束与云雾粒子之间发生碰撞并改变其传播方向。光散射模型主要由米氏散射和瑞利散射来描述。根据粒子大小参数：，（）式（）中，为激光波长；当时，光散射主要是米氏散射；当时，光散射主要是瑞利散射。激光波长与云雾粒径接近，所以这种情况下选用米氏散射模型更合适，瑞利散射可以忽略不计。当激光束照射在单个粒子上时，如式（）所示，激光能量沿散射角的分布情况由米氏相函数描述。（）（）（）（）（），（）式（）中，（）和（）分别为平行散射振幅函数和垂直散射振幅函数，和为米氏散射系数。式（）中的相位函数（）描述了散射角中的光的能量。当粒子复折射率 ，粒径参数 时，散射角和激光能量分布（）之间的关系如图所示。图复折射率 和粒子尺寸参数 时的米氏相函数 光子在云雾中传播，随着散射次数的增加，光子的能量会逐次减少，当光子能量减少到一定程度，可视为光子被云雾粒子吸收。图显示了粒子尺寸参数分别为，和 的情况下，光子后向散射能量强度与散射次数之间的关系。当 时，第次的能量强度占总能量强度的，随着散射次数的增加光子后向散射能量强度逐渐降低。图光子散射次数与后向散射强度的关系图 激光引信探测模型激光引信测距的原理是，引信通过发射探 测 与 控 制 学 报镜头发射具有一定线性调制的激光信号，并将一部分激光信号作为本地振荡信号，再对从激光接收镜头接收到的回波信号与本地振荡器信号进行混频，经过低通滤波器，获得频率单一的差频信号。差频信号的频率与目标距离呈线性关系，目标的距离信息可以从差频信号的频率信息中得到。系统结构图如图所示。图系统结构图 常用的线性频率调制方式有锯齿波调制和三角波调制。锯齿波调制实现目标的测距和测速，需要对差频信号进行二维 运算，计算量较大。但是三角波调制相比锯齿波调制在一个周期内做一维 运算就可以确定目标距离和速度，处理简单，易于实现，所以被广泛应用。对称三角波的调制模型如图所示。图对称三角波调制模型 对称三角波调制上下扫频段差频信号的表达式分别为（）（）（）（），（）（）（）（）（），（）式中，（）为上扫频段的差频信号，（）为下扫频段的差频信号；为某一时刻；为信号幅度；为线性调频的起始频率；为线性调频的最高频率；（）为延迟时间；为调制斜率，为调频带宽，为调频周期。当目标静止时，回波信号可以被认为是发射信号的延时副本。为了简化模型，只进行上扫频周期分析，上扫频周期为，差频信号满足 ，因此目标距离的表达式（）可写为（），（）式（）中，为引信与目标之间的距离，为光速。激光在云雾中传输的蒙特卡洛模型现如今基于 散射理论对于单个粒子的激光散射已有大量的研究，然而激光束在云雾中传输过程中，往往是经过多个粒子散射后到达引信的接收镜头，所以必须考虑多次散射。本文基于米氏散射理论，建立了激光在云雾中传输的蒙特卡洛模型。蒙特卡洛方法是将激光束在云雾的传播过程转化为大量独立的光子与云雾粒子发生碰撞的概率模型。激光在云雾中多重散射过程示意图如图所示。图激光在云雾中多重散射示意图 激光发射模型 激光引信利用线性调频电信号来调制传输激光信号的光强，瞬时发射光强度可以用光子数表示。激光束通过引信发射镜头发射到探测视角，光子的初始发射时刻的位置遵循高斯分布，由式（）表示。，（）刘海明等：调频连续波激光引信精确定距幅度平方相干算法式（）中，为激光束腰半径；，是标准正态分布的随机数。光子在发射瞬间的方向 可表示为，（）式（）中，为光子初始发射方向的天顶角，是光子初始发射方向的方位角。激光传输模型光子在云雾中传输时，位置（，）由上一位 置（，），方 向（，）和散射自由程决定，根据比尔 朗伯定律 计算。所以光子空间位置坐标可以表示为，。（）当光子与云雾粒子发生碰撞后，方向发生变化，如图所示，（，）为碰撞后新的方向向量，（，）为碰撞前光子的方向向量。图光子散射后运动方向示意图 光子与云雾粒子碰撞散射后的移动方向由散射角和天顶角决定。散射角根据式（）相位函数（）得到，天顶角是，之间随机角度。那么光子碰撞后的新方向 可以表示为，（），（），。（）当，时，式（）可以简化为，。（）当光子在传播过程中与目标碰撞，会被目标反射改变方向和能量，则光子经目标反射后的方向为，（），（）式（）中，为目标反射天顶角；为目标射方位角，假设目标时表面粗糙的朗伯体反射体，为，区间上均匀分布的随机数。云雾衰减模型云雾对激光衰减特性由消光系数和散射系数来描述。光子与云雾粒子或者目标碰撞后，沿着散射角继续传播。衰减函数与云雾消光系数和散射系数的关系为。（）建立权重函数来表示光子每次与云雾粒子碰撞后的状态，每个光子的初始权重为，随着碰撞次数增加，光子的权重会越来越小。权重函数为 ，（）式（）中，为预先设定的小值。当 时表示光子继续传播；当 时 表示光子被云雾粒子吸收。激光引信精确定距计算 幅度平方相干算法 算法可以在频域中比较两个信号的相似性。因此，当使用 算法进行固定距离时，计算机根据预期距离计算理想条件下的差频信号频率，得探 测 与 控 制 学 报到信号波形，即参考差频信号的波形。理想情况下，差频信号是单频余弦信号，其频率与引信和目标之间的距离成线性关系。理想条件下的差频信号称为参考差频信号。根据 激光探测原理，参考差频信号（）是具有一定频率 的余弦信号。可通过式（）计算：，（）式（）中，为预先设定的距离。本文利用 算法比较参考差频信号和测量差频信号在频域上的相似性，作为引信固定距离的判断依据。如果频谱分布中的 处出现频谱峰值，则意味着两个信号都包含该频率，且目标的实际距离与预先设定的距离相等。在激光引信工程的实施中，测量差频信号（）和参考差频信号（）是离散的数字信号，可以表示为（）和（），测量差频信号（）和参考差频信号（）信号的功率谱密度分别为，（），（），（），（）式中，和，分别为测量差频信号（）和参考差频信号（）的自相关函数。同时，测量 差 频信 号（）和 参 考 差 频 信 号（）的交叉功率谱密度 为，（），（）式（）中，为测量差频信号（）和参考差频信号（）的互相关函数。因此它们之间的幅度平方相干性表达式 为，（），（），（），（），。（）在没有 云 雾 干 扰 的 情 况 下，测 量 差 频 信 号（）中只有目标回波信号（）和随机噪声。此时，如果目标与引信之间的距离等于预先设定的距离，则参考差频信号的频率与目标信号频率一致，（）和（）只有振幅和相位差。相应地，在幅度平方相干的频率 处出现一个明显的峰值。此时，可以确定与目标的当前距离，以满足固定距离条件。在云雾干扰条件下，目标信息淹没在云雾回波和噪声中难以提取。通过 算法，只保留两个信号中的相关部分，这可以大大降低提取目标信息的难度。幅度平方相干算法实现 激光引信固定距离目标检测基于 的固定距离算法的基本步骤如下：）获取当前测量的差频信号（）；）根据预先设定的距离计算频率，并且生成参考差频信号（）；）通过 算法计算（）和（）的幅度平方相干函数（）；）在 结果中，如果在参考差频频率 处出现大于阈值的峰值，证明目标被检测到；如果小于阈值的峰值，则返回步骤）重新开始检测。仿真与实验基于 激光引信在云雾中传输的蒙特卡洛模型进行回波信号仿真，仿真流程图如图所示。图仿真流程图 刘海明等：调频连续波激光引信精确定距幅度平方相干算法基于 算 法 对“目 标 无 云 雾 干 扰 条 件”、“目标有云雾干扰条件”和“无目标有云雾干扰条件”三种场景进行了仿真。将 算法的处理结果与 算法的处理结果进行比较，验证了 算法在烟雾干扰条件下的优势。云雾对激光的干扰程度采用空气能见度指标来描述，能见度、分 别 对 应 轻 度、中 度、重 度 浓 度 干 扰条件。仿真中的引信参数如表所示。目标无云雾干扰条件仿真图显示了无云雾干扰条件下固定距离为的仿真结果。图中实线表示测量差频信号，虚线表示参考差频信号。从图可以看出，在距离处，算法和 算法的处理结果在参考差频处都有明显的峰值，满足引信起爆条件。在、距离处，两种方法的处理结果在参考差频附近没有峰值，不满足引信起爆条件，满足预先设定的在固定距离处目标识别。表 激光引信仿真参数 参数数值激光波长 扫频带宽 扫频周期 初始频率 参考距离 参考差频频率 图当固定距离为 且无云雾干扰时的仿真波形与频谱 目标有云雾干扰条件仿真图图 为固定距离为时，低浓度（能见度）、中浓度（能见度）和高浓度（能见度）云雾干扰条件下的仿真结果。图中实线表示测量差频信号，虚线表示参考差频信号。从图中可以看出，在低浓度干扰条件下，算法和 算法的峰值集中在与目标距离相关的频率处，并未出现干扰峰值；在中浓度干扰条件下，在与目标距离相关的频率处存在峰值，在其他频率点处也出现峰值，但峰值水平相对较小；在高浓度干扰条件下，既存在目标峰值，也存在与目标峰值水平相当的干扰峰值。随着云雾浓度的增大，云雾对激光的后向散射强度也越大。算法频谱中干扰峰值水平变高，在高浓度云雾条件下存在与目标峰值强度相近的干扰峰值，会造成引信误判，无法准确识别目标。当目标距离等于固定距离时，算法频谱中探 测 与 控 制 学 报参考差频 处都有明显的峰值，满足引信起爆条件。当目标距离不等于固定距离时，干扰峰值水平增高幅度较小，峰值水平依然小于 ，不满足起爆阈值，与预期的效果一致。图当固定距离为且低浓度云雾干扰时的仿真波形与频谱 图 当固定距离为且中浓度云雾干扰时的仿真波形与频谱 刘海明等：调频连续波激光引信精确定距幅度平方相干算法图 当固定距离为且高浓度云雾干扰时的仿真波形与频谱 无目标有云雾干扰条件仿真图 为固定距离为时，在中浓度（能见度）云雾干扰条件下无目标的仿真结果。从图 可以看出，在 处云雾距离与固定 距 离 相 等，算法结果中出现幅值较小的峰值，是由于云雾的后向散射引起的，峰值水平不满足引信的起爆阈值；在、处，算法结果中无明显峰值，和预期结果一致。图 中浓度云雾干扰条件下无目标仿真 综上所述，算法可以有效地抑制烟雾干扰。通过参考信号 和测量信号之间的相似性比较，保留了具有 频率的目标信号，与 算法相比有较好的目标识别精度。为了验证所提出的 算法相较于传统 算法具有良好的抗干扰性能，进行了对比实验，对比两种算法在不同云雾干扰下的识别正确率。通过不同云雾干扰下分别进行 次目标识别实验，结果如表所示。由表可看出，在不同云雾干扰下，提出的 算法具有较高的目标识别正确率。表不同云雾干扰下两种算法的识别正确率 云雾干扰实验次数算法正确率 算法 算法正确次数正确率正确次数 正确率低浓度 中浓度 高浓度 探 测 与 控 制 学 报结论本文基于 激光引信测距原理和 散射理论，采用蒙特卡罗方法建立了 激光在云雾中传输的数学模型。该模型根据光子空间位置和云雾粒径分布实现了光子在云雾中的多次散射和多径传播。通过仿真不同干扰条件下目标回波信号和云雾回波信号的差别，采用固定距离的方法并利用 算法对回波数据进行处理。仿真结果表明 算法在云雾环境中有良好的抗干扰能力，在低、中和高浓度云雾条件下，目标识别正确率均高于 算法，为激光引信抗云雾干扰提供了参考。参 考 文 献：胡俊雄，张艳激光引信抗干扰技术综述制导与引信，（）：陈鹏，赵继广，宋一铄激光近距探测抗气溶胶干扰方法研究进展兵器装备工程学报，（）：，（）：，（）：李华，王彦斌，邹前进，等包容性非球形气溶胶粒子的散射特性 及 其 对 激 光 传 输 的 影 响 研 究 光 子 学 报，（）：王凤杰，刘锡民，陆长平相干激光引信在云雾干扰下的回波特性研究红外与激光工程，（）：，：，：，：，：陈慧敏，王凤杰，冯星泰，等 调频连续波激光引信云雾回波特性仿真红外与激光工程，（）：，（）：宋承天，张垚彦，刘强，等 气溶胶环境影响下连续波激光引信回波特性兵工学报，（）：，（）：，：，（）：，（）：，：，（）：刘海明等：调频连续波激光引信精确定距幅度平方相干算法