一种再入动态等离子鞘套马尔科夫信道建模方法技术

本发明专利技术公开了一种再入动态等离子鞘套马尔科夫信道建模方法，首先输入特定条件下等离子鞘套的幅度衰减序列，建立多状态马尔科夫信道模型，利用可逆跳变马尔科夫蒙特卡洛算法求出模型参数，然后根据马尔科夫信道模型状态转移矩阵生成状态序列，最后根据马尔科夫信道模型每个状态下高斯分布参数以及状态序列，生成信号随机衰减序列。本发明专利技术通过可逆跳变马尔科夫蒙特卡洛算法求出最佳的马尔科夫信道的状态数及相应每个状态下随机过程的分布参数，与现有的马尔科夫等离子信道建模方法相比优势在于不需要预设状态数，可减少人为因素引入的误差。

【技术实现步骤摘要】
一种再入动态等离子鞘套马尔科夫信道建模方法
本专利技术属于测控通信
，具体涉及一种再入动态等离子鞘套信道马尔科夫建模及模拟方法。
技术介绍
再入飞行器在高超声速飞行过程中，会在飞行器周围包覆一层高温热致等离子体层，被称为等离子鞘套。鞘套中以自由电子为主的带电粒子将会吸收、反射和散射电磁波，产生类似金属屏蔽的效应，使得电磁信号发生严重衰减。这些效应将会导致通信质量恶化，严重时将导致通信信号中断(黑障现象)。黑障现象将严重影响地面站对飞行器的捕获、跟踪和实时遥测数据传输，导致不可预料的后果。从无线通信角度而言，等离子鞘套信道环境和通信体制的不匹配导致了通信信号中断。等离子体鞘套物理参数具有明显的动态时变特性，其产生主要是由于飞行器姿态调整、湍流扰动、非均匀烧蚀等因素。物理参数的动态时变必然导致电磁介电参数的时变，进而引起电波传播和信号的畸变。目前通过实验已可明显观察到信号的幅相调制效应。因此为寻找适应于这一特殊复杂信道下的通信体制和方法，再入动态等离子鞘套信道模型的建立至关重要。由于影响等离子电磁参数动态变化的因素不同，造成幅度衰减的波动范围较大，若采用单一的随机过程描述将不够准确，因此可以考虑采用多状态马尔科夫链来进行描述，多状态马尔科夫过程比单个随机过程建模方法精确度更高。已有研究采用等概率划分准则的多状态马尔科夫方法对信道进行建模，但是该方法需要预设信道状态数，信道建模的适应性较差。
技术实现思路
针对现有技术的不足，本专利技术旨在提出一种新的多状态马尔科夫信道建模方法，该方法通过可逆跳变马尔科夫蒙特卡洛算法求出最佳的马尔科夫信道的状态数及每个状态下随机过程的分布参数，与现有的马尔科夫等离子信道建模方法相比优势在于不需要预设状态数，可减少人为因素引入的误差。为了实现上述目的，本专利技术采用如下技术方案：一种再入动态等离子鞘套马尔科夫信道建模方法，包含如下步骤：S1输入特定条件下信号穿过等离子鞘套的衰减序列；S2建立等离子鞘套的多状态马尔科夫信道模型，根据步骤S1输入的等离子鞘套衰减序列，通过可逆跳变马尔科夫蒙特卡洛算法求出信道模型参数，包括马尔科夫信道模型状态数、信道模型状态转移矩阵以及每个状态下高斯分布参数；S3利用步骤S2求出的状态转移矩阵生成马尔科夫状态序列zt，t＝1，2，...n，n为序列长度；S4利用步骤S2求出的每个状态下高斯分布参数和步骤S3下产生的马尔科夫状态序列，产生相对应的随机衰减序列仿真值yt，t＝1，2，...n，n为序列长度。需要说明的是，步骤S2的具体实施如下：2.1)将等离子鞘套的马尔科夫信道的每个状态建模为高斯过程，信道的一阶统计特性用混合高斯模型来表示：其中Y表示信道的概率密度函数，K为马尔科夫过程的状态数即混合高斯模型中高斯过程的个数，ωj表示状态概率，即标号为j的高斯密度函数在混合密度函数中所占的权重，并且有表示标号为j的高斯过程，而表示标号为j的高斯过程的均值和方差，其中j＝1，2，...，K；2.2)确定模型参数(K，ωj，μj，σj)的先验分布：状态数K服从泊松分布，μj服从高斯分布，即μj～N(ξ，κ-1)，而σj服从逆伽马分布，即其中β服从伽马分布，即β～G(g，f)，混合高斯过程权重ωj构成ω向量，而ω服从Dirichlet分布，即ω～D(δ1，δ2，....，δK)，其中(ξ，κ，g，f，δ)为超验分布参数，为定值，δ为(δ1，...，δK)的集合；2.3)初始化模型各参数，包括：初始化状态数K(1)、每个状态下高斯分布参数先验分布超参数(ξ，κ，g，f，δ)以及出生概率bk和死亡概率dk＝1-bk，k为状态数的标号；出生概率和死亡概率服从二项分布，并设定迭代次数N以及假定不平衡迭代次数M，M取初始化迭代计数器h；另外，增加参数隐状态序列z，其初始化由下式确定：z(1)i＝argmax{Pr(zi＝j)}i＝1，2....nj＝1，2，...K；其中，yi，i＝1，2，....n为信号衰减序列，n为序列的长度，zi为参数隐状态序列z中标号为i的值，zi(1)则为其初始化值；2.4)进入第h次迭代的分裂合并过程：产生一个(0，1)均匀分布的随机数U1，若U1＜bk，则进入分裂过程；否则，进入合并过程；参数空间更新为其中分别为状态数K、状态概率ωj、标号为j的高斯过程的均值和方差、参数隐状态序列z完成第h次迭代的分裂合并过程后的更新，超参数β因为在分裂过程中没有更新，所以保留第h-1次迭代完成后更新值β(h)；2.5)进入第h次迭代的出生死亡过程：产生一个(0，1)均匀分布随机数U2，若U2＜bk，则进入出生过程，否则进入死亡过程；参数空间更新为其中分别为状态数K、状态概率ωj、标号为j的高斯过程的均值和方差、参数隐状态序列z完成第h次迭代的出生死亡过程后的更新，超参数β因为在合并过程中没有更新，所以保留第h-1次迭代完成后更新值β(h)；2.6)利用Gibbs采样更新参数空间得到第h次迭代结果K(h+1)、ωj(h+1)、μj(h+1)、σj(h+1)、z(h+1)和β(h+1)分别为状态数K、状态概率ωj、标号为j的高斯过程的均值和方差、参数隐状态序列z和超参数β完成第h次迭代后得到的最终更新值；2.7)h＝h+1，若h＜N，重复2.4)-2.6)；2.8)对于参数空间迭代结果x(1)，x(2)，....，x(N)，去除前M次不平衡迭代，得到x(M+1)，x(M+2)，....，x(N)，对其中参数K(M+1)，K(M+2)，...，K(N)作直方图，求出最大值kop，即为最佳状态数；2.9)在参数空间x(M+1)，x(M+2)，....，x(N)中，对于K(i)＝kop，i＝M+1，M+2，...，N的参数(ω(i)，μ(i)，σ(i))，求出数学期望(ω，μ，σ)，即为马尔科夫模型不同状态下的高斯参数；2.10)利用2.3)中的参数隐状态序列z的似然函数和2.9)中求出的马尔科夫模型不同状态下高斯参数，求出隐状态序列t＝1，2，...n，并进一步求出状态转移矩阵P为：P(i，j)＝Nij/Ni；其中，P(i，j)表示从状态i跳转到状态j的概率，Nij为从状态i转移到状态j的个数，Ni为处在状态i的个数，i，j＝1，2，...，n。需要说明的是，步骤2.4)中，所述分裂过程具体如下：2.4.1.1)对于第h-1次迭代完成后最终得到的参数空间x(h)，从j＝1，2，...，K(h)随机选择一个状态js*，参数为以js1，js2表示分裂后状态的标号，分裂后状态数变为K′＝K(h)+1，K(h)和z(h)为第h-1次迭代完成后得到的状态数K和参数隐状态序列z的最终更新值，而分别为第h-1次迭代完成后状态js*的ωj、μj和σj的更新值；2.4.1.2)产生三个随机数u1，u2，u3，满足：u1～beta(2，2)u2～beta(2，2)u3～beta(1，1)；beta为贝塔分布，分裂后状态权重及相应状态高斯分布参数：和分别表示状态js1和js2所对应的第h-1次迭代完成后得到μj的更新值；2.4.1.3)更新参数隐状态序列，将原本的样本点按照如下概率分配给js1和js2：Palloc是分配给js1的概率，q＝1-Palloc是分配给js2的概率本文档来自技高网...

【技术保护点】
一种再入动态等离子鞘套马尔科夫信道建模方法，其特征在于，所述方法包含如下步骤：S1输入特定条件下信号穿过等离子鞘套的衰减序列；S2建立等离子鞘套的多状态马尔科夫信道模型，根据步骤S1输入的等离子鞘套衰减序列，通过可逆跳变马尔科夫蒙特卡洛算法求出信道模型参数，包括马尔科夫信道模型状态数、信道模型状态转移矩阵以及每个状态下高斯分布参数；S3利用步骤S2求出的状态转移矩阵生成马尔科夫状态序列zt，t＝1，2，...n，n为序列长度；S4利用步骤S2求出的每个状态下高斯分布参数和步骤S3下产生的马尔科夫状态序列，产生相对应的随机衰减序列仿真值yt，t＝1，2，...n，n为序列长度。

【技术特征摘要】
1.一种再入动态等离子鞘套马尔科夫信道建模方法，其特征在于，所述方法包含如下步骤：S1输入特定条件下信号穿过等离子鞘套的衰减序列xt，t＝1,2…n；S2建立等离子鞘套的多状态马尔科夫信道模型，根据步骤S1输入的等离子鞘套衰减序列xt，通过可逆跳变马尔科夫蒙特卡洛算法自适应的估计出求出信道模型参数，包括马尔科夫信道模型状态数、信道模型状态转移矩阵以及每个状态下高斯分布参数；S3利用步骤S2求出的状态转移矩阵生成马尔科夫状态序列zt,t＝1,2,...n，n为序列长度；S4利用步骤S2求出的每个状态下高斯分布参数和步骤S3下产生的马尔科夫状态序列，产生相对应的随机衰减序列仿真值yt，t＝1,2,...n，n为序列长度；步骤S2的具体实施如下：2.1)将等离子鞘套的马尔科夫信道的每个状态建模为高斯过程，信道的一阶统计特性用混合高斯模型来表示：其中Y表示信道的概率密度函数，K为马尔科夫过程的状态数即混合高斯模型中高斯过程的个数，ωj表示状态概率，即标号为j的高斯密度函数在混合密度函数中所占的权重，并且有表示标号为j的高斯过程，而表示标号为j的高斯过程的均值和方差，其中j＝1,2,...,K；2.2)确定模型参数(K,ωj,μj,σj)的先验分布：状态数K服从泊松分布，μj服从高斯分布，即μj～N(ξ,κ-1)，ξ为期望均值＝xmin+R，R＝xmax-xmin,κ＝R-2；xmax和xmin分别为输入序列xt的最大值和最小值；而σj服从逆伽马分布，即其中α为伽马分布的形状参数，α＝2；β为伽马分布的尺度参数服从伽马分布，即β～G(g,f)，其中g＝0.2，f＝10/R2；混合高斯过程权重ωj构成ω向量，而ω服从Dirichlet分布，即ω～D(δ1,δ2,....,δK)，其中(ξ,κ,g,f,δ)为超验分布参数，为定值，δ为(δ1,...,δK)的集合；2.3)初始化模型各参数，包括：初始化状态数K(1)、每个状态下高斯分布参数超验分布参数(ξ,κ,g,f,δ)以及出生概率bk和死亡概率dk＝1-bk，k为状态数的标号；出生概率和死亡概率服从二项分布，并设定迭代次数N以及假定不平衡迭代次数M，M取初始化迭代计数器h；另外，增加参数隐状态序列z，其初始化由下式确定：其中，yi,i＝1,2,....n为信号衰减序列，n为序列的长度，zi为参数隐状态序列z中标号为i的值，zi(1)则为其初始化值；2.4)进入第h次迭代的分裂合并过程：产生一个(0,1)均匀分布的随机数U1，若U1＜bk，则进入分裂过程；否则，进入合并过程；参数空间更新为其中分别为状态数K、状态概率ωj、标号为j的高斯过程的均值和方差、参数隐状态序列z完成第h次迭代的分裂合并过程后的更新，超参数β因为在分裂过程中没有更新，所以保留第h-1次迭代完成后更新值β(h)；2.5)进入第h次迭代的出生死亡过程：产生一个(0,1)均匀分布随机数U2，若U2＜bk，则进入出生过程，否则进入死亡过程；参数空间更新为其中分别为状态数K、状态概率ωj、标号为j的高斯过程的均值和方差、参数隐状态序列z完成第h次迭代的出生死亡过程后的更新，超参数β因为在合并过程中没有更新，所以保留第h-1次迭代完成后更新值β(h)；2.6)利用Gibbs采样更新参数空间得到第h次迭代结果K(h+1)、ωj(h+1)、μj(h+1)、σj(h+1)、z(h+1)和β(h+1)分别为状态数K、状态概率ωj、标号为j的高斯过程的均值和方差、参数隐状态序列z和超参数β完成第h次迭代后得到的最终更新值；2.7)h＝h+1，若h＜N，重复2.4)-2.6)；2.8)对于参数空间迭代结果x(1),x(2),....,x(N)，去除前M次不平衡迭代，得到x(M+1),x(M+2),....,x(N)，对其中参数K(M+1),K(M+2),...,K(N)作直方图，求出最大值kop，即为最佳状态数；2.9)在参数空间x(M+1),x(M+2),....,x(N)中，对于K(i)＝kop,i＝M+1,M+2,...,N的参数(ω(i),μ(i),σ(i))，求出数学期望(ω,μ,σ)，即为马尔科夫模型不同状态下的高斯参数；2.10)利用2.3)中的参数隐状态序列z的似然函数和2.9)中求出的马尔科夫模型不同状态下高斯参数，求出隐状态序列t＝1,2,...n，并进一步求出状态转移矩阵P为：P(i,j)＝Nij/Ni；其中，P(i,j)表示从状态i跳转到状态j的概率，Nij为从状态i转移到状态j的个数，Ni为处在状态i的个数,i,j＝1,2,...,n。2.根据权利要求1所述的再入动态等离子鞘套马尔科夫信道建模方法，其特征在于，步骤2.4)中，所述分裂过程具体如下：2.4.1.1)对于第h-1次迭代完成后最终得到的参数空间x(h)，从j＝1,2,...,K(h)随机选择一个状态js*，参数为以js1,js2表示分裂后状态的标号，分裂后状态数变为K'＝K(h)+1，K(h)和z(h)为第h-1次迭代完成后得到的状态数K和参数隐状态序列z的最终更新值，而分别为第h-1次迭代完成后状态js*的ωj、μj和σj的更新值；2.4.1.2)产生三个随机数u1,u2,u3，满足：u1～beta(2,2)u2～beta(2,2)u3～beta(1,1)；beta为贝塔分布，分裂后状态权重及相应状态高斯分布参数：2和分别表示状态js1和js2所对应的第h-1次迭代完成后得到μj的更新值；2.4.1.3)更新参数隐状态序列，将原本的样本点按照如下概率分配给js1和js2：Palloc是分配给js1的概率，q＝1-Palloc是分配给js2的概率，更新后的参数隐状态序列记为z'；2.4.1.4)计算接受概率split_accept＝min(1,split_A)，其中：

【专利技术属性】
技术研发人员：石磊，方水汛，李小平，姚博，杨敏，刘彦明，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西;61