-
认知无线电(cognitive radios, CRs)主要用于克服无线通信频谱资源的短缺问题,使次用户在不干扰主用户的情况下获得对空闲频谱资源的使用。CRs中谱感知方式可通过非协作或协作的方式实施,本文研究的是一群次用户通过协作的方式实施频谱感知的情况。协作方式可以避免和改善影深和遮掩衰落对检测结果的影响,另外空间分布的用户具有更好的区域覆盖性,有助于提高准确度。
CRs中频谱感知[1]主要有3种方式:匹配滤波理论最优,但需要主用户先验信息;能量检测次优,它简单易实施,对主用户没有太多要求;循环平稳特征检测方法能够检测信噪比(signal to noise ratio, SNR)很低的信号,但它也需要主用户的先验知识[2]。文献[3]提出了分布协作感知,通过分布的次用户协作达到共识并最终判决,特别适应网络中无线衰落严重的用户检测,但是并没有考虑和利用网络拓扑资源。随着无线通信的发展,CRs技术和应用得到广泛关注,文献[4]研究了时变衰落信道下认知无线电技术;文献[5]利用CRs技术改善工业环境中干扰的影响;文献[6]研究了CRs协作信息传递时的信息碰撞问题。近年来,网络拓扑可知理论用于改善网络性能,得到关注和研究[7-9]。连通支配集(connected dominating set, CDS)是网络拓扑的重要方面。文献[9]研究了无线ad hoc网络高效路由机制,提出了连通支配集的简单有效计算方法,以及拓扑改变时的动态更新/重算方法。
本文将文献[3]与文献[9]中的方法相结合,在原算法中考虑了网络拓扑的有利因素,提出基于连通支配集的分布式协作共识频谱感知方法。
-
认知无线电协作共识频谱感知分为两个阶段。第一阶段,次用户进行数据测量,次用户i的测量数据表示为Yi。第二阶段,各次用户与邻居用户建立通信链路并进行信息交互,根据获得的数据计算各自的本地结果,此处理过程通过反复迭代进行。在起始时刻k=0,每个次用户i设置xi(0)=Yi作为本地状态变量的初始值。在时刻k=0, 1, 2, …, 根据实时网络拓扑,次用户邻居间交互它们的状态x(k),各次用户计算并生成状态更新值x(k+1),迭代重复执行直到所有个体的状态x(k)收敛到定值x*。
-
在第一阶段,次用户在每一个时隙的开始测量主用户数据,这里主用户的先验知识未知。图 1为能量检测模块图。为方便采用能量检测频谱感知方法[10],使信号输入中心频率为fs、带宽为W的带通滤波器,再通过求平方器件,然后关于时间区间T求积分。积分器输出Y是次用户接收到的能量,其分布情况依赖于主用户信号是否存在。频谱感知的目的就是在两个假设$ {\mathit{\mathcal{H}}_0} $和$ {\mathit{\mathcal{H}}_1} $间进行判决:
$$ x(t) = \left\{ \begin{array}{l} n(t)\;\;\;\;\;\;\;\;\;\;\;\;\;{\mathcal{H}_0}\\ hs(t) + n(t)\;\;\;\;{\mathcal{H}_1}{\rm{ }} \end{array} \right. $$ (1) 式中,x(t)是次用户接收到的信号;s(t)是主用户传输的信号;n(t)是加性高斯白噪声;h是信道的振幅增益。定义SNR为γ,积分器的输出Y作为判决统计量[10]:
$$ Y = \left\{ \begin{array}{l} \chi _{2TW}^2{\rm{\;\;\;\;\;\;\;\;}}{\mathcal{H}_0}\\ \chi _{2TW}^2(2\gamma ){\rm{\;\;\;\;}}{\mathcal{H}_1}{\rm{ }} \end{array} \right. $$ (2) 式中,χ2TW2和χ2TW2(2γ)分别代表中心和偏心卡方分布的随机数,各自具有2TW维自由度,后者的偏心参数为2γ。为了简化,假设时间带宽积TW为整数,用m表示。
在瑞利衰落信道条件下,增益h是随机量,SNR值γ服从指数分布。在此条件下,输出Y的分布依赖于SNR平均值γ。当主用户不存在时,Y依然服从χ2TW2分布;当主用户出现时,Y可以表示为两个独立随机变量和的形式[11]:
$$ Y = {Y_x} + {Y_e}{\rm{\;\;\;\;\;\;\;\;}}{\mathcal{H}_1} $$ (3) 式中,Yx服从χ2TW2分布;Ye是参数为2(γ+1) 的指数分布。经过T秒后,每个次用户i进行能量检测并得到测量数据$ {Y_i} \in {{\rm{\mathbb{R}}}^ + } $。
-
在第二阶段,次用户与邻居建立通信链路进行信息交互。为简化,次用户组成的网络可以表示成无向图G=(N, ℑ)[12],由节点集合N={i=1, 2, …, n}和边集合ℑ⊂N×N构成。每个边表示为无向对(i, j),两节点单步连接。图G中的路径由节点序列$ {i_1}, {i_2}, \cdots, {i_l}(l \ge 2) $构成,那么$ ({i_m}, {i_{m + 1}}) \in \Im $对于所有1≤m≤l-1均成立。图G中任两个节点通过路径连通,称图G是连通图,本文所述的图指连通图。若$ (i, j) \in \Im $且i≠j,表示次用户j(节点j)和次用户i(节点i)是邻居。节点i的邻居表示为$ {N_i} = \{ j|(j, i) \in \Im \} \subset N $。Ni集合的节点数为|Ni|,称为节点i的度。
考虑协作频谱感知问题为共识问题,在达到共识以前次用户交互各自检测的输出结果信息。对于图G模型中的n个节点,对应分配状态变量xi(i∈N)。xi为共识变量,它本质上代表节点i的能量检测值。当它逐渐趋于稳定时,xi收敛到共识固定值x*,即:
$$ {x_i}(k) \to {x^ * }{\rm{\;\;\;\;}}i \in N, \;\;\;\;\;\;k \to \infty $$ (4) 式中,k=0, 1, 2, …是离散时间;而xi(k)是节点i在k时刻的状态,它是由k-1时刻节点i和它的邻居状态值计算获得的状态更新值。协作频谱感知采用平均共识$ {x^ * } = {\rm{Ave}}(x) = (1/n)\sum\limits_{i = 1}^n {{x_i}(0)} $。
-
1) 对于固定连通网络模型,网络拓扑不发生改变的情况
对于节点i,在k=0时,测量值Yi表示为$ {x_i}(0) = {Y_i} \in {{\rm{\mathbb{R}}}^ + } $。各节点在离散时刻k=0, 1, 2, …的状态更新表示为[3]:
$$ {x_i}(k + 1) = {x_i}(k) + \varepsilon \sum\limits_{j \in {{\rm{N}}_i}} {[{x_j}(k)-{x_i}(k)]} $$ (5) 其中,
$$ 0 < \varepsilon < {(\mathop {\max }\limits_i \left| {{N_i}} \right|)^{ - 1}} = 1/\Delta $$ (6) 式中,∆表示网络的最大自由度。最终分布式协作达成共识会收敛到值x*。最后通过比较共识值x*和预先定义的门限值λ,各节点i实现判决:
$$ \mathcal{H} = \left\{ \begin{array}{l} {\mathcal{H}_0} = 0\;\;\;\;{\rm{其他}}\\ {\mathcal{H}_1} = 1\;\;\;\;\;{x^ * }{\rm{ > }}\lambda \end{array} \right. $$ (7) 2) 对于随机网络模型,网络拓扑发生随机改变的情况
设图G=(N, ℑ)表示无链路错误情况下的最大通信链路集。考虑到随机链路错误,在时刻k用图$ G(k) = (N{\rm{, }}\Im (k)) \subseteq G $来描述随机网络。当k时刻节点j和i连通,则边$ (j, i) \in \Im (k) $,其中$ (j, i) \in \Im $。图G各链路错误相互独立,错误概率p∈(0, 1),k时刻节点i的邻居$ {N_i}(k) = \{ j|(j, i) \in \Im (k)\} $,邻居节点数为|Hi|(k)=∅。随机网络中,每个节点按照如下规则进行状态更新:
$$ {x_i}(k + 1) = {x_i}(k) + \varepsilon \sum\limits_{j \in {N_i}(k)} {[{x_j}(k)-{x_i}(k)]} $$ (8) 式中,如果Ni(k)=∅,则$ {x_i}(k + 1) = {x_i}(k) $。当k→∞时,最终分布式协作达成共识,以概率1收敛到x*[3]。
-
分布式协作共识虽然能够达到一致收敛,但没有借助拓扑结构提升共识速度。实际上可以采用全网络节点分布式协作,并在连通支配集上快速形成稳定共识,然后将结果分享给其他非支配集节点来提升网络共识状态的收敛速度,且最终一致收敛结果不变。
连通支配集由连通的支配节点构成,网络所有节点要么在当前连通支配集中,要么是连通支配集的邻居节点[9]。如图 2所示,12个节点随机分布于长宽各为1的正方形区域,节点间的边用连线表示,图中黑色方块是连通支配集节点{6, 8, 9, 10, 11, 12}。连通支配集经过单步可以到达网络中的任意节点,故网络信息通过它可以高效传递。因此将文献[3]的方法与连通支配集相结合,可以借助网络拓扑结构,达到快速达成共识的效果。对于包含链路误差的随机网络模型,文献[9]提出了相应的更新和重算措施。
-
设$ G'{\rm{ = (}}N{\rm{', }}\Im '{\rm{)}} $为网络G的连通支配集,其节点集合$ N' \subset N $,边集合$ \Im ' \subset \Im $。支配集中节点i的邻居$ {N'_i} = \{ j|(j, i) \in \Im '\} \subset N' $。由于通过支配集G'及其邻居节点可以覆盖全网络,所以支配集状态收敛速度比整个网络状态收敛速度快。由于支配集先于网络收敛到稳定值,可以将支配集状态收敛的结果分享给网络邻居节点,各个节点再按式(7) 进行判决,进而完成网络分布协作频谱感知。
定理 1 无向连通图$ G{\rm{ = (}}N{\rm{, }}\Im {\rm{)}} $固定网络,采用式(5) 进行分布式协作状态更新,当k→∞时,对于任何初始状态$ {x_i}(0)(i \in N) $,最终各个节点逐渐达成共识,以概率1收敛到稳定状态值x*。
证明:图G的拉普拉斯矩阵为$ \mathit{\boldsymbol{L}} = {({l_{i, j}})_{n \times n}} $,其中:
$$ {l_{i, j}} = \left\{ \begin{array}{l} \left| {{N_i}} \right|\;\;\;\;\;\;j = i\\ - 1\;\;\;\;\;\;\;j \in {N_i}\\ 0\;\;\;\;\;\;\;\;\;{\rm{其他}} \end{array} \right. $$ (9) 式(9) 表示矩阵L为半正定。进而,式(5) 的矢量形式为:
$$ \mathit{\boldsymbol{x}}(k + 1) = \mathit{\boldsymbol{Px}}(k) $$ (10) 式中,P=I-εL。设P第i行j列的元素为pi, j,xi(k)是x(k)的第i行元素,易知ε的取值确保P是双随机矩阵(即其各行、各列值的和为1),且各主对角线元素pi, i>0。
式(10) 描述了一个齐次马尔可夫链模型,P为网络各节点状态的一步转移概率矩阵。一步转移概率pi, j只与节点i, j, k有关,与时刻k无关。设图G中节点i, j, k构成一个路径,且$ (i, j) \in \Im $,$ (j, k) \in \Im $,$ (i, k) \notin \Im $,即节点i通过节点j与节点k相连。则$ 0 < {p_{i, i}}, {p_{i, j}}, {p_{j, j}}, {p_{j, k}} < 1 $且pi, i=0。容易得到二步状态转移为:
$$ \mathit{\boldsymbol{x}}(k + 2) = {\mathit{\boldsymbol{P}}^{(2)}}\mathit{\boldsymbol{x}}(k) = \mathit{\boldsymbol{PPx}}(k) = {\mathit{\boldsymbol{P}}^2}\mathit{\boldsymbol{x}}(k) $$ (11) 式中,P(2)代表二步转移概率矩阵,其第i行j列元素为:
$$ {p_{i, k}}^{(2)} = \sum\limits_{m \in N} {{p_{i, m}}{p_{m, k}}} \ge {p_{i, i}}{p_{i, k}} + {p_{i, k}}{p_{k, k}} > 0 $$ (12) 另:
$$ {p_{i, j}}^{(2)} = \sum\limits_{m \in N} {{p_{i, m}}{p_{m, j}}} \ge {p_{i, i}}{p_{i, j}} + {p_{i, j}}{p_{j, j}} > 0 $$ (13) 由于P各行、各列和为1,故对于P(2)各行、各列求和也为1,如下:
$$ \left\{ \begin{array}{l} \sum\limits_{k = 1}^n {{p_{i, k}}^{(2)}} = \sum\limits_{k = 1}^n {\sum\limits_{m = 1}^n {{p_{i, m}}{p_{m, k}}} } = \sum\limits_{m = 1}^n {{p_{i, m}}} \sum\limits_{k = 1}^n {{p_{m, k}}} = 1\\ \sum\limits_{i = 1}^n {{p_{i, k}}^{(2)}} = \sum\limits_{i = 1}^n {\sum\limits_{m = 1}^n {{p_{i, m}}{p_{m, k}}} } = \sum\limits_{m = 1}^n {{p_{m, k}}} \sum\limits_{i = 1}^n {{p_{i, m}}} = 1 \end{array} \right. $$ (14) 同时根据$ 0 < {p_{i, i}}, {p_{i, j}}, {p_{j, j}}, {p_{j, k}} < 1 $,参照式(14) 可推知n>1步转移概率为:
$$ 0 < {p_{i, i}}^{(n)}, {p_{i, j}}^{(n)}, {p_{j, j}}^{(n)}, {p_{j, k}}^{(n)} < 1 $$ (15) 且如果节点k与节点最短路径为l步,则$ {p_{i, k}}^{(l)} > 0 $。如果网络节点间最长的最短路径为l'步,则转移概率矩阵P(l')的所有元素都大于0,同时P(l')各行、各列值求和为1。由此可知,矩阵P(l')符合马尔可夫链遍历定理要求,存在极限分布:
$$ \mathop {\lim }\limits_{k \to \infty } {p_{i, j}}^{(k)} = {p_j} $$ (16) $$ {\mathit{\boldsymbol{P}}^{(k)}}\overrightarrow {\;\;k \to \infty \;\;} \left[{\begin{array}{*{20}{c}} {{p_1}}&{{p_2}}& \cdots &{{p_n}}\\ {{p_1}}&{{p_2}}& \cdots &{{p_n}}\\ \vdots & \vdots &{}& \vdots \\ {{p_1}}&{{p_2}}& \cdots &{{p_n}} \end{array}} \right] $$ (17) 由无向图G可知P是对称矩阵,故而k步转移概率矩阵P(k)也是对称矩阵,且各行、各列和为1。因此可以推出$ {p_1} = {p_2} = \cdots = {p_n} = 1/n $,有:
$$ \mathop {\lim }\limits_{k \to \infty } {\mathit{\boldsymbol{P}}^{(k)}}\mathit{\boldsymbol{x}}(k) = {x^ * } \cdot 1 $$ (18) 不同于文献[3]中的算法,这里各节点按照连通支配集分享的收敛值进行判决,进而有效地利用拓扑资源提升算法性能。连通支配集收敛值作为全网络的收敛值分享给周围非支配集的邻居节点,最后根据式(7) 进行判决。
-
随机网络节点间的链路$ (i, j) \in \Im $以概率p∈(0, 1) 失败,各个链路相互独立。k时刻的随机网络图表示为$ G(k) = (N{\rm{, }}\Im (k)) $。图G(k)的支配集为$ G{\rm{'}}(k) = {\rm{(}}N{\rm{'}}(k){\rm{, }}\Im {\rm{'}}(k){\rm{)}} $,其中N'(k)和$ \Im {\rm{'}}(k) $分别为支配集的节点集合和边集合,文献[9]介绍了支配集节点的更新和重新计算过程。
定理 2 无向随机连通图$ G(k) = (N{\rm{, }}\Im (k)) $,节点间链路$ (i, j) \in \Im $以概率p∈(0, 1) 失败,各个链路相互独立,用式(8) 进行分布式协作状态更新。当k→∞时,对初始状态$ {x_i}(0)(i \in N) $,最终各节点达成共识,以概率1收敛到稳定状态值x*。
证明:随机图G(k)的拉普拉斯随机矩阵$ \mathit{\boldsymbol{L}}(k) = {({l_{i, j}}(k))_{n \times n}} $,其中:
$$ {l_{i, j}}(k) = \left\{ \begin{array}{l} \left| {{N_i}(k)} \right|\;\;\;\;j = i\\ 0{\rm{或}} - 1\;\;\;j \in {N_i}(k);(0, 1){\rm{分布以概率}}p{\rm{取}}0\\ 0\;\;\;\;\;\;\;\;\;\;\;\;{\rm{其他}} \end{array} \right. $$ (19) 由式(19) 可知矩阵L(k)是半正定的随机矩阵。进而,式(8) 的矢量的形式为:
$$ \mathit{\boldsymbol{x}}(k + 1) = \mathit{\boldsymbol{P}}(k)\mathit{\boldsymbol{x}}(k) $$ (20) 式中,$ \mathit{\boldsymbol{P}}(k) = \mathit{\boldsymbol{I}} - \varepsilon \mathit{\boldsymbol{L}}(k) $,设pi, j(k)表示P(k)第i行j列元素,xi(k)是x(k)的第i行元素。易知ε的取值确保了P(k)是双随机矩阵形式的随机过程(即其各行、各列和为1);且保证了其主对角线元素pi, j(k)>0。
由于图G各节点间链路失败与否是独立同分布的,所以P(k)是一个随机过程。进而体现随机链路连接情况的P(k)统计特性不随时间改变而改变,所以P(k)是严平稳随机过程。P(k)在各个时刻的期望相等,可表示为双随机常数矩阵$ {\mathit{\boldsymbol{P}}_\mu }{\rm{ = E[}}\mathit{\boldsymbol{P}}(k)] $。P(k)的各元素pi, j(k)的期望为:
$$ {\rm{E}}[{p_{i, j}}(k)] = \left\{ \begin{array}{l} \left| {{N_i}(k)} \right|(1 - p)/\varepsilon \;\;\;j = i\\ (1 - p)/\varepsilon \;\;\;{\rm{ }}j \in {N_i}(k);(0, 1){\rm{分布以概率}}p{\rm{取}}0\\ 0\;\;\;\;\;\;\;\;\;\;\;\;\;\;{\rm{其他}} \end{array} \right. $$ (21) 且P(k)在各时刻k=1, 2, …时独立同分布,根据大数定律,当k→∞时,P(k)的时间平均以概率1等于统计平均Pμ,即P(k)的均值具有各态历经性。
由此可知,当k→∞时,随机图G(k)的平均一步转移概率P(k)为Pμ以概率1成立。进而当k→∞时,随机图G(k)的状态转移可以表示成固定图G以一步转移概率Pμ进行状态转移的齐次马尔可夫链。参照定理1不难得出,此马尔可夫链是遍历的,结合Pμ的双随机性,各个节点逐渐达成共识,以概率1收敛到稳定状态值x*。这说明当k→∞时,随机图G(k)逐步到达状态x*以概率1成立,也即证明了可达性。再参照式(17) 所示其各个元素值$ {p_1} = {p_2} = \cdots = {p_n} = 1/n $,并结合一步转移概率P(k)的双随机性,可知此状态值x*不再随P(k)的变化而发生改变,即证明了稳定性。因此,当k→∞时各节点以概率1收敛到状态值x*。
按照固定无向网络模型类似的方法,随机支配集G'(k)其状态更新收敛速度会比整个网络收敛速度快,提前收敛到x*。因此,这里用支配集G'(k)收敛值x*分享给周围非支配集的邻居节点,网络各节点再参照式(7) 进行判决。
Distributed Consensus Cooperative Spectrum Sensing Method Based on Connected-Dominating-Set
-
摘要: 针对原有全网络分布式协作共识方法信息交互量大、收敛速度慢、共识收敛结果不稳定的问题,本文在分布式协作共识频谱感知方法的基础上,提出了基于连通支配集的分布式协作共识频谱感知方法。该方法通过网络连通支配子集进行网络频谱感知信息的收集和共识计算,获得稳定的共识结果,再将共识结果分享给网络其他非支配集节点,实现全网络快速共识收敛。与原有分布式协作共识方法相比,降低了网络节点间的信息交互量,并且能快速收敛到稳定、精确的共识结果,仿真结果验证了算法的优良特性。Abstract: Aiming at the problems of large amount of information exchange, slow convergence and unstable convergence result of the original whole-network distributed cooperative consensus method, a distributed consensus cooperative spectrum sensing method based on connected-dominating-set is proposed. In this paper, the connected-dominating-set method collects information of network spectrum sensing and does consensus calculation through the network connected dominating subset, which gets a stable consensus result. Then, the consensus results are shared to other non-dominated nodes in the network, so rapid convergence is achieved. There are two advantages of the proposed method compared with the original network distributed cooperative consensus method. One is that the amount of information exchanged among network nodes reduces, the other is that the network can quickly converge to a stable and accurate consensus result. The proof of convergence theorem of distribution consensus is given in this paper. Lastly, the simulation results show the excellent characteristics of this new algorithm.
-
Key words:
- cognitive radios /
- connected dominating set /
- cooperative spectrum sensing /
- consensus
-
[1] CABRIC D, MISHRA S, BRODERSEN R. Implementation issues in spectrum sensing for cognitive radios[C]//The 38th Asilomar Conference on Signals, Systems, and Computers. California:IEEE, 2004:772-776. [2] SUN C, ZHANG W, LETAIEF K B. Cluster-based cooperative spectrum sensing in cognitive radio systems[C]//The 2007 IEEE International Conference on Communications. Glasgow, Scotland:IEEE, 2007:2511-2515. [3] LI Zhi-qiang, RICHARD Y U F, HUANG Min-yi. A distributed consensus-based cooperative spectrum-sensing scheme in cognitive radios[J]. IEEE Transactions on Vehicular Technology, 2010, 59(1):383-392. doi: 10.1109/TVT.2009.2031181 [4] LI B, SUN M, LI X, et al. Energy detection based spectrum sensing for cognitive radios over time-frequency doubly selective fading channels[J]. IEEE Transactions on Signal Processing, 2015, 63(2):402-417. doi: 10.1109/TSP.2014.2368996 [5] CHIWEWE T M, MBUYA C F, HANCKE G P. Using cognitive radio for interference-resistant industrial wireless sensor networks:an overview[J]. IEEE Transactions on Industrial Informatics, 2015, 11(6):1466-1481. doi: 10.1109/TII.2015.2491267 [6] LUNDÉN J, MOTANI M, POOR H V. Distributed algorithms for sharing spectrum sensing information in cognitive radio networks[J]. IEEE Transactions on Wireless Communications, 2015, 14(8):4667-4678. doi: 10.1109/TWC.2015.2424436 [7] BERTRAND A, MOONEN M. Seeing the bigger picture:How nodes can learn their place within a complex ad hoc network topology[J]. IEEE Signal Processing Magazine, 2013, 30(3):71-82. doi: 10.1109/MSP.2012.2232713 [8] BERTRAND A, MOONEN M. Distributed computation of the fiedler vector with application to topology inference in ad hoc networks[J]. Signal Processing, 2013, 93(5):1106-1117. doi: 10.1016/j.sigpro.2012.12.002 [9] JIE W U, HAILAN L I. A dominating-set-based routing scheme in ad hoc wireless networks[J]. Telecommunication Systems, 2001, 18(1-3):13-36. doi: 10.1023/A:1016783217662 [10] URKOWITZ H. Energy detection of unknown deterministic signals[J]. Proceedings of the IEEE, 1967, 55(4):523-531. doi: 10.1109/PROC.1967.5573 [11] KOSTYLEV V. Energy detection of a signal with random amplitude[C]//The 2002 IEEE International Conference on Communications. New York:IEEE, 2002:1606-1610. [12] HUANG M, MANTON J H. Coordination and consensus of networked agents with noisy measurements:Stochastic algorithms and asymptotic behavior[J]. SIAM Journal on Control and Optimization, 2009, 48(1):134-161. doi: 10.1137/06067359X