一种快速的动态等离子鞘套电波传播计算方法
【技术领域】
[0001] 本发明涉及无线电波传播领域，具体涉及一种动态等离子鞘套电波传播计算方 法。
【背景技术】
[0002] 高超声速飞行过程中会在飞行器周围包覆一层高温热致等离子体层，被称为等离 子鞘套。鞘套中以自由电子为主的带电粒子将会吸收、反射和散射电磁波，使得电磁信号发 生严重衰减，严重时将导致通信信号中断（黑障现象）。等离子鞘套中的电波传播问题的研 究对高超声速飞行器的测控通信体制设计至关重要。
[0003] 传统的等离子鞘套（体）中的电波传播计算针对的多是稳态等离子体，经典的方 法有时域有限差分法、传输矩阵方法以及这些方法的改进方法等。稳态等离子体中的电波 计算获取的是电波透射系数和反射系数，这一结果反映的是特定高度、速度、攻角下的稳态 环境下的电波衰减结果。实际上，飞行器的高度、速度和攻角也是缓变的，且流体中仍存在 湍流扰动、烧蚀剥落过程等不稳定的因素，这些不稳定因素通常是快变随机过程。高速飞行 器飞行条件和流场条件的变化将导致等离子物理参数的动态变化，其中最为主要的是电子 密度的动态变化。前述经典电磁计算方法将不能直接应用于高超声速飞行器的动态等离子 体中的电波传播问题。这种动态时变等离子体中的电波传播问题一直是这一领域的难点问 题。
【发明内容】
[0004] 针对高超声速飞行器动态等离子鞘套电波传播问题，本发明提出一种快速的动态 等离子鞘套电波传播计算方法，该方法从电子密度变化快慢的物理机理出发，建立动态等 离子鞘套电子密度的数学模型，然后在此基础上采取准稳态蒙特卡洛方法和经典电磁计算 方法获取这种复杂随机变化介质的电波传播的时变场强结果。该方法是一种简单有效快速 的电波计算方法，可解决动态等离子体的电波传播计算，避免直接采用随机介质建模方法 计算的高难度和高复杂度。
[0005] 为了实现上述目的，本发明采用如下技术方案：
[0006] -种快速的动态等离子鞘套电波传播计算方法，包含周围包裹有等离子体层的飞 行器，所述方法包括如下步骤：
[0007] Sl输入T。时刻非均匀的电子密度分布函数、攻角变化规律和抖动方差，建立动态 等离子鞘套电子密度广义高斯数学模型；
[0008] S2根据电子密度数学模型的广义高斯随机过程，产生一个攻角变化周期时间t。内 的N个非均匀电子密度样本：
[0009] Nesample (z，t) = [Ne1 (z，t!)，Ne2 (z，t2)，···，NeN(z，tN) ]，tN< 10;
[0010] 其中Ne1(Zd1), i = 1，2，...，NSt1时刻非均匀电子密度样本，Z表示电子距离 飞行器表面的距离；
[0011] S3针对h时刻非均匀电子密度样本Ne i (Lt1), i = 1，2，...，N，建立均匀分层电 磁参数模型；
[0012] S4利用步骤S3中的均匀分层电磁参数模型，通过计算得到t。时间内电波传播N 个时变场强幅度结果Yanip (i)和时变场强相位结果Ypha (i)。
[0013] 需要说明的是，步骤Sl具体如下：
[0014] 1. 1)输入T。时刻等离子高度和速度下的非均匀的稳态电子密度#^,(ζΚΛ妁和碰 撞频率
表示Τ。时刻（錢_方位角度上的与飞行器表面的距离为Z 的电子密度，％表示Τ。时刻（认的方位角度上的与飞行器表面的距离为ζ的碰撞频 率；
[0015] 1. 2)输入攻角变化规律函数f (t)，得到等离子体缓变拟合参数a(t)和b(t)，a(t) 表示幅度变化因子函数，b (t)表示厚度变化因子函数，t表示时间；
[0016] 1. 3)输入电子密度的归一化抖动方差σ，建立动态时变电子密度的数学模型：
[0019] 其中(武州为t时刻(久妁方位角度上的与飞行器表面距离为ζ的电子密 度，u为某时刻的非均匀电子密度分布函数，由#% ,幅度变化因子函数和厚度变化 因子函数确定，Gauss表示高斯函数。
[0020] 进一步需要说明的是，需要说明的是，步骤1.2)中，有
[0022] b(t) = bXf (t)+l, -I < b < 0 ；
[0023] aJP a 2为取值范围为0到1的系数，b为取值范围为-1到0的系数。
[0024] 需要说明的是，步骤S3的具体实施如下：
[0025] 3. 1)输入电磁波角频率ω，以及均匀分层的划分层数N_;
[0026] 3. 2)根据产生的电子密度样本，得到h时刻样本下第m层等离子特征频率ω p,
其中ε c为真空中绝对介电常数，P表示等离子体，Ne ^(z， D为h时刻样本下第m层z处的电子密度，e为自由电子电荷数，Iiif3为自由电子质量；
[0027] 3. 3)建立分层介电常数模型，&时刻样本下第m层的复介电常数：
[0030] <力,）为h时刻样本下第m层相对复介电常数，下标r表示相对，V "为时刻t时等 离子碰撞频率，ω为入射电磁波角频率，/ = 了为复数虚部单位。
[0031] 需要说明的是，步骤S4包括如下步骤：
[0032] 4. 1)计算h时刻样本下第m层等离子电磁波幅度衰减a Jt1)和相移Mt1),其 中m表示等离子层数，me (1，Ν_);
[0033] 4. 2)得到h时刻样本下电磁波穿过整个等离子体的电磁波幅度衰减Yallip U1)和相 位偏移YphaU1):
[0036] 4. 3)重复步骤4. 1)和4. 2)，得到t。时间内N个时变场强幅度结果Y _⑴和时变 场强相位结果Ypha (i)
[0039] 进一步需要说明的是，步骤4. 1)中，a Jti)和PmUi)计算如下：
[0042] 其中，光速c = 3Χ 108m/s，Re[ ·]表示取实部，Im[ ·]表示取虚部。
[0043] 本发明的有益效果在于：
[0044] 1、利用广义高斯过程建立电子密度变化的动态随机过程，克服了无法从飞行实验 中获取数据的困难；
[0045] 2、利用准稳态蒙特卡洛计算方法计算时变等离体中的电波传播问题，大大地降低 直接随机介质建模和计算的困难；
[0046] 3、对于不同条件的动态等离子体，可以用统一的数学形式来进行描述；
[0047] 4、适用于再入动态等离子鞘套电波传播计算，也可适用于临近空间飞行器长期飞 行的动态等离子鞘电波传播计算，所提出的方法可为动态等离体体中信号传输特性、信道 特性及信道模型研究的提供电波理论计算支撑。
【附图说明】
[0048] 图1为本发明的实施流程示意图；
[0049] 图2为本发明一定条件下T。时刻非均匀电子密度分布示意图；
[0050] 图3为本发明正弦规律攻角变化规律示意图；
[0051] 图4为本发明Τ。时刻开始0.1 s内时间内电子密度动态变化示意图；
[0052] 图5为本发明&时刻均匀分层介电常数模型原理示意图；
[0053] 图6为本发明一定条件下计算的得到的时变幅度结果；
[0054] 图7为本发明一定条件下计算的得到的时变相移结果。
【具体实施方式】
[0055] 以下将结合附图对本发明作进一步的描述，需要说明的是，本实施例以本技术方 案为前提，给出了详细的实施方式和具体的操作过程，但本发明的保护范围并不限于本实 施例。
[0056] 如图1所示，一种快速的等离子鞘套电波传播计算方法包括如下步骤：
[0057] SlOl输入T。时刻非均匀的电子密度分布函数、攻角变化规律和抖动方差，建立动 态等离子鞘套电子密度数学模型。
[0058] SI. 1 :输入T。时刻等离子高度和速度下的非均匀的稳态电子密度(z|夂的，如 图2所示，碰撞频率>^(2成的恒定=2G，Α^(ζ|艮的表示τ。时刻说(6))方位角度上的与飞 行器表面距离为z的电子