A Fast Random-valued Impulse Noise Detection Algorithm Based on Deep Belief Network
-
摘要:
为提高现有随机脉冲噪声(RVIN)检测算法的检测准确率和执行效率,该文试图从构建描述能力更强的特征矢量和训练非线性映射更为准确的预测模型两个方面入手,实现一种基于训练策略的快速RVIN检测算法。一方面,提取多个不同阶的对数绝对差值排序统计值并结合一个能够反映图像边缘特性的统计值作为刻画图块中心像素点是否为噪声的特征矢量。在计算量增加极少的情况下,显著提升了特征矢量的描述能力。另一方面,基于深度置信网络(DBN)训练RVIN预测模型(RVIN检测器)将特征矢量映射为噪声类型标签,实现了比浅层预测模型更为准确的映射。大量实验数据表明:与现有的RVIN检测算法相比,所提算法在检测准确率和执行效率两个方面都更有优势。
Abstract:To improve the detection accuracy and execution efficiency of the existing Random-Valued Impulse Noise (RVIN) detectors, a fast training-based RVIN detection algorithm is implemented by constructing a more descriptive feature vector and training a detection model with more accurate nonlinear mapping. On the one hand, multiple Rank-Ordered Logarithmic absolute Deviation (ROLD) statistics are extracted and combined with a statistical value reflecting the edge characteristics in the form of feature vector to describe how RVIN-like the center pixel of a patch is. The description ability of the feature vector is improved significantly while the computational complexity is just increased in small amount. On the other hand, an RVIN prediction model (RVIN detector) is obtained by training a Deep Belief Network (DBN) to map the feature vectors to noise labels, which is more accurate than the shallow prediction model. Extensive experimental results show that, compared with the existing RVIN detectors, the proposed one has better performance in terms of detection accuracy and execution efficiency.
-
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中,接收示孔径位置的
η 代表光轴的偏离方位,γ 代表角向倾斜的角度。文献[19,20]表明,波束的横向矢量远远小于纵向矢量的条件下,当出现孔径失配中的角向倾斜,相当于光束穿过一个相位楔,表示为U(r,ϕ,z)=E(r,ϕ,z)exp(jksin(γ)rcos(ϕ−η)) (1) 其中,
E(r,ϕ,z) 表示在空间传输后的波束,U(r,ϕ,z) 表示接收孔径探测到的波束。2.1 轨道角动量及光束在大气湍流中的传播理论
2.1.1 轨道角动量基本理论
根据经典电动力学理论可知,电磁场具备能量与动量且满足电磁守恒。其中,动量主要包括线动量与角动量,线动量通常与平动或者力的作用有关;角动量可以分为SAM和OAM,SAM与电磁波的极化相联系,OAM与电磁波的相位波前分布有关。粒子态OAM沿着传播方向的量子算符可以写成
ˆLz=−jℏ∂/∂ϕ ,本征值方程为ˆLz|l⟩=L|l⟩ ,其中|l⟩ 在基坐标方位角ϕ 表象下,可以写作⟨ϕ|l⟩= 1/√2πexp(jlϕ) ,表示具有l重螺旋相位[21]。LG光束在旁轴近似条件下的亥姆霍兹方程表示为[22]
1ρ∂∂r(ρ∂E∂r)+1ρ2∂2E∂ϕ2+2jk∂E∂z=0 (2) 通过对式(2)求解,可以得到沿着z轴传播的光场
Elp(r,ϕ,z)=√2p!π(p+|l|)!×1ω(z)×Llp[2r2ω2(z)]×[√2rω(z)]|l|×exp[−r2ω2(z)]×exp[−jkr2z2(z2+Z2R)]×exp[−j(2p+|l|+1)arctan−1(zZR)]×exp(jlϕ) (3) ω(z)=ω0√1+(zZR)2ZR=πω20λ} (4) 其中,r是径向距离,
ϕ 为方位角,z为传输距离,ω0 代表束腰半径,ω(z) 是距离z处的光斑半径,ZR 是瑞利距离,Llp 是广义拉盖尔多项式,exp[−j(2p+|l|+1)arctan−1(zZR)]×exp(jlϕ) 代表LG光束相位,exp(jlϕ) 体现了LG的旋转相位特性。2.1.2 OAM在大气湍流中的传播理论
为了能够将OAM光束更好地应用于光无线通信领域,OAM光束在大气湍流中的传播特性分析至关重要。实验室条件下,常采用功率谱反演法及Zernike多项式法等数值仿真方法产生随机相位屏来模拟大气湍流。本文基于修正的冯卡曼折射率功率谱,采用功率谱反演法模拟大气湍流,折射率功率谱模型表示为[23]
Φn(k)=0.033C2nk−11/3exp(−k2k2m)(k2+k20)−11/6,0≤k<∞ (5) 其中,
C2n 表示折射率结构系数,可以衡量湍流的强度,k=√k2x+k2y ,kx 与ky 代表空间频率,k0= 2π/L0 ,km=5.92/l0 ,l0与L0分别代表湍流涡旋的内外尺寸。根据折射率功率谱可以得到大气湍流相位功率谱Φφ(k)=2πk2ΔzΦn(k) (6) 其中,
Δz 代表两个随机相位屏之间的距离,Φn(k) 是大气湍流折射率功率谱。功率谱反演法的基本思想:首先构造一个N维的0均值、单位方差的复高斯随机矩阵h(kx,ky) ,随后用大气湍流相位功率谱对h(kx,ky) 滤波,最后对滤波后的函数做傅里叶逆变换,可得到大气随机相位φ(x,y)=∑kx∑kyh(kx,ky)√Φφ(k)×exp(j(kxx+kyy))ΔkxΔky (7) 其中,
Δkx ,Δky 代表在空间频率间隔。LG光束E(r,ϕ,z) 在大气湍流中的传播过程,可以等效为光束在真空与随机相位屏之间的迭代,具体过程如表1所示。表 1 OAM在大气湍流中的传播过程初始化:源端光场E(r,ϕ,z),空间传输函数H,相位屏个数N,相位φ(x,y) (1) Forj=1:N (2) 傅里叶变换:U(K,z)=FFT(E(r,ϕ,z)) (3) 真空传播:U′(K,z)=IFFT(U(K,z)×H) (4) 穿过随机相位屏:E(r,ϕ,z)=U′(K,z)×exp(j×φ(x,y)) (5) End For 2.2 基于VBLAST-OAM复用通信系统理论分析
本文采用VBLAST方案对发送数据进行编码,将编码后的数据加载到N列OAM光束,那么一个周期内的数据信息表示为
XVBLAST=[xl11,1xl11,2⋯xl11,Mxl22,1xl22,2⋯xl22,M⋮⋮⋱⋮xlNN,1xlNN,2⋯xlNN,M] (8) 其中,
xlNN,M 表示经过VBLAST编码后的数据,上标lN 表示第N列OAM光束,下标(N,M)代表第M时刻在第N列OAM光束上加载的数据。为讨论方便,只取某一个时刻进行讨论。N列OAM光束复用表示为UMUX=N∑n=1xlnn,T×En(r,ϕ,z) (9) 接收孔径得到的光束表示为
U′MUX=N∑n=1xlnn,T×E′n(r,ϕ,z) (10) 接收端通常利用OAM的正交特性对复用的光束进行分离。当OAM的正交性受到破坏后,此时从第k列OAM上得到的信息为
yk=⟨U′MUX(r,ϕ,z),Ek(r,ϕ,z)⟩=∬U′MUX(r,ϕ,z)×E∗k(r,ϕ,z)rdrdϕ=∬N∑n=1xlnn,T×E′n(r,ϕ,z)×E∗k(r,ϕ,z)rdrdϕ=xlkk,T×∬E′k(r,ϕ,z)×E∗k(r,ϕ,z)rdrdϕ+N∑n=1,n≠kxlnn,T∬E′n(r,ϕ,z)×E∗k(r,ϕ,z)rdrdϕ (11) 其中,
(⋅)∗ 表示共轭。根据式(11)可以看出,当复用的OAM光束受外界影响后,会出现信号间串扰。式(11)推导最后一步,首项代表模态k的期望信号,第2项代表除模态k外,其他模态串扰到模态k上的信号。将式(11)写成矩阵形式YN×1=AN×N×XN×1 (12) 其中,Y代表接收到的信号,A代表OAM模态串扰矩阵,X是发送信号。将式(12)展开
[y1y2 ⋮yN]=[a1,1a1,2⋯a1,Na2,1a2,2⋯a2,N⋮⋮⋱⋮aN,1aN,2⋯aN,N]×[xl11,Txl22,T⋮xlNN,T] (13) 其中,矩阵A的对角线元素
aN,N 对应式(11)推导中的第1项,非对角元素对应式(11)推导中的第2项。根据文献[24]可知,任意模态的OAM光束均可以展开成谐波函数
exp(imϕ) 的线性叠加En(r,ϕ,z)=1√2π∞∑m=−∞am,n(r,z)exp(jmϕ) (14) 利用傅里叶变换可以求得
am,n am,n(r,z)=1√2π∫2π0En(r,ϕ,z)exp(−jmϕ)Cm,n=∫∞0|am,n(r,z)|2rdr} (15) 其中,
Cm,n 代表模态为n的OAM在模态m上的谐波功率,模态n的OAM的纯度表示为pn,n=Cn,n∞∑q=−∞Cq,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层对应的加权向量
Wk ,将Wk 与r相乘Wkr=WkXk+Wkn (17) 其中,
Xk 对应第k层的发送数据,n是高斯白噪声。此时第k层检测后信噪比表示为SNRk=⟨|Xk|2⟩σ2‖Wk‖2 (18) 其中,分子表示第k层符号的发射功率,分母表示噪声功率与加权矢量范数平方的乘积。从式(18)可以看出,每一层检测信噪比与行范数成反比,因此可以从范数小的数据层开始检测,降低大信噪比数据层对其他数据层的干扰。虽然MMSE-OSIC算法通常根据信干噪比进行顺序检测,为了便于分析,本文采用行范数排序规则。
表2给出了OSIC算法在VBLAST-OAM系统中的实施步骤。
表 2 OSIC算法初始化:复用的OAM数目Nt,串扰信道H,接收到的信号r,噪声方差σ2 (1) Forj=1:Nt (2) 加权矢量W: 基于ZF/MMSE准则:W=pinv(H)/pinv(HH+σ2×INt−j+1)×HH (3) 首先对W的每一行求范数,并对范数排序,选取最小范数行k (OSIC) (4) 省略步骤(3) (OAM-OSIC) (5) y(k)=W(k,:)×r 判决统计量(优化前) (6) y(k)=W(1,:)×r 判决统计量(优化后) (7) 根据数据判决得到x(k) (8) r=r−x(k)×H 消除前一次检测的数据 (9) H(:,k)=[] 将信道矩阵的第k列清除(优化前) H(:,1)=[] 将信道矩阵的第1列清除(优化后) (10) 重复步骤(2),更新W (11) EndFor 3.2 基于MRF-BP的VBLAST-OAM复用系统
图2是基于马尔科夫随机场的概率图模型。可以看出,马尔科夫随机场中存在许多相互连接的节点,节点又可以分为观察节点和隐藏节点。如果将上述模型与OAM复用系统相对应,那么图中任意观察节点对应从任意模态的OAM上解调信息
yi ,隐藏节点对应任意模态OAM上加载的信息xi 。此外,观察节点与隐藏节点的关系可以用相容函数ϕi(xi) 描述,同样将相邻隐藏节点的关系用势函数ψi,j(xi,xj) 描述。MRF-BP算法就是利用置信度传播(Belief Propagation, BP)算法对马尔科夫随机场中的各个节点的状态信息不断更新,当算法达到收敛时,即可估计出隐藏节点的状态信息概率。光无线OAM复用系统中,场内隐藏节点的联合概率[25]p(x)=1Z∏iϕi(xi,yi)∏i,jψi,j(xi,xj) (19) 其中,Z称配分函数,保证分布函数标准化。利用BP算法可以将实现信息在场内节点相互传递,从而更新MRF中各个节点状态。将节点j传递到节点i的信息表示为
mj,i(xi) ,节点i的置信度用bi(xi) ,那么bi(xi) 可以表示为bi(xi)∝ϕi(xi)∏j∈N(i)mj,i(xi) (20) N(i) 表示节点i的相邻节点集合,从节点j到节点i的传递信息近似表示mj,i(xi)∝∑xjϕj(xj)ψi,j(xi,xj)∏k∈N(j)∖imk,j(xj) (21) 利用势函数及相容函数,可以将信息的传递与更新构建成一个迭代过程。由于多模态OAM在复杂环境传播后会出现模态间的相互串扰,因此基于OAM的马尔科夫随机场是全连接图,此时算法的收敛性较差,最终得到的信息存在较大误差。本文采用信息阻尼的方式加速算法收敛,具体操作就是把上一次迭代得到的估计值与新的估计值各自分配一个加权系数,每次得到的估计值是新旧估计值的加权组合。
mti,j(xj) 是t次迭代时,节点i到j的新信息,那么更新的信息为[26]¯mti,j(xj)∝∑xiϕi(xi)ψi,j(xi,xj){∏k∈N(i)∖jmt−1k,i(xi)}mti,j(xj)=ammt−1i,j(xj)+{(1−am)¯mti,j(xj)} (22) 其中,
αm∈[0,1] 是信息阻尼因子,取0.2。基于马尔科夫随机场图模型的VBLAST-OAM系统信息后验概率表示为p(x|y,H)=(∏i≠jexp(−xiℜ(Ri,j)xj))×(∏iexp(xiℜ(zi)+lnp(xi))) (23) 其中,
R=(1/δ2)HHH ,z=(1/δ2)HHy 。将式(19)与式(23)比较,可以得到两者的关系ψi,j(xi,xj)=exp(−xiℜ(Ri,j)xj) (24) ϕi(xi)=exp(xiℜ(zi)+lnp(xi)) (25) 将MRF-BP算法应用于VBLAST-OAM系统,详细的实施步骤如表3所示。
表 3 MRF-BP算法初始化:m0i,j=b0i, z,R,p(xi=1)=p(xi=−1),∀i,j∈(1,2,⋯,N), M是信息迭代次数 (1) Fori=1:N 势函数 (2) Forj=1:N i≠j (3) 根据式(24)计算ψi,j; (4) EndFor (5) EndFor (6) Fori=1:N 相容函数 (7) 根据式(25)计算ϕi; (8) EndFor (9) Fort=1:M 迭代更新 (10)Fori=1:N (11) Forj=1:N i≠j (12) 第t次迭代得到更新的信息¯mti,j,利用信息阻尼方式得到新信息mti,j; (13) EndFor (14) EndFor (15) EndFor (16) Fori=1:N 置信度计算 (17) 根据式(20)计算置信度bi; (18) EndFor 4. 仿真参数设置及结果分析
仿真参数:光束波长
λ=1.55×10−6m ,传播距离z=1km ,l0=0.001m ,L0=50m ,栅格数为300,调制方式为二进制相移键控(Binary Phase Shift Keying, BPSK)。4.1 大气湍流、孔径失配场景下的OAM模态谱分析
本节首先对复杂场景下OAM的模态纯度进行分析,之后分析多模态OAM之间的串扰特征。
图3是大气湍流、孔径失配条件下的OAM模态纯度3维图。仿真选择模态为3的OAM,从图中可以看出,随着倾斜角、湍流强度的增大,接收端检测到的OAM模态纯度逐渐降低。表明OAM模态纯度会受到多种外界条件的影响,大气湍流、孔径失配均会对OAM的模态纯度带来一定程度的影响。
为了直观地描述式(11)中的模态串扰,图4给出了单模态OAM及多模态OAM的能量串扰分布图谱。仿真参数:折射率结构常数
C2n=1×10−14m−2/3 ,角向倾斜角γ=5×10−5rad 。图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.2 基于信号检测的VBLAST-OAM性能分析
图5是VBLAST-OAM通信系统的误码率随着信噪比变化曲线。仿真参数:角向倾斜角
γ=5× 10−5rad ,折射率结构常数C2n=1×10−14m−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虽然属于次优类算法,却在算法计算开销方面有明显优势。图7是角向倾斜角
γ=5×10−5rad 时,系统误码率随着折射率结构常数的变化曲线。仿真参数:OAM模态选择1和3,信噪比为18 dB。从图中可以看出,随着折射率结构常数的增加,虽然系统误码率处于不断上升状态,但采用本文算法的系统性能得到了明显改善。例如,折射率结构常数C2n=1× 10−14m−2/3 时,基于MRF-BP的系统误码率相对采用本文算法前,大约降低了2个数量级;而基于OSIC的系统误码率下降了1.4~1.8个数量级。虽然本文算法可以在一定程度上改善系统性能,但是随着湍流强度的增大,系统的性能依旧在不断下降。图8是折射率结构常数为
C2n=1×10−14m−2/3 时,系统误码率随着角向倾斜角的变化曲线。仿真参数:OAM模态选择1和3,信噪比为18 dB。从图中可以看出,随着倾斜角的增大,系统误码率在逐渐上升,表明孔径失配带给OAM通信系统不可忽视的影响。此外,孔径失配条件下,基于OAM-OSIC的系统误码率依旧能够与传统的OSIC算法性能保持一致。可以推断出,孔径失配造成的模态串扰现象与大气湍流造成的串扰现象类似,基于大气湍流特性的OAM-OSIC算法同样适应于孔径失配角向倾斜的OAM系统。由图8看出,基于MRF-BP算法的系统依旧保持最佳的性能,而OAM-OSIC仍然属于次优算法。例如,在倾斜角γ=3×10−5rad 时,基于MRF-BP的系统相对未采用本文算法的系统,误码率下降了2.6个数量级;而基于OAM-ZF-OSIC和OAM-MMSE-OSIC的系统误码率分别降低了1.7个和2个数量级。5. 结论
本文首先建立了大气湍流、孔径失配场景下的VBLAST-OAM复用通信系统,并对该系统进行相应的理论推导;之后,分析了OAM在复杂环境下模态纯度变化特征;最后,对比分析了OSIC, OAM-OSIC, MRF-BP应用于VBLAST-OAM系统中的性能。结果表明,OAM的模态纯度受多个外界条件的影响,且本文所提算法均能提升VBLAST-OAM系统在复杂场景的性能。当系统开销没有受到严格限制时,MRF-BP算法可以作为一个可靠的应用方案;OAM-OSIC虽然属于次优算法,但在系统开销受限的条件下,可以作为一个优秀的备选方案。
-
表 1 图1中a1和b1图块上所提取的前m阶ROLD值比较
图块 阶数m 1 2 3 4 5 6 7 8 9 10 11 12 a1 1.63 3.29 4.99 6.68 8.37 10.09 11.80 13.52 15.24 16.98 18.72 20.46 b1 1.00 2.20 3.66 5.12 6.84 8.62 10.50 12.44 14.40 16.40 18.42 20.45 表 2 各噪声检测算法在常用图像集合上的各项性能指标的平均值比较
方法 含噪20% 含噪40% 含噪60% 漏检数 误检数 错检总数 MEMH 漏检数 误检数 错检总数 MEMH 漏检数 误检数 错检总数 MEMH ASWM 3462 10687 14149 14.23 7478 10005 17483 16.40 14720 9804 24524 23.17 PSMF 10695 3585 14279 15.14 23038 3603 26641 30.27 39096 5634 44730 45.81 ROLD-EPR 6567 5106 11673 18.77 9462 8956 18419 15.95 10417 11616 22034 14.04 ROR-NLM 5068 9354 14421 14.91 11906 8873 20779 17.13 22553 12856 35408 23.60 MLP-EPR 8505 2081 10586 22.78 13244 5759 19003 18.09 15017 10113 25130 16.18 本文方法 4084 5909 9992 11.77 7975 8586 16561 12.17 10076 12594 22670 11.53 表 3 各检测算法统一用相同修复算法降噪后在PSNR指标上的比较(dB)
方法 含噪20% 含噪40% 含噪50% 含噪60% Lena House Bridge Lena House Bridge Lena House Bridge Lena House Bridge ASWM 39.06 33.31 25.76 34.27 31.21 24.33 30.66 28.81 23.25 26.04 26.13 21.61 PSMF 30.24 27.82 23.25 29.26 26.03 22.77 26.03 24.09 21.91 22.04 21.98 20.00 ROLD-EPR 34.77 33.31 26.75 31.77 31.21 24.25 30.54 28.81 23.12 28.78 26.13 22.20 ROR-NLM 36.94 28.92 25.28 31.58 28.91 23.59 27.61 27.39 22.35 22.92 24.68 20.39 MLP--EPR 36.45 39.36 27.71 33.83 37.49 24.40 31.86 36.48 23.33 29.42 33.96 22.25 本文方法 40.35 42.95 26.97 35.89 40.31 24.91 33.04 38.27 23.36 29.66 36.12 22.18 表 4 各噪声检测算法平均执行时间的比较(s)
方法 ASWM PSMF ROLD-EPR ROR-NLM MLP-EPR 本文方法 时间 102.72 0.86 10.40 77.19 0.79 0.70 -
LIU Licheng, Chen C L P, and ZHOU Yicong. A new weighted mean filter with a two-phase detector for removing impulse noise[J]. Information Sciences, 2015, 315(September): 4052–4057. doi: 10.1016/j.ins.2015.03.067 SINGH N, THILAGAVATHY T, LAKSHMIPRIYA R T, et al. Some studies on detection and filtering algorithms for the removal of random valued impulse noise[J]. IET Image Processing, 2017, 11(11): 953–963. doi: 10.1049/iet-ipr. 2017.0346 DONG Yiqiu and XU Shufang. A new directional weighted median filter for removal of random-valued impulse noise[J]. IEEE Signal Processing Letters, 2007, 14(3): 193–196. doi: 10.1109/LSP.2006.884014 刘万军, 梁雪剑, 曲海成. 自适应增强卷积神经网络图像识别[J]. 中国图象图形学报, 2017, 22(12): 1723–1736. doi: 10.11834/jig.170079LIU Wanjun, LIANG Xuejian, and QU Haicheng. Adaptively enhanced convolutional neural network algorithm for image recognition[J]. Journal of Image and Graphics, 2017, 22(12): 1723–1736. doi: 10.11834/jig.170079 于海平, 何发智, 潘一腾, 等. 一种基于多特征的距离正则化水平集快速分割方法[J]. 电子学报, 2017, 45(3): 534–539. doi: 10.3969/j.issn.0372-2112.2017.003.004YU Haiping, HE Fazhi, PAN Yiteng, et al. A fast distance regularized level set method for segmentation based on multi-features[J]. Acta Electronica Sinica, 2017, 45(3): 534–539. doi: 10.3969/j.issn.0372-2112.2017.003.004 王森, 王余, 王易川, 等. 水下高速目标声谱图特征提取及分类设计[J]. 电子与信息学报, 2017, 39(11): 2684–2689. doi: 10.11999/JEIT170283WANG Sen, WANG Yu, WANG Yichuan, et al. Feature extraction and classification of spectrum of radiated noise of underwater high speed vehicle[J]. Journal of Electronics &Information Technology, 2017, 39(11): 2684–2689. doi: 10.11999/JEIT170283 周箩鱼, 汤佳欣. 基于图像稀疏表达的模拟退火图像复原[J]. 光电子·激光, 2018, 29(2): 218–223. doi: 10.16136/j.joel.2018.02.0088ZHOU Luoyu and TANG Jiaxin. Image restoration based on simulated annealing with image sparse representation[J]. Journal of Optoelectronics · Laser, 2018, 29(2): 218–223. doi: 10.16136/j.joel.2018.02.0088 AKKOUL S, LEDEE R, LECONGE R, et al. A new adaptive switching median filter[J]. IEEE Signal Processing Letters, 2010, 17(6): 587–590. doi: 10.1109/LSP.2010.2048646 YU Hancheng, ZHAO Li, and WANG Haixian. An efficient procedure for removing random-valued impulse noise in images[J]. IEEE Signal Processing Letters, 2008, 15(12): 922–925. doi: 10.1109/LSP.2008.2005051 GARNETT R, HUEGERICH T, CHUI C, et al. A universal noise removal algorithm with an impulse detector[J]. IEEE Transactions on Image Processing, 2005, 14(11): 1747–1754. doi: 10.1109/TIP.2005.857261 DONG Yiqiu, RAYMOND H C, and XU Shufang. A detection statistic for random-valued impulse noise[J]. IEEE Transactions on Image Processing, 2007, 16(4): 1112–1120. doi: 10.1109/TIP.2006.891348 XIONG Bo and YIN Zhouping. A universal denoising framework with a new impulse detector and nonlocal means[J]. IEEE Transactions on Image Processing, 2012, 21(4): 1663–1675. doi: 10.1109/TIP.2011.2172804 ROY A, SINGHA J, DEVI S S, et al. Impulse noise removal using SVM classification based fuzzy filter from gray scale images[J]. Signal Processing, 2016, 128(11): 262–273. doi: 10.1016/j.sigpro.2016.04.007 TURKMEN I. The ANN based detector to remove random-valued impulse noise in images[J]. Journal of Visual Communication and Image Representation, 2015, 34(10): 28–36. doi: 10.1016/j.jvcir.2015.10.011 SOLEIMANY S and HAMGHALAM M. A novel random-valued impulse noise detector based on MLP neural network classifier[C]. Artificial Intelligence and Robotics (IRANOPEN), Qazvin, Iran, 2017: 165–169. HINTON G, DENG Li, YU Dong, et al. Deep neural networks for acoustic modeling in speech recognition: The shared views of four research groups[J]. IEEE Signal Processing Magazine, 2012, 29(6): 82–97. doi: 10.1109/MSP.2012.2205597 WANG Zhou and ZHANG D. Progressive switching median filter for the removal of impulse noise highly corrupted images[J]. IEEE Transactions on Circuits System II: Analog Digital Signal Processing, 1999, 46(1): 78–80. doi: 10.1109/82.749102 CHAN S H, WANG Xiran, and ELGENDY O A. Plug-and-play ADMM for image restoration: Fixed-point convergence and applications[J]. IEEE Transactions on Computational Imaging, 2017, 3(1): 84–98. doi: 10.1109/TCI.2016.2629286 期刊类型引用(2)
1. 王洋 ,修艳磊 ,胡韬 ,施盼盼 ,廖希 . 基于相位补偿的非理想无线轨道角动量复用通信系统研究. 电子与信息学报. 2022(09): 3212-3219 . 本站查看
2. 卜洋,杨志,赵丽娟,徐志钮. 大气湍流对基于轨道角动量的自由空间光通信影响及解决方案综述. 半导体光电. 2022(06): 1099-1108 . 百度学术
其他类型引用(2)
-