基于体目标雷达截面积的毫米波雷达道路场景仿真​

在智能驾驶测试领域， 毫米波雷达射频仿真测试以效率高、重复性强等特点被广泛应用，雷达目标的模拟方法随之成为射频仿真的研究重点。

目前， 通过雷达模拟器生成一个点目标来实现目标模拟是行业内的主流方法， 众多研究人员通过获取雷达截面积并进行修正， 以及控制优化点目标回波，提升点目标模拟方法的仿真精度。

在道路交通场景中， 处于高频区的复杂目标由于自身表面特征复杂多样而产生多散射中心， 通过点目标模拟无法完整表征复杂目标。 对此， 体目标模拟是一种解决方案。 韩晓东等［４］阐述了多散射点模型在雷达仿真中的必要性， 分析了体目 标仿真的相关参数，但没有提出一套完整的体目标仿真测试方案。

笔者通过分析复杂目标雷达截面积特性的差异及多散射中心特性， 提出基于体目 标雷达截面积的毫米波雷达道路场景仿真方法， 通过获取体目标雷达截面积和位姿状态， 求解道路场景中的多散射点位姿序列。 笔者通过实例验证了这一仿真方法的可

行性。

复杂目标的不同表面特征及其所处场区造成雷达截面积特性的差异， 这也是射频仿真中的主要误差来源。 分析复杂目标雷达截面积特性， 是实现高精度道路场景仿真的重要前提。

通过经典远场条件公式，可以将目标的散射场区近似划分为菲涅耳区和夫琅禾费区 ［５］ ， 即近场和远场散射区域。 对于复杂目 标而言， 表面结构的差异会造成远场条件下雷达截面积的测量误差。 Ｋｎｏｔｔ等 ［６］ 研究发现， 在经典远场条件下， 具有强棱边结构特征的目 标雷达截面积存在至少 ６ ｄＢ 的误差。雷达与目标之间的远场距离 Ｒ 为：

式中：Ｄ 为目标有效散射最大尺寸；λ 为雷达工作波长；Ｎ 为修正因数，取 ４ ～ １０。由此可知，当雷达的工作频率较高时，复杂目标的远场距离较大。在道路场景中， 车载毫米波雷达的作用范围通常处于近场区， 对近场条件下的雷达截面积特性进行分析十分有必要。

近场条件下目标雷达截面积不仅是目标自身特性的函数，而且与入射和散射条件有关。广义雷达散射截面积的提出对传统雷达截面积定义公式进行了推广，考虑近场条件下的目标散射性质，并将远场雷达截面积公式推论为从 Ｒ 到∞ 时的一种特殊情

况 ［７-８］ 。广义雷达散射截面积 σ 为：

式中：为目标散射场强；为目标入射场强； 为雷达接收天线极化方向。

此时， 广义雷达散射截面积为一个包含幅度和相位信息的复数量，将其转化为复数二次方根，有：

式中：ｋ 为波数。

当目 标处于高频区时， 目标表面特征成为雷达截面积的主要影响因素。目 标由于局部性原理， 存在多散射中心 ［９］ 。车辆目 标不同结构对应的主要散射情况如图 １ 所示。

图 １ 中， 车辆目标在不同入射角度下形成镜面散射、边缘绕射、多次散射等各类散射中心，在对车辆目标进行场景仿真时，应考虑多散射点的幅度和坐标随时间序列的变化。通过体目标模拟，可以进一步分析雷达的目标跟踪能力，以及检测概率、虚警概率、信噪比等检测性能。

为了便于求解复杂目标雷达截面积， 在计算多散射点的雷达截面积时通常采用分别预估局部等效散射中心的散射场，再进行矢量叠加的方法 ［１０］ 。 复杂目标总雷达截面积为：

式中：为第 ｋ 个等效散射中心的雷达截面积；ｎ 为等效散射中心数量。

为了实现体目标射频仿真， 需要分别对被测目标的特征参数和位姿状态进行模拟。 特征参数中，目标雷达截面积计算是仿真的重点。 位姿状态模拟时，主要将目标的运动轨迹参数化， 输出目标位置、姿态、速度、加速度等信息。 通过获取参数， 可以对仿真场景进行搭建。

在体目标射频仿真的过程中， 对道路场景进行参数化分析，通过目标特征参数获取、目标运动特性建模、多散射点序列构建等环节搭建仿真道路场景。道路场景仿真流程如图 ２ 所示。

仿真道路场景搭建的具体实施步骤如下：

（１） 对某一特定道路场景进行参数化分析， 主要对交通参与者的类别、数量，以及各个目标与雷达的相对位置关系、姿态角、相对速度、相对加速度等参数进行信息提取；

（２） 对交通参与者建立用于电磁散射特性计算的面元模型，并对目标雷达截面积进行仿真计算，或通过实际测量方法获取目标雷达截面积；

（３） 通过雷达坐标系与目 标坐标系的坐标变换，将从场景中提取的目 标运动轨迹信息转化为随时间变化的目标位置、姿态序列；

（４） 计算目 标在不同姿态角 下的散射中心场强，结合雷达天线方向图与波束指向矢量对各散射点能量进行排序， 确定对目 标雷达截面积贡献较大的主散射中心；

（５） 通过坐标变换方法求解多散射点在目标坐标系中的位置信息， 将对应雷达截面积注入目 标位姿参数序列，形成随时间变化的多散射点分布序列；

（６） 通过三元组阵列模拟体目标多散射点分布序列，实现道路场景的目标射频仿真。在场景搭建及目 标仿真过程中， 目 标雷达截面积的获取及雷达坐标系与目标坐标系坐标变换是两个关键步骤，笔者将进行重点讨论。

目标雷达截面积主要通过理论计算或实际测量方法获取。 对于处于高频区的目 标， 可通过高频近似方法对雷达截面积进行预估计算 ［１１］ 。 实际测量时，主要通过动态测量、缩比模型测量、全尺寸静态测量等方法获取目标雷达截面积 ［１２］ 。

对于复杂目标， 采用单一理论方法计算雷达截面积存在局限性。 例如几何光学法和物理光学法能够准确预估镜面反射场，但是无法考虑绕射效应，而等效电磁流等理论则可较好地预估边缘绕射场。 结合多种理论方法， 可以实现高精度的雷达截面积预

估计算。

弹跳射线法是较为常用的高频方法， 这一方法结合几何光学法和物理光学法，通过路径追踪、场强追踪等步骤对经多次反射后的反射波射线路径坐标及射线反射场强进行求解 ［１３-１４］ 。 根据 Ｓｎｅｌｌ 定律确定反射方向，反射波射线路径坐标 ｒ（ ｘ，ｙ，ｚ） 为：

对于边缘绕射场， 可以基于等效电流物理绕射理论，运用等效边缘电磁流法求解，并叠加复杂目标的多个面元和棱边散射场计算总雷达截面积。

场景仿真时， 在雷达坐标系下获取的目标车辆运动轨迹信息需要通过坐标变换转换为目标坐标系下的位姿序列，进一步生成多散射点位姿序列。坐标变换如图 ３ 所示。图 ３ 中， ＯＸＹＺ 为以雷达为中心建立的雷达坐标系，为车辆目标所在位

置，对应坐标为（ ｘ １ （ ｔ） ， ｙ １ （ ｔ） ， ｚ １ （ ｔ） ） 。目标在动过程中同时存在姿态的变化，即存在三个欧拉角：横滚角 η（ ｔ） 、俯仰角 γ（ ｔ） 、方位角 α（ ｔ） 。通过以上六自由度参数可以对目标的位姿关系进行完整描述。

以 Ｏ ｔ 为中心建立目 标坐标系 Ｏ ｔ Ｘ ｔ Ｙ ｔ Ｚ ｔ ， 目标在目标坐标系中的坐标表示为（ ｘ ｔ （ ｔ） ， ｙ ｔ （ ｔ） ， ｚ ｔ （ ｔ） ） ，则雷达坐标系与目标坐标系的坐标变换关系为：

随着目标运动， 多散射点在目标坐标系中的位置序列可以通过式（８） 求得。 同样， 雷达相对目 标的位置为：

两个坐标系的三阶转换矩阵 Ｐ 中， 各元素表示分别为：

选取相邻车道车辆超越主车这一道路场景实例对笔者所述仿真方法进行可行性验证，如图 ４ 所示。设置场景参数，主车车速 ｖｆ 为 ６ ｍ／ｓ， 车辆目标速度ｖｔ 为 １０ ｍ／ｓ，两车几何中心侧向距离为 ４ ｍ， 场景时间序列 ０ 时刻两车为平行状态，运行时间为 ７５ ｓ。

结合弹跳射线法与等效边缘电磁流法， 对车辆的面元模型进行雷达截面积仿真计算。 仿真参数中，频率为 ２３.７５ ＧＨｚ ～ ２４.２５ ＧＨｚ， 步长为 ０.０１２ ５ＧＨｚ，俯仰角为 ９０°， 方位角为 ０° ～ １８０°。 目标雷达截面积仿真结果如图 ５ 所示。

对每条射线的雷达截面积贡献在二维图像域中进行累加，可以得到不同方位角的雷达截面积成像图，如图 ６ 所示。 成像参数中，方位角中心值间隔为５°，扫角步长为 0.0３°， 扫角宽度为 １.２°， 极化方式为 ＶＶ。

对各散射点能量进行排序， 选取三个强散射点作为主散射中心，生成多散射点位姿及散射场序列。主散射中心幅度序列如图 ７ 所示， 主散射中心坐标序列如图 ８ 所示。

根据三元组等效相位中心理论及三元组幅度质心公式， 可以通过三元组天线阵列系统对某时刻目标的三个主散射中心进行幅相参数设置， 并发出回波信号，从而实现道路场景的射频仿真。

　 结束语

针对复杂目标雷达截面积具有的多散射中心特性，笔者提出基于体目 标雷达截面积的毫米波雷达道路场景仿真方法， 对特定道路场景中体目标多散射点进行模拟，并对这一仿真方法进行理论分析和应用实例验证。

基于体目标雷达截面积的射频仿真测试对分析雷达检测能力具有重要意义， 基于体目标雷达截面积的毫米波雷达道路场景仿真方法为实现高精度道路场景射频仿真提供了思路， 在仿真测试能力及仿真精度提升方面具有参考价值。