基于多层平面波分解的一维单站rcs近远场转换方法
【技术领域】
[0001] 本发明所涉及的是目标电磁散射与逆散射领域，尤其涉及电磁散射的一维单站 RCS近远场转换方法。
【背景技术】
[0002] 目标RCS测试需要满足远场条件，目前较为成熟的方法有远场地平场测试和紧缩 场测试。然而，对于电大尺寸目标RCS测试，远场地平场需要巨大的测试场地支持，紧缩场 需要昂贵的反射面设施建设。近年来发展活跃的RCS近场测试技术，在不满足远场条件的 近场进行测试，再通过近远场转换得到目标RCS，具有低廉便捷的特点。近远场转换方法是 上述近场测试的关键。
[0003] 在检索到的国内外公开及有限范围发表的文献中，有论文和专利 (CN201410432104)介绍基于近场成像的近远场转换方法，未见基于多层平面波分解的近远 场转换方法。另有国外论文介绍多层平面波分解近远场转换方法，但该方法不能用于一维 单站RCS近远场转换。综上所述，现有技术未公开基于多层平面波分解的一维单站RCS近 远场转换方法。
【发明内容】
[0004] 本发明的目的在于提出一种基于多层平面波分解的一维单站RCS近远场转换方 法，达到非余量采样，具有快速近远场转换、自动天线方向图补偿、自动采样位置补偿、误差 可控等特点，可指导近场采样设置，从而为获取目标RCS(雷达散射截面）提供一种快速便 捷的方法。
[0005] 本发明提供一种基于多层平面波分解的一维单站RCS近远场转换方法，其包含：
[0006] 步骤1:确定近场散射数据采样平面，在该平面上以目标中心为圆心，在最近半径 和最远半径所确定的圆环内，使用任意天线在该圆环的任意位置进行采样，记录每个采样 点的天线接收电压及采样点位置；
[0007] 步骤2:基于近场散射的多层平面波分解公式，即根据近场采样数据的参数和系 统所要达到的精度，计算方程常数系数和矩阵大小；
[0008] 步骤3:根据步骤2计算所得矩阵大小，将单位圆周上的平面波谱均匀离散化，逐 一填充近场转移矩阵元素，完成矩阵方程的建立；
[0009] 步骤4 :用最陡下降的共辄梯度法求解矩阵方程，当迭代收敛时得到目标反射率 方向图；
[0010] 步骤5 :根据RCS与目标反射率方向图之间的关系获得目标RCS。
[0011] 本发明通过一种多层平面波分解快速算法，将近场一维单站散射数据转换为远场 散射数据，从而获得目标RCS。基于优选的实施例可知，本发明所述基于多层平面波分解的 一维单站RCS近远场转换方法，根据加法定理将目标散射近场用多层平面波展开，写成转 移算子与目标反射率方向图函数在单位角谱球上的积分；根据目标扁平体特征将球面积分 简化为圆周积分，依据算法所需要达到的精度截断转移算子，并使用可快速收敛的共辄梯 度法进行矩阵求逆。
[0012] 本发明带来以下有益效果：
[0013] 本发明提出一种在水平面上进行一维单站近场采样的近远场转换算法，其优点在 于近场采样数与模式截断阶数一致，达到了非余量采样，并与算法所要达到的精度具有明 确的关系，使得误差可控；天线方向图和采样位置引起的影响在迭代求逆过程中被完全补 偿了，使得近场测试时可选用任意天线，可在平面内随机位置采样，避免了转台或扫描架等 定位设备的使用，极大化简了近场测试系统。
【附图说明】
[0014] 图1是本发明中的基于多层平面波分解的一维单站RCS近远场转换算法的流程 图。
【具体实施方式】
[0015] 以下结合【附图说明】本发明的较佳实施例。
[0016] 如图1所示为本发明中基于多层平面波分解的一维单站RCS近远场转换算法流程 图。本发明的计算方法原理如下：
[0017] 基于一阶Born近似和单站设定，建立入射场与散射场的辐射反射模型；根据加法 定理和目标的扁平体特征，将上述辐射反射模型用多层平面波展开，转换成单位圆周上平 面波谱的积分；使用共辄梯度法进行矩阵方程求解。
[0018] 计算具体步骤如下：
[0019] 步骤1 :确定近场散射数据采样平面，在该平面上以目标中心为圆心，在最近和最 远半径所确定的圆环内，使用任意天线在该圆环的任意位置进行采样，记录每个采样点的 天线接收电压及采样点位置；设采样点数为M。
[0020] 步骤2 :基于近场散射的多层平面波分解公式，即根据近场采样数据的各项参数 和系统所要达到的精度，计算方程常数系数和矩阵大小；近场散射的多层平面波分解公式 为：
[0021]
[0022] f为入射波矢，g,k和|分别是波数和波矢方向，Z是自由空间中的波阻抗，A 是入射电压，g是采样点位置，匕?Υ是转移算子，ris是与距离无关的因子，表达式为
[0023]
(9)
[0024] 是天线方向图，<2^是与目标有关的函数。
[0025] 在预处理阶段需要计算式（8)的系数
-维圆周积分的积分点 数KL
[0026] KL= 2(L+1) (10)
[0027] L是转移算子?Υ的截断阶数，其选取需要满足如下条件
[0028] L=kd+alog(π+kd) (11)
[0029] d是测试天线与目标的距离，α=-l〇ge，ε是算法所要达到的精度。结合步骤 1，得到转移算子矩阵的大小为ΜΧ&。
[0030] 步骤3 :根据步骤2计算所得矩阵大小，将单位圆周上的平面波谱均匀离散化，逐 一填充近场转移算子矩阵元素，完成矩阵方程的建立；所述矩阵方程是式（8)的离散形式
[0031]
:(1?
[0032] 下标表示矩阵的大小。转移算子矩阵C的元素为
[0033]
[0034] 是第二类球汉克尔函数，乃@ 是勒让德多项式。
[0035] 步骤4 :用最陡下降共辄梯度法求解矩阵方程（12)，当迭代收敛时得到包含目标 反射率方向图的Η矩阵；
[0036] 步骤5 :根据RCS与目标反射率方向图之间的关系获得目标RCS。
[0037]
[0038] 综上所述，本发明中基于一阶Born近似和单站设定建立辐射散射模型，使用加法 定理将近场散射用多层平面波分解，并针对扁平体目标将单位球面上的平面波积分缩减到 了单位圆周上，实现了任意天线任意位置非余量采样下一维单站RCS快速近远场转换方 法，解决了在目标散射近场进行低廉便捷测试并获取误差可控的RCS的需要。
[0039] 尽管本发明的内容已经通过上述优选实施例作了详细介绍，但应当认识到上述的 描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后，对于本发明的 多种修改和替代都将是显而易见的。因此，本发明的保护范围应由所附的权利要求来限定。
【主权项】
1. 一种基于多层平面波分解的一维单站RCS近远场转换方法，其特征在于， 所述转换方法，包含： 步骤1 :确定近场散射数据采样平面，在该平面上W目标中屯、为圆屯、，在最近半径和最 远半径所确定的圆环内，使用任意天线在该圆环的任意位置进行采样，记录每个采样点的 天线接收电压及采样点位置； 步骤2 :基于近场散射的多层平面波分解公式，即根据近场采样数据的参数和系统所 要达到的精度，计算方程常数系数和矩阵大小； 步骤3 :根据步骤2计算所得矩阵大小，将单位圆周上的平面波谱均匀离散化，逐一填 充近场转移矩阵元素，完成矩阵方程的建立； 步骤4 :用最睹下降的共辆梯度法求解矩阵方程，当迭代收敛时得到目标反射率方向 图； 步骤5 :根据RCS与目标反射率方向图之间的关系获得目标RCS。2. 如权利要求1所述的转换方法，其特征在于， 步骤1中的采样点数为M; 步骤2中近场散射的多层平面波分解公式为：(1) 其中，I为入射波矢，系=，k和A:分别是波数和波矢方向，Z是自由空间中的波阻抗，Ui是入射电压，是采样点位置，弓=7:记，1；是转移算子，nS是与距离无关的因子，表达式 为傑 是天线方向图，是与目标有关的函数； 在预处理阶段需要计算式（1)的系数:巧日一维圆周积分的积分点数町 Kl= 2 化+1) (3) L是转移算子的截断阶数，其选取需要满足如下条件L=kd+al〇g(jT+kd) (4) d是测试天线与目标的距离，a= -Ioge，e是算法所要达到的精度； 结合步骤1，得到转移算子矩阵的大小为MX町。3. 如权利要求2所述的转换方法，其特征在于， 步骤3中建立的矩阵方程是式（1)的离散形式 下标表示矩阵的大小；转移算子矩阵C的元素为 (5)片户(2Ay；i)是第二类球汉克尔函数，与(每慮勒让德多项式。4.如权利要求3所述的转换方法，其特征在于， 步骤4中用最睹下降共辆梯度法求解矩阵方程巧），当迭代收敛时得到包含目标反射 率方向图的H矩阵； 步骤5 :根据RCS与目标反射率方向图之间的关系获得目标RCS:
【专利摘要】本发明涉及一种多层平面波分解的一维单站RCS近远场转换方法，根据加法定理将目标散射近场用多层平面波展开，写成转移算子与目标反射率方向图函数在单位角谱球上的积分；根据目标扁平体特征将球面积分简化为圆周积分，依据算法所需要达到的精度截断转移算子，并使用可快速收敛的共轭梯度法进行矩阵求逆。本发明通过多层平面波分解快速算法，将近场一维单站散射数据转换为远场散射数据，达到非余量采样，具有快速近远场转换、自动天线方向图补偿、自动采样位置补偿、误差可控等特点，可指导近场采样设置，从而为获取目标RCS提供一种快速便捷的方法。
【IPC分类】G01S7/41
【公开号】CN105372640
【申请号】CN201510802086
【发明人】贺新毅, 童广德, 徐秀丽, 王晓冰
【申请人】上海无线电设备研究所
【公开日】2016年3月2日
【申请日】2015年11月19日