1.一种基于粘性隐马尔可夫模型的协同频谱感知方法，其特征在于包括以下步骤：步骤1：在认知无线电系统中设定存在一个主用户和Num个次级用户，从时域的角度分析，仅当主用户空闲时，次级用户才能动态地接入主用户的授权频段；从空域的角度分析，位于主用户的通信范围以外的次级用户使用主用户的授权频段不会影响主用户的正常通信；然后结合时域和空域两方面，构建频谱感知模型，描述为：其中，Num表示认知无线电系统中的次级用户的个数，Num≥2，频谱感知模型实际为一个二元假设检验问题，S0代表在一个感知时隙内次级用户能够接入主用户的授权频段的情况，ri(n)表示第i个次级用户接收的信号的第n个样本值，1≤i≤Num，1≤n≤N，N表示对每个次级用户接收的信号进行采样得到的样本数，N≥25，wi(n)表示ri(n)中包含的噪声，wi(n)服从均值为0、方差为 的复高斯分布，ri(n)＝wi(n)代表主用户空闲时第i个次级用户动态地接入主用户的授权频段的情况，gi表示第i个次级用户的信道增益，gi服从自由空间损耗模型，hi(n)表示第i个次级用户的信道系数，hi(n)服从瑞利分布，si(n)表示ri(n)中包含的主用户发送的信号，si(n)服从均值为零、方差为 的复高斯分布，di表示第i个次级用户到主用户的欧氏距离，rpu表示主用户的通信范围，di＞rpu表示第i个次级用户位于主用户的通信范围以外，di＞rpu代表位于主用户的通信范围以外的第i个次级用户在不影响主用户的正常通信的情形下使用主用户的授权频段的情况，S1代表在一个感知时隙内主用户发送信号的情况下次级用户不能够接入主用户的授权频段的情况，di≤rpu表示第i个次级用户位于主用户的通信范围以内， di≤rpu代表位于主用户的通信范围以内的第i个次级用户不能够接入主用户的授权频段的情况；

步骤2：针对连续的Tslot个感知时隙，计算在每个感知时隙内各个次级用户接收的信号的测量功率，将在第t个感知时隙内第i个次级用户接收的信号的测量功率记为yt,i，然后将在第t个感知时隙内Num个次级用户接收的信号的测量功率构成的T

功率向量记为yt，yt＝[yt,1,…,yt,i,…,yt,Num] ；再将在Tslot个感知时隙内Num个次级用户接收的信号的测量功率构成的功率矩阵记为Y， 其中，Tslot≥4，1≤t≤Tslot，符号“||”为求模运算符号，rt,i(n)表示在第t个感知时隙内第i个次级用户接收的信号的第n个样本值，yt,1表示在第t个感知时隙内第1个次级用户接收的信号的测量功率，yt,Num表示在第t个感知时隙内第Num个次级用户接收的信号的测量功率，上标“T”表示向量或矩阵的转置，y1表示在第1个感知时隙内Num个次级用户接收的信号的测量功率构成的功率向量，yTslot表示在第Tslot个感知时隙内Num个次级用户接收的信号的测量功率构成的功率向量，y1、yt、yTslot的维数为Num×1，Y的维数为Num×Tslot；

步骤3：将粘性隐马尔可夫模型的三个模型参数即状态转移概率矩阵、初始概率分布、状态参数集合对应记为A、π、θ，θ＝{θ1,θ2}；然后对A、π、θ1,θ2进行初始化，A的初始化值为π的初始化值为[0.5 0.5]，θ1,θ2的初始化值为随机数；再将Y视为粘性隐马尔可夫模型的观测序列；其中，θ1表示第1类状态的参数，θ2表示第2类状态的参数，μ1表示第1类状态的均值， 表示第1类状态的方差，μ2表示第2类状态的均值， 表示第2类状态的方差；

步骤4：在Y和A、π、θ已知的情况下，依次求取各个感知时隙的隐藏状态，将第t个感知时隙的隐藏状态记为zt， 再将Tslot个感知时隙的隐 藏 状态构 成隐藏 状态向 量 ，记为Z ， 其中 ，表示取使得

1:t

的值最大时的zt的值，p(zt,Y |A,π,θ)表示在A、π、θ

1:t

已知的情况下zt、Y 的概率， 表示在A、π、θ已知的情况下zt、 的概

1:t

率，Y 表示在第1个感知时隙至第t个感知时隙内Num个次级用户接收的信号的测量功率构

1:t

成的功率矩阵，Y 亦表示粘性隐马尔可夫模型的观测序列Y中的第1个至第t个观测值，表示在第t+1个感知时隙至第Tslot个感知时隙内Num个次级用户接收的信号的测量功率构成的功率矩阵， 亦表示粘性隐马尔可夫模型的观测序列Y中的第t+1个至第Tslot个观测值，zt的值为1或2，z1表示第1个感知时隙的隐藏状态， 表示第Tslot个感知时隙的隐藏状态，Z的维数为1×Tslot；

* * *

步骤5：在Y和Z已知的情况下，更新A和θ，将A和θ更新后的值对应记为A 和θ，将A中第j*

行第j'列的元素的值记为 将θ中的第j*

类状态的均值和第j类状态的方差对应记为μj和其中，1≤j≤2，1≤j'≤2，κ为粘性隐马尔可夫模型的粘性因子，δ()为克罗内克函数，当j＝j'时δ(j,j')的值为1，当j≠j'时δ(j,j')的值为0，当j＝1时δ(j,1)的值为1，当j≠1时δ(j,1)的值为0，当j＝2时δ(j,2)的值为1，当j≠

2时δ(j,2)的值为0，Nj,j'表示Tslot个感知时隙中相邻的两个感知时隙的前一个感知时隙的隐藏状态为j而下一个感知时隙的隐藏状态为j'的发生次数，Nj,1表示Tslot个感知时隙中相邻的两个感知时隙的前一个感知时隙的隐藏状态为j而下一个感知时隙的隐藏状态为1的发生次数，Nj,2表示Tslot个感知时隙中相邻的两个感知时隙的前一个感知时隙的隐藏状态为j而下一个感知时隙的隐藏状态为2的发生次数， 表示 中元素的个数， 表示 中的所有元素的值的平均值， 表示由Y中位置与Z中值为j的隐藏状态的位置对应的所有的观测值构成的集合；

* *

步骤6：令A＝A 、θ＝θ；然后返回步骤4继续执行，在迭代执行至少10次后将得到的隐藏final

状态向量作为最终的隐藏状态向量，记为Z ， 其中，A＝* * final

A和θ＝θ中的“＝”为赋值符号， 表示第1个感知时隙的最终隐藏状态，zt 表示第t个感知时隙的最终隐藏状态， 表示第Tslot个感知时隙的最终隐藏状态；

步骤7：计算每个次级用户的检验统计量，将第i个次级用户的检验统计量记为然后将每个次级用户的检验统计量与门限值进行比较来判断该次级用户能否接入主用户的授权频段，对于第i个次级用户，若 大于或等于门限值，则表示第i个次级用户能够接入主用户的授权频段；若 小于门限值，则表示第i个次级用户不能够接入主用户的授权频段；其中，1≤b≤Num，ωb表示第b个次级用户的权重， e为自然基数，β0为融合参数，β0可调， 表示向量 中的第b个元素， 表示向量 中的第i个final

元素，表示获得Z 时对应的θ中的第 类状态的均值，门限值是根据给定的虚警概率计算得到的。

2.根据权利要求1所述的一种基于粘性隐马尔可夫模型的协同频谱感知方法，其特征在于所述的步骤5中的κ取值为0.1。

3.根据权利要求1或2所述的一种基于粘性隐马尔可夫模型的协同频谱感知方法，其特征在于所述的步骤7中β0取值为2.0。