CN106230529A - 一种基于部分中继机制的协同频谱检测方法 - Google Patents

Abstract

本发明公开了一种基于部分中继机制的协同频谱检测方法，包括：选择中继策略并确定中继时间T_R；分别获取同时按能量检测方法进行检测的两个次级用户在时隙的前T_d‑T_R时间内接收信号的能量；获取低信噪比次级用户在时隙的T_d‑T_R到T_d时间内接收信号的能量；计算低信噪比次级用户检测到的总能量；分别将两个次级用户检测到的能量与各自的判决门限做比较，根据比较结果确定是否可以利用主用户的频段传输数据。本发明技术方案实现了高信噪比次级用户将分配到的检测时间部分用于检测，另一部分充当中继节点，协助低信噪比次级用户提高主用户检测精度，从而提高整个次级网络的性能。

Description

一种基于部分中继机制的协同频谱检测方法

技术领域

本发明属于通信技术中认知无线电领域，涉及认知无线电网络中多个次级用户协同检测频谱空洞的方法。

背景技术

频谱是无线通信的宝贵资源。在传统的无线通信系统中，频谱分配是固定的。一旦一段频谱授权给某个用户，这个用户(授权用户，主用户)便对其具有独占性，即使自身不使用也禁止其它用户(非授权用户，次级用户)占用。此方案实现简单，有利于保证主用户的服务质量，在无线通信发展伊始就得到了广泛应用，并一直沿用至今。然而，近来的研究表明，主用户在独享授权频段的同时，却没有充分利用该频段。根据美国联邦通信委员会的调查报告，授权频段的利用率在15％到85％之间波动；很多授权频段在不同的时间段内、不同的地理区域里都未被利用，形成了大量的频谱空洞。在频谱需求日益增大的今天，合理回收利用这些频谱空洞，可以极大地缓解当前频谱资源紧张的问题。

认知无线电是一种有效利用频谱空洞的新技术。这一技术允许次级用户在主用户未使用授权频段(频谱空洞存在)的前提下使用该频段；一旦主用户重新使用授权频段(频谱空洞消失)，次级用户必须马上停止使用，以免对主用户产生干扰。为了达到这一目标，次级用户首先必须检测频谱空洞是否存在。根据检测是否由多个次级用户协同完成，主用户检测可分为非协同检测和协同检测两大类。由于和非协同检测相比，协同检测能够避免单用户检测的不确定性、大幅提高检测精度，有效抑制隐藏终端和暴露终端的危害，因此逐渐成为了当前研究的热点。协同检测的本质是利用次级网络中的多个次级用户协同工作，一起完成主用户检测。在分布式次级网络中，所有次级用户地位均等，它们之间的协同可借助于中继机制实现。

中继协同提供了另一种获得多用户增益的途径，这种方法能够有效降低主用户检测所需的检测时间，提高检测精度、保护主用户不被干扰，因而得到了广泛关注。中继协同检测的核心思想是让空闲的次级用户为正在进行检测的次级用户充当中继节点。如图1所示，次级用户在空闲时可以扮演中继节点的角色，将自己接收到的主用户信号转发给其他次级用户。此时，对于其他正在进主用户检测的次级用户而言，它们不仅会接收到原有主用户信号，还会接收到中继次级用户转发过来的主用户信号；这些次级用户利用以上多个信号作出判决，就实现了与中继次级用户协同完成主用户检测。

目前，大部分的中继协同检测都假定次级用户基于时分多址(TDMA,TimeDivision Multiple Access)模式感知和接入授权频段，如图2所示。即每个时隙内，只有一个次级用户实施主用户检测，如果检测到频谱空洞，它就接入授权频段。但是当次级用户工作在频分多址(FDMA,Frequency Division Multiple Access)或者码分多址(CDMA,CodeDivision Multiple Access)模式下，如图3所示，即每个时隙内，次级用户同时进行检测，如果没有检测到主用户，它们就利用不同的频分或者码分子信道完成数据传输。在此情形下，可能并不存在真正意义上空闲的次级用户，传统的中继协同检测方法也不一定适用。

发明内容

本发明的目的在于克服现有技术的不足，提供一种基于部分中继机制的协同检测方法来检测频谱空洞，实现了高信噪比次级用户将分配到的检测时间部分用于检测，另一部分充当中继节点，协助低信噪比次级用户提高主用户检测精度，从而提高整个次级网络的性能。

本发明解决其技术问题所采用的技术方案是：

为了弥补传统的协同检测方法中次级用户必须工作在TDMA模式下的不足，以适应FDMA/CDMA型次级网络的需要，本发明提出了一种新的部分中继模型，如图3所示，在每个时隙T_d内，两个次级用户SU₁和SU₂均首先检测主用户，然后在T_c时间内进行数据传输；其中，高信噪比次级用户SU₂并不将T_d时间全部用来完成主用户检测，而只是用到了其中T_d-T_R的一部分；在剩下的T_R时间内，SU₂将充当中继节点，协助低信噪比次级用户SU₁提高检测精度。

一种基于部分中继机制的协同频谱检测方法，具体步骤包括：

步骤1，选择中继策略并确定中继时间T_R；

步骤2，分别获取同时按能量检测方法进行检测的两个次级用户在时隙的前T_d-T_R时间内接收信号的能量；所述T_d表示一次检测的时隙宽度；

步骤3，获取低信噪比次级用户在时隙的T_d-T_R到T_d时间内接收信号的能量；

步骤4，计算低信噪比次级用户检测到的总能量；

步骤5，分别将两个次级用户检测到的能量与各自的判决门限做比较，根据比较结果确定是否可以利用主用户的频段传输数据。

进一步的，步骤1包括采用平衡检测精度准则或采用最大化检测敏捷性准则确定中继策略。

进一步的，采用平衡检测精度准则时，所述中继策略η由如下公式确定：

$\mathrm{\η}\left(\mathrm{\ρ}\right)=\frac{5\mathrm{\ρ}-2+\sqrt{53{\mathrm{\ρ}}^2-20\mathrm{\ρ}-24}}{14\mathrm{\ρ}}$

其中，γ₁和γ₂为低信噪比次级用户和高信噪比次级用户的接收信噪比。

进一步的，采用最大化检测敏捷性准则时，所述中继策略η由如下公式确定：

$min\left\{E\right\{T_{\mathrm{\η}}\}=(\frac{1}{1-P_{m1}}+\frac{1}{1-P_{m2}}-\frac{1}{1-P_{m1}{P}_{m2}})(T_d+T_c)\}$

其中，P_m1和P_m2分别表示低信噪比次级用户和高信噪比次级用户的漏检概率，(T_d+T_c)表示时隙宽度。

进一步的，所述漏检概率P_m1用如下公式表示：

$P_{m1}=1-Q\left(\frac{Q^{-1}\left(\mathrm{\α}\right)\sqrt{2(1+6\mathrm{\η})}-\sqrt{n}({\mathrm{\γ}}_1+{\mathrm{\η\γ}}_2)}{\sqrt{2{(1+{\mathrm{\γ}}_1)}^2+2\mathrm{\η}\[{(2+{\mathrm{\γ}}_1+{\mathrm{\γ}}_2)}^2-{(1+{\mathrm{\γ}}_1)}^2+3{(1+{\mathrm{\γ}}_2)}^2\]}}\right)$

所述漏检概率P_m2用如下公式表示：

$P_{m2}=1-Q\left(\frac{Q^{-1}\left(\mathrm{\α}\right)\sqrt{2}-\sqrt{n(1-\mathrm{\η})}{\mathrm{\γ}}_2}{\sqrt{2{(1+{\mathrm{\γ}}_2)}^2}}\right)$

其中，α表示低信噪比次级用户和高信噪比次级用户的虚警概率，γ₁和γ₂为低信噪比次级用户和高信噪比次级用户的接收信噪比，n＝T_d f_s，f_s为接收信号采样周期。

进一步的，所述中继时间T_R用如下公式表示：

T_R＝η·T_d。

进一步的，时隙前T_d-T_R时间内，低信噪比次级用户接收信号的能量服从如下高斯分布：

${v_1}^{NR}~\left\{\begin{array}{cc}N(n-n_R,2(n-n_R\left)\right)& H_0\\ N\left(\right(n-n_R\left)\right(1+{\mathrm{\γ}}_1),2(n-n_R\left){\left(1+{\mathrm{\γ}}_1\right)}^2\right)& H_1\end{array}\right.$

高信噪比次级用户接收信号的能量服从如下高斯分布：

$v_2~\left\{\begin{array}{cc}N(n-n_R,2(n-n_R\left)\right)& H_0\\ N\left(\right(n-n_R\left)\right(1+{\mathrm{\γ}}_2),2(n-n_R\left){\left(1+{\mathrm{\γ}}_2\right)}^2\right)& H_1\end{array}\right.$

其中，n_R＝T_R f_s。

进一步的，时隙的T_d-T_R到T_d时间内，低信噪比次级用户接收信号的能量服从如下高斯分布：

${v_1}^R~\left\{\begin{array}{cc}N(2n_R,14n_R)& H_0\\ N\left(n_R\right(2+{\mathrm{\γ}}_1+{\mathrm{\γ}}_2),2n_R\[{(2+{\mathrm{\γ}}_1+{\mathrm{\γ}}_2)}^2+3{(1+{\mathrm{\γ}}_2)}^2\])& H_1\end{array}\right..$

进一步的，低信噪比次级用户检测到的总能量服从如下高斯分布：

$v_1~\left\{\begin{array}{cc}N(n+n_R,2(n+6n_R\left)\right)& H_0\\ N\left(n\right(1+{\mathrm{\γ}}_1)+n_R(1+{\mathrm{\γ}}_2),2n{\left(1+{\mathrm{\γ}}_1\right)}^2+2n_R\[{(2+{\mathrm{\γ}}_1+{\mathrm{\γ}}_2)}^2-{(1+{\mathrm{\γ}}_1)}^2+3{(1+{\mathrm{\γ}}_2)}^2\])& H_1\end{array}\right..$

进一步的，步骤5具体包括：

将低信噪比次级用户检测到的能量v₁与低信噪比次级用户的判决门限λ₁做比较，如果v₁大于λ₁，表示低信噪比次级用户不能利用该频段传输数据，否则可以利用该频段传输数据；

将高信噪比次级用户检测到的能量v₂与高信噪比次级用户的判决门限λ₂做比较，如果v₂大于λ₂，表示高信噪比次级用户不能利用该频段传输数据，否则可以利用该频段传输数据；

所述判决门限λ₁和λ₂由如下公式确定：

$\begin{array}{c}{\mathrm{\λ}}_1=Q^{-1}\left(\mathrm{\α}\right)\sqrt{2n(1+6\mathrm{\η})}+n(1+\mathrm{\η})\\ {\mathrm{\λ}}_2=Q^{-1}\left(\mathrm{\α}\right)\sqrt{2n(1-\mathrm{\η})}+n(1-\mathrm{\η})\end{array}.$

本发明具有如下有益效果：

(1)条件简单。本发明并不要求中继次级用户完全处于空闲状态，因此即使当网络中所有次级用户同时进入到授权频段时，也能正常工作。

(2)有效提高检测精度。接收信噪比较高的次级用户并不将分配到的检测时间全部用于检测，而是抽出部分时间充当中继节点，尽管这样做会给自己的检测性能带来损害，但它能够协助其他次级用户提高主用户检测的检测精度，从而提高整个次级网络的性能。

(3)大幅改善检测敏捷性。检测敏捷性的目的是为了描述次级用户迅速检测出主用户状态变化的能力。如果所需时间较短，则说明检测敏捷性较高。本发明中，如果SU₂按照搜索到的最优中继策略转发自己的接收信号，就能保证SU₁和SU₂在最短的时间内检测出主用户的存在，从而最大化整个次级系统的检测敏捷性。

以下结合附图及实施例对本发明作进一步详细说明，但本发明的一种基于部分中继机制的协同频谱检测方法不局限于实施例。

附图说明

图1为本发明方法的两个次级用户组成的分布式次级网络示意图；

图2为TDMA工作模式下次级用户中继示意图；

图3本发明方法的FDMA和CDMA模式下次级用户中继示意图；

图4为本发明方法的主流程图。

具体实施方式

如图4所示，本发明提供的一种基于部分中继机制的协同频谱检测方法，具体包括以下步骤：

步骤1，选择中继策略并确定中继时间T_R。

具体的，包括如下步骤：

(1)SU₁和SU₂接收信号计算

假定主用户信号为s(k)，SU₁和SU₂与主用户之间的信道分别为h₁(k)和h₂(k),对应的信道噪声分别为n₁(k)和n₂(k)。H₀表示主用户发射机信号不存在，授权频段为频谱空洞,H₁表示主用户发射机信号存在，授权频段不是频谱空洞。在每个时隙的前T_d-T_R时间内，没有中继发生，SU₁和SU₂同时进行检测,它们的接收信号可分别表示为：

${y_1}^{NR}\left(k\right)=\left\{\begin{array}{cc}{n}_1\left(k\right)& H_0\\ h_1\left(k\right)\·s\left(k\right)+n_1\left(k\right)& H_1\end{array}\right.---\left(1\right)$

$y_2\left(k\right)=\left\{\begin{array}{cc}{n}_2\left(k\right)& H_0\\ h_2\left(k\right)\·s\left(k\right)+n_2\left(k\right)& H_1\end{array}\right.---\left(2\right)$

在剩余的T_R时间内，SU₁继续实施检测，SU₂停止检测，并将自己的接收信号中继给SU₁。假定中继基于放大转发的方式进行，此时SU₁的接收信号可以表示如下

$\begin{array}{c}{y_1}^R\left(k\right)={y_1}^{NR}\left(k\right)+\sqrt{k}\·h_{12}\left(k\right)\·y_2\left(k\right)\\ =\left\{\begin{array}{cc}{n}_1\left(k\right)+\sqrt{k}\·h_{12}\left(k\right)\·n_2\left(k\right)& H_0\\ \begin{array}{c}\[h_1\left(k\right)+\sqrt{k}\·h_{12}\left(k\right)\·h_2\left(k\right)\]\·s\left(k\right)\\ +\[n_1\left(k\right)+\sqrt{k}\·h_{12}\left(k\right)\·n_2\left(k\right)\]\end{array}& H_1\end{array}\right.\end{array},---\left(3\right)$

其中，k为中继增益因子，h ₁₂(k)表示SU₂与SU₁之间的信道影响。

(2)中继协同能量检测

由于能量检测具有计算复杂度低、对主用户信号的相关信息要求低等特点，因此比较常用。假设S₁(k),S₂(k)以及h₁₂(k)是零均值实高斯独立随机过程，n₁(k)和n₂(k)为零均值单位方差的实高斯白噪声。在此情形下，SU₁和SU₂接收信号的能量可以分别用如下公式计算：

$\begin{array}{c}{v}_1={v_1}^{NR}+{v_1}^R\\ =\underset{k=1}{\overset{n-n_R}{\Σ}}|{y_1}^{NR}\left(k\right){|}^2+\underset{k=n-n_R+1}{\overset{n}{\Σ}}|{y_1}^R\left(k\right){|}^2\end{array},---\left(4\right)$

$v_2=\underset{k=1}{\overset{n-n_R}{\Σ}}|y_2\left(k\right){|}^2,$

其中，v₁ ^NR和v₁ ^R分别表示SU₂不中继和中继时SU₁接收到的能量，n＝T_d f_s，n_R＝T_Rf_s，f_s为接收信号采样周期。v₁ ^NR和v₂分别近似服从如下高斯分布：

$\begin{array}{c}{v_1}^{NR}=\underset{k=1}{\overset{n-n_R}{\Σ}}|{y_1}^{NR}\left(k\right){|}^2\\ ~\left\{\begin{array}{cc}N(n-n_R,2(n-n_R\left)\right)& H_0\\ N\left(\right(n-n_R\left)\right(1+{\mathrm{\γ}}_1),2(n-n_R\left){\left(1+{\mathrm{\γ}}_1\right)}^2\right)& H_1\end{array}\right.\end{array},---\left(5\right)$

$v_2~\left\{\begin{array}{cc}N(n-n_R,2(n-n_R\left)\right)& H_0\\ N\left(\right(n-n_R\left)\right(1+{\mathrm{\γ}}_2),2(n-n_R\left){\left(1+{\mathrm{\γ}}_2\right)}^2\right)& H_1\end{array}\right.,---\left(6\right)$

其中，γ₁＝E{|h₁(k)·s(k)|²}/E{|n₁(k)|²}和γ₂＝E{|h₂(k)·s(k)²}/E{|n₂(k)|²}分别代表不发生中继时SU₁和SU₂的接收信噪比。

考虑到中继时间T_R一般不会很小(过小的T_R几乎等同于没有发生中继)，因此n_R通常较大。根据中心极限定理，v₁ ^R同样近似服从高斯分布：

${v_1}^R~\left\{\begin{array}{cc}N(2n_R,14n_R)& H_0\\ \begin{array}{c}N\left(n_R\right(2+{\mathrm{\γ}}_1+{\mathrm{\γ}}_2),\\ 2n_R\[{(2+{\mathrm{\γ}}_1+{\mathrm{\γ}}_2)}^2+3{(1+{\mathrm{\γ}}_2)}^2\])\end{array}& H_1\end{array}\right.,---\left(7\right)$

不考虑SU₂的发射功率限制，令中继增益正好补偿SU₂与SU₁之间信道的功率衰减，可知v₁＝v₁ ^NR+v₁ ^R服从如下高斯分布：

$v_1~\left\{\begin{array}{cc}N(n+n_R,2(n+6n_R\left)\right)& H_0\\ \begin{array}{c}N\left(n\right(1+{\mathrm{\γ}}_1)+n_R(1+{\mathrm{\γ}}_2),2n{(1+{\mathrm{\γ}}_1)}^2+\\ 2n_R\[{(2+{\mathrm{\γ}}_1+{\mathrm{\γ}}_2)}^2-{(1+{\mathrm{\γ}}_1)}^2+3{(1+{\mathrm{\γ}}_2)}^2\])\end{array}& H_1\end{array}\right.,---\left(8\right)$

推导出了v₁和v₂的概率分布之后，SU₁和SU₂的虚警概率和漏检概率用如下公式表示：

$P_{f1}=Q\left(\frac{{\mathrm{\λ}}_1-(n+n_R)}{\sqrt{2(n+6n_R)}}\right),---\left(9\right)$

$P_{f2}=Q\left(\frac{{\mathrm{\λ}}_2-(n-n_R)}{\sqrt{2(n-n_R)}}\right),---\left(10\right)$

$P_{m1}=1-Q\left(\frac{\mathrm{\λ}\[-\[n(1+{\mathrm{\γ}}_1)+n_R(1+{\mathrm{\γ}}_2)\]}{\sqrt{2n{(1+{\mathrm{\γ}}_1)}^2+2n_R\[{(2+{\mathrm{\γ}}_1+{\mathrm{\γ}}_2)}^2-{(1+{\mathrm{\γ}}_1)}^2+3{(1+{\mathrm{\γ}}_2)}^2\]}}\right),---\left(11\right)$

$P_{m2}=1-Q\left(\frac{{\mathrm{\λ}}_2-(n-n_R)(1+{\mathrm{\γ}}_2)}{\sqrt{2(n-n_R){(1+{\mathrm{\γ}}_2)}^2}}\right),---\left(12\right)$

其中，λ₁和λ₂分别为SU₁和SU₂的判决门限，P_f1表示SU₁的虚警概率，P_f2表示SU₂的虚警概率，P_m1表示SU₁的漏检概率，P_m2表示SU₂的漏检概率。

(3)中继协同能量检测中继策略的选择

在上述中继协同能量检测中，SU₂从T_d时间中抽出T_R时间用以中继转发接收信号，帮助SU₁提高接收信噪比、改善检测性能；而自己真正用于检测的时间只有T_d-T_R。很明显，如果采取保守的中继策略，选择较小T_R，给予SU₁的帮助就会比较有限；反之，如果采取激进的中继策略，选择较大T_R，SU₂自身的检测性能就会受到较大影响。因此，有必要选择合理的中继策略，确定最佳的T_R。

为了量化中继策略的保守和激进程度，定义：

$\mathrm{\η}\stackrel{\mathrm{\Δ}}{=}\frac{T_R}{T_d}=\frac{n_R}{n},---\left(13\right)$

其中，η(0<η<1)越接近0，则说明中继策略越保守；η越接近1，则说明中继策略越激进。为了选择最佳中继策略，可以从平衡检测精度和最大化检测敏捷性两个角度确定。

①平衡检测精度

SU₁和SU₂组成的次级网络中，SU₂的接收信噪比相对较高，因而其检测精度也相对较高。平衡检测精度，实际上就是要借助部分中继协同使得SU₁和SU₂实现相同的虚警概率和漏检概率。

采用恒虚警准则，令SU₁和SU₂的虚警概率均为α，结合公式(9)和(10)给出的虚警概率表达式，SU₁和SU₂的判决门限可以分别表示为：

${\mathrm{\&lambda;}}_1=Q^{-1}\left(\mathrm{\&alpha;}\right)\sqrt{2n(1+6\mathrm{\&eta;})}+n(1+\mathrm{\&eta;}),---\left(14\right)$

${\mathrm{\&lambda;}}_2=Q^{-1}\left(\mathrm{\&alpha;}\right)\sqrt{2n(1-\mathrm{\&eta;})}+n(1-\mathrm{\&eta;}),---\left(15\right)$

将上式中λ₁和λ₂分别代入公式(11)和(12)中，可得SU₁和SU₂的漏检概率：

$P_{m1}=1-Q\left(\frac{Q^{-1}\left(\mathrm{\&alpha;}\right)\sqrt{2(1+6\mathrm{\&eta;})}-\sqrt{n}({\mathrm{\&gamma;}}_1+{\mathrm{\&eta;\&gamma;}}_2)}{\sqrt{2{(1+{\mathrm{\&gamma;}}_1)}^2+2\mathrm{\&eta;}\&lsqb;{(2+{\mathrm{\&gamma;}}_1+{\mathrm{\&gamma;}}_2)}^2-{(1+{\mathrm{\&gamma;}}_1)}^2+3{(1+{\mathrm{\&gamma;}}_2)}^2\&rsqb;}}\right),---\left(16\right)$

$P_{m2}=1-Q\left(\frac{Q^{-1}\left(\mathrm{\&alpha;}\right)\sqrt{2}-\sqrt{n(1-\mathrm{\&eta;})}{\mathrm{\&gamma;}}_2}{\sqrt{2{(1+{\mathrm{\&gamma;}}_2)}^2}}\right),---\left(17\right)$

再令P_m1＝P_m2，即可求解出所需的η。考虑到认知无线电场景中次级用户的接收信噪比通常都比较小，将SU₂与SU₁接收信噪比之比定义为为了平衡SU₂和SU₁的检测精度，应当采用的中继策略η只与ρ关，具体可以写成

$\mathrm{\&eta;}\left(\mathrm{\&rho;}\right)=\frac{5\mathrm{\&rho;}-2+\sqrt{53{\mathrm{\&rho;}}^2-20\mathrm{\&rho;}-24}}{14\mathrm{\&rho;}},---\left(18\right)$

②最大化检测敏捷性

检测敏捷性是指次级用户迅速检测出主用户状态变化的能力。主用户出现之后次级用户检测出主用户存在所需的时间较短，说明检测敏捷性较高；反之，则说明检测敏捷性较低。最大化检测敏捷性，实际上就是要最小化该检测时间。该检测时间有别于前面的单次检测所需的时间T_d，具体可以计算如下

T_η＝k_η(T_d+T_c), (19)

其中，k_η表示次级用户检测出主用户存在所需的检测次数，T_d+T_c表示时隙宽度(每个时隙内进行一次检测)。

在中继协同能量检测中，只有全部次级用户均检测到主用户存在，才能算是真正检测到了主用户的存在。然而，主用户出现后，SU₁和SU₂经过多少次检测才能全部检测出主用户存在并不确定，因此k_η是一个随机变量。当k_η＝l时，对应的概率为

P_r(k_η＝l)＝P_r(SU₁在第l次，SU₂在前l次检测出主用户存在)+P_r(SU₁在前l次，SU₂在第l次检测出主用户存在)-P_r(SU₁，SU₂均在第l次检测出主用户存在)

基于k_η＝l时的概率分布，k_η的均值可表示为

$\begin{array}{c}E\left\{k_{\mathrm{\&eta;}}\right\}=\underset{l=1}{\overset{\mathrm{\&infin;}}{\&Sigma;}}l\&CenterDot;P_r(k_{\mathrm{\&eta;}}=l)\\ =\frac{1}{1-P_{m1}}+\frac{1}{1-P_{m2}}-\frac{1}{1-P_{m1}{P}_{m2}}\end{array},---\left(21\right)$

相应地，检测出主用户存在所需的平均时间为

$E\left\{T_{\mathrm{\&eta;}}\right\}=(\frac{1}{1-P_{m1}}+\frac{1}{1-P_{m2}}-\frac{1}{1-P_{m1}{P}_{m2}})(T_d+T_c),---\left(22\right)$

给定α、γ₁、γ₂和n时，将公式(16)和(17)代入公式(22)中，并基于该公式展开搜索，可以很容易地找出令检测平均时间E{T_η}达到最小的η。在执行部分中继协同能量检测方案时，如果SU₂按照搜索到的最优中继策略转发自己的接收信号，就能保证SU₁和SU₂在最短的时间内检测出主用户的存在，从而最大化整个次级系统的检测敏捷性。

求得η后，根据公式(13)可求出合适的中继时间T_R和n_R。

步骤2，分别获取同时按能量检测方法进行检测的两个次级用户在时隙的前T_d-T_R时间内接收信号的能量；所述T_d表示一次检测的时隙宽度。

具体的，根据公式(5)获得SU₁接收到的能量v₁ ^NR,根据公式(6)获得SU₂接收到的能量v₂。

步骤3，获取低信噪比次级用户在时隙的T_d-T_R到T_d时间内接收信号的能量。

具体的，根据公式(7)，获得SU₁接收到的能量v₁ ^R。

步骤4，计算低信噪比次级用户SU₁检测到的总能量。

具体的，根据公式(8)，获得SU₁接收到的总能量v₁。

步骤5，分别将两个次级用户检测到的能量与各自的判决门限做比较，根据比较结果确定是否可以利用主用户的频段传输数据。

具体的，将低信噪比次级用户检测到的能量v₁与低信噪比次级用户的判决门限λ₁做比较，如果v₁大于λ₁，表示低信噪比次级用户不能利用该频段传输数据，否则可以利用该频段传输数据。

将高信噪比次级用户检测到的能量v₂与高信噪比次级用户的判决门限λ₂做比较，如果v₂大于λ₂，表示高信噪比次级用户不能利用该频段传输数据，否则可以利用该频段传输数据。

所述判决门限λ₁和λ₂由公式(14)和公式(15)求出。

步骤6，返回步骤2继续进行下一次检测。

以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种基于部分中继机制的协同频谱检测方法，其特征在于：两个次级用户同时对主用户进行检测，高信噪比次级用户将分配到的检测时间部分用于中继，转发接收信号给低信噪比次级用户，具体步骤包括：

步骤1，选择中继策略并确定中继时间T_R；

步骤2，分别获取同时按能量检测方法进行检测的两个次级用户在时隙的前T_d-T_R时间内接收信号的能量；所述T_d表示一次检测的时隙宽度；

步骤3，获取低信噪比次级用户在时隙的T_d-T_R到T_d时间内接收信号的能量；

步骤4，计算低信噪比次级用户检测到的总能量；

步骤5，分别将两个次级用户检测到的能量与各自的判决门限做比较，根据比较结果确定是否可以利用主用户的频段传输数据。

2.根据权利要求1所述的基于部分中继机制的协同频谱检测方法，其特征在于：步骤1包括采用平衡检测精度准则或采用最大化检测敏捷性准则确定中继策略。

3.根据权利要求2所述的基于部分中继机制的协同频谱检测方法，其特征在于：采用平衡检测精度准则时，所述中继策略η由如下公式确定：

$\mathrm{\&eta;}\left(\mathrm{\&rho;}\right)=\frac{5\mathrm{\&rho;}-2+\sqrt{53{\mathrm{\&rho;}}^2-20\mathrm{\&rho;}-24}}{14\mathrm{\&rho;}}$

其中，γ₁和γ₂为低信噪比次级用户和高信噪比次级用户的接收信噪比。

4.根据权利要求2所述的基于部分中继机制的协同频谱检测方法，其特征在于：采用最大化检测敏捷性准则时，所述中继策略η由如下公式确定：

$min\left\{E\right\{T_{\mathrm{\&eta;}}\}=(\frac{1}{1-P_{m1}}+\frac{1}{1-P_{m2}}-\frac{1}{1-P_{m1}{P}_{m2}})(T_d+T_c)\}$

其中，P_m1和P_m2分别表示低信噪比次级用户和高信噪比次级用户的漏检概率，(T_d+T_c)表示时隙宽度。

5.根据权利要求4所述的基于部分中继机制的协同频谱检测方法，其特征在于：所述漏检概率P_m1用如下公式表示：

$P_{m1}=1-Q\left(\frac{Q^{-1}\left(\mathrm{\&alpha;}\right)\sqrt{2(1+6\mathrm{\&eta;})}-\sqrt{n}({\mathrm{\&gamma;}}_1+{\mathrm{\&eta;\&gamma;}}_2)}{\sqrt{2{(1+{\mathrm{\&gamma;}}_1)}^2+2\mathrm{\&eta;}\&lsqb;{(2+{\mathrm{\&gamma;}}_1+{\mathrm{\&gamma;}}_2)}^2-{(1+{\mathrm{\&gamma;}}_1)}^2+3{(1+{\mathrm{\&gamma;}}_2)}^2\&rsqb;}}\right)$

所述漏检概率P_m2用如下公式表示：

$P_{m2}=1-Q\left(\frac{Q^{-1}\left(\mathrm{\&alpha;}\right)\sqrt{2}-\sqrt{n(1-\mathrm{\&eta;})}{\mathrm{\&gamma;}}_2}{\sqrt{2{(1+{\mathrm{\&gamma;}}_2)}^2}}\right)$

其中，α表示低信噪比次级用户和高信噪比次级用户的虚警概率，γ₁和γ₂为低信噪比次级用户和高信噪比次级用户的接收信噪比，n＝T_d f_s，f_s为接收信号采样周期。

6.根据权利要求5所述的基于部分中继机制的协同频谱检测方法，其特征在于：所述中继时间T_R用如下公式表示：

T_R＝η·T_d。

7.根据权利要求6所述的基于部分中继机制的协同频谱检测方法，其特征在于：时隙前T_d-T_R时间内，低信噪比次级用户接收信号的能量服从如下高斯分布：

${v_1}^{NR}~\left\{\begin{array}{cc}N\left(n-n_R,2\left(n-n_R\right)\right)& H_0\\ N\left(\left(n-n_R\right)\left(1+{\mathrm{\&gamma;}}_1\right),2\left(n-n_R\right){\left(1+{\mathrm{\&gamma;}}_1\right)}^2\right)& H_1\end{array}\right.$

高信噪比次级用户接收信号的能量服从如下高斯分布：

$v_2~\left\{\begin{array}{cc}N(n-n_R,2(n-n_R\left)\right)& H_0\\ N\left(\right(n-n_R\left)\right(1+{\mathrm{\&gamma;}}_2),2(n-n_R\left){\left(1+{\mathrm{\&gamma;}}_2\right)}^2\right)& H_1\end{array}\right.$

其中，n_R＝T_R f_s。

8.根据权利要求7所述的基于部分中继机制的协同频谱检测方法，其特征在于：时隙的T_d-T_R到T_d时间内，低信噪比次级用户接收信号的能量服从如下高斯分布：

${v_1}^R~\left\{\begin{array}{cc}N(2n_R,14n_R)& H_0\\ N\left(n_R\right(2+{\mathrm{\&gamma;}}_1+{\mathrm{\&gamma;}}_2),2n_R\&lsqb;{(2+{\mathrm{\&gamma;}}_1+{\mathrm{\&gamma;}}_2)}^2+3{(1+{\mathrm{\&gamma;}}_2)}^2\&rsqb;)& H_1\end{array}\right..$

9.根据权利要求8所述的基于部分中继机制的协同频谱检测方法，其特征在于：低信噪比次级用户检测到的总能量服从如下高斯分布：

$v_1~\left\{\begin{array}{cc}N\left(n+n_R,2\left(n+6n_R\right)\right)& H_0\\ N\left(n\left(1+{\mathrm{\&gamma;}}_1\right)+n_R\left(1+{\mathrm{\&gamma;}}_2\right),2n{\left(1+{\mathrm{\&gamma;}}_1\right)}^2+2n_R\&lsqb;{\left(2+{\mathrm{\&gamma;}}_1+{\mathrm{\&gamma;}}_2\right)}^2-{\left(1+{\mathrm{\&gamma;}}_1\right)}^2+3{\left(1+{\mathrm{\&gamma;}}_2\right)}^2\&rsqb;\right)& H_1\end{array}\right..$

10.根据权利要求9所述的基于部分中继机制的协同频谱检测方法，其特征在于：步骤5包括：

将低信噪比次级用户检测到的能量v₁与低信噪比次级用户的判决门限λ₁做比较，如果v₁大于λ₁，表示低信噪比次级用户不能利用该频段传输数据，否则可以利用该频段传输数据；

将高信噪比次级用户检测到的能量v₂与高信噪比次级用户的判决门限λ₂做比较，如果v₂大于λ₂，表示高信噪比次级用户不能利用该频段传输数据，否则可以利用该频段传输数据；

所述判决门限λ₁和λ₂由如下公式确定：

$\begin{array}{c}{\mathrm{\&lambda;}}_1=Q^{-1}\left(\mathrm{\&alpha;}\right)\sqrt{2n(1+6\mathrm{\&eta;})}+n(1+\mathrm{\&eta;})\\ {\mathrm{\&lambda;}}_2=Q^{-1}\left(\mathrm{\&alpha;}\right)\sqrt{2n(1-\mathrm{\&eta;})}+n(1-\mathrm{\&eta;})\end{array}.$