一种大规模天线阵中基于人工噪声的窃听编码方法与流程

本发明属于无线通信技术领域，具体涉及一种大规模天线阵中基于人工噪声的窃听编码方法。

背景技术：

无线通信的流量经历了前所未有的增长，预计将从2016年到2021年以47％的复合年增长率(cagr)增长，到2020年，每月全球移动数据流量将达到30.6艾字节。越来越多的隐私数据将通过无线网络传输，无线传输介质的电磁特性和开放特性会给这些数据的安全通信带来威胁。

据申请人了解，传统的安全通信方法是通过复杂的数学运算来加密信号，寄希望于窃听者没有足够的运算能力。但随着智能硬件以及物联网的发展，这种需要密钥交换的加密方式已无法保证通信的安全。

物理层安全是一种不需要上层加密算法就能获得安全通信的信息理论方法，它利用了无线信道的物理特性，比如衰落、噪声、信道干扰等。常见的信道编码技术如维纳窃听编码、低密度奇偶校验码、极性码等在保障通信可靠性的同时也保证了安全性。

常见的物理层安全方法有波束成形法，通过调节天线阵列因子使得发射波束聚焦到合法接收者的方向保证了通信安全；天线选择技术，通过在众多的发射天线中选择部分天线发射，增加了随机性，造成了窃听方向用户的星座图的扰乱。

前述的方法都是主动增加合法接收者的信道质量，除此之外，还可以通过分配一定的功率给噪声，联合有用的信号一起发射出去，达到故意干扰窃听者的目的。而在众多的人工噪声方法中，有些方法在确定传输参数时需要过于复杂的计算。有些方法的性能比较低。

技术实现要素：

本发明的目的在于：提供一种大规模天线阵中基于人工噪声的窃听编码方法，利用严格的数学公式推导出安全吞吐量与冗余速率以及功率分配因子的关系，在安全中断概率的约束下，利用黄金分割搜索法求出最优的冗余速率以及功率分配因子，进而使得安全吞吐量达到最大。

为了达到以上目的，提供一种大规模天线阵中基于人工噪声的窃听编码方法，包括如下步骤，

发送者与接受者之间建立主信道，并与窃听者之间建立窃听信道；

根据信道估计获取主信道csi以及主信道和窃听信道的发射信噪比；

根据主信道和窃听信道的发射信噪比计算出对应的信噪比的累计分布函数；

根据累计分布函数计算出安全吞吐量和安全中断概率的函数，证明对于任意的功率分配因子，都存在唯一的冗余速率使安全吞吐量打到最优；

通过黄金分割搜索法获取最优的冗余速率和分配因子。并通过安全吞吐量的函数求得最大的安全吞吐量。

本发明的优选方法是：发送者的发射信号矢量x由信息信号xi和(n-1)×1的噪声信号xn组成，xn中每个元素的方差为χn，发送者的发射信号功率pt中发射给主信道的功率占比为η(0≤η≤1)，则主信道的信号功率为χi＝ηpt，当同时传输信息信号xi和噪声信号xn时，可得波束赋形矩阵为：

w＝[wiwn]；

其中，wi表示传输的信息信号xi，wn表示传输的噪声信号xn，对矩阵f进行特征值分解f＝h^hh，其中h^h是h的共轭转置，h为主信道信道状态信息，取最大特征值对应的特征向量作为wi，并对wi归一化处理使||wi||²＝1，剩余的n-1个特征向量组成wn,wn的列向量形成h的零空间的正交基，即:hwn＝0，从而可得发送者的发射信号为：

接受者收到的信号表示为：

y＝hx+nb＝hwixi+nb；

其中，nb是加性高斯白噪声，满足e[nbnb^h]＝σb²，则主信道的发射信噪比为：因此，接收者的接收信噪比为||h||为欧几里得范数；

同理，窃听者的接收信号为：

z＝gx+ne＝gwixi+gwnxn+ne；

其中，ne是加性高斯白噪声，满足e[nene^h]＝σe²，窃听信道的发射信噪比为

发射者将人工噪声均匀分配到wn的每个向量上，故而可得：χn＝(1-η)pt/(n-1)；

由上述可得，窃听者的接收信噪比为：

优选地，由于||h||²是n个独立高斯随机变量的平方和，γb服从卡方分布，所以可得接收者的接收信噪比γb的概率密度函数(pdf)为：

则累计分布函数为：

其中，为不完整伽玛函数，γ(n)为伽玛函数；

窃听者的接收信噪比γe的累计分布函数为：

优选地，在窃听信道中，可获得的最大保密速率为：

其中，cb＝log2(1+γb)为主信道的瞬时信道容量，ce＝log2(1+γe)为窃听信道的瞬时信道容量，γb和γe分别表示接收者和窃听者的接收信噪比；窃听编码的速率参数包括总的编码速率rb和保密速率rs，冗余速率re＝rb-rs；

设置rb＝cb，则传输信号的概率为：

根据安全吞吐量的定义：t＝e[rs]；

可得：

其中，η为功率分配因子；

当信噪比高时，cb＝log2(1+γb)可以转化为cb＝log2(γb)，则可得：

根据meijer函数g-function的定义：

转化为：

当ce＝log2(1+γe)大于冗余速率re时发生安全中断时，则：

其中，

优选地，对安全吞吐量t(η,re)计算其关于冗余速率的一阶导数，根据莱布尼兹定理可得：

已知[cb-re]⁺＞0且所以进而得到吞吐量t(η,re)为关于re的单调递减函数。接着对安全中断概率关于re求一阶导数，可得：

在n→∞和两个渐进过程中，可以化简出最优冗余速率的两个闭合表达式，

当n→∞时，

关于η的最优冗余速率表达式为：

当时，

关于η的最优冗余速率表达式为：

本发明有益效果为：该方法将信道估计后的矩阵进行特征值分解，求出各个特征值对应的特征向量，形成波束成形矩阵，接着求出合法接收者和窃听者的信噪比的累积分布函数。再求出安全吞吐量和安全中断概率关于冗余速率的函数，证明对于任意的功率分配因子，都存在唯一的冗余速率使得安全吞吐量达到最大，通过黄金分割搜索法求得最优的冗余速率和分配因子，在传输的过程中保持不变，达到最优的安全吞吐量。

附图说明

下面结合附图对本发明作进一步的说明。

图1为本发明的示意图。

图2为本发明的流程图。

图3为本发明的方法流程示意图。

图4为本发明的系统结构示意图。

具体实施方式

实施例一

请参阅图1-3，本实施例的一种大规模天线阵中基于人工噪声的窃听编码方法，包括如下步骤，

发送者与接受者之间建立主信道，并与窃听者之间建立窃听信道；

根据信道估计获取主信道csi以及主信道和窃听信道的发射信噪比；

根据主信道和窃听信道的发射信噪比计算出对应的信噪比的累计分布函数；

根据累计分布函数计算出安全吞吐量和安全中断概率的函数，证明对于任意的功率分配因子，都存在唯一的冗余速率使安全吞吐量达到最优；

通过黄金分割搜索法获取最优的冗余速率和分配因子。并通过安全吞吐量的函数求得最大的安全吞吐量。

如图4所示的系统模型。发射端配置多根天线的系统，发送者alice配置n根天线，合法接收者bob以及窃听者eve均配置一根天线。当alice向bob发送信号时，由于无线信号的电磁特性，也会把一部分信息传给eve，eve尝试恢复出原始信息。

alice与bob之间建立主信道h，而alice与eve之间建立窃听信道g，都用一个1×n的复向量表示，是平坦瑞丽慢衰落信道。h和g中的每一个元素都是独立同分布的高斯随机变量。考虑到被动窃听的场景，alice不能获得g的瞬时信息但是可以知道h的瞬时值。同时，假设alice可以得到g的一些统计信息，可以把eve看作系统前一个时隙的合法的多用户之一，因此也会反馈它的一些信道信息以及热噪声的大小，假设eve可以估计出bob到alice之间的信道反馈，索性直接认为eve知道的所有信息。因此，alice需要对信号进行编码处理之后再发送，本发明的做法是在原始信号中加入人工噪声一起发送。eve只是被动地窃听alice发往bob的信号，不会对主信道造成干扰。在miso窃听信道中，可获得的最大保密速率为：

其中cb＝log2(1+γb)是主信道的瞬时信道容量，ce＝log2(1+γe)是窃听信道的瞬时信道容量，γb和γe分别表示bob和eve的接收信噪比(snr)。窃听编码方法中有两个主要的速率参数，一个是总的编码速率rb，还有一个是保密速率rs，他们的差值re＝rb-rs就是用于干扰窃听者的冗余速率。当主信道的信道容量cb＞rb时，可以解码出原始信号；当窃听信道容量ce＜re时，窃听者不能解码出加密信息。因此，一个理想的做法是，保证re比ce稍微大一点，rb比cb稍微小一点。这样的话加密信息就可以最大保密速率cs传输，同时可以保证合法接收者可以译码，窃听者不能解码。因此，本传输方法是：提前计算好一个最优冗余速率，并且在传输的过程中保持恒定不变，而且保证冗余速率和保密速率的总和稍稍小于信道容量cb。如果最后一个条件不能满足则发射端停止发射信息。

人工噪声的添加中，把发射信号矢量用x表示，它由信息信号xi和(n-1)×1的噪声信号xn组成，信号的方差为χi，xn中每个元素的方差为χn。假设alice总的发射功率为pt，其中发射给主信道的功率占比为η(0≤η≤1)，所以有信号的功率χi＝ηpt。为了同时传输信息信号和噪声信号，设计一个n×n的波束赋形矩阵如下：

w＝[wiwn]

其中wi用来传输信息信号xi，wn用来传输人工噪声xn。设计wn的目的是用来干扰窃听者的，同时要求不对合法接收者造成干扰。因此，对矩阵f＝h^hh进行特征值分解(h^h是h的共轭转置)，取其最大的特征值对应的特征向量作为wi，而且对wi做归一化处理使得||wi||²＝1。剩下的n-1个特征向量组成wn，因此wn的列向量形成h的零空间的正交基，也即hwn＝0，而且w是一个酉矩阵。因此，alice的发射信号可以表示为：

所以接收端bob收到的信号可以表示为：

y＝hx+nb＝hwixi+nb

nb是加性高斯白噪声，满足e[nbnb^h]＝σb²。定义主信道的发射信噪比为所以bob的接收信噪比为||h||是欧几里得范数。

同理，eve的接收信号可以表示为：

z＝gx+ne＝gwixi+gwnxn+ne

ne是加性高斯白噪声，满足e[nene^h]＝σe²，窃听信道的发射信噪比为由上述讨论可以，和在发射端都是已知的。由于alice不知道g的瞬时值，所以他把人工噪声均匀地分配到wn的每个向量上，因此有χn＝(1-η)pt/(n-1)。

由上面的讨论可知：尽管eve知道h的信息，但是仍然不能消除wnxn带来的干扰，所以eve的接收信噪比为：

接下来求出γb和γe的累积分布函数(cdf)便于后续的分析。因为||h||²是n个独立高斯随机变量的平方和，所以γb服从卡方分布，所以可得其概率密度函数(pdf)为：

累积分布函数(cdf)为：

其中是上不完整伽玛函数，γ(n)是伽玛函数。

因为h和g分别代表两条相互独立的信道状态信息，所以它们之间是相互独立的，而由上面的分析可知，波束赋形矩阵w完全是依赖于h的值产生的，所以g和w也是相互独立的。又因为g中的元素是均值为零的独立同分布复高斯随机变量，w是酉矩阵，所以gw与g的有着同样的分布，其中的每一个元素都是独立同分布的零均值独立复高斯随机变量。因此γe的累积分布函数可以推导如下：

本传输方法中，rb被设置为rb＝cb，re的值根据接收信号信噪比提前设置好以满足安全传输的条件，并且在整个传输过程中保持不变。如果不满足cb＞re的条件则停止传输。所以本方法中，传输信号的概率是：

安全吞吐量的定义：

t＝e[rs]

所以该方法的安全吞吐量为：

其中，[x]⁺＝max{0,x}。

信噪比较高时，cb＝log2(1+γb)可以写成cb＝log2(γb)，可以得到：

由meijer函数g-function的定义：

整理后可得：

对于安全中断概率。根据熟知的维纳窃听编码理论，当窃听信道的信道容量ce＝log2(1+γe)大于冗余速率re时发生安全中断，可得：

其中：

本传输方法中，只有rb是依据主信道瞬时信道状态改变的，re和η都是根据主信道和窃听信道的信噪比提前计算好并保持恒定的。它们的最优值在满足一定的安全性限制的条件下，安全吞吐量最大时取得。也即：(η^*,re^*)使得：

maxt(η,re)

s.t.pso(η,re)≤p0

下面求取最优传输参数：

对安全吞吐量求其关于冗余速率的一阶导数，对求导，根据莱布尼兹定理可知：

已知[cb-re]⁺＞0且所以进而得到吞吐量t(η,re)是关于re的单调递减函数。接着对安全中断概率关于re求一阶导数，可得：

在大规模miso窃听信道中，n＞1，是大于零的，因此是小于零的，所以中断概率也是关于re的单调递减函数。因此，对于给定的一个功率分配因子η，存在唯一的re使得在满足安全条件约束下，吞吐量达到最大值，而且在安全约束边界处达到，也即pso(η,re)＝p0。因此，对于一个给定的η，最优冗余速率re^*的取值可以这样得到：从零开始慢慢增加re的值，直至最终使得pso(η,re)＝p0。因此可以用黄金分割搜索法获得最优解。将最优传输参数代回到可以得到最大的安全吞吐量。特别地，n→∞和两个渐进过程中，可以化简出最优冗余速率的两个闭合表达式。如下：

当n→∞时，可以写成：

因此关于η的最优冗余速率表达式为：

当时，可以写成：

因此关于η的最优冗余速率表达式为：

除上述实施例外，本发明还可以有其他实施方式。凡采用等同替换或等效变换形成的技术方法，均落在本发明要求的保护范围。