Robust Secure Transmit Methods for Multibeam Satellite Communication with Imperfect Channel State Information
-
摘要: 考虑多波束卫星通信中存在多天线窃听用户的场景以及非理想信道状态信息带来的系统安全性能恶化,针对合法信道与窃听信道均不理想的多用户MISOME(Multi Input, Single Output, Multi Eavesdropper)系统,该文基于两种类型的信道误差提出了两种鲁棒安全传输方法。在确定误差模型下,利用发送波束成形以及人工噪声协方差的联合优化设计解决系统最差情况安全速率最大化问题(Worst-Case Secrecy Rate Maximization, WC-SRM)问题,通过推导分析将非凸的WC-SRM问题等价为一个1维搜索问题,而这个等价问题的最优解可通过一系列半定规划(SemiDefinite Program, SDP)问题的求解而得到。在统计误差模型下,为求解中断概率约束安全速率最大化(Outage Probability Constrained Secrecy Rate Maximization, OP-SRM)问题,该文提出一种基于WC-SRM的鲁棒安全设计的可靠近似方案。仿真分析验证了该文所提方法的鲁棒性与有效性。Abstract: Considering a general eavesdropping case where the eavesdropper is equipped with multiple antennas and two uncertainty models for the imperfect Channel State Information (CSI) on the main and eavesdroppers channels at the transmitter. Two robust transmission methods are proposed for multiuser Multiple-In, Single-Out, Multiple-antenna Eavesdropper (MISOME) systems to provide SINR guarantee for the legitimate user by the optimization of beamforming and artificial noise covariance matrix. For the deterministic uncertainty model, the Worst-Case Secrecy Rate Maximization (WC-SRM) problem is derived an equivalent problem of the WC-SRM problem through analysis. The equivalent problem can be recast as a single-variable optimization problem, which can be handled by solving a sequence of convex SemiDefinite Programs (SDPs). For the stochastic uncertainty model, a suboptimal scheme to solve the Outage Probability constrained Secrecy Rate Maximization (OP-SRM) problem is proposed based on the robust design for the WC-SRM problem. Finally, the effectiveness and the robustness of the proposed algorithms are validated by the simulation results.
-
Key words:
- Multibeam satellite communication /
- Physical layer security /
- Secrecy rate /
- Robustness
-
1. 引言
在阵列信号处理中,互质阵列[1]的研究是近年来备受关注的课题。利用互素数的性质和互素数组的互质阵,可以得到一个更大的虚拟阵列。众所周知,更大的阵列意味着更大的阵列孔径,这将有助于提高分辨率和干扰抑制能力。互质采样的优点是有助于提高DOA估计的自由度。此外,一些研究将互质采样用于自适应波束形成研究。在文献[2]中,基于压缩感知理论[3],提出了一种适用于互质阵列的自适应波束形成算法,该算法充分利用了压缩感知和互质阵列的优点,结果表明该方法比其他方法具有更好的性能。文献[4-6]提出了几种稳健自适应波束形成算法,但这些算法均是在理想条件下,未考虑幅相误差的影响。
然而,在实际阵列中,从阵列接收的信号包含幅相误差。如果幅相误差存在但被忽视,则自适应波束形成的性能会受到影响。因此,如何校准幅相误差是一个非常重要的问题。为了解决该问题。许多算法被提出。大多数提出的算法,例如基于特征结构的自校准方法[7,8],可以精确校准幅相误差。然而,在低信噪比和快照数较少的情况下,这些算法的性能会急剧下降。文献[9,10]提出了一种基于压缩感知的幅相误差校准算法,但该算法主要用于直线阵列。文献[11]提出了一种幅相误差条件下的稳健自适应波束形成算法,该算法主要用于均匀直线阵。文献[12]提出了一种扰动情况下的阵列信号处理模型,该方法针对均匀直线阵,将阵列扰动视为一个与信号相关的误差,建立了一个扰动存在条件下的信号接收模型,将扰动的校准问题转化为了对误差参数的估计问题,并利用文献[13,14]中的最小二乘算法实现参数的有效估计。
为了在低信噪比和少量快照的情况下准确地校准幅相误差并重构信号的协方差矩阵,本文提出了一种基于最小二乘(TLS)的方法。首先,建立了具有幅相误差的协方差矩阵重构的基本模型。然后将问题转化为变量误差(EIV)模型,此时问题就转化为估计与幅相误差相关的误差矩阵,而不是直接校准幅相误差,并提出一种交替下降算法来估计误差矩阵和重构协方差矩阵。该方法有望消除幅相误差的影响,并优于部分其他算法。最后利用重构的信号协方差矩阵进行自适应波束形成,从而降低幅相误差的影响。
2. 信号模型
考虑两个均匀直线阵列,第1个阵列包含2M个阵元,阵元之间的间距为Nd,第2个阵列包含N个阵元,阵元之间的间距为Md。其中,d=λ/2,λ表示波长,M和N满足M<N。根据质数的性质,这两个均匀直线阵列可以等效为一个包含2M+N−1个阵元的互质阵列。此时,该互质阵的信号接收模型可表示为
x(t)=xs(t)+xi(t)+n(t) (1) 其中,t表示第t次快拍。xs(t)=a(θs)s(t)表示阵列输出的目标信号,a(θs)表示导向矢量,θs表示信号的实际方位信息,xi(t)表示干扰信号,n(t)表示均值为0方差为1的高斯白噪声。
为了自适应地抑制干扰,实现对实际目标的有效探测,常常会选择合适的自适应波束形成算法对干扰进行抑制。其中,最常用的一种自适应波束形成算法就是线性约束最小方差(Linearly Constrained Minimum Variance, LCMV)算法。LCMV算法中最重要的步骤是获得优化的权值向量。该权值向量可以表示为
w=R−1a(θs)aH(θs)R−1a(θs) (2) 其中,R表示信号的协方差矩阵,其表达式为
R=E[xxH]=K∑i=1E[|si|2]a(θi)aH(θi)+σ2I (3) 其中,K表示信号的个数,si表示信号的幅度。I表示单位阵。
由压缩感知理论可知,压缩感知重构算法能够精确重构信号的一个前提是信号在空域上满足稀疏性,为了满足信号的稀疏性,将整个信号空间[−90∘,90∘]划分成Q等份,同时将信号的协方差矩阵向量化,此时可得
r=vec(R)=Bp+σ2nvec(I) (4) 其中,B=[a∗(θ1)⊗a(θ1),a∗(θ2)⊗a(θ2),⋯,a∗(θQ)⊗a(θQ)], p=[0,⋯,s2s,⋯,s22,⋯,s2K,⋯,0]∈CQ×1, ⊗表示Kronecker积,s2s表示目标信号的功率,s22,⋯,s2K表示K–1个干扰信号的功率,σ2n表示噪声信号的功率。
以上模型建立在阵列精确已知的理想情况下,但在实际系统中,由于雷达系统器件本身的非理想性,比如滤波器特性不一致、放大器的幅相特性不一致等因素,会导致阵列幅相误差的出现,且这些误差在实际系统中无法避免,幅相误差的存在,会导致雷达参数估计性能的下降。考虑幅相误差,令
Γ=diag[1,Γ2,⋯,Γ2M+N−1]T=diag[1,ρ2ejϕ2,⋯,ρ2M+N−1ejϕ2M+N−1]T (5) 表示幅相误差矩阵。其中,Γi表示第i个阵元的幅相误差,ρi和ϕi分别表示第i个阵元的幅度误差和相位误差。此时,向量化的信号协方差矩阵可表示为
r′=vec(R′)=B′p+σ2nvec(I) (6) 其中,B′=[Γ∗a∗(θ1)⊗Γa(θ1),⋯,Γ∗a∗(θQ)⊗Γa(θQ)], “∗”表示共轭。
根据Kronecker积的运算性质,
Γ∗a∗(θi)⊗Γa(θi)=(Γ∗⊗Γ)[a∗(θi)⊗a(θi)] (7) 因此有
B′=(Γ∗⊗Γ)B=Γ′B (8) 其中,Γ′=Γ∗⊗Γ。
此时,信号的协方差矩阵的向量化形式可表示为
r′=Γ′Bp+n′ (9) 其中,n′=σ2nvec(I)。
从上述模型中可以看出,由于B跟阵列相关,因此可认为是已知的,如果能够准确估计出Γ′和p,则可以在存在幅相误差的情况下准确构建信号协方差矩阵。
3. 提出的方法
为了得到该问题的求解,我们首先将式(9)中的模型转换为一个EIV模型
r′=Γ′Bp+n′=[B+(Γ′−I)B]p+n′=[B+E]p+n′ (10) 其中,E=(Γ′−I)B表示和幅相误差矩阵相关的误差矩阵。此时,对幅相误差的估计问题就转化为了对误差矩阵E的求解问题。通过文献[12-14],式(10)中的问题可转换为对以下问题的求解
argminp,E,n′‖[En′]‖2F+λ‖p‖1s.t.r′=[B+E]p+n′} (11) 其中‖⋅‖F表示F范数,‖⋅‖1表示1范数。
式(11)的问题中,主要目标是利用有效的参数估计算法,实现对未知参数的有效估计。但是由于和误差矩阵相关的Ep的存在,式(11)中的求解问题是一个非凸优化的问题,难以直接通过凸优化算法进行参数的求解。为了解决该问题,将式(11)中的优化问题转换为一个无约束的优化问题,可转化为
minp,E‖E‖2F+‖r′−[B+E]S‖22+λ‖p‖1,2 (12) 通过式(12)可以得到E和p的求解,但该问题仍然是一个非凸优化问题,仍然无法直接求解,为了得到问题的求解,本文中采用梯度下降算法,通过迭代的方式进行求解。该算法主要通过迭代对多个未知参数进行求解,首先假设其中一个待求解的参数是已知的,然后利用该参数求解另一个未知参数,再利用求解的参数估计另外的未知参数,直到迭代收敛,从而得到多个未知参数的求解。
首先,假设误差矩阵E是已知的,在第i次迭代过程中,稀疏目标信号和干扰信号向量p可通过式(13)进行求解:
minpi‖pi‖1,2s.t.‖r′−[B+Ei−1]pi‖F≤ε (13) 其中,ε为常数。
式(13)中的问题可通过压缩感知信号重构方法进行求解,例如l1/l2混合范数优化方法、贪婪算法等,本文选择贪婪算法进行问题的求解,主要选择OMP(Orthogonal Matching Pursuit)算法,OMP算法的最大特点就是操作简单、计算复杂度较低。得到了稀疏目标信号和干扰信号向量p的估计,将估计得到的稀疏信号向量p代入式(12),可以得到误差矩阵Ei的估计如式(14)所示
Ei=minE‖E‖22+‖r′−[B+E]pi‖2F (14) 对式(14)关于E求偏导数并等于0,可得到误差矩阵E的求解
Ei=[r′−Bpi](pi)T[I+pi(pi)T]−1 (15) 本文提出的算法主要是将贪婪算法和最小二乘算法相结合,进行参数的估计。其中误差矩阵E的初始化为E0=0。该算法主要包括两个关键步骤:式(13)中对稀疏目标信号和干扰信号向量p的估计以及式(14)中对和幅相误差相关的误差矩阵E的估计。当满足迭代终止条件时,迭代终止。通过估计的稀疏目标信号和干扰信号向量p,以及已知的B,可以通过式(4)得到信号的协方差矩阵向量化的有效估计:
r=Bp (16) 从而得到重构的信号协方差矩阵,并将其代入式(2)中,得到自适应波束形成的加权向量的有效估计,从而得到互质阵列的自适应波束形成结果
y(t)=wHx(t) (17) 建立了信号模型,并提出了估计幅相误差和重建信号协方差矩阵的方法。第4节将进行仿真分析,以验证所提方法的有效性。
4. 仿真分析
对于所有的仿真试验,假设两个均匀直线阵的阵元数目分别为2M=10以及N=11,假设1个目标信号和1个干扰信号,信号角度分别为θs=15∘, θ1=45∘。对于幅相误差,假设幅度误差和相位误差均满足均匀分布,其中,幅度误差为U=[−0.2,0.2],相位误差为U=[−20∘,20∘]。
首先利用重构得到的信号协方差矩阵,得到自适应波束形成的最优权系数,得到自适应波束形成之后的方向图,并和含有幅相误差情况下直接进行自适应波束形成以及理想情况下的自适应波束形成结果比较。
从图1可以看出,在含有幅相误差的高信噪比条件下,利用本文提出的算法对信号的协方差矩阵进行重构之后再做自适应波束形成,得到的波束方向图基本和理想的均匀直线阵的结果相近。和含有幅相误差的结果相比,提出的算法在很大程度上降低了幅相误差的影响,得到的方向图旁瓣更低,干扰方向的零陷更深。
为了进一步验证提出的方法的有效性,进行蒙多卡罗仿真试验,分别仿真了相干信号和非相干信号的蒙多卡罗仿真试验,通过输出信干噪比进行不同方法的衡量比较
OutputSINR=σ2s|wHa(θs)|2wHRw (18) 其中,σ2s是目标信号的功率,a(θs)是目标信号的导向矢量。蒙特卡罗分析次数为500次。我们将本文方法与文献[15]中基于稀疏表示技术(SR technique)的鲁棒自适应波束形成方法和文献[2]中的互质阵列自适应波束形成方法(CA-ABF)进行比较。
首先仿真分析不同信噪比条件下,不同方法的输出信干噪比,信噪比从–20 dB到20 dB间变化,结果如图2所示。
从图2可以看出,当SNR较低时,输出SINR将显著降低。由于信噪比较低,信号协方差矩阵的估计偏差较大。当信噪比较高时,信号协方差矩阵被准确估计,输出信噪比显著提高。与其他方法相比,本文提出的方法的输出信噪比更高。当考虑到信号相干问题时,基于压缩感知的方法基本不受信号相干性的影响。
接着仿真分析不同快拍数条件下,不同方法的输出信干噪比,快拍数的变化范围为10~100,结果如图3所示。
图3中,当快拍数量较少时,信号的协方差矩阵估计有一定的偏差,导致输出SINR低于快照较多的情况。当快拍数量增加时,输出SINR增加并趋于稳定。与其他方法相比,本文方法的输出信噪比也更高。
5. 结论
在本文中,针对存在幅相误差的互质阵列,提出了一种基于协方差矩阵重构的鲁棒自适应波束形成方法。首先,建立了协方差矩阵重构的基本模型。然后我们要解决的问题变成了EIV问题。然后,将幅相误差校准问题转换为与幅相误差相关的误差矩阵的估计。为了解决这一问题,提出了一种基于TLS的交替下降算法。从仿真结果可以看出,相对于不同的信噪比和快照,该方法可以获得比其他方法更好的性能。但目前我们的研究还存在一些局限性,仍然面临着失配问题。这影响了它在实际系统中的应用。
-
ZHANG Xi, MCKAY M R , ZHOU Xiangyun, et al. Artificial-noise-aided secure multi-antenna transmission with limited feedback[J]. IEEE Transactions on Wireless Communications, 2015, 14(5): 2742-2754. doi: 10.1109/TWC. 2015.2391261. LI Na, TAO Xiaofeng, and XU Jin. Artificial noise assisted communication in the multiuser downlink: optimal power allocation[J]. IEEE Wireless Communations Letters, 2015, 19(2): 295-298. doi: 10.1109/LCOMM.2014.2385779. WYNER A D. The wire-tap channel[J]. Bell System Technical Journal, 1975, 54(8): 1355-1387. doi: 10.1002/j. 1538-7305.1975.tb02040.x. CSISZAR I and KOMER J. Broadcast channels with confidential messages[J]. IEEE Transactions on Information Theory, 1978, 24(3): 339-348. doi: 10.1109/TIT.1978. 1055892. GOEL S and NEGI R. Guaranteeing secrecy using artificial noise[J]. IEEE Transactions on Wireless Communications, 2008, 7(6): 2180-2189. doi: 10.1109/TWC.2008.060848. WANG Huiming, ZHENG Tongxin, and XIA Xianggen. Secure MISO wiretap channels with multiantenna passive eavesdropper artificial noise vs. artificial fast fading[J]. IEEE Transactions on Wireless Communications, 2015, 14(1): 94-106. doi: 10.1109/TWC.2014.2332164. LIU Shuiyin, HONG Yi, and VITERBO Emanuele. Guaranteeing positive secrecy capacity for MIMOME wiretap channels with finite-rate feedback using artificial noise[J]. IEEE Transactions on Wireless Communications, 2015, 14(8): 4193-4203. doi: 10.1109/ TWC.2015.2417886. TANG Yanqun, XIONG Jun, MA Dongtang, et al. Robust artificial noise aided transmit design for MISO wiretap channels with channel uncertainty[J]. IEEE Communications Letters, 2013, 17(11): 2096-2099. doi: 10.1109/LCOMM.2013. 100713.131673. LI Qiang and MA Wingkin. Spatially selective artificial-noise aided transmit optimization for MISO multi-eves secrecy rate maximization[J]. IEEE Transactions on Signal Processing, 2013, 61(10): 2704-2717. doi: 10.1109/TSP.2013. 2253771. JIANG Lei, ZHU Han, Marangeles Vazquez-Castro, et al. Secure satellite communication systems design with individual secrecy rate constraints[J]. IEEE Transactions on Information Forensics and Security, 2011, 6(3): 661671. doi: 10.1109/TIFS.2011.2148716. GAN Zheng, ATHANASIOS D, and OTTERSTEN Bjorn. Physical layer security in multibeam satellite systems[J]. IEEE Transactions on Wireless Communications, 2012, 11(2): 852863. doi: 10.1109/TWC.2011.120911.111460. DESTOUNIS A and ATHANASIOS D. Dynamic power allocation for broadband multi-beam satellite communication networks[J]. IEEE Communications Letters, 2011, 15(4): 380382. doi: 10.1109/LCOMM.2011.020111.102201. 赵家杰, 彭建华, 黄开枝, 等. 基于人工噪声的多用户 MIMO 系统加密算法[J]. 电子与信息学报, 2012, 34(8): 19391943. doi: 10.3724/SP.J.1146.2011.01068. ZHAO Jiajie, PENG Jianhua, HUANG Kaizhi, et al. A multi- user MIMO system encryption algorithm based on artificial noise[J]. Journal of Electronics Information Technology, 2012, 34(8): 19391943. doi: 10.3724/SP.J.1146.2011.01068. ZHU Jun, SCHOBER Robert, and BHARGAVA V K. Linear precoding of data and artificial noise in secure massive MIMO systems[J]. IEEE Transactions on Wireless Communications, 2016, 15(3): 22452261. doi: 10.1109/TWC.2015.2500578. LI Qiang, MA Wingkin, and SO Anthony Mancho. A safe approximation approach to secrecy outage design for MIMO wiretap channels[J]. IEEE Signal Processing Letters, 2014, 21(1): 118-121. doi: 10.1109/LSP.2013.2293884. LI Qiang and MA Wingkin. Optimal and robust transmit designs for MISO channel secrecy by semidefinite programming[J]. IEEE Transactions on Signal Processing, 2011, 59(8): 3799-3812. doi: 10.1109/TSP.2011.2146775. LUO Zhiquan, MA Wingkin, SO Anthony Mancho, et al. Semidefinite relaxation of quadratic optimization problems [J]. IEEE Signal Processing Magazine, 2010, 27(3): 20-34. doi: 10.1109/MSP.2010.936019. DATTORRO J. Convex Optimization and Euclidean Distance Geometry[M]. Meboo Publishing USA, 2005: 404-410. LUO Zhiquan, STURM F, and ZHANG Shuzhong. Multivariate nonnegative quadratic mappings[J]. SIAM Journal on Control and Optimization, 2004, 14(4): 1140-1162. doi: 10.2139/ssrn.545582. KOLDA T, LEWIS R, and TORCZON V. Optimization by direct search: new perspectives on some classical and modern methods[J]. SIAM Review, 2003, 45(3): 385-482. doi: 10.1137 /S003614450242889. XIONG Jun, MA Dongtang, WONG Kaikit, et al. Robust masked beamforming for MISO cognitive radio networks with unknown eavesdroppers[J]. IEEE Transactions on Vehicular Technology, 2016, 65(2): 744-755. doi: 10.1109/TVT.2015. 2400452. 期刊类型引用(2)
1. 张文青,李胤辰,陈胜垚,何成,田巳睿. 基于深度展开ADMM网络的稳健自适应波束形成. 现代雷达. 2024(06): 43-49 . 百度学术
2. 佟昊阳,刘玉财,易文胜,李水. 基于波束形成的水下声源精确测量. 计量科学与技术. 2023(12): 27-33+66 . 百度学术
其他类型引用(4)
-
计量
- 文章访问数: 1329
- HTML全文浏览量: 129
- PDF下载量: 443
- 被引次数: 6