利用功率谱极值和几何平均的频谱感知算法

韩仕鹏 赵知劲,2 毛翊君

(1. 杭州电子科技大学通信工程学院，浙江杭州 310018；2. 中国电子科技集团第36所研究所通信系统信息控制技术国家级重点实验室，浙江嘉兴 314001)

1 引言

频谱感知是认知无线电(Cognitive Radio，CR)的关键技术之一，根据所采用的信号处理技术可以分为：能量检测法(Energy detection，ED)、匹配滤波检测法(Matched-Filtering，MF)、循环特征检测法(Cyclostationary Feature Detection，CFD)、基于随机矩阵理论(Random Matrix Theory, RMT)的盲检测法、功率谱检测法等方法，其中，ED、MF和CFD是经典频谱感知算法[1]。

由于经典频谱感知算法容易受到阴影衰落、噪声不确定性和隐藏终端的影响，文献[2]提出了一种基于RMT的盲频谱感知方法，它不需要主用户(Primary User, PU)的任何先验知识。但是基于RMT的频谱感知方法需要进行复杂的特征值分解，且只能给出渐进而非准确的判决门限，因此难以控制算法的检测概率和虚警概率。而基于功率谱的频谱感知算法不需要特征值分解，复杂度低，且能准确给出算法的判决门限。文献[3-5]采用功率谱分段对消构造检测统计量，降低了信号功率谱成份，对功率谱最小值估计精度较低。为了提高频谱感知性能和降低信号载波频偏[6-7]对功率谱极值估计的影响，文献[8]提出了基于功率谱的平均比值算法(Power Spectral Ratio Averaging，PSRA)，文献[9]提出了将功率谱最大值的平均值与最小值的平均值之比作为检测统计量的检测算法(Power Spectral Maximum and Minimum Averaging Ratio，PSMMAR)，PSMMAR算法抗频偏性能优于PSRA算法，但是其功率谱最大值平均估计中包含噪声功率分量，算法的频谱感知性能有待进一步提高。对此，本文提出了一种利用功率谱的最大值与最小值的平均值之差与功率谱的几何平均之比作为判决统计量的频谱感知算法(Power Spectral Extremum Geometric Averaging Radio，PSEGAR)。通过功率谱差值可消除噪声功率分量，而几何平均受极值的影响比算术平均小，因此该频谱感知算法具有更好的检测性能。

2 算法提出

2.1 算法理论分析

用H0和H1分别表示PU信号不存在和存在的假设，则认知无线电中，次用户(Secondary User, SU)感知PU存在与否问题可表示成如下二元假设检验：

(1)

由于x(n)和ν(n)相互独立，则式(1)的二元假设检验模型可用功率谱表示如下：

(2)

(3)

由文献[10]可知，当采样点数N较大时，H0假设下，Ry的最大特征值λmax(=λ1)为：

(4)

(5)

本文讨论单个CU的频谱感知，因此B=1。将式(4)代入式(5)，当采样样本数N远大1时，可得：

(6)

(7)

2.2 判决规则和门限推导

yi(n)(i=0～L-1,n=0～M-1)

则平滑处理后信号y(n)的离散功率谱估计为：

(8)

(9)

(10)

(11)

因此可得本文提出的频谱感知算法判决统计量计算如式(12)所示，判决规则如式(13)所示。

(12)

(13)

(14)

(15)

(16)

(17)

则有

(18)

根据式(13)、式(16)和式(18)可得虚警概率Pf表达式为：

(19)

根据式(13)、式(17)和式(18)，类似式(19)的推导，可得检测概率Pd表达式为：

(20)

其中，erfc(·)为误差函数。给定虚警概率Pf，可得判决门限为：

(21)

从式(21)可以看出，判决门限γ与噪声ν(n)的方差无关，因此理论上本文算法对噪声不确定性是鲁棒的。

3 算法仿真与性能分析

本节仿真分析本文算法、PSMMAR算法[9]、PSRA算法[8]和PSC算法[4]的抗噪声不确定性和抗载波频偏性能。仿真中采用QPSK信号作为PU信号，过采样倍数即每个码元的采样点数p=8，码元数为500，虚警概率Pf=0.01。

(1)K0的选取

信噪比SNR=-15 dB，频偏Δf=0.04、0.05和0.06，K1=4，M=16，K0=1,...,8时，得到四种算法的检测概率与K0的关系曲线如图1所示。由图1可知，K0值取2时，本文算法和PSMMAR算法的性能较好；而PSC算法和PSRA算法的性能不受K0值影响。因此，下文中取K0=2。

图1 不同K0值时算法检测概率曲线Fig.1 Algorithms detection probability curves with different K0 values

(2)K1的选取

频偏系数Δf=0，信噪比SNR=-15 dB，K0=2和M=16时，四种算法的检测概率与K1的关系曲线如图2所示。由图2可知，K1在一定范围内取值越大，算法性能越好；当K1取值足够大时，算法性能就不再提高。这是因为当K0和K2确定后，K1越大，用于最小值估计的样本数越多，估计越准确；当样本数足够多时，算法性能就不再受K1值影响。因此，下文中取K1=4。

图2 不同K1值时算法检测概率曲线Fig.2 Algorithms detection probability curves with different K1 values

(3)M的选取

频偏系数Δf=0，信噪比SNR=-15 dB，K0=2和K1=4，M=16,24,...,48时四种算法的检测概率与帧长M的关系曲线图如图3所示。由图3可知，随着帧长M值增大，四种算法的检测性能逐渐降低。这是因为帧长增加，样本帧数减少，功率谱平滑效果降低。因此，下文中取M=16。

图3 不同M值时算法检测概率曲线Fig.3 Algorithms detection probability curves with different M values

(4)算法抗频偏性能

频偏Δf=0.04、0.05和0.06，K0=2，K1=4，四种算法的检测概率与信噪比的关系曲线如图4所示。

图4 四种算法抗频偏性能Fig.4 Anti-frequency offset performance of algorithms

由图可见，本文算法性能最好, PSMMAR算法次之， PSRA算法第三，PSC算法最差；当信噪比为-15 dB、载波频偏系数由0.04增大到0.06时，本文算法的检测概率由0.95下降到0.91，PSMMAR算法的检测概率由0.89下降到0.84，PSRA算法的检测概率由0.54下降到0.19，PSC算法的检测概率由0.32下降到0.10。由图可见，随着载波频偏的增大，本文算法性能恶化相对其他几种算法是最小的， PSMMAR算法次之，而PSRA算法和PSC算法的性能均有较明显的下降。因此本文算法具有更优的抗频偏性能。

(5)算法抗噪声不确定性性能

由图可见，在两种噪声不确定性情况下，本文算法性能都要优于其他三种算法，其中本文算法和PSMMAR算法都具有较强的抗噪声不确定性能力。当信噪比较小时，Rayleigh衰落信道中四种算法的检测概率要大于其在AWGN信道中的检测概率，因为此时Rayleigh衰落信道的起伏增大了信号能量，从而提高了信噪比。

从图4和图5(a)可见，当信噪比低于-20 dB时本文算法性能略优于其他三种算法，且四种算法的检测概率曲线逐渐趋于重叠。这是因为信噪比太低，主用户信号已淹没在噪声中，不同算法的检测统计量反映的主用户功率谱信息差别缩小，最后趋于都是噪声功率谱信息。

图5 两种信道下四种算法的抗噪声不确定性性能Fig.5 Anti-noise uncertainty performance of four algorithms under two channels

(6)ROC曲线

信噪比SNR=-12 dB，其他参数同(5)，四种算法在AWGN信道和Rayleigh衰落信道中的ROC(Receiver Operating Characteristics)曲线分别如图6(a)和图6(b)所示。

由图可见，在AWGN和Rayleigh衰落信道中，本文算法的检测概率值高于其他算法，当ρ=1.1时本文算法的检测概率虚线的起伏相对其他算法更小。因此本文算法在AWGN和Rayleigh衰落信道下具有更优的检测性能和更强的抗噪声不确定性。

图6 两种信道下不同算法的ROC曲线Fig.6 ROC curves of different algorithms under two channels

(7)实际信号检测

利用基于AD9361的软件无线电平台采集通过空中传输的QPSK信号。实验中在52 MHz带宽内共划分42个信道，随机选择4个信道产生主用户信号，单个信道带宽1.2 MHz。每个主用户信号码速率为0.2 MHz，升余弦成形滤波器滚降因子为0.5。针对产生的实际信号，算法参数取L=15，K0=20和K1=20。实验中产生的主用户信号所在信道编号及中心频率如表1所示，接收的宽带基带信号功率谱如图7所示。每个信道应用PSMMAR算法和本文算法计算得到的检测统计量与各自门限的比值，分别记为RPSMMAR和RPSEGAR，如图8和图9所示。

表1 主用户信道编号及信道中心频点分布

图7 接收信号功率谱图Fig.7 Received signal power spectrum

图8 PSMMAR算法判决统计量与门限比值Fig.8 The ratio of PSMMAR algorithm decision statistic to threshold

图9 本文算法判决统计量与门限比值Fig.9 The ratio of PSEGAR algorithm decision statistics and threshold ratio

由图可见，两种算法主用户所在信道的RPSEGAR和RPSMMAR都高于1，可检测到主用户存在，且RPSEGAR远大于1，而RPSMMAR略大于1；主用户不存在信道的RPSEGAR值明显小于主用户所在信道，而RPSMMAR的二者差别很小，本文算法容易检测主用户。由此可知本文算法检测实际信号性能优于PSMMAR算法。

4 结论

提出了一种利用功率谱极值和几何平均的频谱感知算法，在AWGN和Rayleigh衰落信道下，该算法频谱感知性能优于已有的基于功率谱的频谱感知算法，抗未知载波频偏和抗噪声不确定性能力强，并且能够较好地检测实际接收信号中主用户信号。