Loading [MathJax]/jax/element/mml/optable/GreekAndCoptic.js
高级搜索

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

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

无人机中继协助的有限码长隐蔽通信

胡锦松 吴林梅 束锋 陈由甲 郑海峰

胡锦松, 吴林梅, 束锋, 陈由甲, 郑海峰. 无人机中继协助的有限码长隐蔽通信[J]. 电子与信息学报, 2022, 44(3): 1006-1013. doi: 10.11999/JEIT211372
引用本文: 胡锦松, 吴林梅, 束锋, 陈由甲, 郑海峰. 无人机中继协助的有限码长隐蔽通信[J]. 电子与信息学报, 2022, 44(3): 1006-1013. doi: 10.11999/JEIT211372
HU Jinsong, WU Linmei, SHU Feng, CHEN Youjia, ZHENG Haifeng. UAV-relay Assisted Covert Communication with Finite Block-length[J]. Journal of Electronics & Information Technology, 2022, 44(3): 1006-1013. doi: 10.11999/JEIT211372
Citation: HU Jinsong, WU Linmei, SHU Feng, CHEN Youjia, ZHENG Haifeng. UAV-relay Assisted Covert Communication with Finite Block-length[J]. Journal of Electronics & Information Technology, 2022, 44(3): 1006-1013. doi: 10.11999/JEIT211372

无人机中继协助的有限码长隐蔽通信

doi: 10.11999/JEIT211372
基金项目: 国家自然科学基金(62001116, 62071234, 61971139, 61771244)
详细信息
    作者简介:

    胡锦松:男,1989年生,博士,讲师,研究方向为隐蔽通信、物理层安全、天线阵列信号处理、无人机通信等

    吴林梅:女,1997年生,硕士生,研究方向为隐蔽通信、无人机通信等

    束锋:男,1973年生,博士,博士生导师,教授,研究方向为无线信息安全传输、大规模MIMO、无人机通信、无线定位技术等

    陈由甲:女,1983年生,博士,博士生导师,教授,研究方向为移动通信、边缘计算、工业物联网、深度学习等

    郑海峰:男,1978年生,博士,博士生导师,教授,研究方向为智能物联网、车联网、无线感知、张量理论及其应用等

    通讯作者:

    束锋 shufeng@hainanu.edu.cn

  • 中图分类号: TN92

UAV-relay Assisted Covert Communication with Finite Block-length

Funds: The National Natural Science Foundation of China (62001116, 62071234, 61971139, 61771244)
  • 摘要: 由于中长距离的无线隐蔽通信中信号发射功率较大,容易被监测者检测到,针对此问题该文提出无人机(UAV)中继协助的有限码长隐蔽传输方案。首先,根据无人机通信的信道模型,推导从发射机经过无人机中继到合法接收机的信噪比(SNR),得到隐蔽传输中断概率及吞吐量的表达式;然后,分析监测者的检测性能, 即KL(Kullback-Leibler)散度,并将其作为隐蔽通信的约束条件;最后,在此约束的条件下,联合优化发射机和中继的发射功率以及无人机的飞行高度,从而最大化隐蔽通信的吞吐量。仿真结果验证了系统的隐蔽性能,并且通过与缺少无人机协助的传统方案对比,揭示所提出的无人机中继隐蔽传输方案可以有效地降低监测者的KL散度。
  • 随着智能无线终端的快速发展,大量机密和私人信息将通过无线的方式进行传输,通信的安全性得到日益广泛的关注。由于无线通信的广播性质,没有明确的物理边界,这使得传输的机密信息容易受到外部攻击和窃听。传统加密算法或是物理层安全技术,都只能保障信息在通信过程中最大限度不被监测者进行正确解码,而无法对信息传输过程起到保护作用,因此如何使监测者无法准确检测到信息的传输情况是十分重要的,于是为了降低机密信息的传输行为被监测者检测的概率,引入了“隐蔽通信”这个概念。当监测者的检测错误概率足够大时,隐蔽通信可以实现[1]。近年来,众多专家学者根据不同的应用场景,对隐蔽通信进行了广泛的研究。文献[2]研究了监测者通过统计不同时隙的接收功率对接收机噪声进行估计来提高发现隐蔽通信的概率。有些研究关注了多天线对隐蔽通信的影响。比如,文献[3]研究了在监测者随机分布环境下,集中式和分布式多天线协助的隐蔽通信。文献[4]发现监测者的检测设备为多天线阵列时,隐蔽传输速率会明显下降。这是因为监测者的检测性能会随着天线数目的增加而提高。文献[5]提出的隐蔽通信方案能够最大限度地提高多天线发射机和全双工干扰接收器之间的连接吞吐量。

    在无线通信的过程中,随着通信距离的进一步增加,直接的点对点传输需要很大的发射功率,极易在监测者处暴露隐蔽传输的行为,采用中继技术是一种切实可行的方案。文献[6]研究了衰落信道中继网络中的隐蔽通信,在此基础上,文献[7]研究了中继进行多跳信息传输有助于提高性能,因此,增加中继节点可以减小传输的距离,节点之间可以采用较小的信号发射功率,这样能有效地降低传输被监测者检测到的概率。由于无人机辅助通信低成本和灵活部署的特点,常被应用于隐蔽通信的场景中[8]。其中,可以把无人机部署为移动中继,连接两个孤立的地面用户,它将场景从静态变为动态,灵活地避免了监测者的检测,提高了传输的隐蔽性。文献[9]提出了基于无人机网络的隐蔽通信方案,通过联合优化无人机的轨迹和发射功率来实现安全隐蔽的传输。当地面存在监测者的位置不确定时,文献[10]通过最坏隐蔽性约束的资源配置,最大限度地提高平均隐蔽速率。文献[11,12]提出了由多天线干扰器辅助无人机的隐蔽通信方案,最大限度地提高地面发射机和无人机接收机针对几个随机分布的监测者之间的传输速率。值得注意的是,在空对地信道进行通信时,视距通信(Line of Sight, LoS)的概率很高,无人机中继极有可能遭到监测者的攻击,给隐蔽通信的实现增加了难度,同时,目前对于隐蔽通信的研究大多数基于传输符号数为无限长,但是在实际应用中,符号数通常是有限的(即低时延)。

    因此,针对以上问题,本文考虑了一个长距离通信的场景,提出基于无人机中继有限码长隐蔽传输方案,隐蔽信息从源节点出发,通过中继节点的放大转发,最终到达目的节点。首先,构建了基于无人机中继形成的隐蔽通信系统模型,分析了信息传输过程的性能;然后,分析了监测者最优检测性能;最后,通过优化发射功率及无人机的高度,得到隐蔽吞吐量的最大值。仿真结果证明了所提出的方案能使监测者的检测性能降低,同时得到无人机中继的最优高度和最佳发射功率,为实际隐蔽通信系统设计及参数配置提供了理论指导。

    本文考虑了一个隐蔽通信的场景。如图1所示,包含发射机Alice、监测者Willie、合法接收机Bob,其中发射机在监测者的检测下,借助无人机(Unmanned Aerial Vehicle, UAV)向合法接收机发送机密信息。定义Alice-Willie和Willie-Bob之间的距离分别表示为dawdwb,UAV-relay高度为h。本文假设UAV-Relay在半双工模式下工作,因此,从Alice到Bob的传输分两个阶段进行:第1阶段(Alice向UAV-Relay传输),第2阶段(UAV-Relay向Bob传输)。与此同时,假设所有的节点只配备了1根天线。

    图 1  无人机-中继协助的无线隐蔽通信

    考虑将UAV到Bob, Alice和Willie之间的信道模拟成一个空对地的无线信道,它具有非视距(Non Line of Sight, NLoS)和概率视距(Line of Sight, LoS)成分。具体地说,对于无人机和地面用户之间的LoS和NLoS的路径损耗为[13,14]

    ρj={dξLj,LoS dξNj,NLoS  (1)

    其中,dj代表无人机到地面各用户的距离,j{ar,rb,rw},且ar,rb,rw分别代表路径Alice-Relay,Relay-Bob,Relay-Willie,ξL是LoS信道部分的路径损耗指数,ξN是NLoS信道部分的路径损耗指数。一般来说,ξL<ξN。在空对地信道中,LoS成分并不总是存在,而在这些信道中出现LoS分量的概率取决于环境参数以及从UAV-Relay到Alice和Bob的角度[15]。在本文所考虑的情况,从Alice或Bob到UAV-Relay的角度并不总是90°,那么信道中出现LoS分量的概率为

    PLos (θ)=11+aexp[b(θa)] (2)

    其中,ab是S型曲线参数,取决于通信环境。θ是UAV-Relay到地面用户的仰角(单位为度)。由于θ为90°时,信道可以简化为LoS信道,本文考虑最坏的情况,Willie位于UAV-relay下方的垂直方向,即UAV-relay到Willie的仰角为90°。因此,可以得到从UAV-Relay到Willie信道中出现LoS分量的概率为PLos =1。此时,路径损耗为ˉρrw=dξLrw

    在第1阶段,UAV-Relay的接收信号可以表示为

    yr=Paˉρarxa+nr (3)

    其中,xa是Alice发送的机密信息,并且需要满足E[xaxHa]=1Pa是Alice的发射功率,ˉρar表示从Alice到UAV-Relay信道的路径损耗,nr是UAV-Relay方差为σ2r的加性高斯白噪声,即nrCN(0,σ2r)。因此,ˉρar可以表示为

    ˉρar=PLos(θar)dξLar+(1PLos(θar))dξNar (4)

    其中,θar=180πarctan(h/daw)

    根据在UAV-Relay接收到的信号,可以得到在中继的信噪比(Signal-to-Noise Ratio, SNR)为

    γr=Paˉρarσ2r (5)

    在第2阶段,Bob从UAV-Relay接收的信号可以表示为

    yb=PrˉρrbGryr+nb=PrˉρrbPaˉρar+σ2r[Paˉρarxa+nr]+nb (6)

    其中,Pr是UAV-Relay的发射功率,ˉρrb表示从UAV-Relay到Bob信道的路径损耗,nb是Bob方差为σ2b的加性高斯白噪声,即nbCN(0,σ2b)。因此,ˉρrb可以表示为

    ˉρrb=PLos(θrb)dξLrb+(1PLos(θrb))dξNrb (7)

    其中,θrb=180πarctan(h/dwb)

    Gr是UAV-Relay的放大系数,用于保证功率约束。根据文献[16],其可以表示为

    Gr=1Paˉρar+σr2 (8)

    根据式(6),可以得到在合法接收机Bob的信噪比为

    γb=PaˉρarPrˉρrbPrˉρrbσ2r+(Paˉρar+σ2r)σ2b (9)

    引理1 :在合法接收机Bob的信噪比γb是关于发射功率PaPr的单调递增函数。

    证明:将式(9)的分子和分母同时除以PaPr,可将γb转换为

    γb=ˉρarˉρrbσ2rˉρrbPa+σ2b(ˉρarPr+σ2rPaPr) (10)

    可知,Bob信噪比γb随着PaPr的增加而增加。证毕

    本文主要使用隐蔽吞吐量η来评估从发射机Alice到合法接收机Bob的通信性能。对于一个给定的传输速率R,它被定义为η=R(1δ)。其中,δ为中断概率。那么,此时对于一个固定的R,这个中断概率可以表示为

    δ=Q(Nln2(1+γb)[log2(1+γb)R]γb(γb+2)) (11)

    首先将Q函数括号中的整个项表示为μ。于是,μ关于γb的1阶导数为

    dμdγb=N(γ2b+2γb)Nln2[ln(1+γb)ln2R][γb(2+γb)]3/2 (12)

    γb非常小时,可以采用近似相等即ln(1+γb)γb。根据这个近似值,可以得到化简后的μ关于γb的1阶导数为

    dμdγb=N(γ2b+γb+Rln2)[γb(2+γb)]3/2>0 (13)

    从式(13)可以得到dμdγb>0,并且Q函数关于μ的偏导有Qμ<0,那么根据链式法则,可得Q函数关于γb的偏导为Qγb=Qμdμdγb<0。因此,Q函数会随着Bob的信噪比γb增大而减小,即中断概率δγb的单调递减函数。由于隐蔽传输吞吐量ηδ成反比,那么可以知道η关于γb是一个单调递增的函数。

    遵循信号检测中的假设检验理论,首先分析监测者Willie的检测方式,根据Willie在第1阶段和第2阶段接收到的所有信号做出判决。因此,在Willie接收到的信号可以表示为

    yw=Paˉρawxa+PrˉρrwGryr+nw=Paˉρawxa+Prˉρrw1Paˉρar+σ2r(Paˉρarxa+nr)+nw (14)

    其中,nw是Willie方差为σ2w的加性高斯白噪声,即nwCN(0,σ2w)。从Alice和UAV-Relay到Willie的路径损耗分别为ˉρaw=dξLawˉρrw=dξLrw

    根据式(14),在Willie的信噪比可以表示为

    γw=Paˉρaw(Paˉρar+σ2r)+PrˉρrwPaˉρarPrˉρrwσ2r+(Paˉρar+σ2r)σ2w (15)

    在零假设H0下,Alice不向Bob发送信号,反之备选假设H1意味着Alice借助UAV中继向Bob发送信号。基于式(14)中接收到的信号,Willie进行二元假设检验对于接收到的信号是来自H0还是H1。其中,D01是在H0下的似然函数对于H1下似然函数的KL(Kullback-Leibler)散度,可以表达为

    D01=N2[γwln(1+γw)] (16)

    其中,N为无人机隐蔽通信中时延约束对应的最大码长。在本文研究的过程采用D012ϵ2作为隐蔽约束要求,其ϵ是一个任意小的值[17]

    在本节中,通过联合优化Alice的发射功率Pa,UAV-Relay的发射功率Pr及高度h,以便在考虑的情况下最大化有效隐蔽吞吐量η。具体来说,优化问题可以写为下式

    η=maxPa,Pr,hη  (17)
    s.tD012ϵ2
    PaPmaxa
    PrPmaxr
    h_h¯h

    其中,PmaxaPmaxr为Alice和UAV中继的最大发射功率,h_¯h分别为UAV的最小高度和最大高度,隐蔽性约束由式(17a)给出。

    引理2 当隐蔽约束满足 {D}_{01}=2{ϵ}^{2} 时,可以在监测者Willie获得最大的信噪比 {\gamma _{\rm w}} ,即 \gamma _{\rm w}^{\ddag}

    证明 根据式(16),求 {D_{01}} 关于 {\gamma _{\rm w}} 的1阶导数,可以得到

    \frac{{{\rm{d}}{D_{01}}}}{{{\rm{d}}{\gamma _{\rm w}}}} = \frac{{N{\gamma _{\rm w}}}}{{2\left( {1 + {\gamma _{\rm w}}} \right)}} > 0 (18)

    因此, {D_{01}} 是关于 {\gamma _{\rm w}} 的单调递增函数。那么,当 {D}_{01}({\gamma }_{w}^{‡})=2{ϵ}^{2} 时,在监测者Willie可以得到最大的SNR,即 \gamma _{\rm w}^\ddag 。 证毕

    引理3  {D_{01}} 是发射机Alice的发射功率 {P_{\rm a}} 的单调递增函数。

    证明 根据式(15),首先通过求 {\gamma _{\rm w}} 关于 {P_{\rm a}} 的偏导,可以得到

    \frac{{\partial {\gamma _{\rm w}}}}{{\partial {P_{\rm a}}}} = \frac{{{P_{\rm a}}{{\bar \rho }_{{\rm{aw}}}}{{\bar \rho }_{{\rm{ar}}}}B + {P_{\rm r}}{{\bar \rho }_{{\rm{rw}}}}\sigma _{\rm r}^2C + {{\bar \rho }_{{\rm{aw}}}}\sigma _{\rm w}^2\sigma _{\rm r}^4}}{{{A^2}}} > 0 (19)

    其中,

    \begin{split} & A = {P_{\rm r}}{\bar \rho _{{\rm{rw}}}}\sigma _{\rm r}^2 + {P_{\rm a}}{\bar \rho _{{\rm{ar}}}}\sigma _{\rm w}^2 + \sigma _{\rm r}^2\sigma _{\rm w}^2\\ & B = 2{P_{\rm r}}{\bar \rho _{{\rm{rw}}}}\sigma _{\rm r}^2 + {P_{\rm a}}{\bar \rho _{{\rm{ar}}}}\sigma _{\rm w}^2 + 2\sigma _{\rm w}^2\sigma _{\rm r}^2\\ & C = {\bar \rho _{{\rm{aw}}}}\sigma _{\rm r}^2 + {P_{\rm r}}{\bar \rho _{{\rm{ar}}}}\bar \rho _{{\rm{rw}}}^2 + {\bar \rho _{{\rm{ar}}}}\sigma _{\rm w}^2 \end{split}

    基于式(18)、式(19),利用链式法则的推导得到了 {D_{01}} 关于 {P_{\rm a}} 的偏导为

    \frac{{\partial {D_{01}}}}{{\partial {P_{\rm a}}}} = \frac{{{\rm{d}}{D_{01}}}}{{{\rm{d}}{\gamma _{\rm w}}}}\frac{{\partial {\gamma _{\rm w}}}}{{\partial {P_{\rm a}}}} > 0 (20)

    因此,根据式(20), {D_{01}}({P_{\rm a}}) 关于 {P_{\rm a}} 是一个递增的函数。 证毕

    引理4 当{\bar \rho _{{\rm{ar}}}}\sigma _{\rm w}^2 > {\bar \rho _{{\rm{aw}}}}\sigma _{\rm r}^2时, {D_{01}} 是关于UAV-Relay发射功率 {P_{\rm r}} 的单调递增函数。反之, {D_{01}} 是关于 {P_{\rm r}} 的单调递减函数。

    证明 根据式(15), {\gamma _{\rm w}} 关于 {P_{\rm r}} 的偏导可以表示为

    \frac{{\partial {\gamma _{\rm w}}}}{{\partial {P_{\rm r}}}} = \frac{{\left( {P_{\rm a}^2{{\bar \rho }_{{\rm{rw}}}}{{\bar \rho }_{{\rm{ar}}}} + {P_{\rm a}}{{\bar \rho }_{{\rm{rw}}}}\sigma _{\rm r}^2} \right)\left( {{{\bar \rho }_{{\rm{ar}}}}\sigma _{\rm w}^2 - {{\bar \rho }_{{\rm{aw}}}}\sigma _{\rm r}^2} \right)}}{{{A^2}}} (21)

    当令式(21)中的偏导\dfrac{{\partial {\gamma _{\rm w}}}}{{\partial {P_{\rm r}}}} = 0时,可以得到

    K = {\bar \rho _{{\rm{ar}}}}\sigma _{\rm w}^2 - {\bar \rho _{{\rm{aw}}}}\sigma _{\rm r}^2 = 0 (22)

    那么, {\gamma _{\rm w}} 关于 {P_{\rm r}} 的单调性质取决于{\bar \rho _{{\rm{ar}}}}\sigma _{\rm w}^2{\bar \rho _{{\rm{aw}}}}\sigma _{\rm r}^2之间的关系,并进行如下讨论

    情况1:若有{\bar \rho _{{\rm{ar}}}}\sigma _{\rm w}^2 > {\bar \rho _{{\rm{aw}}}}\sigma _{\rm r}^2

    在第1种情况下,由式(21),可得\dfrac{{\partial {\gamma _{\rm w}}}}{{\partial {P_{\rm r}}}} > 0,那么使用链式法则的推导得到

    \frac{{\partial {D_{01}}}}{{\partial {P_{\rm r}}}} = \frac{{{\rm{d}}{D_{01}}}}{{{\rm{d}}{\gamma _{\rm w}}}}\frac{{\partial {\gamma _{\rm w}}}}{{\partial {P_{\rm r}}}} > 0 (23)

    这表明了此时, {D_{01}} 是关于 {P_{\rm r}} 单调递增的函数。

    情况2:若有{\bar \rho _{{\rm{ar}}}}\sigma _{\rm w}^2 < {\bar \rho _{{\rm{aw}}}}\sigma _{\rm r}^2

    在第2种情况下,由式(21),可以得到\dfrac{{\partial {\gamma _{\rm w}}}}{{\partial {P_{\rm r}}}} < 0,这意味着 {\gamma _{\rm w}} 关于 {P_{\rm r}} 是单调递减的,由求导的链式法则得到 {D_{01}} 关于 {P_{\rm r}} 的偏导\dfrac{{\partial {D_{01}}}}{{\partial {P_{\rm r}}}} = \dfrac{{{\rm{d}}{D_{01}}}}{{{\rm{d}}{\gamma _{\rm w}}}}\dfrac{{\partial {\gamma _{\rm w}}}}{{\partial {P_{\rm r}}}} < 0,即此时 {D_{01}} 是关于 {P_{\rm r}} 单调递减的函数。 证毕

    考虑到发射功率 {P_{\rm a}} {P_{\rm r}} 之间的关系会产生不同的结果,本文将优化问题式(17)转化为3个子优化问题。

    定理1 给定在UAV-Relay的发射功率 {P_{\rm r}} ,根据式(17)中的优化问题会得到在发射机Alice存在一个最优发射功率(即 P_{\rm a}^{\text{*}} ),它表示为

    P_{\rm a}^* = \min \{ P_{\rm a}^\ddag ,P_{\rm a}^{\max }\} (24)

    其中, P_{\rm a}^\ddag 是能够满足Willie的信噪比取到最大值的解,即 {\gamma _{\rm w}}(P_{\rm a}^\ddag ) = \gamma _{\rm w}^\ddag

    证明 若给定 {P_{\rm r}} ,那么优化问题式(17)可以进一步写为

    \qquad\qquad\qquad\quad {\eta ^*} = \mathop {\max }\limits_{_{{P_{\rm a}}}} \eta {\text{ }} (25)
    \qquad\qquad \text{s}\text{.t}\text{.     }{D}_{01}\le 2{\epsilon}^{2}\tag{25a}
    \qquad\qquad\qquad\quad {P_{\rm a}} \le P_{\rm a}^{\max } \tag{25b}

    根据引理2,当隐蔽约束满足{D}_{01}=2{\epsilon}^{2}时,可以在Willie获得最大的信噪比 {\gamma _{\rm w}} ,即 \gamma _{\rm w}^\ddag 。此时由引理3可知, {\gamma _{\rm w}} 是关于 {P_{\rm a}} 的单调递增函数,那么当 {\gamma _{\rm w}} 取到最大值时,可以得到一个最大的发射功率 P_{\rm a}^\ddag ,即 {\gamma _{\rm w}}(P_{\rm a}^\ddag ) = \gamma _{\rm w}^\ddag

    于是,将最大的Willie信噪比 \gamma _{\rm w}^\ddag 代入式(15),可以反解 {P_{\rm a}} ,得到如下两个值

    P_{\rm a}^ \pm = \frac{{ - {b_1} \pm \sqrt {b_1^2 + 4{{\bar \rho }_{{\rm{ar}}}}{{\bar \rho }_{{\rm{aw}}}}{c_1}} }}{{2{{\bar \rho }_{{\rm{ar}}}}{{\bar \rho }_{{\rm{aw}}}}}} (26)

    其中,

    \begin{split} & {b_1} = - \gamma _{\rm w}^\ddag {\bar \rho _{{\rm{ar}}}}\sigma _{\rm w}^2 + {\bar \rho _{{\rm{aw}}}}\sigma _{\rm r}^2 + {P_{\rm r}}{\bar \rho _{{\rm{rw}}}}{\bar \rho _{{\rm{ar}}}}\\ & {c_1} = \gamma _{\rm w}^\ddag \sigma _{\rm r}^2(\sigma _{\rm w}^2 + {P_{\rm r}}{\bar \rho _{{\rm{rw}}}}) \end{split}

    由于在式(26)中有- \sqrt {b_1^2 + 4{{\bar \rho }_{{\rm{ar}}}}{{\bar \rho }_{{\rm{aw}}}}{c_1}} < 0,从而导致 P_{\rm a}^ - < 0 ,不满足实际发射功率为正值的意义,故舍去该值。同理可知, P_{\rm a}^ + > 0 ,满足条件,因此可以得到最大发射功率为

    P_{\rm a}^\ddag = P_{\rm a}^ + = \frac{{ - {b_1} + \sqrt {b_1^2 + 4{{\bar \rho }_{{\rm{ar}}}}{{\bar \rho }_{{\rm{aw}}}}{c_1}} }}{{2{{\bar \rho }_{{\rm{ar}}}}{{\bar \rho }_{{\rm{aw}}}}}} (27)

    考虑到式(25b)中的限制约束条件,可以得到在Alice的最优发射功率 P_{\rm a}^* 表示为 P_{\rm a}^* = \min \{ P_{\rm a}^\ddag ,P_{\rm a}^{\max }\}

    证毕

    定理2 当满足{\bar \rho _{{\rm{ar}}}}\sigma _{\rm w}^2 > {\bar \rho _{{\rm{aw}}}}\sigma _{\rm r}^2情况时,对于给定任何在Alice的发射功率 {P_{\rm a}} 都存在一个在UAV-Relay的最优发射功率 P_{\rm r}^* ,表示为

    P_{\rm r}^* = \min \{ P_{\rm r}^\ddag ,P_{\rm r}^{\max }\} (28)

    证明 若给定 {P_{\rm a}} ,那么优化问题式(17)可以进一步写为

    \qquad\qquad\qquad\quad {\eta ^*} = \mathop {\max }\limits_{_{{P_{\rm a}}}} \eta {\text{ }} (29)
    \qquad\qquad \text{s}\text{.t}\text{.     }{D}_{01}\le 2{\epsilon}^{2} \tag{29a}
    \qquad\qquad\qquad\quad {P_{\rm r}} \le P_{\rm r}^{\max } \tag{29b}

    由引理4知,在{\bar \rho _{{\rm{ar}}}}\sigma _{\rm w}^2 > {\bar \rho _{{\rm{aw}}}}\sigma _{\rm r}^2时, {D_{01}} 是关于 {P_{\rm r}} 的单调递增函数,于是将Willie的最大信噪比 \gamma _{\rm w}^\ddag 代入式(15),可以得到满足约束条件式(29a)在UAV-Relay的最大发射功率 P_{\rm r}^\ddag ,表示为

    P_{\rm r}^\ddag = \frac{{\left( { - \gamma _{\rm w}^\ddag \sigma _{\rm w}^2 + {P_{\rm a}}{{\bar \rho }_{{\rm{aw}}}}} \right)\left( {{P_{\rm a}}{{\bar \rho }_{{\rm{ar}}}} + \sigma _{\rm r}^2} \right)}}{{{{\bar \rho }_{{\rm{rw}}}}\left( {\gamma _{\rm w}^\ddag \sigma _{\rm r}^2 - {P_{\rm a}}{{\bar \rho }_{{\rm{ar}}}}} \right)}} (30)

    与此同时,考虑到式(29b)的条件,可以得到在UAV-Relay的最优发射功率表示为 P_{\rm r}^* = \min \{ P_{\rm r}^\ddag , P_{\rm r}^{\max }\} 。证毕

    定理3 当满足{\bar \rho _{{\rm{ar}}}}\sigma _{\rm w}^2 < {\bar \rho _{{\rm{aw}}}}\sigma _{\rm r}^2的情况,并且{P_{\rm r}} = P_{\rm r}^{\max }时,在Alice的最优发射功率 P_{\rm a}^* 可以表示为

    P_{\rm a}^* = \min \{ P_{\rm a}^\ddag ,P_{\rm a}^{\max }\} (31)

    证明 由引理4知,在{\bar \rho _{{\rm{ar}}}}\sigma _{\rm w}^2 < {\bar \rho _{{\rm{aw}}}}\sigma _{\rm r}^2的情况下, {D_{01}} 是关于 {P_{\rm r}} 的单调递减函数,于是当 {D_{01}} 将取到最大值时, {P_{\rm r}} 会得到一个最小发射功率,此时令 {P_{\rm r}} = P_{\rm r}^{\max } ,那么可以得到一个新的满足约束(17b)的在Alice的最大发射功率,表示为

    P_{\rm a}^\ddag = P_{\rm a}^ + = \frac{{ - {b_2} + \sqrt {b_2^2 + 4{{\bar \rho }_{{\rm{ar}}}}{{\bar \rho }_{{\rm{aw}}}}{c_2}} }}{{2{{\bar \rho }_{{\rm{ar}}}}{{\bar \rho }_{{\rm{aw}}}}}} (32)

    其中,

    \begin{split} & {b_2} = - \gamma _{\rm w}^\ddag {\bar \rho _{{\rm{ar}}}}\sigma _{\rm w}^2 + {\bar \rho _{{\rm{aw}}}}\sigma _{\rm r}^2 + {P_{\rm r}}^{\max }{\bar \rho _{{\rm{rw}}}}{\bar \rho _{{\rm{ar}}}}\\ & {c_1} = \gamma _{\rm w}^\ddag \sigma _{\rm r}^2(\sigma _{\rm w}^2 + {P_{\rm r}}^{\max }{\bar \rho _{{\rm{rw}}}}) \end{split}

    于是,同定理1的证明相似地,可以在满足限制约束的条件下得到在Alice的最优发射功率表示为 P_{\rm a}^* = \min \{ P_{\rm a}^\ddag ,P_{\rm a}^{\max }\} 。所提算法的流程如表1 所示。证毕

    表 1  关于最大隐蔽吞吐量及最优发射功率算法
     输入:初始化 {D_{01}} \epsilon \delta {P}_{a}^{\mathrm{max}} {P}_{r}^{\mathrm{max}} {\gamma _w} {\gamma _b} ;给定传输速率 R
     输出: {\eta ^*} P_a^{\text{*}} P_r^{\text{*}}
     (1)将优化问题式(17)转化为两个子优化问题式(25)、式(29);
      (a)根据式(16),求得 {D_{01}} 关于 {\gamma _w} 的1阶导数,将隐蔽性约束转化为一个等式,通过 {D}_{01}({\gamma }_{w}^{‡})=2{\epsilon}^{2} 得到 \gamma _w^\ddag
      (b)通过求 {\gamma _w} 关于 {P_a} 的偏导,利用链式法则,得到 {D_{01}} 关于 {P_a} 的偏导 \dfrac{{\partial {D_{01}}}}{{\partial {P_a}}} = \dfrac{{{\text{d}}{D_{01}}}}{{{\text{d}}{\gamma _w}}}\dfrac{{\partial {\gamma _w}}}{{\partial {P_a}}}
      (c)通过求 {\gamma _w} 关于 {P_r} 的偏导,讨论{\bar \rho _{{\rm{ar}}} }\sigma _w^2{\bar \rho _{{\rm{aw}}} }\sigma _r^2的关系,利用链式法则,得到 {D_{01}} 关于 {P_r} 的偏导 \dfrac{{\partial {D_{01}}}}{{\partial {P_r}}} = \dfrac{{{\text{d}}{D_{01}}}}{{{\text{d}}{\gamma _w}}}\dfrac{{\partial {\gamma _w}}}{{\partial {P_r}}}
     (2)根据式(15),由 {\gamma _w}(P_a^\ddag ) = \gamma _w^\ddag ,可以反解出 P_a^\ddag ,通过给定 {P_r} ,得到最优发射功率 P_a^* = \min \{ P_a^\ddag ,P_a^{\max }\} 来求解优化问题式(25);
     (3)根据式(15),由 {\gamma _w}(P_r^\ddag ) = \gamma _w^\ddag ,得到 P_r^\ddag
      (a)当{\bar \rho _{{\rm{ar}}} }\sigma _w^2 > {\bar \rho _{{\rm{aw}}} }\sigma _r^2时,得到最优发射功率 P_r^* = \min \{ P_r^\ddag ,P_r^{\max }\} 求解优化问题式(29);
      (b)当{\bar \rho _{{\rm{ar}}} }\sigma _w^2 < {\bar \rho _{{\rm{aw}}} }\sigma _r^2时,给定 {P_r} = P_r^{\max } ,得到 P_a^* = \min \{ P_a^\ddag ,P_a^{\max }\}
     (4)根据 \eta = R(1 - \delta ) ,得到最大隐蔽吞吐量 {\eta ^*}
      (a)根据式(11),利用{\rm{Q}}函数,求得 \delta 关于 {\gamma _b} 的1阶导数;
      (b)由于 \eta \delta 成反比,当 \delta 取最小值时,得到最大的隐蔽吞吐量 {\eta ^*}
    下载: 导出CSV 
    | 显示表格

    在本节中,主要使用仿真分析得出的数值结果来评估本文所提出关于无人机中继系统模型的性能。假设发射机Alice与合法接收机Bob之间的距离为一固定值 c ,即{d_{{\rm{aw}}}} + {d_{{\rm{bw}}}} = c。仿真参数分别为,S型曲线参数 a = 4.88 b = 0.1 ,Alice-Willie的距离{d_{{\rm{aw}}}} = 500 m,LoS信道的路径损耗指数 {\xi _L} = - 2 ,NLoS信道的路径损耗指数 {\xi _N} = - 3 ,Alice发射功率 {P_{\rm a}} = 30 dBm,Willie-Bob距离{d_{{\rm{wb}}}} = 500 m, \sigma _{\rm r}^2 = \sigma _{\rm b}^2 = –20 dBm。

    图2揭示了在不同Alice-Willie间的距离{d_{{\rm{aw}}}}以及不同无人机发射功率 {P_{\rm r}} 的情况下,隐蔽吞吐量 \eta 与Alice发射功率 {P_{\rm a}} 之间的关系。可以发现,随着 {P_{\rm a}} 的增大, \eta 从0开始增加最后达到一个最大值,这是因为当发射功率增大后,Bob的信噪比也会随之增大,通过引理1可以知道,隐蔽吞吐量与Bob的信噪比是正相关的关系。但是隐蔽吞吐量不会无限制的不断增大,是因为对于一个固定的传输速率 R ,隐蔽吞吐量的最大值取决于中断概率 \delta 。并且还可以很明确地发现当{d_{{\rm{aw}}}}增大时,必须增大 {P_{\rm a}} 才可以更好地提高隐蔽吞吐量。

    图 2  发射机Alice的发射功率与隐蔽吞吐量的变化关系

    图3展示了在不同隐蔽系数 ϵ 的情况下发射机Alice的最大发射功率 P_{\rm a}^{\max } 与最佳发射功率 P_{\rm a}^* 间的关系。在本图中,可以观察到在 P_{\rm a}^{\max } 相对较小的时候, P_{\rm a}^* = P_{\rm a}^{\max } ,随着 P_{\rm a}^{\max } 的增加,根据定理1,可以得到 P_{\rm a}^* = P_{\rm a}^\ddag ,这是因为当{D}_{01}=2{\epsilon}^{2}时, P_{\rm a}^\ddag 是其对应的解。另外,还可以发现 P_{\rm a}^* 随着 \epsilon 的增加而增加,这是因为当无人机飞行的隐蔽约束不那么严格时,可以得到更高的传输功率。

    图 3  发射机Alice的最大发射功率与最佳发射功率的变化关系

    图4说明了在无人机中继发射功率 {P_{\rm r}} 不同的情况下,无人机垂直高度 h 对合法接收机Bob信噪比 {\gamma _{\rm b}} 的影响,其中,UAV的最小高度 \underline h = 100 m,最大高度 \overline h = 1000 m。由于 {\gamma _{\rm b}} 关于 h 的表达式过于复杂,只能通过1维搜索。由图4可以看出,随着 h 的增加, {\gamma _{\rm b}} 先增加再减小,这是因为随着无人机飞行高度的增加,视距链路的概率会增加,此时,路径损耗也会增大。但是在无人机飞行高度小于最优高度时,Bob的信噪比主要由视距链路决定,于是信噪比增大;然而,随着 h 大于最优高度,路径损耗在信噪比中起决定性作用,信号衰减增大,信噪比减小。与此同时,还可以看出,在相同的高度下,中继发射功率越大,在Bob的信噪比也会相对较大。

    图 4  合法接收机Bob的信噪比与无人机垂直高度的变化关系

    图5通过对比是否存在无人机作为中继的情况下,得出KL散度 {D_{01}} 。可以很清晰地观察到在无人机作为中继的情况下, {D_{01}} 比较低,即 {\gamma _{\rm w}} 比较低,说明此时运用本文提出的系统模型进行长距离通信的性能比较好,隐蔽通信很大程度上能够避免被监测者检测。并且随着Alice发射功率 {P_{\rm a}} 的增大,可以得到 {D_{01}} 增大,那么 {\gamma _{\rm w}} 也逐渐增大,因此,在长距离通信的过程中,要控制发射机Alice的发射功率,使其不超过定理1所提出的最佳发射功率 P_{\rm a}^*

    图 5  Alice-Willie距离与KL散度的变化关系

    本文研究了在长距离通信时,基于无人机中继有限码长的隐蔽传输方案,降低机密信息的传输行为被监测者检测的概率。首先,构建基于无人机中继形成的隐蔽通信场景,分析从发射机经过无人机中继到合法接收机传输过程的性能;然后,研究监测者最优检测性能,根据KL(Kullback-Leibler)散度得出了系统的隐蔽约束;最后,针对考虑的通信场景进行问题的优化,实现最大化隐蔽吞吐量。通过实验仿真得到无人机中继的最优高度和最佳发射功率,表明了所提出的无人机中继方法能使监测者的检测性能降低,得到最大的隐蔽吞吐量,验证了系统的隐蔽性能。在未来的工作中,将进一步研究UAV在噪声环境中以及UAV装备多天线的条件下如何进行隐蔽传输。

  • 图  1  无人机-中继协助的无线隐蔽通信

    图  2  发射机Alice的发射功率与隐蔽吞吐量的变化关系

    图  3  发射机Alice的最大发射功率与最佳发射功率的变化关系

    图  4  合法接收机Bob的信噪比与无人机垂直高度的变化关系

    图  5  Alice-Willie距离与KL散度的变化关系

    表  1  关于最大隐蔽吞吐量及最优发射功率算法

     输入:初始化 {D_{01}} \epsilon \delta {P}_{a}^{\mathrm{max}} {P}_{r}^{\mathrm{max}} {\gamma _w} {\gamma _b} ;给定传输速率 R
     输出: {\eta ^*} P_a^{\text{*}} P_r^{\text{*}}
     (1)将优化问题式(17)转化为两个子优化问题式(25)、式(29);
      (a)根据式(16),求得 {D_{01}} 关于 {\gamma _w} 的1阶导数,将隐蔽性约束转化为一个等式,通过 {D}_{01}({\gamma }_{w}^{‡})=2{\epsilon}^{2} 得到 \gamma _w^\ddag
      (b)通过求 {\gamma _w} 关于 {P_a} 的偏导,利用链式法则,得到 {D_{01}} 关于 {P_a} 的偏导 \dfrac{{\partial {D_{01}}}}{{\partial {P_a}}} = \dfrac{{{\text{d}}{D_{01}}}}{{{\text{d}}{\gamma _w}}}\dfrac{{\partial {\gamma _w}}}{{\partial {P_a}}}
      (c)通过求 {\gamma _w} 关于 {P_r} 的偏导,讨论{\bar \rho _{{\rm{ar}}} }\sigma _w^2{\bar \rho _{{\rm{aw}}} }\sigma _r^2的关系,利用链式法则,得到 {D_{01}} 关于 {P_r} 的偏导 \dfrac{{\partial {D_{01}}}}{{\partial {P_r}}} = \dfrac{{{\text{d}}{D_{01}}}}{{{\text{d}}{\gamma _w}}}\dfrac{{\partial {\gamma _w}}}{{\partial {P_r}}}
     (2)根据式(15),由 {\gamma _w}(P_a^\ddag ) = \gamma _w^\ddag ,可以反解出 P_a^\ddag ,通过给定 {P_r} ,得到最优发射功率 P_a^* = \min \{ P_a^\ddag ,P_a^{\max }\} 来求解优化问题式(25);
     (3)根据式(15),由 {\gamma _w}(P_r^\ddag ) = \gamma _w^\ddag ,得到 P_r^\ddag
      (a)当{\bar \rho _{{\rm{ar}}} }\sigma _w^2 > {\bar \rho _{{\rm{aw}}} }\sigma _r^2时,得到最优发射功率 P_r^* = \min \{ P_r^\ddag ,P_r^{\max }\} 求解优化问题式(29);
      (b)当{\bar \rho _{{\rm{ar}}} }\sigma _w^2 < {\bar \rho _{{\rm{aw}}} }\sigma _r^2时,给定 {P_r} = P_r^{\max } ,得到 P_a^* = \min \{ P_a^\ddag ,P_a^{\max }\}
     (4)根据 \eta = R(1 - \delta ) ,得到最大隐蔽吞吐量 {\eta ^*}
      (a)根据式(11),利用{\rm{Q}}函数,求得 \delta 关于 {\gamma _b} 的1阶导数;
      (b)由于 \eta \delta 成反比,当 \delta 取最小值时,得到最大的隐蔽吞吐量 {\eta ^*}
    下载: 导出CSV
  • [1] YAN Shihao, ZHOU Xiangyun, HU Jinsong, et al. Low probability of detection communication: Opportunities and challenges[J]. IEEE Wireless Communications, 2019, 26(5): 19–25. doi: 10.1109/MWC.001.1900057
    [2] GOECKEL D, BASH B, GUHA S, et al. Covert communications when the warden does not know the background noise power[J]. IEEE Communications Letters, 2016, 20(2): 236–239. doi: 10.1109/LCOMM.2015.2507594
    [3] ZHENG Tongxing, Wang Huiming, NG D W K, et al. Multi-antenna covert communications in random wireless networks[J]. IEEE Transactions on Wireless Communications, 2019, 18(3): 1974–1987. doi: 10.1109/TWC.2019.2900915
    [4] SHAHZAD K, ZHOU Xiangyun, and YAN Shihao. Covert wireless communication in presence of a multi-antenna adversary and delay constraints[J]. IEEE Transactions on Vehicular Technology, 2019, 68(12): 12432–12436. doi: 10.1109/TVT.2019.2948608
    [5] CHEN Xinying, SUN Wen, XING Chengwen, et al. Multi-antenna covert communication via full-duplex jamming against a warden with uncertain locations[J]. IEEE Transactions on Wireless Communication, 2021, 20(8): 5467–5480. doi: 10.1109/TWC.2021.3068096
    [6] HU Jinsong, YAN Shihao, ZHOU Xiangyun, et al. Covert communication achieved by a greedy relay in wireless networks[J]. IEEE Transactions on Wireless Communications, 2018, 17(7): 4766–4779. doi: 10.1109/TWC.2018.2831217
    [7] SHEIKHOLESLAMI A, GHADERI M, TOWSLEY D, et al. Multi-hop routing in covert wireless networks[J]. IEEE Transactions on Wireless Communications, 2018, 17(6): 3656–3669. doi: 10.1109/TWC.2018.2812881
    [8] JIANG Xu, CHEN Xinying, TANG Jie, et al. Covert communication in UAV-assisted air-ground networks[J]. IEEE Wireless Communications, 2021, 28(4): 190–197. doi: 10.1109/MWC.001.2000454
    [9] ZHOU Xiaobo, YAN Shihao, HU Jinsong, et al. Joint optimization of a UAV's trajectory and transmit power for covert communications[J]. IEEE Transactions on Signal Processing, 2019, 67(16): 4276–4290. doi: 10.1109/TSP.2019.2928949
    [10] JIANG Xu, YANG Zhutian, ZHAO Nan, et al. Resource allocation and trajectory optimization for UAV-enabled multi-user covert communications[J]. IEEE Transactions on Vehicular Technology, 2021, 70(2): 1989–1994. doi: 10.1109/TVT.2021.3053936
    [11] CHEN Xinying, CHANG Zheng, TANG Jie, et al. UAV-aided multi-antenna covert communication against multiple wardens[C]. The 2021 IEEE International Conference on Communications, Montreal, Canada, 2021: 1–6.
    [12] CHEN Xinying, ZHANG Ning, TANG Jie, et al. UAV-aided covert communication with a multi-antenna jammer[J]. IEEE Transactions on Vehicular Technology, 2021, 70(11): 11619–11631. doi: 10.1109/TVT.2021.3112121
    [13] ZHOU Yi, YEOH P L, CHEN He, et al. Improving physical layer security via a UAV friendly jammer for unknown eavesdropper location[J]. IEEE Transactions on Vehicular Technology, 2018, 67(11): 11280–11284. doi: 10.1109/TVT.2018.2868944
    [14] ALZENAD M, EL-KEYI A, and YANIKOMEROGLU H. 3-D placement of an unmanned aerial vehicle base station for maximum coverage of users with different QoS requirements[J]. IEEE Wireless Communications Letters, 2018, 7(1): 38–41. doi: 10.1109/LWC.2017.2752161
    [15] AL-HOURANI A, KANDEEPAN S, and LARDNER S. Optimal LAP altitude for maximum coverage[J]. IEEE Wireless Communications Letters, 2014, 3(6): 569–572. doi: 10.1109/LWC.2014.2342736
    [16] KRIKIDIS I, THOMPSON J, MCLAUGHLIN S, et al. Amplify-and-forward with partial relay selection[J]. IEEE Communications Letters, 2008, 12(4): 235–237. doi: 10.1109/LCOMM.2008.071987
    [17] SHU Feng, XU Tingzhen, HU Jinsong, et al. Delay-constrained covert communications with a full-duplex receiver[J]. IEEE Wireless Communications Letters, 2019, 8(3): 813–816. doi: 10.1109/LWC.2019.2894617
  • 期刊类型引用(10)

    1. 吉小鹏,张秀再. 基于极化脏星座的无线隐蔽通信教学实验设计. 实验室研究与探索. 2025(01): 30-34+49 . 百度学术
    2. 李杰,张培航. 远距离无人机信号传输技术研究. 物联网技术. 2025(05): 68-70 . 百度学术
    3. 石火财,张虎,孟绘. 基于甲虫搜索的无人机中继网络部署算法设计. 计算机仿真. 2024(10): 54-58+100 . 百度学术
    4. 郝振超,贾向东,陈智,许晋. 非线性能量采集的WSN状态更新系统信息年龄研究. 计算机工程. 2023(05): 198-205+214 . 百度学术
    5. 章礼玮,李国鑫,陈瑾,王海超,贺文辉,黄育侦. 主动干扰辅助下的无人机隐蔽通信功率与位置联合优化. 信息对抗技术. 2023(02): 28-38 . 百度学术
    6. 白恒志,王海超,李国鑫,龚玉萍. 无人机隐蔽通信网络研究综述. 电信科学. 2023(08): 1-16 . 百度学术
    7. 于季弘,蔡雨庭,刘家豪,王帅,马永锋. 空天隐蔽通信技术综述. 重庆邮电大学学报(自然科学版). 2023(05): 775-788 . 百度学术
    8. 廖晓闽,韩双利,朱璇,林初善,王海鹏. 无人机干扰辅助认知隐蔽通信资源优化算法. 西安电子科技大学学报. 2023(06): 75-83 . 百度学术
    9. 黄英,雷菁,万泽含,安宁. 基于Spinal码的噪声式隐蔽传输. 信号处理. 2023(12): 2153-2164 . 百度学术
    10. 李世党,魏明生,赵娟,刘加跃,唐守锋. 蜂窝物联网中短包域能量效率最大化波束成形优化与设计. 电子与信息学报. 2022(09): 3075-3082 . 本站查看

    其他类型引用(3)

  • 加载中
图(5) / 表(1)
计量
  • 文章访问数:  1401
  • HTML全文浏览量:  813
  • PDF下载量:  178
  • 被引次数: 13
出版历程
  • 收稿日期:  2021-11-30
  • 修回日期:  2022-02-19
  • 录用日期:  2022-02-21
  • 网络出版日期:  2022-03-01
  • 刊出日期:  2022-03-28

目录

/

返回文章
返回