一种临近空间动态等离子鞘套信道建模及模拟的方法
【专利摘要】本发明公开了一种临近空间动态等离子鞘套信道建模及模拟的方法，包括根据大尺度空时频特征划分信道资源，构建大尺度空时频信道资源矩阵；其中高度对应大尺度空域特征，速度对应大尺度时域特征，频点对应大尺度频域特征；针对每个信道资源矩阵，获取动态小尺度信道参数；根据给定的小尺度信道参数，实现小尺度信道模型的模拟；重复近空间飞行轨迹上多个状态的小尺度信道模型，即实现大尺度信道模型的模拟。本发明可适用于不同动态规律下的信道建模，对等离子鞘套下所可能产生的任意幅度相位概率密度分布进行模拟，保证信道建模的可行性；适用于临近空间飞行器等离子鞘套信道建模和航天再入返回飞行器等离子鞘套信道建模。
【专利说明】一种临近空间动态等离子鞘套信道建模及模拟的方法
【技术领域】
[0001] 本发明属于测控通信【技术领域】，尤其涉及一种临近空间动态等离子鞘套信道建模 及模拟的方法。
【背景技术】
[0002] 高超声速飞行器在临近空间高速飞行过程中，其周围将包覆一层高温热致等离子 体（等离子鞘套），这一鞘套内部以自由电子为主的带电粒子将会吸收、反射和散射电磁 波，产生类似金属屏蔽的效应，使得电磁信号发生严重衰减，同时使天线的阻抗特性发生改 变、方向图畸变。这些效应将会导致通信质量恶化，严重时导致通信信号中断，产生黑障现 象。黑障现象将严重影响地面站对飞行器的捕获、跟踪和实时遥测数据传输，导致不可预料 的后果。
[0003] 为了对抗或者削弱缓解等离子鞘套对通信的影响，研究者们提出了很多物理化学 方法试图削弱等离子体，提高通信的可能性。这些措施理论上可一定程度地降低衰减的程 度。对于通信而言，关注的更多是等离子鞘套的特殊性给传输信道带来的影响。众多研究 将这一特殊信道归结为衰减和相位偏移，这一特性通常可在稳态的非均匀等离子中通过电 磁计算方法获取。随着对等离子鞘套认知的深入，发现等离子鞘套不仅具有稳态特性，而且 具有动态特性，等离子鞘套的参数动态分布特性与高度、速度、轨迹、姿态变化、湍流、压力 脉动、烧蚀剥落等众多随机因素及多因素的紧密耦合有关。等离子鞘套的动态性导致了复 杂的随机介质电磁特性，进而对信号产生更为复杂的幅相调制效应。即使在信号穿过等离 子体时，其动态特性引起的幅相调制效应也会严重恶化通信质量。因此为设计适应于这一 复杂信道下的通信体制服务，动态等离子鞘套信道建模的研究至关重要。对动态等离子体 信道的研究及其模型建模的研究几乎处于空白。
[0004] 实际上，从信道的物理特征角度而言，可将等离子的动态性分为大尺度动态性和 小尺度动态性。针对飞行高度（空）、速度（时间）、频点（频）等物理状态引入"大尺度" 动态特性，可进行大尺度空时频信道资源的划分，对大尺度动态特征建模，为通信系统需求 更多的信道资源窗口协同资源优化利用和拓展信道容量服务。针对烧蚀湍流等物理"小尺 度"动态特性，体现微观动态性对信号幅度、相位和功率谱弥散的影响，进行小尺度信道建 模，为通信系统设计更有针对性的编码调制从而提升性能逼近容量界服务。
【发明内容】
[0005] 本发明实施例的目的在于提供一种临近空间动态等离子鞘套信道建模及模拟的 方法，由于目前信道探测在实际在飞实验中无法实现，以及对动态等离子体对信号和信道 特性影响认知不足问题，旨在从理论建模方法和信道参数的理论计算两方面解决动态等离 子体信道模型建模问题，为设计缓解或者适应性的通信技术提供指导服务。
[0006] 本发明实施例是这样实现的，一种临近空间动态等离子鞘套信道建模及模拟的方 法，该临近空间动态等离子鞘套信道建模及模拟的方法包括以下步骤：
[0007] 步骤一，针对每个信道资源矩阵，获取动态小尺度信道参数；通过电磁计算方法获 取远场接收点三维场强数据，获取小尺度信道的幅度和相位概率密度函数，获取小尺度信 道多普勒功率谱函数；
[0008] 步骤二，根据给定的小尺度信道参数，实现小尺度信道模型的模拟；具体包括：
[0009] 第一步，采用均匀帽子舍弃法模拟满足所提取幅度概率分布的随机序列R(k)和 相位概率分布随机序列9 (k);
[0010] 第二步，采用IIR数字滤波器直接设计法设计归一化离散的多普勒成型滤波器 S(k)；
[0011] 第三步，叠加幅度随机序列和相位随机序列生成衰落信道包络复随机序列z(k) = R(k)eJ0(k)；
[0012] 第四步，信道包络复随机序列通过多普勒成型滤波器序列，生成动态等离子鞘套 信道随机序列C p(k) = Z(k) ? S(k);
[0013] 步骤三，针对每个信道资源矩阵，重复步骤二，得到不同接收极化下的动态等离子 鞘套信道随机序列。
[0014] 进一步，在步骤一之前需要根据大尺度空时频特征划分信道资源，构建大尺度空 时频信道资源矩阵；其中高度对应大尺度空域特征，速度对应大尺度时域特征，频点对应大 尺度频域特征；
[0015] 进一步，在步骤一中，通过电磁计算方法获取远场接收点三维场强数据；具体包 括：
[0016] 第一步、根据大尺度特征下对应的N组等离子体电子密度分布及其动态性规律， 构建等离子鞘套电磁参数模型，其中在每组非均匀电子分布下，对非均匀电子密度分层均 匀化处理，每层电磁介质参数模型为：
【权利要求】
1. 一种临近空间动态等离子鞘套信道建模及模拟的方法，其特征在于，该临近空间动 态等离子鞘套信道建模及模拟的方法包括以下步骤： 步骤一，针对每个信道资源矩阵，获取动态小尺度信道参数；通过电磁计算方法获取远 场接收点三维场强数据，获取小尺度信道的幅度和相位概率密度函数，获取小尺度信道多 普勒功率谱函数； 步骤二，根据给定的小尺度信道参数，实现小尺度信道模型的模拟；具体包括： 第一步，采用均匀帽子舍弃法模拟满足所提取幅度概率分布的随机序列R(k)和相位 概率分布随机序列θ (k); 第二步，采用IIR数字滤波器直接设计法设计归一化离散的多普勒成型滤波器S(k); 第三步，叠加幅度随机序列和相位随机序列生成衰落信道包络复随机序列Z(k)= R(k)eJ0(k)； 第四步，信道包络复随机序列通过多普勒成型滤波器序列，生成动态等离子鞘套信道 随机序列 Cp (k) = Z (k) · S (k); 步骤三，针对每个信道资源矩阵，重复步骤二，得到不同接收极化下的动态等离子鞘套 信道随机序列。
2. 如权利要求1所述的临近空间动态等离子鞘套信道建模及模拟的方法，其特征在 于，在步骤一之前需要根据大尺度空时频特征划分信道资源，构建大尺度空时频信道资源 矩阵；其中高度对应大尺度空域特征，速度对应大尺度时域特征，频点对应大尺度频域特 征，其中高度以5km间隔划分，速度以0. 5马赫间隔划分，频点以500MHz间隔划分。
3. 如权利要求1所述的临近空间动态等离子鞘套信道建模及模拟的方法，其特征在 于，在步骤一中，通过电磁计算方法获取远场接收点三维场强数据；具体包括： 第一步、根据大尺度特征下对应的N组等离子体电子密度分布及其动态性规律，构建 等离子鞘套电磁参数模型，其中在每组非均匀电子分布下，对非均匀电子密度分层均匀化 处理，每层电磁介质参数模型为：
<为第m层等离子体相对介电常数，^为自由空间介电常数，ω为入射电波频率，Um 为第m层等离子体碰撞频率，ωρ π为第m层等离子体特征频率，m为等离子体层数； 第二步、采用三维FDTD电磁计算算法获取每组准静态下接收点处的场强Ep，Ep = Exax+Eyay+Eza z，其中 Ex = (a+bi)，Ey = (c+di)，Ez = (e+fi)，ax、ay、az 为 X，y，z 方向的单位 向量； 第三步、重复第二步仿真计算，得到在大尺度动态特征下的N组复场强数据。
4. 如权利要求1所述的临近空间动态等离子鞘套信道建模及模拟的方法，其特征在 于，在步骤一中，获取小尺度信道的幅度和相位概率密度函数，具体包括： 第一步、根据天线接收极化形式，选择相应的极化分量并获取N个接收信号幅度数据 和和相位数据古⑴； 左旋接收时々⑴=丨=丨.."V丨，6⑴=丨九=丨···#}， 右旋接收时々(，)，卜U}，知)=[?,&H};
第二步、对接收信号幅度数据和龟O和相位数据4(0进行统计分析和多项式拟合，得 到幅度相位概率密度函数fK(x)和相位概率密度函数fe (X);得到不同接收极化下的幅度 概率密度函数；得到不同接收极化下的相位概率密度函数。
5. 如权利要求1所述的临近空间动态等离子鞘套信道建模及模拟的方法，其特征在 于，在步骤一中，获取小尺度信道多普勒功率谱函数，具体包括： 第一步、由幅度数据和力(〇和相位数据命?)得到接收信号复数据 /XO = P⑴+ ./如)H，即通过N组复场强数据得到近似的连续时间数据； 第二步、求解接收信号MO的其自相关函数价「) = J· + ，进行傅里叶变换 得到多普勒功率谱函数SM = ;得到不同接收极化下的多普勒功率谱。
6. 如权利要求1所述的临近空间动态等离子鞘套信道建模及模拟的方法，其特征在 于，在步骤二的第一步中均匀帽子法步骤如下，产生f K(x)的随机序列X ; 第一步、产生在（0, ax)上均匀分布的Vix，其中ax是X的最大值； 第二步、产生在（〇，bx)上均匀分布的V2x，其中bx不小于的匕〇〇最大值； 第三步、如果Vix彡fx (Vix)，令X = Vix ;如果不等式不满足，则丢弃Vix和V2x，从第一步 开始重复以上过程； 第四步、得到服从幅度概率分布函数&00的随机序列R(k) =X;不同接收极化下的幅 度模拟随机序列。
7. 如权利要求1所述的临近空间动态等离子鞘套信道建模及模拟的方法，其特征在 于，在步骤二的第一步中均匀帽子法步骤如下，产生f e (X)的随机序列Y ; 第一步、产生在（0, aY)上均匀分布的Viy，其中aY是Y的最大值； 第二步、产生在（0, bY)上均匀分布的V2y，其中bY不小于的fx(x)最大值； 第三步、如果Viy彡fx (Viy)，令Y = Viy ;如果不等式不满足，则丢弃Viy和V2y，从第一步 开始重复以上过程； 第四步、得到服从相位概率分布函数^00的随机序列Θ (k) =Y;不同接收极化下的 相位模拟随机序列。
8. 如权利要求1所述的临近空间动态等离子鞘套信道建模及模拟的方法，其特征在 于，在步骤二中的第三步采用IIR数字滤波器直接设计法设计归一化离散的多普勒成型滤 波器S (k)，具体包括： 第一步、由K个二阶网络级联而成多普勒成型滤波器，系统函数表示为S(Z)为：
式中，A是常数，a」，bp Cp 4是待求第j个滤波器系数，K为级联二阶滤波器个数； 第二步、采用Yule-Walker方程法求解系数求解滤波器系数4K+1个滤波器系数，在（0， 区间内取M点数字频率ωρ? = 1，2……M，M为频点个数，在这些频率上，使得Sd (e>) 和期望S (e>)之间幅度平方误差E最小；其中
^是期 望滤波器的频率响应，WZb) = 是设计滤波器的频率响应； 第三步、离散化Sd(ejw)波形，得到多普勒成型滤波器序列S (k)，其中k的个数与第一步 和第二步中相同。
9.如权利要求1所述的临近空间动态等离子鞘套信道建模及模拟的方法，其特征在 于，在步骤二之后重复步骤一和步骤二近空间飞行轨迹上多个状态的小尺度信道模型，即 实现大尺度信道模型的模拟。
【文档编号】H04B17/391GK104378170SQ201410583216
【公开日】2015年2月25日 申请日期:2014年10月27日 优先权日:2014年10月27日
【发明者】石磊, 李小平, 刘彦明, 周辉, 姚博, 方水汛 申请人:西安电子科技大学