Construction of Optimal Zero Correlation Zone Aperiodic Complementary Sequence Sets

1.   引言

随着智能无线终端的快速发展，大量机密和私人信息将通过无线的方式进行传输，通信的安全性得到日益广泛的关注。由于无线通信的广播性质，没有明确的物理边界，这使得传输的机密信息容易受到外部攻击和窃听。传统加密算法或是物理层安全技术，都只能保障信息在通信过程中最大限度不被监测者进行正确解码，而无法对信息传输过程起到保护作用，因此如何使监测者无法准确检测到信息的传输情况是十分重要的，于是为了降低机密信息的传输行为被监测者检测的概率，引入了“隐蔽通信”这个概念。当监测者的检测错误概率足够大时，隐蔽通信可以实现^[1]。近年来，众多专家学者根据不同的应用场景，对隐蔽通信进行了广泛的研究。文献[2]研究了监测者通过统计不同时隙的接收功率对接收机噪声进行估计来提高发现隐蔽通信的概率。有些研究关注了多天线对隐蔽通信的影响。比如，文献[3]研究了在监测者随机分布环境下，集中式和分布式多天线协助的隐蔽通信。文献[4]发现监测者的检测设备为多天线阵列时，隐蔽传输速率会明显下降。这是因为监测者的检测性能会随着天线数目的增加而提高。文献[5]提出的隐蔽通信方案能够最大限度地提高多天线发射机和全双工干扰接收器之间的连接吞吐量。

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

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

2.   系统模型

2.1   考虑的方案和采用的假设

本文考虑了一个隐蔽通信的场景。如图1所示，包含发射机Alice、监测者Willie、合法接收机Bob，其中发射机在监测者的检测下，借助无人机(Unmanned Aerial Vehicle, UAV)向合法接收机发送机密信息。定义Alice-Willie和Willie-Bob之间的距离分别表示为${d_{{\rm{aw}}}}$和${d_{{\rm{wb}}}}$，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]

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

其中，$a$和$b$是S型曲线参数，取决于通信环境。$\theta$是UAV-Relay到地面用户的仰角(单位为度)。由于$\theta$为90°时，信道可以简化为LoS信道，本文考虑最坏的情况，Willie位于UAV-relay下方的垂直方向，即UAV-relay到Willie的仰角为90°。因此，可以得到从UAV-Relay到Willie信道中出现LoS分量的概率为${\mathcal{P}_{{\rm{Los}}{\text{ }}}} = 1$。此时，路径损耗为${\bar \rho _{{\rm{rw}}}} = d_{{\rm{rw}}}^{{\xi _{\rm{L}}}}$。

2.2   经过UAV-Relay从Alice传输到Bob

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

其中，${{\boldsymbol{x}}_{\rm{a}}}$是Alice发送的机密信息，并且需要满足$\mathbb{E}\left[ {{{\boldsymbol{x}}_{\rm{a}}}{\boldsymbol{x}}_{\rm{a}}^{\rm{H}}} \right] = 1$。${P_{\rm{a}}}$是Alice的发射功率，${\bar \rho _{{\rm{ar}}}}$表示从Alice到UAV-Relay信道的路径损耗，${{\boldsymbol{n}}_{\rm{r}}}$是UAV-Relay方差为$\sigma _{\rm{r}}^2$的加性高斯白噪声，即${{\boldsymbol{n}}_{\rm{r}}}\sim\mathcal{C}\mathcal{N}\left( {0,\sigma _{\rm{r}}^2} \right)$。因此，${\bar \rho _{{\rm{ar}}}}$可以表示为

其中，${\theta _{{\rm{ar}}}} = \dfrac{{180}}{\pi }\arctan \left( {h/{d_{{\rm{aw}}}}} \right)$。

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

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

其中，${P_{\rm{r}}}$是UAV-Relay的发射功率，${\bar \rho _{{\rm{rb}}}}$表示从UAV-Relay到Bob信道的路径损耗，${{\boldsymbol{n}}_{\rm{b}}}$是Bob方差为$\sigma _{\rm{b}}^2$的加性高斯白噪声，即${{\boldsymbol{n}}_{\rm{b}}}\sim\mathcal{C}\mathcal{N}\left( {0,\sigma _{\rm{b}}^2} \right)$。因此，${\bar \rho _{{\rm{rb}}}}$可以表示为

其中，${\theta _{{\rm{rb}}}} = \dfrac{{180}}{\pi }\arctan \left( {h/{d_{{\rm{wb}}}}} \right)$。

${G_{\rm{r}}}$是UAV-Relay的放大系数，用于保证功率约束。根据文献[16]，其可以表示为

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

引理1 ：在合法接收机Bob的信噪比${\gamma _{\rm{b}}}$是关于发射功率${P_{\rm{a}}}$和${P_{\rm{r}}}$的单调递增函数。

证明：将式(9)的分子和分母同时除以${P_{\rm{a}}}{P_{\rm{r}}}$，可将${\gamma _{\rm{b}}}$转换为

可知，Bob信噪比${\gamma _{\rm b}}$随着${P_{\rm a}}$或${P_{\rm r}}$的增加而增加。证毕

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

首先将$Q$函数括号中的整个项表示为$\mu$。于是，$\mu$关于${\gamma _{\rm b}}$的1阶导数为

当${\gamma _{\rm b}}$非常小时，可以采用近似相等即$\ln \left( {1 + {\gamma _{\rm b}}} \right)\sim{\gamma _{\rm b}}$。根据这个近似值，可以得到化简后的$\mu$关于${\gamma _{\rm b}}$的1阶导数为

从式(13)可以得到$\dfrac{{{\text{d}}\mu }}{{{\text{d}}{\gamma _{\rm b}}}} > 0$，并且$Q$函数关于$\mu$的偏导有$\dfrac{{\partial Q}}{{\partial \mu }} \lt 0$，那么根据链式法则，可得$Q$函数关于${\gamma _{\rm b}}$的偏导为$\dfrac{{\partial Q}}{{\partial {\gamma _{\rm b}}}} = \dfrac{{\partial Q}}{{\partial \mu }}\dfrac{{{\text{d}}\mu }}{{{\text{d}}{\gamma _{\rm b}}}} < 0$。因此，$Q$函数会随着Bob的信噪比${\gamma _{\rm b}}$增大而减小，即中断概率$\delta$是${\gamma _{\rm b}}$的单调递减函数。由于隐蔽传输吞吐量$\eta$与$\delta$成反比，那么可以知道$\eta$关于${\gamma _{\rm b}}$是一个单调递增的函数。

2.3   监测者Willie检测性能分析

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

其中，${{\boldsymbol{n}}_{\rm w}}$是Willie方差为$\sigma _{\rm w}^2$的加性高斯白噪声，即${{\boldsymbol{n}}_{\rm w}}\sim\mathcal{C}\mathcal{N}\left( {0,\sigma _{\rm w}^2} \right)$。从Alice和UAV-Relay到Willie的路径损耗分别为${\bar \rho _{{\rm{aw}}}} = d_{{\rm{aw}}}^{{\xi _{\text{L}}}}$及${\bar \rho _{{\rm{rw}}}} = d_{{\rm{rw}}}^{{\xi _{\text{L}}}}$。

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

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

其中，$N$为无人机隐蔽通信中时延约束对应的最大码长。在本文研究的过程采用${D}_{01}\le 2{\epsilon}^{2}$作为隐蔽约束要求，其$\epsilon$是一个任意小的值^[17]。

3.   无人机隐蔽通信的优化问题

在本节中，通过联合优化Alice的发射功率${P_{\rm a}}$，UAV-Relay的发射功率${P_{\rm r}}$及高度$h$，以便在考虑的情况下最大化有效隐蔽吞吐量${\eta ^*}$。具体来说，优化问题可以写为下式

其中，${P}_{{\rm{a}}}^{\mathrm{max}}$和${P}_{{\rm{r}}}^{\mathrm{max}}$为Alice和UAV中继的最大发射功率，$\underline h$和$\overline h$分别为UAV的最小高度和最大高度，隐蔽性约束由式(17a)给出。

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

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

因此，${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}}$的偏导，可以得到

其中，

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

因此，根据式(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}}$的偏导可以表示为

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

那么，${\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$，那么使用链式法则的推导得到

这表明了此时，${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}^\ddag$是能够满足Willie的信噪比取到最大值的解，即${\gamma _{\rm w}}(P_{\rm a}^\ddag ) = \gamma _{\rm w}^\ddag$。

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

根据引理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}}$，得到如下两个值

其中，

由于在式(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$，满足条件，因此可以得到最大发射功率为

考虑到式(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 a}}$，那么优化问题式(17)可以进一步写为

由引理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$，表示为

与此同时，考虑到式(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}^*$可以表示为

证明　由引理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的最大发射功率，表示为

其中，

于是，同定理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 

| 显示表格

4.   仿真结果及讨论

在本节中，主要使用仿真分析得出的数值结果来评估本文所提出关于无人机中继系统模型的性能。假设发射机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散度的变化关系

下载: 全尺寸图片 幻灯片

5.   结论

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