一种基于时变等离子体鞘套识别技术仿真方法

1.本发明属于等离子体鞘套与目标识别问题的分析领域，特别是一种基于时变等离子体鞘套识别技术仿真方法。


背景技术：

2.据已有文献国内外对包覆等离子体目标识别研究较少，因此开展包覆等离子体鞘套目标的isar成像和目标识别的研究是在国防建设与对空天信息安全的需求背景下进行的，通过获取等离子体对目标远区电场的影响并与isar成像可以得到等离子体鞘套对 isar成像影响。
3.目前，对于包覆等离子体目标的识别，国外鲜有文章发表，国内的研究也处于起步阶段，停留在等离子体对成像的影响并未深入研究，到目前为止，没有对包覆等离子体目标识别有个明确定义和方法，在该背景下，研究包覆等离子体鞘套的高超声速目标的识别技术具有重要的意义。


技术实现要素：

4.本发明的目的在于提供一种基于时变等离子体鞘套识别技术仿真方法，实现包覆等离子体目标的识别，相干性能得到了显著提高，可以获得目标所在位置的信息。
5.实现本发明目的的技术解决方案为：一种基于时变等离子体鞘套识别技术仿真方法，包括以下步骤：
6.建立空间电磁场仿真模型和流场仿真模型，并采用四面体和六面体混合分别对两模型进行网格离散，得到两模型的结构信息，包括四面体、六面体的单元信息以及节点信息；
7.将飞行器对应流场模型区域的节点信息映射到电磁场模型区域相应位置，求出电磁场模型每个单元对应的等离子体参数；
8.利用时域不连续伽辽金方法对散射问题进行求解，获取等离子体影响下的回波信号；
9.通过对等离子体密度的测量构建补偿因子，对等离子体影响下的回波信号进行补偿；
10.通过上述步骤的仿真结果确定目标位置信息。
11.进一步地，所述流场模型的节点信息包括原流场四面体网格上的温度、电子密度、压强及各气体质量分数。
12.进一步地，所述网格离散采用ansys对模型进行网格划分，通过赋不同材料号将等离子体区域和其他区域进行区分，根据材料号找到所属的飞行器网格。
13.进一步地，所述每个单元对应的等离子体参数为：
[0014][0015]
νc(i,j,k)＝5.2e18
·
ne(i,j,k)
·
kk
·
t(i,j,k)
ꢀꢀꢀ
(1.2)
[0016][0017]
式中，i,j,k为x,y,z方向的流场点坐标，t代表流场点处的温度，kk为波尔兹曼常数，ω
p
和νc分别是等离子体频率与碰撞频率，ne为电子密度，εr为等离子体相对介电参数。
[0018]
进一步地，所述利用时域不连续伽辽金方法对散射问题进行求解，获取等离子体影响下的回波信号获具体包括：
[0019]
步骤3-1，构建等离子体散射问题方程组，对方程组进行伽辽金法测试、未知量展开和矢量恒等变换，获取新的方程组；
[0020]
步骤3-2，将迎风通量并代入新的方程组中，求取近场值；
[0021]
步骤3-3，算出近场值后，在散射场设置外推面，通过求外推面上的面电流与面磁流，再由面电流与面磁流计算得出电流矩与磁流矩，获得时域远场的电场值，将时域远场值进行傅里叶变换求取物体rcs；
[0022]
步骤3-4，求取等离子体鞘套影响下的雷达回波。
[0023]
进一步地，所述步骤3-1具体包括：
[0024]
构建等离子体散射问题方程组为:
[0025][0026][0027][0028]
其中e为电场强度，单位为v/m，h为磁场强度，单位为a/m，ε为介电常数，单位是f/m，μ为磁导率，单位是ω/m，ε0为真空介电常数，ε
∞
为无穷远处介电常数，值为1，ε
drude
为等离子体介电常数，是nabla算符；
[0029]
将介电常数代入(1.4)第一式中：
[0030][0031][0032][0033]
令：方程组写成：
[0034]
[0035]
对公式(1.8)进行伽辽金测试，测试基函数为ni，对公式(1.8)采用矢量恒等式和高斯定理得到新的方程组为：
[0036][0037][0038][0039]
其中，s为求解区域外表面，en为该外表面的法向分量，ω为求解区域，δt为时间步进。
[0040]
进一步地，所述近场值为：
[0041][0042]
其中，参数参数
[0043]
矩阵元素表达式如下：
[0044][0045][0046][0047][0048]
[0049][0050][0051]
式中，y＝1/z代表本征导纳，上标“+”表示的是相邻体的信息，没有上标表示的是本体的信息，代表本征阻抗；
[0052]
式(1.12)写成dgtd半离散形式为：
[0053][0054]
式中，分别表示电场的测试基函数与展开基函数；分别表示磁场的测试基函数与展开基函数。
[0055]
进一步地，所述步骤3-3具体包括：
[0056]
由等效原理求外推面上的面电流与面磁流分别为：
[0057][0058]
式中，e(r,t)为近场值；
[0059]
由面电流与面磁流计算得出电流矩与磁流矩为：
[0060][0061]
式中，r表示场点的位置矢量，r
′
表示等效源的位置矢量，α表示两个矢量之间的夹角；r
′
cosα为：
[0062][0063]
时域远场电场分量e
θ
和为：
[0064][0065][0066]
物体rcs为：
[0067][0068]
式中，f{
·
}表示的傅里叶变换，和表示的是远场电场值的两个分量，e
inc
(r,t)表示的是入射电场。
[0069]
进一步地，所述步骤3-4中等离子体鞘套影响下的雷达回波为：
[0070][0071]
其中，t1和t2分别是雷达信号在发射和接收周期的传输系数，t1(ne(t-τ),fr(t))表示时变等离子体鞘层在传输周期所产生的影响，它是ne(t)和fr(t)的函数， t2(ne(t-τ),fr(t))是接收周期时变等离子体鞘层所产生的效应，τ为透射系数t1和t2的电子密度时间间隔，fr为发射周期和接收周期的瞬时频率。
[0072]
进一步地，所述补偿后的雷达回波为：
[0073][0074]
其中是透射系数t1的相位，是透射系数t2的相位。
[0075]
本发明与现有技术相比，其显著优点为：1)引入了由等离子体鞘套的电子密度构造了相应的补偿因子，通过对相位的补偿，使回波得到了很好的补偿，相干性能得到了显著提高，可以获得目标所在位置的信息；2)对高超声速飞行器目标识别结果更加准确；3)通过将等离子体包覆目标与逆合成孔径雷达成像技术结合可以填补目前成像技术的空缺；4)对空中隐身目标的监测和识别也是一项迫切和极具意义的工作，开展等离子体目标识别技术的研究也必然是反隐身技术研究的重要组成部分，具有重要的意义和价值；5)反过来也可以为工程上实现目标物的等离子体隐身提供相关经验。
[0076]
下面结合附图对本发明作进一步详细描述。
附图说明
[0077]
图1为获取目标远区电场的流程图。
[0078]
图2为等离子体鞘套对信号影响图。
[0079]
图3为时变等离子体鞘套电子密度分布图。
[0080]
图4为等离子体鞘套影响下的相干结果示意图。
[0081]
图5为补偿后的相干积分结果图。图6为基于时域不连续伽辽金法分析的散射区域划分图。
具体实施方式
[0082]
为了使本技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本技术进行进一步详细说明。应当理解，此处描述的具体实施例仅仅用以解释本技术，并不用于限定本技术。
[0083]
在一个实施例中，结合图1，提供了一种基于时变等离子体鞘套识别技术仿真方法，包括以下步骤：
[0084]
步骤1，建立空间电磁场仿真模型和流场仿真模型，并采用四面体和六面体混合对两者模型进行网格离散，得到模型的结构信息，包括四面体、六面体的单元信息以及节点信息；
[0085]
步骤2，建立包覆等离子体鞘套的飞行器的二维与三维模型图，得到流场的节点信息后，将流场区域的节点信息映射到电磁场模型区域相应位置，求出每个单元对应的等离子体参数；
[0086]
步骤3，利用时域不连续伽辽金方法对散射问题进行求解，获取等离子体影响下的回波信号；
[0087]
步骤4，通过对等离子密度的测量构建补偿因子，对等离子体影响下的回波信号进行补偿；
[0088]
在没有任何补偿的情况下，雷达回波严重失真，相干积累没有峰值。如果电子密度得到了很好的测量，那么通过补偿，回波得到了很好的补偿，相干积累结果在目标所在的位置有一个峰值。
[0089]
步骤2所描述的流场的节点信息包括原流场四面体网格上的温度、电子密度、压强、各气体质量分数等信息，求出每个单元对应的等离子体参数，其方程为：
[0090][0091]
νc(i,j,k)＝5.2e18
·
ne(i,j,k)
·
kk
·
t(i,j,k)
ꢀꢀꢀ
(1.2)
[0092][0093]
式中i,j,k为x，y，z方向的流场点编号，t代表流场点处的温度，kk为波尔兹曼常数，ω
p
和νc分别是等离子体频率与碰撞频率，ω为电子密度，ne为电子密度，εr为等离子体相对介电参数，以上参量都是流场点对应坐标i,j,k的函数。
[0094]
得到流场的节点信息后，在ansys对模型进行网格划分，通过赋不同材料号将等离子体区域和其他区域进行区分，根据材料号找到所属的飞行器网格，此时采用四面体网格离散，将流场区域的节点信息映射到电磁场模型区域相应位置，求出每个单元对应的等离子体参数ω
p
和νc，代入式(1.3)可以求出每个单元的介电参数。
[0095]
所描述的利用时域不连续伽辽金方法对散射问题进行求解具体为：
[0096][0097]
[0098][0099]
其中e为电场强度，单位为v/m，h为磁场强度，单位为a/m，ε为介电常数，单位是f/m，μ为磁导率，单位是ω/m，ε0为真空介电常数，ε
∞
为无穷远处介电常数，值为1，ε
drude
为等离子体介电常数，是nabla算符。
[0100]
将介电常数代入(1.4)第一式中：
[0101][0102][0103][0104]
令：方程组可以写成：
[0105][0106]
对公式(1.8)进行伽辽金测试，测试基函数为ni，对公式(1.8)采用矢量恒等式和高斯定理方程组变为：
[0107][0108][0109][0110]
其中，s为求解区域外表面，en为该外表面的法向分量，ω为求解区域，δt为时间步进，上式的面积分是为了一阶吸收边界条件的引入，不连续是指相邻单元不共用基函数的信息，在各自单元上定义未知量。dgtd方法通过引入数值通量来保证相邻单元间场的连续性，常见的数值通量有中心通量，迎风通量以及惩罚通量。
[0111]
步骤3-1，对方程组(1.8)进行伽辽金法测试、未知量展开、矢量恒等变换；
[0112]
步骤3-2，将下式所示的迎风通量(1.12)代入步骤2-1的(1.9)(1.10)中；
[0113]
一般而言，迎风通量的计算精度较高，本发明的相邻单元间场的连续性由迎风通量条件给出，具体定义如下：
[0114][0115]
式中代表本征阻抗，y＝1/z代表本征导纳，上标“+”表示的是相邻体的信息，没有上标表示的是本体的信息。
[0116]
使用迎风通量来保证相邻单元的连续性，代入上式(1.12)可得：
[0117][0118]
并用六面体叠层矢量基函数展开电场和磁场：
[0119][0120]
将(1.14)代入(1.13)，可得：
[0121][0122]
对上述方程在n时刻采用中心差分展开：
[0123]
[0124][0125][0126]
最终得到：
[0127][0128]
其中，
[0129]
矩阵元素表达式如下：
[0130][0131][0132][0133][0134][0135][0136]
[0137]
最终写成dgtd半离散形式：
[0138][0139]
式中，分别表示电场的测试基函数与展开基函数；分别表示磁场的测试基函数与展开基函数。
[0140]
由于本专利数值方法采用的是dgtd，单元与单元之间是不连续且独立的，t
ee
和t
hh
表示本体基函数和本体基函数的相互作用，当形成等号左边全局矩阵时，该全局矩阵为块对角矩阵，因为本体基函数和相邻体基函数不存在相互作用。
[0141]
步骤3-3，将上式所得的近场值(1.20)的e
n+1
代入下式结果中，
[0142]
算出近场值后，在散射场设置外推面，通过求外推面上的面电流与面磁流，再由面电流与面磁流计算得出电流矩与磁流矩即可获得远场的电场值。
[0143]
由等效原理推导出外推面上的等效电流与磁流：
[0144][0145]
e(r,t)为近场电场值，计算雷达散射截面时，需要设置球坐标系下的俯仰角θ和方位角的值，这两个角度决定了远场观察方向。
[0146]
电流矩n与磁流矩l通过下面的公式求得：
[0147][0148]
在上面的式子中，r表示场点的位置矢量，r
′
表示等效源的位置矢量，α表示两个矢量之间的夹角。则上式中r
′
cosα为：
[0149][0150]
时域远场电场分量e
θ
和可以根据以下公式得到：
[0151][0152][0153]
如图6所示，时域不连续伽辽金法分析散射问题将计算区域分为总场区和散射场区，最外层由abc所包围，各区域划分范围如图所示：最中间部分是散射体，tf代表总场区，sf代表散射场区，cb代表总场-散射场边界，ob代表外推面，最外层abc是吸收边界条件用于
截断整个计算区域。
[0154]
rcs的公式为：
[0155][0156]
上式f{
·
}表示的傅里叶变换，和表示的是远场电场值的两个分量，e
inc
(r,t)表示的是入射电场。
[0157]
得到的时域远场值进行傅里叶变换再代入上式即可得到物体rcs。
[0158]
结合图2，步骤3-4获取等离子体鞘套影响下的雷达回波模型：
[0159]
应该注意的是，等离子体鞘层是时变的，随空间位置变化的。电子密度从飞行器表面到自由空间的分布与从自由空间到飞行器表面的分布不同。此外，对于时变等离子体鞘层，电子密度值随时间而变化。当雷达发射信号通过时变等离子体鞘层时，其电子密度不同于当雷达接收信号通过时变等离子体鞘层时等离子体鞘层的时变电子密度。因此，发送周期和接收周期中的传输系数是不同的。时变等离子体鞘层对雷达信号发射周期和接收周期的影响应该用两种不同的函数来描述。推导了时变等离子体鞘层下的单脉冲雷达回波模型：
[0160][0161]
式中t1(ne(t-τ),fr(t))表示时变等离子体鞘层在传输周期所产生的影响，它是ne(t) 和fr(t)的函数。t2(ne(t-τ),fr(t))是接收周期时变等离子体鞘层所产生的效应。t1和t2分别是雷达信号在发射和接收周期的传输系数。fr(t)是计算t1和t2所需的瞬时频率。t1和t2是时变等离子体鞘层下单脉冲雷达回波模型的关键，其推导如下。
[0162]
电子密度是计算透射系数t的重要参数。由于雷达信号发射和接收的时间间隔为τ，因此计算透射系数t1和t2的电子密度时间间隔为τ。因此，在(1.28)中，发射周期透射周期t1的电子密度为：
[0163][0164]
电子密度的时间跨度为：
[0165][0166]
接收周期透射系数t2的电子密度为
[0167][0168]
电子密度的时间跨度为：
[0169]
[0170]
发射周期和接收周期的瞬时频率fr是计算透射系数的另一个重要参数，时间t雷达信号的瞬时频率为：
[0171][0172]
其中是没有等离子体影响的雷达回波s(t)的相位。
[0173][0174][0175]
获得电子密度与瞬时频率后，可以计算透射系数t1和t2，那么就可以得到受到等离子体鞘套影响下的雷达回波：
[0176][0177]
所述步骤4通过对等离子密度的测量构建补偿因子，对等离子体影响下的回波信号进行补偿；包括：
[0178]
步骤4-1，图3是时变等离子体鞘套电子密度示意图，通过获取电子密度，可构建补偿因子，对等离子体影响后的回波进行补偿。
[0179]
当出现时变等离子体鞘层时，相干积累结果发生变化，目标检测性能下降。该部分分析了随时间变化的等离子体鞘层的相干积累结果，为研究时变等离子体鞘层对目标探测性能的影响打下基础，推导时变等离子体鞘套下的多脉冲雷达回波：对上式(1.36)
[0180]
解调可得：
[0181][0182]
其中tm＝(m-1)t
p
是慢时间
[0183]
将慢时间代入式(1.37)：
[0184][0185]
经过匹配滤波器后为：
[0186]
[0187]
脉冲压缩结果为：
[0188][0189]
上式中代表卷积，对于高超声速飞行器，由于其高超特性，会出现距离偏移，需要进行距离校正，忽略加速度影响，进行一阶keystone变换，将代入上式可得：
[0190][0191]
步骤4-2获取补偿后的雷达回波。
[0192]
雷达回波中等离子体鞘层的补偿方法：时变等离子体鞘层会引起严重的脉内调制和脉间调制，导致脉冲压缩失败和相干性能下降。这是一个棘手的问题，因为时变等离子体鞘层不仅影响多普勒维度(慢时间)，而且影响距离维度(快时间)。分析表明，时变等离子体鞘层引起的相位变化是导致高超声速机载目标检测性能下降的主要原因。根据测量的电子密度，补偿后的雷达回波可以写成：
[0193][0194]
其中是透射系数t1的相位，是透射系数t2的相位。
[0195]
图4是电磁波通过等离子体鞘套的相干积分结果，当等离子体鞘层电子密度是时变时，它会引起时变振幅衰减和相移。时变幅度衰减和相移会使雷达信号失真，引起脉内调制。失真的回波与匹配的滤波器不匹配。峰值消失，副瓣异常升高，导致距离色散和峰值消失。当信号通过时变等离子体鞘层时，电子密度的变化不同，由于的传输系数不同，振幅和相位也不同，这将导致脉间调制。目标包覆等离子体会导致识别性能的下降，在没有任何补偿的情况下，相干积累没有峰值。
[0196]
图5是对时变等离子体密度的补偿后的相干积分结果，从图中可以看出回波得到了很好的补偿，仿真结果表示相干累积结果在目标所在的位置有一个峰值，相干性能得到了显著提高。
[0197]
本发明结合电磁波在等离子体中传播理论，提出等离子体对空间分辨率的影响，接着对包覆等离子体目标进行识别研究，分析有无等离子体对高超目标识别的影响，最后通过估计等离子体鞘套的相位误差并将其补偿，最终通过仿真结果验证了该方法的正确性。
[0198]
需要指出的本发明中涉及的方法未做详细阐述的均采用公知方法，例如建立空间电磁场仿真模型和流场仿真模型、对两者模型进行网格离散、建立包覆等离子体鞘套的飞行器的二维与三维模型等，在此不再累述。
[0199]
以上显示和描述了本发明的基本原理、主要特征及优点。本行业的技术人员应该
了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。