GPS接收机中晶振误差的模拟方法

寇艳红; 张其善

摘要: GPS(Global Positioning System)接收机中参考频率源的性能对接收机整体性能有着诸多影响。该文对在GPS接收机中使用最多的晶振的误差特性和性能影响进行了分析，在此基础上提出了一种模拟晶振频率误差的方法，即首先确定晶振在给定取样时间间隔上各主要随机误差成分的Allan(阿仑)方差参数，由此反演产生相应的各项晶振频率误差模拟序列并进行叠加。文中给出了计算机数值模拟的方法，并对仿真实验结果进行了分析，从而为GPS接收机设计、性能评估和测试平台如GPS信号模拟器的研制提供有关的理论支持和设计方法。

关键词:

GPS; 频率稳定度; 晶振; Monte Carlo模拟; Allan方差

Abstract: The frequency source in Global Positioning System (GPS) receiver has a great many influence on the whole performance of the receiver. Based on the analysis of error characteristics and performance effect of the mostly used crystal oscillator, a method to simulate the frequency error of crystal oscillator is proposed. The major idea is to determine the Allan variance parameter of each main stochastic error component at the given sampling interval at first, then to deduct a set of corresponding frequency error series and superpose them together. The computer numerical simulation approach and the experimental results are given. This research provides related both theoretical reference and practical design for the development of an integrated design, test and evaluation platform such as a GPS signal simulator.

1. 引言

轨道角动量(Orbital Angular Momentum, OAM)作为量子力学及经典力学中的基本物理量，是独立于光的波长、偏振、相位、振幅外的一个空间维度^[1]。近年来，将OAM作为新的通信自由度应用到通信领域得到了人们广泛的关注。

1992年，Allen等人^[2]证实拉盖尔高斯(Laguerre-Gaussian, LG)光束除了拥有自旋角动量(Spin Angular Momentum, SAM)，还携带OAM，且不同模式的OAM光束空间中彼此正交。2004年，Gibson等人^[3]将携带OAM的波束作为通信载波，首次实现单模态OAM通信。2011年，Wang等人^[4]首次完成了多模态OAM复用光通信，实现了4个模态OAM复用通信传输实验。2013年，Huang等人^[5]完成了32路独立OAM复用通信传输，达到了2.56 Tbit/s的传输速率。2014年，Wang等人^[6]将OAM模式与光的偏振结合，达到了230 bit/(s·Hz)的频谱效率。

上述研究表明，OAM应用于光无线通信领域可以提升频谱利用率、信息传输速率。然而现实中大气湍流、多径效应、孔径失配等因素均降低OAM通信系统的性能^[7]。为了推进光无线OAM通信系统在实际场景中的应用，部分学者对大气湍流场景下的OAM通信系统性能展开了广泛的研究。2015年，邹丽等人^[8,9]利用GS(Gerchberg-Saxton)算法对OAM波束进行相位校正，降低2～3个数量级的系统误码率。2016年，Zou等人^[10]将OAM模式串扰等效为多用户间的相关信息，借鉴多用户检测的解相关技术，有效地降低了大气湍流对OAM通信系统的影响。同年，Zhao等人^[11]联合信道编码技术与波前校正技术，提升了OAM系统在大气湍流下的性能。2017年，Zou等人^[12]基于分集思想，发送端采用多孔径进行数据传输，接收端采用等增益合并技术，提升了OAM在大气湍流中的性能。同年，Zhang等人^[13]基于空时编码准则，建立大气湍流场景下的OAM复用通信系统，对比分析垂直分层空时码(Vertical Bell LAyered Space Time, VBLAST)及空时分组码(Space Time Block Code, STBC)方案对OAM系统性能的提升效果。2018年，Yousif等人^[14]将MIMO均衡的思想应用于OAM通信系统，降低数据间的串扰，系统性能得到明显提升。2018年，Wang的等人^[15]将恒模均衡算法与模态选择策略相结合，降低模态串扰对OAM复用系统的影响。2019年，Dedo等人^[16]通过改变随机相位屏的个数、湍流强度等条件，分析GS算法下OAM波束修复的极限。同年，Amhoud等人^[17]根据OAM不同模态间的串扰差异性，采取OAM模态选择结合STBC编码的方案，提升了大气湍流下OAM通信性能。2020年，Song等人^[18]根据OAM串扰信道信息，对OAM光束进行预处理和后处理，在一定程度上降低了OAM模态间的串扰，提升了OAM通信系统在大气湍流中的性能。

综上所述，光无线OAM通信系统主要从光学自适应及数字信号处理两个角度应对大气湍流。光学自适应虽然有良好的波形修复效果，但采用该方式时需要考虑设备的复杂性及成本；数字信号处理多数是借鉴多输入多输出(Multiple Input Multiple Output, MIMO)通信中的均衡、分集及一些改进的技术，在应对大气湍流时同样有不俗的表现。然而，通信系统中的数字信号处理包含了预编码、均衡及信号检测等多个方面技术，以往研究很少从信号检测的角度分析复杂场景下的光无线OAM复用系统。另外，目前多数光无线OAM复用系统默认收发端孔径适配，然而实际中激光通信很容易出现孔径失配。针对于以上两个方面，本文建立了大气湍流、孔径失配场景下，基于VBLAST编码准则的OAM(VBLAST-OAM)复用通信系统，从信号检测的角度对所建系统进行分析。主要的研究工作如下：

(1) 建立了大气湍流、孔径失配场景下基于VBLAST准则的OAM复用通信系统，并进行相应的理论推导。

(2) 基于(1)中所建系统，首先分析了大气湍流、孔径失配场景下OAM的模态串扰的特征；之后，将排序连续干扰消除(Ordered Successive Interference Cancellation, OSIC)、马尔科夫随机场置信度传播算法(Markov Random Field Belief Propagation, MRF-BP)应用于VBLAST-OAM复用通信系统；此外，根据OAM模态间的串扰差异性，给出了基于轨道角动量的排序连续干扰消除算法(Orbital Angular Momentum with Ordered Successive Interference Cancellation, OAM-OSIC)。最后，对比了文中所给算法的性能，分析各个算法的优缺点。

2. OAM空间传播理论及复杂场景下OAM复用系统理论推导

光束在大气湍流中传播时，会被附加随机相位因子，导致接收到的波形出现失真。除了大气湍流外，孔径失配也会造成波形失真，孔径失配通常包含孔径偏离、角向倾斜两种情况，前者多数是设备抖动引起的；后者是光束长距离传播后，波束传播轴与接收系统参考轴的非重合现象。本文主要分析OAM的角向倾斜情形，图1是大气湍流、孔径失配场景下的VBLAST-OAM复用系统模型。首先通过串并变换操作将串行数据转变为N路并行数据，并将其送入信号调制器，随后将调制过的数据送入空时编码器，将编码后的数据加载到N列高斯光束，利用OAM转换器将高斯光束转换成一系列不同模态的OAM光束，最后利用单孔径望远镜将复用的波束发射到空间。假设接收孔径将复用的OAM波束接收后，首先利用OAM解复用器将复用的OAM光束分别检测出来，之后采用OAM转换器将OAM光束转为高斯波束，最后对每一列高斯波束上面的信息进行译码和解调。图1中，接收示孔径位置的 $\eta$ 代表光轴的偏离方位， $\gamma$ 代表角向倾斜的角度。文献[19,20]表明，波束的横向矢量远远小于纵向矢量的条件下，当出现孔径失配中的角向倾斜，相当于光束穿过一个相位楔，表示为

图 1 大气湍流、孔径失配场景下的VBLAST-OAM复用系统

下载: 全尺寸图片幻灯片

$U(r,\phi ,z) = E(r,\phi ,z)\exp ({\rm{j}}k\sin (\gamma )r\cos (\phi - \eta ))$

(1)

其中， $E(r,\phi ,z)$ 表示在空间传输后的波束， $U(r,\phi ,z)$ 表示接收孔径探测到的波束。

2.1 轨道角动量及光束在大气湍流中的传播理论

2.1.1 轨道角动量基本理论

根据经典电动力学理论可知，电磁场具备能量与动量且满足电磁守恒。其中，动量主要包括线动量与角动量，线动量通常与平动或者力的作用有关；角动量可以分为SAM和OAM，SAM与电磁波的极化相联系，OAM与电磁波的相位波前分布有关。粒子态OAM沿着传播方向的量子算符可以写成 ${\widehat L_z} = - {\rm{j}}\hbar {\partial }/{{\partial \phi }}$ ，本征值方程为 ${\widehat L_z}\left| l \right\rangle = L\left| l \right\rangle$ ，其中 $\left| l \right\rangle$ 在基坐标方位角 $\phi$ 表象下，可以写作 $\left\langle {\phi |l} \right\rangle =$ ${1 / {\sqrt {2\pi } \exp ({\rm{j}}l\phi )}}$ ，表示具有l重螺旋相位^[21]。

LG光束在旁轴近似条件下的亥姆霍兹方程表示为^[22]

$\frac{1}{\rho }\frac{\partial }{{\partial r}}\left(\rho \frac{{\partial E}}{{\partial r}}\right) + \frac{1}{{{\rho ^2}}}\frac{{{\partial ^2}E}}{{\partial {\phi ^2}}} + 2{\rm{j}}k\frac{{\partial E}}{{\partial z}} = 0$

(2)

通过对式(2)求解，可以得到沿着z轴传播的光场

$\;\begin{split} E_p^l(r,\phi ,z) =& \sqrt {\frac{{2p!}}{{\pi (p + \left| l \right|)!}}} \times \frac{1}{{\omega (z)}} \times L_p^l\left[ {\frac{{2{r^2}}}{{{\omega ^2}(z)}}} \right] \\ & \times {\left[ {\frac{{\sqrt 2 r}}{{\omega (z)}}} \right]^{\left| l \right|}} \times \exp \left[ { - \frac{{{r^2}}}{{{\omega ^2}(z)}}} \right] \\ & \times {\rm{exp}}\left[ { - \frac{{{\rm{j}}k{r^2}z}}{{2({z^2} + Z_R^2)}}} \right] \\ & \times \exp \left[ { - {\rm{j}}(2p + \left| l \right| + 1){{{\rm{arctan}} }^{ - 1}}\left(\frac{z}{{{Z_R}}}\right)} \right] \\ & \times \exp ({\rm{j}}l\phi ) \\[-10pt] \end{split}$

(3)

$\left.\begin{aligned} &\omega (z) = {\omega _0}\sqrt {1 + {{\left( {\frac{z}{{{Z_R}}}} \right)}^2}} \\ &{Z_R} = \frac{{\pi \omega _0^2}}{\lambda } \end{aligned}\right\}$

(4)

其中，r是径向距离， $\phi$ 为方位角，z为传输距离， ${\omega _0}$ 代表束腰半径， $\omega (z)$ 是距离z处的光斑半径， ${Z_R}$ 是瑞利距离， $L_p^l$ 是广义拉盖尔多项式， $\exp \left[ { - {\rm{j}}(2p + \left| l \right| + 1){{{\rm{arctan}} }^{ - 1}}\left( {\frac{z}{{{Z_R}}}} \right)} \right] \times \exp ({\rm{j}}l\phi )$ 代表LG光束相位， $\exp ({\rm{j}}l\phi )$ 体现了LG的旋转相位特性。

2.1.2 OAM在大气湍流中的传播理论

为了能够将OAM光束更好地应用于光无线通信领域，OAM光束在大气湍流中的传播特性分析至关重要。实验室条件下，常采用功率谱反演法及Zernike多项式法等数值仿真方法产生随机相位屏来模拟大气湍流。本文基于修正的冯卡曼折射率功率谱，采用功率谱反演法模拟大气湍流，折射率功率谱模型表示为^[23]

$\begin{split} {\varPhi _n}(k) =& 0.033C_n^2{k^{ - 11/3}}\exp \left( - \frac{{{k^2}}}{{k_m^2}}\right){({k^2} + k_0^2)^{ - 11/6}},\\ & 0 \le k < \infty\\[-10pt] \end{split}$

(5)

其中， $C_n^2$ 表示折射率结构系数，可以衡量湍流的强度， $k = \sqrt {k_x^2 + k_y^2}$ , ${k_x}$ 与 ${k_y}$ 代表空间频率， ${k_0} =$ $2\pi /{L_0}$ , ${k_m} = 5.92/{l_0}$ ,l₀与L₀分别代表湍流涡旋的内外尺寸。根据折射率功率谱可以得到大气湍流相位功率谱

${\varPhi _\varphi }(k) = 2\pi {k^2}\Delta z{\varPhi _n}(k)$

(6)

其中， $\Delta z$ 代表两个随机相位屏之间的距离， ${\varPhi _n}(k)$ 是大气湍流折射率功率谱。功率谱反演法的基本思想：首先构造一个N维的0均值、单位方差的复高斯随机矩阵 $h({k_x},{k_y})$ ，随后用大气湍流相位功率谱对 $h({k_x},{k_y})$ 滤波，最后对滤波后的函数做傅里叶逆变换，可得到大气随机相位

$\begin{split} \varphi (x,y) =& \sum\limits_{{k_x}} {\sum\limits_{{k_y}} {{\boldsymbol{h}}({k_x},{k_y})} \sqrt {{\varPhi _\varphi }(k)} } \\ & \times \exp ({\rm{j}}({k_x}x + {k_y}y))\Delta {k_x}\Delta {k_y} \end{split}$

(7)

其中， $\Delta {k_x}$ , $\Delta {k_y}$ 代表在空间频率间隔。LG光束 $E(r,\phi ,z)$ 在大气湍流中的传播过程，可以等效为光束在真空与随机相位屏之间的迭代，具体过程如表1所示。

表 1 OAM在大气湍流中的传播过程

初始化：源端光场 $E(r,\phi ,z)$ ，空间传输函数 $H$ ，相位屏个数N，相位 $\varphi (x,y)$
(1) ${\rm{For }}j{\rm{ = 1:}}N$
(2) 傅里叶变换： $U(K,z) = {\rm{FFT}}\left( {E(r,\phi ,z)} \right)$
(3) 真空传播： ${U'}(K,z) = {\rm{IFFT}}\left( {U(K,z) \times H} \right)$
(4) 穿过随机相位屏： $E(r,\phi ,z) = {U'}(K,z) \times \exp \left( {{\rm{j}} \times \varphi (x,y)} \right)$
(5) End For

下载: 导出CSV

| 显示表格

2.2 基于VBLAST-OAM复用通信系统理论分析

本文采用VBLAST方案对发送数据进行编码，将编码后的数据加载到N列OAM光束，那么一个周期内的数据信息表示为

${{\boldsymbol{X}}_{{\rm{VBLAST}}}} = \left[ {\begin{array}{*{20}{l}} {x_{1,1}^{{l_1}}}&{x_{1,2}^{{l_1}}}& \cdots &{x_{1,M}^{{l_1}}}\\ {x_{2,1}^{{l_2}}}&{x_{2,2}^{{l_2}}}& \cdots &{x_{2,M}^{{l_2}}}\\ \quad \vdots &\quad \vdots & \ddots &\quad \vdots \\ {x_{N,1}^{{l_N}}}&{x_{N,2}^{{l_N}}}& \cdots &{x_{N,M}^{{l_N}}} \end{array}} \right]$

(8)

其中， $x_{N,M}^{{l_N}}$ 表示经过VBLAST编码后的数据，上标 ${l_N}$ 表示第N列OAM光束，下标(N,M)代表第M时刻在第N列OAM光束上加载的数据。为讨论方便，只取某一个时刻进行讨论。N列OAM光束复用表示为

${U_{{\rm{MUX}}}} = \sum\limits_{n = 1}^N {x_{n,T}^{{l_n}} \times {E_n}(r,\phi ,z)}$

(9)

接收孔径得到的光束表示为

$U_{{\rm{MUX}}}' = \sum\limits_{n = 1}^N {x_{n,T}^{{l_n}} \times E_n'} (r,\phi ,z)$

(10)

接收端通常利用OAM的正交特性对复用的光束进行分离。当OAM的正交性受到破坏后，此时从第k列OAM上得到的信息为

$\begin{split} {y_k} =& \langle U_{{\rm{MUX}}}'(r,\phi ,z),{E_k}(r,\phi ,z)\rangle {\rm{ }} \\ =& {\rm{ }}\iint {U_{{\rm{MUX}}}'(r,\phi ,z) \times E_k^*(r,\phi ,z)r{\rm{d}}r{\rm{d}}\phi } \\ =& \iint {\sum\limits_{n = 1}^N {x_{n,T}^{{l_n}} \times E_n'} (r,\phi ,z) \times E_k^*(r,\phi ,z)r{\rm{d}}r{\rm{d}}\phi } \\ =& x_{k,T}^{{l_k}} \times \iint {E_k'(r,\phi ,z)} \times E_k^*(r,\phi ,z)r{\rm{d}}rd\phi \\ & + \sum\limits_{n = 1,n \ne k}^N {x_{n,T}^{{l_n}}\iint {E_n'}} (r,\phi ,z) \times E_k^*(r,\phi ,z)r{\rm{d}}r{\rm{d}}\phi \end{split}$

(11)

其中， ${( \cdot )^*}$ 表示共轭。根据式(11)可以看出，当复用的OAM光束受外界影响后，会出现信号间串扰。式(11)推导最后一步，首项代表模态k的期望信号，第2项代表除模态k外，其他模态串扰到模态k上的信号。将式(11)写成矩阵形式

${{\boldsymbol{Y}}_{N \times 1}} = {{\boldsymbol{A}}_{N \times N}} \times {{\boldsymbol{X}}_{N \times 1}}$

(12)

其中，Y代表接收到的信号，A代表OAM模态串扰矩阵，X是发送信号。将式(12)展开

$\left[ {\begin{array}{*{20}{l}} {{y_1}}\\ {{y_2}}\\ \ \vdots \\ {{y_N}} \end{array}} \right] = \left[ {\begin{array}{*{20}{l}} {{a_{1,1}}}&{{a_{1,2}}}& \cdots &{{a_{1,N}}}\\ {{a_{2,1}}}&{{a_{2,2}}}& \cdots &{{a_{2,N}}}\\ \quad\vdots & \quad\vdots & \ddots & \quad\vdots \\ {{a_{N,1}}}&{{a_{N,2}}}& \cdots &{{a_{N,N}}} \end{array}} \right] \times \left[ {\begin{array}{*{20}{l}} {x_{1,T}^{{l_1}}}\\ {x_{2,T}^{{l_2}}}\\ \quad \vdots \\ {x_{N,T}^{{l_N}}} \end{array}} \right]$

(13)

其中，矩阵A的对角线元素 ${a_{N,N}}$ 对应式(11)推导中的第1项，非对角元素对应式(11)推导中的第2项。

根据文献[24]可知，任意模态的OAM光束均可以展开成谐波函数 $\exp (im\phi )$ 的线性叠加

${E_n}(r,\phi ,z) = \frac{1}{{\sqrt {2\pi } }}\sum\limits_{m = - \infty }^\infty {{a_{m,n}}(r,z)\exp ({\rm{j}}m\phi )}$

(14)

利用傅里叶变换可以求得 ${a_{m,n}}$

$\left.\begin{split} & {a_{m,n}}(r,z) = \frac{1}{{\sqrt {2\pi } }}\int_0^{2\pi } {{E_n}(r,\phi ,z)} \exp ( - {\rm{j}}m\phi ) \\ & {C_{m,n}} = \int_0^\infty {{{\left| {{a_{m,n}}(r,z)} \right|}^2}} r{\rm{d}}r \\[-15pt] \end{split}\right\}$

(15)

其中， ${C_{m,n}}$ 代表模态为n的OAM在模态m上的谐波功率，模态n的OAM的纯度表示为

${p_{n,n}} = \frac{{{C_{n,n}}}}{{\displaystyle\sum\limits_{q = - \infty }^\infty {{C_{q,n}}} }}$

(16)

3. 基于信号检测算法的VBLAST-OAM复用通信系统

3.1 基于OSIC的VBLAST-OAM复用系统

对基于空时分层编码方案的通信系统进行数据恢复时，通常采用分层处理方式。具体来说，就是将所要处理的数据层作为期望数据，那么除了期望层外的其他层数据均视为干扰，该方式需要生成加权矢量对干扰进行抑制，利用加权向量可以将所需数据提取出来。当数据层间的差异性较大时，随机地对数据层进行检测，最终的检测效果可能会不理想。因此有人提出了连续干扰消除的方法，这种方法是每检测完一层就将该层数据从接收数据中消除，避免了高信噪比数据层对低信噪比数据层造成的干扰。

根据文献[17]可知，同等大气湍流强度下，不同模态的OAM串扰具有明显的差异性。OAM的模态绝对值越小，主模态功率就保留得越多。那么串扰到相邻模态的功率就越少；反之模态绝对值越大，主功率就越低。因此，可以认为低模态数据层信噪比就越大。本文基于不同模态OAM间的能量串扰特点，优化了传统OSIC算法中的数据层排序步骤，给出一种低阶模态数据层优先检测的OSIC算法。

抑制干扰的加权矢量可以基于迫零、最小均方误差准则获取，因此将OSIC算法分为基于迫零准则的排序连续干扰消除(Zero Forcing criterion with Ordered Successive Interference Cancellation, ZF-OSIC)和基于最小均方误差准则的排序连续干扰消除(Minimum Mean Squared Error criterion with Ordered Successive Interference Cancellation, MMSE-OSIC)；而本文OAM-OSIC算法同样可以分为基于迫零准则的排序连续干扰消除(Orbital Angular Momentum with Zero Forcing criterion with Ordered Successive Interference Cancellation, OAM-ZF-OSIC)和基于最小均方误差准则排序连续干扰消除(Orbital Angular Momentum with Minimum Mean Squared Error criterion with Ordered Successive Interference Cancellation, OAM-MMSE-OSIC)算法。

以ZF-OSIC为例，假设接收到的信号r，产生一个对第k层对应的加权向量 ${{\boldsymbol{W}}_k}$ ，将 ${{\boldsymbol{W}}_k}$ 与r相乘

${{\boldsymbol{W}}_k}{\boldsymbol{r}} = {{\boldsymbol{W}}_k}{{\boldsymbol{X}}_k} + {{\boldsymbol{W}}_k}{\boldsymbol{n}}$

(17)

其中， ${X_k}$ 对应第k层的发送数据，n是高斯白噪声。此时第k层检测后信噪比表示为

${\rm{SN}}{{\rm{R}}_k} = \frac{{\left\langle {{{\left| {{{\boldsymbol{X}}_k}} \right|}^2}} \right\rangle }}{{{\sigma ^2}{{\left\| {{{\boldsymbol{W}}_k}} \right\|}^2}}}$

(18)

其中，分子表示第k层符号的发射功率，分母表示噪声功率与加权矢量范数平方的乘积。从式(18)可以看出，每一层检测信噪比与行范数成反比，因此可以从范数小的数据层开始检测，降低大信噪比数据层对其他数据层的干扰。虽然MMSE-OSIC算法通常根据信干噪比进行顺序检测，为了便于分析，本文采用行范数排序规则。

表2给出了OSIC算法在VBLAST-OAM系统中的实施步骤。

表 2 OSIC算法

初始化：复用的OAM数目 ${ {{N} }_{{t} } }$ ，串扰信道 $H$ ，接收到的信号r，噪声方差 ${\sigma ^2}$
(1) ${\rm{For } }j{\rm{ = 1:} }{ {{N} }_{{t} } }$
(2) 加权矢量W: 基于ZF/MMSE准则： $W = {\rm{pinv(} }H{\rm{)/pinv(} }{H^{\rm{H} } }{\rm{ + } }{\sigma ^2} \times {I_{ {N_{\rm{t} } } - j + 1} }{\rm{)} } \times {H^{\rm{H} } }$
(3) 首先对W的每一行求范数，并对范数排序，选取最小范数行k (OSIC)
(4) 省略步骤(3) (OAM-OSIC)
(5) $y(k) = W(k,:) \times r$ 判决统计量(优化前)
(6) $y(k) = W(1,:) \times r$ 判决统计量(优化后)
(7) 根据数据判决得到x(k)
(8) $r = r{\rm{ - }}x(k) \times H$ 消除前一次检测的数据
(9) $H(:,k) = [\;]$ 将信道矩阵的第k列清除(优化前)
$H(:,1) = [\;]$ 将信道矩阵的第1列清除(优化后)
(10) 重复步骤(2)，更新W
(11) ${\rm{End\; For}}$

下载: 导出CSV

| 显示表格

3.2 基于MRF-BP的VBLAST-OAM复用系统

图2是基于马尔科夫随机场的概率图模型。可以看出，马尔科夫随机场中存在许多相互连接的节点，节点又可以分为观察节点和隐藏节点。如果将上述模型与OAM复用系统相对应，那么图中任意观察节点对应从任意模态的OAM上解调信息 ${y_i}$ ，隐藏节点对应任意模态OAM上加载的信息 ${x_i}$ 。此外，观察节点与隐藏节点的关系可以用相容函数 ${\phi _i}({x_i})$ 描述，同样将相邻隐藏节点的关系用势函数 ${\psi _{i,j}}({x_i},{x_j})$ 描述。MRF-BP算法就是利用置信度传播(Belief Propagation, BP)算法对马尔科夫随机场中的各个节点的状态信息不断更新，当算法达到收敛时，即可估计出隐藏节点的状态信息概率。光无线OAM复用系统中，场内隐藏节点的联合概率^[25]

图 2 基于MRF的BP信号检测模型

下载: 全尺寸图片幻灯片

$p(x) = \frac{1}{Z}\prod\limits_i {{\phi _i}({x_i},{y_i})\prod\limits_{i,j} {{\psi _{i,j}}({x_i},{x_j})} }$

(19)

其中，Z称配分函数，保证分布函数标准化。利用BP算法可以将实现信息在场内节点相互传递，从而更新MRF中各个节点状态。将节点j传递到节点i的信息表示为 ${m_{j,i}}({x_i})$ ，节点i的置信度用 ${b_i}({x_i})$ ，那么 ${b_i}({x_i})$ 可以表示为

${b_i}({x_i}) \propto {\phi _i}({x_i})\prod\limits_{j \in N(i)} {{m_{j,i}}({x_i})}$

(20)

$N(i)$ 表示节点i的相邻节点集合，从节点j到节点i的传递信息近似表示

${m_{j,i}}({x_i}) \propto \sum\limits_{{x_j}} {{\phi _j}({x_j}){\psi _{i,j}}({x_i},{x_j})} \prod\limits_{k \in N(j)\backslash i} {{m_{k,j}}({x_j})}$

(21)

利用势函数及相容函数，可以将信息的传递与更新构建成一个迭代过程。由于多模态OAM在复杂环境传播后会出现模态间的相互串扰，因此基于OAM的马尔科夫随机场是全连接图，此时算法的收敛性较差，最终得到的信息存在较大误差。本文采用信息阻尼的方式加速算法收敛，具体操作就是把上一次迭代得到的估计值与新的估计值各自分配一个加权系数，每次得到的估计值是新旧估计值的加权组合。 $m_{i,j}^t({x_j})$ 是t次迭代时，节点i到j的新信息，那么更新的信息为^[26]

$\begin{split} & \overline m _{i,j}^t({x_j}) \propto \sum\limits_{{x_i}} {{\phi _i}({x_i}){\psi _{i,j}}({x_i},{x_j})} \left\{ {\prod\limits_{k \in N(i)\backslash j} {m_{k,i}^{t - 1}({x_i})} } \right\} \\ & m_{i,j}^t({x_j}) = {a_m}m_{i,j}^{t - 1}({x_j}) + \left\{ {(1 - {a_m})\overline m _{i,j}^t({x_j})} \right\} \end{split}$

(22)

其中， ${\alpha _m} \in \left[ {0,1} \right]$ 是信息阻尼因子，取0.2。基于马尔科夫随机场图模型的VBLAST-OAM系统信息后验概率表示为

$\begin{split} p(x|y,{\boldsymbol{H}}) =& \left( {\prod\limits_{i \ne j} {\exp ( - {x_i}\Re ({R_{i,j}}){x_j})} } \right)\\ & \times \left( {\prod\limits_i {\exp ({x_i}\Re ({z_i}) + \ln p({x_i}))} } \right) \end{split}$

(23)

其中， $R = (1/{\delta ^2}){\boldsymbol{H}^{\rm{H}}}\boldsymbol{H}$ , $z = (1/{\delta ^2}){\boldsymbol{H}^{\rm{H}}}y$ 。将式(19)与式(23)比较，可以得到两者的关系

${\psi _{i,j}}({x_i},{x_j}) = \exp ( - {x_i}\Re ({R_{i,j}}){x_j})$

(24)

${\phi _i}({x_i}) = \exp ({x_i}\Re ({z_i}) + \ln p({x_i}))\;\;$

(25)

将MRF-BP算法应用于VBLAST-OAM系统，详细的实施步骤如表3所示。

表 3 MRF-BP算法

初始化： $m_{i,j}^0 = b_i^0$ , $z,R, p({x_i} = 1) = p({x_i} = - 1), \forall i,j \in$ $(1,2,\cdots,N)$ , M是信息迭代次数
(1) ${\rm{For }}\;i{\rm{ = 1:}}N$ 势函数
(2) ${\rm{For }}\;j{\rm{ = 1:}}N$ ${\rm{ }}i \ne j{\rm{ }}$
(3) 根据式(24)计算 ${\psi _{i,j}}$ ；
(4) ${\rm{End\; For}}$
(5) ${\rm{End\; For}}$
(6) ${\rm{For }}\;i{\rm{ = 1:}}N$ 相容函数
(7) 根据式(25)计算 ${\phi _i}$ ；
(8) ${\rm{End\; For}}$
(9) ${\rm{For} }\; t = 1:M$ 迭代更新
(10) ${\rm{For }}\;i{\rm{ = 1:}}N$
(11) ${\rm{For\; }}j{\rm{ = 1:}}N$ ${\rm{ }}i \ne j{\rm{ }}$
(12) 第t次迭代得到更新的信息 $\overline m _{i,j}^t$ ，利用信息阻尼方式得到新信息 $m_{i,j}^t$ ；
(13) ${\rm{End\; For}}$
(14) ${\rm{End\; For}}$
(15) ${\rm{End\; For}}$
(16) ${\rm{For }}\;i{\rm{ = 1:}}N$ 置信度计算
(17) 根据式(20)计算置信度 ${b_i}$ ；
(18) ${\rm{End\; For}}$

下载: 导出CSV

| 显示表格

4. 仿真参数设置及结果分析

仿真参数：光束波长 $\lambda = 1.55 \times {10^{ - 6}}\;{\rm{ m}}$ ，传播距离 $z = 1\;{\rm{ km}}$ , ${l_0} = 0.001\;{\rm{ m}}$ , ${L_0} = 50\;{\rm{ m}}$ ，栅格数为300，调制方式为二进制相移键控(Binary Phase Shift Keying, BPSK)。

4.1 大气湍流、孔径失配场景下的OAM模态谱分析

本节首先对复杂场景下OAM的模态纯度进行分析，之后分析多模态OAM之间的串扰特征。

图3是大气湍流、孔径失配条件下的OAM模态纯度3维图。仿真选择模态为3的OAM，从图中可以看出，随着倾斜角、湍流强度的增大，接收端检测到的OAM模态纯度逐渐降低。表明OAM模态纯度会受到多种外界条件的影响，大气湍流、孔径失配均会对OAM的模态纯度带来一定程度的影响。

图 3 不同大气湍流、倾斜角条件下的模态纯度

下载: 全尺寸图片幻灯片

为了直观地描述式(11)中的模态串扰，图4给出了单模态OAM及多模态OAM的能量串扰分布图谱。仿真参数：折射率结构常数 $C_n^2 = 1 \times {10^{ - 14}}\;{{\rm{m}}^{ - 2/3}}$ ，角向倾斜角 $\gamma = 5 \times {10^{ - 5}}\;{\rm{ rad}}$ 。图4(a)是l = 1的OAM模态串扰分布图谱，虽然仅发送了l = 1的OAM，但检测时，却得到了除l = 1外其他模态的能量。图4(b)是与图4(a)同等仿真条件下，l = 3的OAM的模态串扰图谱，两个模态的OAM具有类似的串扰规律，唯一不同的是l = 3的OAM主模态能量低于l = 1的主模态能量，可以看出，OAM模态值越大，发散性越严重。图4(c)是l = 1与l = 3的OAM复用后的串扰图谱，结合图4(a)及图4(b)可以看出l = 1和l = 3的OAM复用后主模态能量明显高于单模态条件下的能量，由于不同OAM间的能量发生了串扰。接下来，分析模态串扰对多模OAM通信系统带来的影响，并从信号检测的角度降低复杂环境下OAM通信系统误码率。

图 4 单值OAM及复用OAM串扰图谱

下载: 全尺寸图片幻灯片

4.2 基于信号检测的VBLAST-OAM性能分析

图5是VBLAST-OAM通信系统的误码率随着信噪比变化曲线。仿真参数：角向倾斜角 $\gamma = 5 \times$ ${10^{ - 5}}\;{\rm{ rad}}$ ，折射率结构常数 $C_n^2 = 1 \times {10^{ - 14}}\;{{\rm{m}}^{ - 2/3}}$ ，选择模态l = 1和l = 3的OAM。从图中可以看出，未采用本文算法时，随着信噪比的增加，系统误码率下降非常微弱，且信噪比达到一定数值，误码率几乎不再下降；而采用本文算法后的系统误码率呈现快速下降趋势，在很大程度上改善了原有系统性能。例如，在信噪比为18 dB时，基于MRF-BP的系统相较采用本文算法前的系统，误码率大约下降了2.2个数量级；而OAM-OSIC算法与传统的OSIC算法性能几乎保持一致，系统误码率大约下降了1.6～1.9个数量级。另外，根据表2中OSIC算法可知，OAM-OSIC算法省去传统OSIC算法中的求解加权矩阵行范数及范数排序两个步骤，结合图5及图6可知，OAM-OSIC算法在保证性能不变的前提下，能够降低程序运行的计算开销。虽然基于MRF-BP的系统表现出最佳的性能，但程序的运行开销较大；而OAM-OSIC虽然属于次优类算法，却在算法计算开销方面有明显优势。

图 5 不同信噪比下的系统误码率变化曲线

下载: 全尺寸图片幻灯片

图 6 算法运行时间对比

下载: 全尺寸图片幻灯片

图7是角向倾斜角 $\gamma = 5 \times {10^{ - 5}}\;{\rm{ rad}}$ 时，系统误码率随着折射率结构常数的变化曲线。仿真参数：OAM模态选择1和3，信噪比为18 dB。从图中可以看出，随着折射率结构常数的增加，虽然系统误码率处于不断上升状态，但采用本文算法的系统性能得到了明显改善。例如，折射率结构常数 $C_n^2 = 1 \times$ ${10^{ - 14}}\;{{\rm{m}}^{ - 2/3}}$ 时，基于MRF-BP的系统误码率相对采用本文算法前，大约降低了2个数量级；而基于OSIC的系统误码率下降了1.4～1.8个数量级。虽然本文算法可以在一定程度上改善系统性能，但是随着湍流强度的增大，系统的性能依旧在不断下降。

图 7 不同湍流强度下的系统误码率曲线

下载: 全尺寸图片幻灯片

图8是折射率结构常数为 $C_n^2 = 1 \times {10^{ - 14}}\;{{\rm{m}}^{ - 2/3}}$ 时，系统误码率随着角向倾斜角的变化曲线。仿真参数：OAM模态选择1和3，信噪比为18 dB。从图中可以看出，随着倾斜角的增大，系统误码率在逐渐上升，表明孔径失配带给OAM通信系统不可忽视的影响。此外，孔径失配条件下，基于OAM-OSIC的系统误码率依旧能够与传统的OSIC算法性能保持一致。可以推断出，孔径失配造成的模态串扰现象与大气湍流造成的串扰现象类似，基于大气湍流特性的OAM-OSIC算法同样适应于孔径失配角向倾斜的OAM系统。由图8看出，基于MRF-BP算法的系统依旧保持最佳的性能，而OAM-OSIC仍然属于次优算法。例如，在倾斜角 $\gamma = 3 \times {10^{ - 5}}\;{\rm{ rad}}$ 时，基于MRF-BP的系统相对未采用本文算法的系统，误码率下降了2.6个数量级；而基于OAM-ZF-OSIC和OAM-MMSE-OSIC的系统误码率分别降低了1.7个和2个数量级。

图 8 不同倾斜角下的系统误码率曲线

下载: 全尺寸图片幻灯片

5. 结论

本文首先建立了大气湍流、孔径失配场景下的VBLAST-OAM复用通信系统，并对该系统进行相应的理论推导；之后，分析了OAM在复杂环境下模态纯度变化特征；最后，对比分析了OSIC, OAM-OSIC, MRF-BP应用于VBLAST-OAM系统中的性能。结果表明，OAM的模态纯度受多个外界条件的影响，且本文所提算法均能提升VBLAST-OAM系统在复杂场景的性能。当系统开销没有受到严格限制时，MRF-BP算法可以作为一个可靠的应用方案；OAM-OSIC虽然属于次优算法，但在系统开销受限的条件下，可以作为一个优秀的备选方案。

期刊类型引用(2)

1.	王洋，修艳磊，胡韬，施盼盼，廖希 . 基于相位补偿的非理想无线轨道角动量复用通信系统研究. 电子与信息学报. 2022(09): 3212-3219 . 本站查看
2.	卜洋，杨志，赵丽娟，徐志钮. 大气湍流对基于轨道角动量的自由空间光通信影响及解决方案综述. 半导体光电. 2022(06): 1099-1108 . 百度学术

其他类型引用(2)

姓名
邮箱
手机号码
标题
留言内容
验证码

留言板