一种目标近场多普勒回波快速仿真方法技术

本发明专利技术公开了一种目标近场多普勒回波快速仿真方法，包括：步骤S1建立弹目交会坐标系，设置弹目相对速度、弹目相对姿态、脱靶量和脱靶方位；步骤S2、进行分级矩阵的构造；步骤S3、填充当前采样时刻激励矩阵，采用迭代方法或直接法求解矩阵方程，获取当前采样时刻后向雷达截面；步骤S4、根据雷达方程计算近场多普勒回波功率；步骤S5、重复步骤S3和步骤S4，计算目标近场多普勒回波时间序列，本发明专利技术实现了回波仿真各采样时刻阻抗矩阵逆矩阵的重用，采用分级矩阵大幅提升了阻抗矩阵求逆的计算效率，对于弹目交会回波仿真这类多右端项求解问题，大幅提高了仿真速度。提高了仿真速度。提高了仿真速度。

【技术实现步骤摘要】
一种目标近场多普勒回波快速仿真方法


[0001]本专利技术涉及雷达目标电磁散射特性仿真
，具体涉及一种目标近场多普勒回波快速仿真方法。

技术介绍

[0002]目标近场多普勒回波在引信设计、研制和评估阶段均发挥着重要作用。常用的目标近场多普勒回波仿真方法包括基于散射中心模型、基于高频近似方法和基于数值方法三类。
[0003]现有技术中基于散射中心模型和基于高频近似方法的近场多普勒回波仿真速度较快，但对于复杂材质和腔体等特殊部件仿真精度有限。基于数值方法的仿真精度较高，但仿真速度难以满足要求。

技术实现思路

[0004]本专利技术的目的是为了提供一种目标近场多普勒回波快速仿真方法。此方法旨在解决传统雷达目标电磁散射特性仿真方法中对于复杂材质和腔体等特殊部件仿真精度有限，仿真速度难以满足要求的问题。
[0005]为达到上述目的，本专利技术通过以下技术方案实现：
[0006]一种目标近场多普勒回波快速仿真方法，包括：
[0007]步骤S1、建立弹目交会坐标系，设置弹目相对速度、弹目相对姿态、脱靶量和脱靶方位；
[0008]步骤S2、建立目标模型的几何多级分区，根据所述多级分区，进行分级矩阵的构造；
[0009]步骤S3、填充当前采样时刻激励矩阵，采用迭代方法或直接法求解矩阵方程，获取当前采样时刻后向雷达截面；
[0010]步骤S4、根据雷达方程，计算近场多普勒回波功率；
[0011]步骤S5、重复步骤S3和步骤S4，计算目标近场多普勒回波时间序列，
[0012]优选的，在步骤S2中，所述分级矩阵的构造包括：
[0013]步骤S201：建立目标模型的基函数多级分区结构；
[0014]步骤S202：建立系统矩阵多级分块结构；
[0015]步骤S203：压缩填充系统矩阵。
[0016]优选的，在步骤S201中，所述建立目标模型的基函数多级分区结构包括：通过三角形网格离散所述目标表面，得到总数为N的局域基函数，将所述目标拆分为两个子区域，每个所述子区域包含近似相同的N/2个未知量，采用递归拆分所述目标，直至每个所述子区域包含的所述未知量不超过预设阈值n，根据基函数的空间分布，构成均分二叉树。
[0017]优选的，在步骤S202中，所述建立系统矩阵多级分块结构包括：对于所述基函数多级分区结构中位于同一层数的两个基函数组，根据所述同一层数的两个基函数组组中心的
相对距离判断远区耦合或近区耦合。
[0018]优选的，判断远区耦合或近区耦合采用相容性条件，所述相容性条件的表达式为：
[0019]max{D(Ω
i
)D(Ω
j
)}＜η
·
dist(Ω
i
Ω
j
)
[0020]其中，Ω表示基函数组；
[0021]Ω
i
表示第i个基函数组；
[0022]Ω
j
表示第j个基函数组；
[0023]D表示几何直径；
[0024]D(Ω
i
)表示第i个基函数组的几何直径；
[0025]D(Ω
j
)表示第j个基函数组的几何直径；
[0026]dist表示两基函数组之间的几何距离；
[0027]η为一预设正实数。
[0028]优选的，若所述同一层数的两个基函数组满足相容性条件，则为远区耦合；若所述同一层数的两个所述基函数组不满足相容性条件，则为近区耦合。
[0029]优选的，在步骤S203中，所述压缩填充系统矩阵包括：若两个所述基函数组为远区耦合，得到对应矩阵块有低秩特性，采用快速低秩压缩方法填充，所述对应矩阵块表示成两个低秩矩阵相乘的形式。若两个所述基函数组为近区耦合，则对下一层子组的两个所述基函数组进行判断，依次类推，直到最底层子组的两个所述基函数组，所有近区耦合的两个所述基函数组采用矩量法直接填充。
[0030]优选的，在步骤S3中，所述填充当前采样时刻激励矩阵包括：根据当前采样时刻的弹目相对姿态和天线方向图，确定局部照射区域，填充所述激励矩阵。
[0031]优选的，在步骤S4中，所述近场多普勒回波功率的计算表达式为：
[0032][0033]其中，P
r
为接收天线处的接收功率；
[0034]P
t
为发射天线的辐射功率；
[0035]R为天线到目标的距离；
[0036]G
t
为雷达发射天线的增益；
[0037]G
r
为雷达接受天线的增益；
[0038]λ为入射波波长；
[0039]σ为后向雷达截面。
[0040]与现有技术相比，本专利技术具有以下有益效果：
[0041]本专利技术通过积分方程框架下的直接求解法，实现了回波仿真各采样时刻阻抗矩阵逆矩阵的重用，采用分级矩阵大幅提升了阻抗矩阵求逆的计算效率，对于弹目交会回波仿真这类多右端项求解问题，计算效率显著优于现有技术中多层快速多极子迭代类方法，使得仿真速度满足要求。
附图说明
[0042]为了更清楚地说明本专利技术的技术方案，下面将对描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本专利技术的一个实施例，对于本领域普通技术人员
来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图：
[0043]图1为本专利技术一实施例提供的一种目标近场多普勒回波快速仿真方法的流程示意图；
[0044]图2为本专利技术一实施例提供的目标近场多普勒回波快速仿真场景示意图；
[0045]图3为本专利技术一实施例提供的目标近场多普勒仿真效果图。
具体实施方式
[0046]以下结合附图1
‑
3和具体实施方式对本专利技术提出的一种目标近场多普勒回波快速仿真方法作进一步详细说明。根据下面说明，本专利技术的优点和特征将更清楚。需要说明的是，附图采用非常简化的形式且均使用非精准的比例，仅用以方便、明晰地辅助说明本专利技术实施方式的目的。为了使本专利技术的目的、特征和优点能够更加明显易懂，请参阅附图。须知，本说明书所附图式所绘示的结构、比例、大小等，均仅用以配合说明书所揭示的内容，以供熟悉此技术的人士了解与阅读，并非用以限定本专利技术实施的限定条件，故不具技术上的实质意义，任何结构的修饰、比例关系的改变或大小的调整，在不影响本专利技术所能产生的功效及所能达成的目的下，均应仍落在本专利技术所揭示的
技术实现思路
能涵盖的范围内。
[0047]鉴于传统雷达目标电磁散射特性仿真方法存在的不足，为了保证对于复杂材质和腔体等特殊部件仿真精度，提高仿真效率，本实施例提供了一种目标近场多普勒回波快速仿真方法，包括以下步骤：
[0048]步骤S1、建立弹目交会坐标系，设置弹目相对速度、弹目相对姿态、脱靶量和脱靶方位；
[0049]步骤S2、建立目标模型的几何多级分区，根据所述本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种目标近场多普勒回波快速仿真方法，其特征在于，包括：步骤S1、建立弹目交会坐标系，设置弹目相对速度、弹目相对姿态、脱靶量和脱靶方位；步骤S2、建立目标模型的几何多级分区，根据所述多级分区，进行分级矩阵的构造；步骤S3、填充当前采样时刻激励矩阵，采用迭代方法或直接法求解矩阵方程，获取当前采样时刻后向雷达截面；步骤S4、根据雷达方程，计算近场多普勒回波功率；步骤S5、重复步骤S3和步骤S4，计算目标近场多普勒回波时间序列。2.如权利要求1所述的目标近场多普勒回波快速仿真方法，其特征在于，在步骤S2中，所述分级矩阵的构造包括：步骤S201：建立目标模型的基函数多级分区结构；步骤S202：建立系统矩阵多级分块结构；步骤S203：压缩填充系统矩阵。3.如权利要求2所述的目标近场多普勒回波快速仿真方法，其特征在于，在所述步骤S201中，所述建立目标模型的基函数多级分区结构包括：通过三角形网格离散所述目标表面，得到总数为N的局域基函数，将所述目标拆分为两个子区域，每个所述子区域包含近似相同的N/2个未知量，采用递归拆分所述目标，直至每个所述子区域包含的所述未知量不超过预设阈值n，根据基函数的空间分布，构成均分二叉树。4.如权利要求3所述的目标近场多普勒回波快速仿真方法，其特征在于，在所述步骤S202中，所述建立系统矩阵多级分块结构包括：对于所述基函数多级分区结构中位于同一层数的两个基函数组，根据所述同一层数的两个基函数组组中心的相对距离判断远区耦合或近区耦合。5.如权利要求4所述的目标近场多普勒回波快速仿真方法，其特征在于，判断远区耦合或近区耦合采用相容性条件，所述相容性条件的表达式为：max{D(Ω
i
) D(Ω
j
)}＜η
·
dist(Ω
i Ω
j

【专利技术属性】
技术研发人员：高鹏程，黄志勇，安锐，都妍，张润俊，
申请(专利权)人：上海无线电设备研究所，
类型：发明
国别省市：