CN102404063B - 一种基于过采样的glrt检测方法 - Google Patents

Abstract

本发明公开了一种基于过采样的GLRT检测方法，针对现有GLRT检测算法未能利用信道相关性和不能在短时间累积大量信号样本的问题而提出的，具体包括：获取过采样信号样本矩阵；计算接收信号的平均能量；计算相关统计量；计算检测统计量；将得到的检测统计量与预先设置的判决门限比较作检测判断。本发明通过对接收信号过采样得到接收信号的平均能量，在计算统计量时，利用MIMO信道的时间相关性，使得本发明的方法较现有频谱检测算法具有更好的检测性能，并且不受到噪声方差估计误差的影响，对噪声方差误差具有很强的鲁棒性。

Description

一种基于过采样的GLRT检测方法

技术领域

本发明属于通信技术领域，涉及认知无线电（Cognitive Radio）中的频谱检测算法。

背景技术

认知无线电技术是为了解决目前面临的频谱资源稀缺提出来的，而频谱检测算法是认知无线电的关键技术之一。为了避免认知无线电系统对授权用户产生有害的干扰，要求频谱检测方法能够在低信噪比下可靠地检测出授权用户信号来。

现有的MIMO(Multiple-Input Multiple-Output)频谱检测算法有匹配滤波器检测算法（Matched Filter Detection,MF）、能量检测算法（Energy Detection,ED）、循环谱检测算法（Cyclostationary Feature Detection）等。然而这些算法都有自己显著的优缺点。循环谱检测算法需要知道主用户信号的频率周期特性，这对于实际的认知无线电用户来说很难得到，并且该算法具有很高的复杂度。当认知无线电用户知道主用户的调制方式、载波频率等信号信息时，匹配滤波器检测算法被认为是最优的检测算法。但这些信息同样在实际系统中很难得到，并且该算法需要认知无线电用户与主用户准确同步，这对于认知无线电用户同样是很难的。能量检测算法不需要主用户信号的先验信息和认知无线电用户与主用户之间的信道信息，并且当主用户信号采样满足独立同分布特性时可以获得很好的检测性能，但是能量检测具有一个致命的缺点，其性能极易受到噪声功率估计误差的影响。

以上检测算法均或多或少需要一些先验信息，这对于实际的认知无线电用户都是很难做到的。为了克服上述缺陷，文献：T.J.Lim,R.Zhang,Y.C.Liang,and Y.Zeng,GLRT-basedspectrum sensing for cognitive radio,in Proc.IEEE Global Telecommun.Conf.,New Orleans,LA,Nov.30–Dec.4,2008,pp.1–5，提出了一种基于广义似然比（Generalized Log-likelihood RatioTest，GLRT）的检测算法—AGM(Arithmetic-to-Geometric Mean)算法；文献：Pu Wang,Jun Fang,Ning Han,and Hongbin Li,Multiantenna-Assisted Spectrum Sensing for Cognitive Radio,IEEETrans.Veh.Technol.,vol.59,no.4,pp.1791-1800,may2010，在AGM算法的基础上假设在感知时间内信道信息不变，得到另一种新的GLRT检测算法—MSE(Maximum-to-Sum ofEigenvalues)算法。以上这些GLRT检测算法不需要任何先验信息，但是都是假设接收信号样本为独立同分布的，所以这就导致接收机必须对信号按符号速率进行采样，这样就带了两个缺点：①无法在短时间内累积大量样本数，从而很难在较短的感知时间内获得很好的检测性能；②未能利用信道的时间相关特性。

发明内容

本发明的目的是为了解决现有GLRT检测算法未能利用信道相关性和不能在短时间累积大量信号样本的问题，提出了一种基于过采样的GLRT检测方法。

本发明技术方案是：一种基于过采样的GLRT检测方法，具体包括如下步骤：

步骤1：认知无线电用户分别对M根接收天线的接收信号进行过采样率为L的过采样，然后将第k个符号周期内过采样得到的样本信号组成过采样信号样本矩阵Y_k，1≤k≤K，其中，K表示过采样的符号周期数；

步骤2：认知无线电用户根据步骤1得到的过采样信号样本矩阵Y_k计算接收信号的平均能量

其中，tr表示矩阵迹运算；

步骤3：认知无线电用户根据步骤1得到的过采样信号样本矩阵Y_k计算矩阵B_k：B_k=Y_kU,1≤k≤K，其中，U通过对归一化信道时间相关矩阵

特征值分解得到，

Σ=diag[λ_h,1,…,λ_h,L],特征值λ_h,l,1≤l≤L按降序排列；

步骤4：根据步骤3的结果B_k,1≤k≤K，计算相关统计量T₂，

其中，B_k(:,l)表示矩阵B_k的第l列，‖B_k(:,l)‖表示对B_k(:,l)二范数运算；

步骤5：根据步骤2和步骤4的结果，计算检测统计量T：

$T=\frac{T_1}{T_2}=\frac{\frac{1}{\mathrm{KL}}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\mathrm{tr}\left(Y_k{Y}_k^H\right)}{{\left(\underset{l=1}{\overset{L}{\mathrm{\Π}}}\left(\frac{1}{K}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}{\left|\right|B_k(:,l)\left|\right|}^2\right)\right)}^{1/L}};$

步骤6：将得到的检测统计量T与预先设置的判决门限γ比较：

61：当检测统计量大于等于判决门限时，即T≥γ，认为当前频谱资源被主用户占用，认知无线电用户不能利用该频谱资源；

62：当检测统计量小于判决门限时，即T<γ，认为当前频谱资源空闲，认知无线电用户可以利用该频谱资源。

本发明的有益效果：本发明的方法通过对接收信号过采样得到接收信号的平均能量T₁；在计算统计量T₂时，利用MIMO信道的时间相关性，即利用归一化信道时间相关矩阵特征值分解得到U计算相关统计量T₂，使得本发明的方法较现有频谱检测算法具有更好的检测性能（尤其在低信噪比环境下），并且不受到噪声方差估计误差的影响，对噪声方差误差具有很强的鲁棒性。另外，本发明的方法在相同的系统要求检测性能下，需要的感知时间更短，可以实现快速、准确、高效地感知主用户信号，适用于实时性要求高的无线通信系统。

附图说明

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

图2为1ms感知时间下的本发明的方法与现有方法检测性能对比示意图。

图3为本发明的方法与现有方法时间消耗对比示意图。

具体实施方式

下面对本发明实施例中的技术方案进行清楚、完整地描述。可以理解的是，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本发明的方法的流程示意图如图1所示，包括：获取过采样信号样本矩阵；计算接收信号的平均能量；计算相关统计量；计算检测统计量；将得到的检测统计量与预先设置的判决门限比较作检测判断。

具体实施步骤如下：

步骤1：认知无线电用户分别对M根接收天线的接收信号进行过采样率为L的过采样，然后将第k个符号周期内过采样得到的样本信号组成为过采样信号样本矩阵Y_k,1≤k≤K，K表示过采样的符号周期数。这里，Y_k可以表示为 $\left[\begin{array}{c}{y}_{1,k}\\ ...\\ y_{M,k}\end{array}\right],$ y_m,k,m=1,…,M，表示第m根接收天线在第k个符号周期内过采样得到的样本信号组成的1×L的行向量。

步骤2：认知无线电用户根据步骤1得到的过采样信号样本矩阵Y_k,1≤k≤K计算接收信号的平均能量T₁：

$T_1=\frac{1}{\mathrm{KL}}\underset{k=1}{\overset{K}{\mathrm{\&Sigma;}}}\mathrm{tr}\left(Y_k{Y}_k^H\right)$

步骤3：认知无线电用户根据步骤1得到的过采样信号样本矩阵Y_k,1≤k≤K计算矩阵B_k：

B_k=Y_kU,1≤k≤K

U通过对归一化信道时间相关矩阵

特征值分解得到，可以提前计算好储存在认知无线电用户的储存器中直接调用。

Σ=diag[λ_h,1,…,λ_h,L],特征值λ_h,l,1≤l≤L按降序排列。这里，根据信道的时间相关模型获取，可以预先得到，属于本领域的现有技术，在这里不再详细描述。

步骤4：根据步骤3的结果B_k,1≤k≤K，计算相关统计量T₂：

$T_2={\left(\underset{l=1}{\overset{L}{\mathrm{\&Pi;}}}\left(\frac{1}{K}\underset{k=1}{\overset{K}{\mathrm{\&Sigma;}}}{\left|\right|B_k(:,l)\left|\right|}^2\right)\right)}^{1/L}$

步骤5：根据步骤2和步骤4的结果，计算检测统计量T：

$T=\frac{T_1}{T_2}=\frac{\frac{1}{\mathrm{KL}}\underset{k=1}{\overset{K}{\mathrm{\&Sigma;}}}\mathrm{tr}\left(Y_k{Y}_k^H\right)}{{\left(\underset{l=1}{\overset{L}{\mathrm{\&Pi;}}}\left(\frac{1}{K}\underset{k=1}{\overset{K}{\mathrm{\&Sigma;}}}{\left|\right|B_k(:,l)\left|\right|}^2\right)\right)}^{1/L}};$

步骤6：将步骤5的结果T与预先设置的判决门限γ比较，判决门限γ一般根据系统的具体虚警概率来设置。

61：当检测统计量大于等于判决门限时，即T≥γ，认为当前频谱资源被主用户占用，认知无线电用户不能利用该频谱资源。

62：当检测统计量小于判决门限时，即T<γ，认为当前频谱资源空闲，认知无线电用户可以利用该频谱资源。

根据本发明提供的基于MIMO过采样的GLRT检测方法，不仅可以在较短的时间内累积大量的信号样本，而且能充分利用信道的时间相关特性，获得最优的检测性能。

采用本发明实施例所述方法进行仿真测试，1ms感知时间下的本发明的方法与现有方法检测性能对比示意图如图2所示，仿真结果表明：在Jakes功率谱平坦瑞利衰落信道模型下，主用户采用BPSK调制，当信噪比为-20dB，过采样率L=16,主用户发射天线数P=1,认知无线电用户接收天线数M=4，频谱感知时间为1ms时，采用本发明的基于MIMO过采样的新型GLRT检测方法的检测概率能达到0.8，而现有的无噪声方差误差ED检测算法只有0.5，有1dB噪声方差误差的ED检测算法检测概率趋于零，而现有的非过采样GLRT检测算法（如AGM）检测概率只有0.38。

本发明的方法与现有方法时间消耗对比示意图如图3所示，在相同条件下达到系统要求的检测概率P_d=0.9，本发明的方法感知时间最短，尤其在低信噪比情况下。在-20dB情况下，由于AGM和MSE均需要较长的感知时间，MSE算法假设在感知时间内信道不变的条件不再满足。因此，MSE检测算法检测性能下降，在达到系统要求的检测概率P_d=0.9时，需要的感知时间超过AGM算法。

综上所述，本发明实施的基于MIMO过采样的新型GLRT检测方法，通过对接收信号过采样，利用MIMO信道的时间相关性，能在短频谱感知时间内，在低信噪比情况下获得较好的检测性能，可以实现快速、准确、高效地感知主用户信号，适用于实时性要求高的无线通信系统；并且，本发明在实施中不受到噪声方差估计误差的影响，对噪声方差误差具有很强的鲁棒性，所以具有稳定的检测性能。

本领域普通技术人员可以理解，实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件完成，所述的程序可以存储于计算机可读存储介质中，例如只读存储器（简称ROM）、随机存取存储器（简称RAM）、磁盘、光盘等。

以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。

Claims (1)

1.一种基于过采样的GLRT检测方法，所述GLRT具体为：广义似然比；其特征在于，具体包括如下步骤：

步骤1：认知无线电用户分别对M根接收天线的接收信号进行过采样率L的过采样，然后将第k个符号周期内过采样得到的样本信号组成过采样信号样本矩阵Y_k，1≤k≤K，其中，K表示过采样的符号周期数；

步骤2：认知无线电用户根据步骤1得到的过采样信号样本矩阵Y_k计算接收信号的平均能量

其中，tr表示矩阵迹运算；

步骤3：认知无线电用户根据步骤1得到的过采样信号样本矩阵Y_k计算矩阵B_k：B_k=Y_kU,1≤k≤K，其中，U通过对归一化信道时间相关矩阵

特征值分解得到，Σ=diag[λ_h,1,…,λ_h,L],特征值λ_h,l,1≤l≤L按降序排列；

步骤4：根据步骤3的矩阵B_k，计算相关统计量

其中，B_k(:,l)表示矩阵B_k的第l列，‖B_k(:,l)‖表示对B_k(:,l)二范数运算；

步骤5：根据步骤2和步骤4的结果，计算检测统计量T：

$T=\frac{T_1}{T_2}=\frac{\frac{1}{\mathrm{KL}}\underset{k=1}{\overset{K}{\mathrm{\&Sigma;}}}\mathrm{tr}\left(Y_k{Y}_k^H\right)}{{\left(\underset{l=1}{\overset{L}{\mathrm{\&Pi;}}}\left(\frac{1}{K}\underset{k=1}{\overset{K}{\mathrm{\&Sigma;}}}{\left|\right|B_k(:,l)\left|\right|}^2\right)\right)}^{1/L}};$

步骤6：将得到的检测统计量T与预先设置的判决门限γ比较：

61：当检测统计量大于等于判决门限时，即T≥γ，认为当前频谱资源被主用户占用，认知无线电用户不能利用该频谱资源；

62：当检测统计量小于判决门限时，即T<γ，认为当前频谱资源空闲，认知无线电用户可以利用该频谱资源。

Images