具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

针对现有技术存在的问题，本发明提供了一种遥测信道非相干大规模SIMO 处理方法、系统及应用，下面结合附图对本发明作详细的描述。

如图1所示，本发明提供的遥测信道非相干大规模SIMO处理方法包括以下步骤：

S101：高超声速飞行器等离子体鞘套信道的建模；

S102：构建高超声速飞行器下行大规模SIMO系统模型；

S103：自适应功率调制；

S104：基于RRJ-MCMC算法的非相干检测。

本发明提供的遥测信道非相干大规模SIMO处理方法业内的普通技术人员还可以采用其他的步骤实施，图1的本发明提供的遥测信道非相干大规模SIMO 处理方法仅仅是一个具体实施例而已。

如图2所示，本发明提供的遥测信道非相干大规模SIMO处理系统包括：

建模模块1，用于实现高超声速飞行器等离子体鞘套信道的建模；

系统模型构建模块2，用于构建高超声速飞行器下行大规模SIMO系统模型；

功率调制模块3，用于自适应功率调制；

非相干检测模块4，用于基于RRJ-MCMC算法的非相干检测。

下面结合附图对本发明的技术方案作进一步的描述。

本发明可以实现等离子体鞘套信道下高可靠性和高效率的信息传输，是缓解黑障的有效手段。

如图3所示，该方案包含有如下步骤：

步骤一：高超声速飞行器等离子体鞘套信道的建模：

1)建立时变等离子体鞘套的分层电子密度模型：输入高超声速飞行器天线窗口附近非均匀等离子体鞘套厚度Z、等离子体鞘套分层总数N与等离子体鞘套各分层厚度d_m，m为等离子体鞘套分层序号m＝1，2，…，N，等离子体鞘套时变时间T，根据等离子体鞘套电子密度沿垂直飞行器天线窗口方向(设为z轴) 服从双高斯分布，第一高斯函数影响参数c₁和第二高斯函数影响参数c₂，峰值电子密度Ne_peak，峰值电子密度在z轴坐标为z_peak。Ne_peak沿时间方向服从正弦分布的分布规律，确定时变等离子体鞘套的电子密度分布函数，建立一个时变等离子体鞘套电子密度模型Ne_dyn(z,t)；

2)计算等离子体鞘套透射系数T_dyn(t,f₀)：输入通信信号的载波频率f₀，电子质量m_e，时变等离子体鞘套的电子碰撞频率v_en，等离子体鞘套各层的本征波阻抗z_m(m＝1，2，…，N)、真空的本征波阻抗z₀，利用时变等离子体鞘套电子密度模型Ne_dyn(z,t)，通过等效传输线法计算得到等离子体鞘套的时变透射系数 T_dyn(t,f₀)；将计算透射系数的幅度，并将结果赋值给等离子体鞘套信道

步骤二：构建高超声速飞行器大规模下行SIMO系统模型：

如图4所示，由高超声速飞行器端单天线和一个多天线接收站组成，其中接收站配备n(n→∞)根天线。该系统在一个时隙内的接收站端等效基带信号接收模型为：

y＝hx+v；

其中y是M×1维接收信号向量，x为发送信号，具体设计将在下节给出，h 是n×1维信道矩阵，v是n×1接收端复高斯白噪声。v的每一个元素均服从CSCG 分布，且均值为零，方差为σ²，即v_i～CN(0，σ²)，i＝1，2，…，n。

为高超声速飞行器下行链路等离子体鞘套-莱斯信道，其中T表示等离子体鞘套信道，由T_dyn(t,f₀)得到；表示空间莱斯信道衰落，即矩阵元素的包络服从莱斯分布，i＝1，…，n。假定噪声方差σ²、T和为未知信息，且已知等离子体鞘套信道的平均功率衰减q²。

步骤三：优化的功率调制：

发送符号其中第k个符号的功率，是编码本，L是的基数。假设p_k(k＝1，2，···，L)等概率发送，且满足平均功率约束依据等离子体鞘套信道的平均功率衰减q²和下行信道的渐进特性设计优化的

1)对于p_k，令p_k的解调区间其中d_l，k>0，d_r，k>0，a_k表示解调边界r(p_k)+d_r，k。以L＝6举例，待设计的星座如图5所示，其中为判决边界，有d_r,k+d_l,k+1＝p_k+1-p_k，k∈[L-1]。令由0均值独立随机变量U_k实现，U_k的矩母函数为令对于优化星座P^opt,有下面讲解如何得到优化t^opt，和优化星座 P^opt。

(3.2)利用算法1和算法2搜索到满足功率约束的最高t^opt，和相应的最优星座P^opt。算法1用于初始化和更新t值，算法1得出t值代入算法2；算法2搜索当前t值下的发送星座；算法1依据算法2所得星座是否满足功率约束条件，更新t值；循环算法1和2，直至t值满足设定精度10^-3，得到输出满足功率约束的最高t^opt，和相应的最优星座P^opt。算法1和2具体如下：

算法1：初始化t_l＝0，t_u＝∞，t代入算法2；若算法2星座满足功率约束条件，令t_l＝t；若算法2星座不满足功率约束条件，令t_u＝t；得到的t值继续代入算法2，再次得到t值，循环算法1，直至t值满足设定精度10^-3。

算法2：令p₁＝0，选择p₂的最小值，使得J(p₂)为t，且d_l， 此过程按顺序执行，直到找到p_L；检查是否满足功率约束条件；如果满足，保留这个星座，利用算法1增加t并重复这个过程；如果不满足，应该丢弃该星座，利用算法1减小t并重复这个过程。

步骤四：基于RRJ-MCMC算法的快速非相干检测：

为了方便接下来的信号检测及性能分析，首先计算接收向量y的平均功率：设计一种简化的可变逆跳蒙特卡洛(RRJ-MCMC)算法来估计的条件概率密度函数(CPDF)则最大似然(ML)解调器(或k 的估计值)为：

其中为似然函数/条件概率密度函数。

本发明为了保证马尔科夫模型能适应多种环境下等离子鞘套的动态性，可以将其建模为高斯过程，这样整个模型的概率密度函数可以用混合高斯过程c 来描述，其优势在于混合高斯过程可以用来拟合任意概率密度函数，实现数学形式上的统一，ω_k表示状态k在混合高斯过程中的状态概率，即第k个高斯过程在混合高斯过程中所占比率，表示均值和方差：

本发明考虑一个二维时变序列(c_t，z_t)描述接收统计量c的非平稳随机过程，其中(c_t)表示接收端的统计序列，(z_t)表示隐藏状态的序列，z_t∈{1，...，L}，L 表示马尔科夫模型中的状态数。需要估计的参数是λ＝{μ，σ，z，ω}其中μ＝(μ_k)，σ＝(σ_k)，ω＝(ω_k)。参数先验分布：ω～D(δ₁，δ₂，...，δ_L)，δ₁，δ₂，...，δ_L初始化为1。μ_k～N(ξ，κ^-1)，其中β～G(g，h)，α＝2，κ＝R^-2，ξ＝y_min＝R/2，g＝0.2， h＝10/R²，R＝c_max-c_min；G(g，h)表示具有参数g和h的伽马分布，D(δ₁，δ2，...，δ_L)代表Dirichlet分布)。

在初始时刻，状态数可以设置为L，剩余变量都基于它们的先验分布初始化；。最大迭代次数设置为N_iteration。迭代N时刻所有相关参数的状态空间 x＝(η，ω，z，μ，σ，c)在迭代N+1时刻更新为新的状态空间，即 x′＝(η′，ω′，z′，μ′，σ′，c′)。每个迭代可以总结为以下四个步骤：

1)更新ω：

其中为N时刻隐藏状态z_t＝k的个数；

2)更新(μ，σ)：

3)更新z：

4)更新β：

下面结合仿真对本发明的技术效果作详细的描述。

1、仿真1，等离子体鞘套信道幅相特性仿真分析

仿真条件：

时变等离子体鞘套的分层电子密度模型：高超声速飞行器天线窗口附近非均匀等离子体鞘套厚度Z＝0.06m、等离子体鞘套分层总数N＝100，与等离子体鞘套各分层厚度d_m＝0.06cm，m为等离子体鞘套分层序号m＝1，2，…，N，等离子体鞘套时变时间T＝0.1s，根据等离子体鞘套电子密度沿垂直飞行器天线窗口方向(设为z轴)服从双高斯分布，第一高斯函数影响参数c₁＝2.56×10³和第二高斯函数影响参数c₂＝3.55×10³，峰值电子密度Ne_peak＝1×10¹⁸m^-3～1×10¹⁹m^-3，峰值电子密度在z轴坐标为z_peak＝0.036m。确定时变等离子体鞘套的电子密度分布函数，建立了一个时变等离子体鞘套电子密度模型Ne_dyn(z,t)。

计算等离子体鞘套透射系数T_dyn(t,f₀)：输入通信信号的载波频率f₀＝28GHz，载波角频率ω＝2πf₀电子质量m_e＝9.10953^-31kg，时变等离子体鞘套的电子碰撞频率 v_en＝1GHz，计算等离子体频率与复介电常数给出等离子体鞘套各层的本征波阻抗 真空的本征波阻抗利用时变等离子体鞘套电子密度模型Ne_dyn(z,t)，通过等效传输线法计算得到等离子体鞘套的时变透射系数T_dyn(t,f₀)。

2、仿真结果和分析：

在图6中，描绘了等离子体鞘套电子密度的时变图。图中看出等离子体鞘套电子密度受到流体扰动影响的出现了小尺度抖动。

在图7中，描绘了载波在28GHz时电波穿过等离子体的透射系数的幅度结果，电波穿过等离子体的透射系数的幅度受到电子密度扰动等因素的影响的出现了高动态扰动。透射系数的幅度结果将赋值给等离子体鞘套信道T。

3、仿真2，仿真分析在高超声速飞行器下行SIMO信道优化的功率调制和非相干检测方案的误码率和信道容量随信噪比或接收天线个数变化的变化情况。

仿真条件：

等离子体鞘套参数与电磁波时变透射系数采用仿真1中的参数与结果。莱斯斯因子K＝{0，1，∞}。

4、仿真结果和分析：

图8(a)和图8(b)中描绘了在信噪比γ＝10dB，莱斯因子K＝{0，1，∞}，电子密度峰值Ne_peak分别为2×10¹⁸m^-3和6×10¹⁸m^-3，误码率随接收站天线数的变化情况；图8(c)中描绘了在接收站天线数100，K＝{0，1，∞}，Ne_peak为6×10¹⁸m^-3时，误码率随信噪比的变化情况。可以看出随着接收站天线数的增加，系统误码率呈现指数型下降；随着信噪比的增加，系统误码率快速下降；莱斯因子K 越大，系统检测性能越好。

在图8(a)-图8(c)仿真条件下，图9(a)-图9(c)描绘相应的信道容量性能。仿真发现系统信道容量损耗和误码率随信噪比或接收天线个数变化有相同的变化趋势。信道容量随着接收站天线数的增加，系统误码率呈现指数型下降；随着信噪比的增加，信道容量快速下降；莱斯因子K越大，系统信道容量性能越好。

本发明公开了一种高超声速飞行器下行遥测信道下非相干大规模SIMO系统。方案包含两部分，一是超声速飞行器下行SIMO信道的功率调制方案，二是针对这种调制方案设计相应的自适应非相干信号检测方法。方案首先建立高超声速飞行器等离子体鞘套信道模型，并建立高超声速飞行器大规模下行SIMO系统模型；在系统模型基础上提出一种优化功率调制方案；针对这种功率调制方案，提出了一种RRJ-MCMC算法，并在算法基础上提出自适应非相干检测。仿真表明本发明方法可通过增加接收站天线的数量和提高莱斯因子和信噪比来提升通信的质量。本发明公开的高超声速飞行器下行SIMO信道下的自适应功率调制方案以及非相干信号检测方案可显著提高等离子体信道下通信的误码率性能，改善通信质量、缓解通信黑障。

应当注意，本发明的实施方式可以通过硬件、软件或者软件和硬件的结合来实现。硬件部分可以利用专用逻辑来实现；软件部分可以存储在存储器中，由适当的指令执行系统，例如微处理器或者专用设计硬件来执行。本领域的普通技术人员可以理解上述的设备和方法可以使用计算机可执行指令和/或包含在处理器控制代码中来实现，例如在诸如磁盘、CD或DVD-ROM的载体介质、诸如只读存储器(固件)的可编程的存储器或者诸如光学或电子信号载体的数据载体上提供了这样的代码。本发明的设备及其模块可以由诸如超大规模集成电路或门阵列、诸如逻辑芯片、晶体管等的半导体、或者诸如现场可编程门阵列、可编程逻辑设备等的可编程硬件设备的硬件电路实现，也可以用由各种类型的处理器执行的软件实现，也可以由上述硬件电路和软件的结合例如固件来实现。

以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，都应涵盖在本发明的保护范围之内。