CN105429717A - 一种卫星通信感知频谱技术中的能量检测方法 - Google Patents

Abstract

本发明涉及一种卫星通信感知频谱技术中的能量检测方法，C谱未授权的通信卫星首先根据接收天线的极化状态，获得星地通信使用的C波段下，接收信号对干扰加噪声功率比(SINR)，然后根据频谱能量检测法，次用户得到接收信号能量的决策统计数值分布，最后结合接收信号对干扰加噪声功率比的最大值，输入到接收信号的决策统计数值分布参数中，在一定时间T内和带宽B，通过决策统计数值大于事先确定的某个阈值概率，确定频谱授权用户信号是否存在。利用这种技术，次用户不仅可以提高接收到的信号对干扰和噪声的信噪比，与现有技术比，还可以提高了信号检测概率，进而改善了频谱感知的效率。

Description

一种卫星通信感知频谱技术中的能量检测方法

技术领域

本发明涉及一种卫星通信感知频谱技术中的能量检测方法，属于卫星通信技术领域。

背景技术

未来通信卫星网络和地面网络聚合在一个标准的平台上，共同提供移动性、固定性和广播服务的可能性已经成为了一个重要的研究热点，这也要求通信卫星和地面网络使用相同的频段。高速宽带和多媒体无线服务的迅速发展，导致宽带无线频谱需求增加，同时由于当前频段分割和静态分配的策略，可用的频段资源变得越来越稀少。

在这种背景下，认知无线电通信技术是一种在共存的异构无线网络里高效地使用频谱效率的技术。使用相同频段的无线网络有多种表现形式，如两种地面网络组合，两种卫星网络组合，卫星和地面的异构网络组合等。这种技术允许主用户和次用户系统共同使用相同的频段，从而不用阻碍授权给主用户系统频段的正常工作。

在卫星通信领域，固定通信卫星服务通常使用C波段和K波段；移动通信卫星服务更多使用L和S波段，原因在于可以用板上天线，更好的穿透力和更小的大气影响。地面服务(如3G、LTE、WiMax和WiFi服务)的发展对通信卫星的L和C波段造成了持续的压力。宽带服务的高需求和L、S波段频率资源的限制使用，使得固定卫星服务已经开始使用更多频率的Ku和Ka波段。因此，探索高效的频率共享技术提高频谱效率，同时保证服务质量是重要的具有挑战的研究问题。这种技术也催生了认知卫星通信的概念，用于在两种通信卫星网络系统或者卫星和地面网络系统中开发频谱共享技术。

频谱感知是认知无线通信中的重要手段，通过多种维度(时间、空间、频率和码)用来获得频段使用特征，决定哪种类型的信号占用频段，以及包括不同的特征，如调制方式、波形、带宽、载波频率等。现有的频谱感知研究主要考虑信号的幅度特征、频率特征或者相位特征，可以认为是标量信号处理过程。在可用的频谱感知技术中，能量检测技术是最常用的技术，原因在于计算和实现复杂度低。空间传播的电磁信号可以被某种极化状态所极化，一种信号的极化状态，也是一种电磁信号向量特征，是电场向量在振荡期间在空间某一点的方向。由于传统的能量检测技术只考虑到信号的幅度部分信息，排斥其正交部分的信息，电磁信号向量正交特征的这部分信息没有充分得到开发。研究信号极化状态特征方面，也就是研究信号极化域，可看作探索频谱共享的一个额外的自由度和认知卫星通信的一种传输机制。通过正交极化，2个独立信号在相同的频带上可以进行传输，使得频段使用效率翻倍，2个正交极化信号传输可认为是2个不同的数据流。从认知的角度讲，如果主用户系统被设计为以一种极化状态传输，那么次用户可被设计为以另外一种极化状态传输。此外，使用极化域可以减少干扰。

异构卫星和地面网络系统，如图1所示，系统由卫星、卫星终端和地面基站组成。在场景中，地面网络使用WiMax技术，其使用C波段并作为主系统；卫星网络使用C波段，并作为次系统。在地面WiMax系统中，把基站之间的信号回传链路作为主链路，卫星终端到卫星之间的上行链路作为次链路。这种实际场景有助于卫星操作者利用地面授权的C波段作为次要的选择。在卫星传输和地面传输中，发射和接收均使用双极化天线。由于相距距离远，地面WiMax基站对卫星的干扰可忽略，但是卫星终端对地面接收机基站的干扰是强的，需要考虑并满足地面主用户的干扰约束。

如何实现卫星与地面网络的频段共用，且克服卫星终端对地面接收机基站的干扰是本领域亟待解决的技术问题。

发明内容

本发明的目的在于克服现有技术的不足，提供一种基于最优极化接收机制的通信卫星感知频谱的能量检测技术，未授权的通信卫星利用这种技术，不仅可以提高接收到的信号对干扰和噪声的信噪比，还可以利用最优信噪比，通过频谱能量检测法，确定授权用户信号是否存在。

本发明目的通过如下技术方案予以实现：

提供一种基于最优化极化接收机制的频谱感知能量检测方法，包括下列步骤：

(1)地面基站为某波段的授权主用户，卫星为所述波段的非授权用户，为次用户，卫星、卫星用户终端、地面基站的发送和接收均采用极化天线；次用户根据接收天线的极化状态，计算次用户的信噪比SINR，其中，S为次用户接收功率，P为主用户接收功率，N₀为噪声平均功率，得出SINR的最大值maxSINR；

(2)推导得出主用户信号存在时，主用户发射信号的平均检测概率关于信噪比SINR的公式；

(3)使用步骤(1)计算的maxSINR代替SINR，通过主用户发射信号的平均检测概率与第一阈值比较，确定主用户信号是否存在。

优选的，得出SINR的最大值maxSINR的具体方法为：次用户调整接收极化天线与接收的信号角度α_s，使得θ_p＝θ_s+α_s，其中θ_s为分别代表极化状态s_r和s_s的在极化球上的两个点的圆弧的对向角，θ_p为分别代表极化状态s_r和s_p的在极化球上的两个点的圆弧的对向角，s_r为次用户接收机的极化状态，s_s表示次用户发射信号的极化状态，s_p表示主用户极化状态；

在所述波段干扰对噪声比INR确定情况下，利用如下公式得到信噪比SINR的最大值maxSINR；

$\max SINR=SNR*\frac{1+INRsin2\left({\mathrm{\α}}_s\right)}{INR+1}$

其中，SNR为信噪比P_s/N₀，P_s是次用户信号的功率，INR为干扰噪声比P_p/N₀，P_p是主用户信号的功率，α_s,α_p分别为接收极化天线与水平极化信号的角度和发射极化天线与水平极化信号的角度。

优选的，步骤(2)中计算主用户发射信号的平均检测概率关于信噪比γ的公式如下：

${\stackrel{\‾}{P}}_{dr}=e^{\frac{-{\mathrm{\λ}}_{th}}{2}}\underset{n=0}{\overset{u-2}{\Σ}}\frac{1}{n!}{\left(\frac{{\mathrm{\λ}}_{th}}{2}\right)}^n+{\left(\frac{1+\stackrel{\‾}{\mathrm{\γ}}}{\stackrel{\‾}{\mathrm{\γ}}}\right)}^{u-1}\×\[e^{-\frac{{\mathrm{\λ}}_{th}}{2(1+\stackrel{\‾}{\mathrm{\γ}})}}-e^{\frac{{\mathrm{\λ}}_{th}}{2}}\underset{n=0}{\overset{u-2}{\Σ}}\frac{1}{n!}{\left(\frac{{\mathrm{\λ}}_{th}\stackrel{\‾}{\mathrm{\γ}}}{2(1+\stackrel{\‾}{\mathrm{\γ}})}\right)}^n\]$

其中λ_th为决策阈值，γ表示信噪比SINR，为平均信噪比，u为时间T秒和带宽B的乘积，u＝TB，n表示第一形式修改贝叶斯函数的阶数。

优选的，第一阈值选择90％，当时，表示主用户存在，停止卫星与卫星用户终端通讯；如果表示主用户不存在，卫星与卫星用户终端能够使用该频段通讯。

优选的，步骤(2)还包括计算，主用户没有发射信号时，次用户探测到发射信号的概率P_f，公式如下：

$P_f=P_r(Y_0>{\mathrm{\λ}}_{th}|H_0)=\frac{\mathrm{\Γ}(u,\frac{{\mathrm{\λ}}_{th}}{2})}{\mathrm{\Γ}\left(u\right)}$

其中 $\mathrm{\Γ}\left(u\right)={\∫}_0^{\mathrm{\∞}}{e}^{-t}{t}^{u-1}dt,\mathrm{\Γ}(u,\frac{{\mathrm{\λ}}_{th}}{2})=\frac{1}{\mathrm{\Γ}\left(\frac{{\mathrm{\λ}}_{th}}{2}\right)}{\∫}_0^u{e}^{-t}{t}^{\frac{{\mathrm{\λ}}_{th}}{2}-1}dt,$ t为积分参数；

步骤(3)中通过主用户发射信号的平均检测概率与第一阈值比较的同时，还比较次用户探测到发射信号的概率P_f与第二阈值，比较方法如下：

如果P_f≥40％，则表示的可信度较低，卫星不使用该频段，停止卫星与卫星用户终端通讯；如果表示主用户存在，停止卫星与卫星用户终端通讯；如果表示主用户不存在，卫星与卫星用户终端能够使用该频段通讯。

优选的，步骤(2)中推导得出主用户信号存在时，主用户发射信号的平均检测概率关于信噪比SINR的公式的方法为，采用信号能量检测法。

优选的，所述波段为C波段、L波段或K波段。

本发明与现有技术相比具有如下优点：

(1)本发明次用户利用最优极化接收的机制，可提高次用户的接收信号对干扰加噪声功率比SINR约为30dB；

(2)本发明次用户采用最优信噪比，进行能量检测，与现有技术比，提高了信号检测概率，进而改善了频谱感知的效率；

(3)本发明实现了卫星与地面基站通讯频率的共用，提高了频谱使用效率，频段资源有限的情况下，检测成功的概率高达90％。

附图说明

图1为异构卫星和地面网络系统结构示意图；

图2为Poincare极化球结构示意图；

图3为本发明的双极化天线下极化接收机制的能量检测方法示意图。

具体实施方式

下面结合附图对本发明作进一步详细说明。

基于最优极化接收机制次用户感知频谱的能量检测技术，包括下列步骤：

步骤1：卫星通信采用发射接收极化天线，根据接收天线的极化状态，频谱非授权用户，即次用户卫星，获得星地通信使用的C波段下，接收信号对干扰加噪声功率比(SINR)，得出最大值。

基于认知通信技术异构的星地场景中包含主用户链路和次用户链路，主用户发射的信号和次用户发射的信号均为全极化后的信号，分别用和表达，用JONES向量，以次用户发射的信号为例，可具体表示如下：

${\stackrel{\→}{E}}_s={\[E_x,E_y\]}^T=\left|{\stackrel{\→}{E}}_s\right|{\[{\mathrm{cos\α}}_s,{\mathrm{sin\α}}_s{e}^{{\mathrm{j\φ}}_s}\]}^T---\left(1\right)$

其中E_x,E_y是互为正交信号，分别为水平极化信号和垂直极化信号。

是极化相位描述符，每个极化状态描述符对应Poincare极化球的一个点，Poincare极化球如图2所示。其中弧PH是球上点P到点H的圆弧，弧度为2α_s，点H是弧PH与极化球赤道的交叉点，极化球赤道表示水平极化状态，弧PH与赤道的夹角为φ_s。

可用s_s:(α_s,φ_s)表示次用户发射信号的极化状态，该状态已知，主用户极化状态用s_p:(α_p,φ_p)表示，该状态未知。次用户接收机的极化状态表示为s_r:(α_r,φ_r)。

次用户接收信号的接收功率S：

$\begin{array}{c}S=P_s\·m_p(s_r,s_s)\\ =P_s\·{\cos }^2\left(\frac{{\mathrm{\θ}}_s}{2}\right)\end{array}---\left(2\right)$

其中P_s是次用户信号的功率，m_p(s_r,s_s)(∈[0,1])为次用户接的信号与次用户接收机极化不匹配引起的功率损失，可表示为：

$m_p(s_r,s_s)={\cos }^2\left(\frac{{\mathrm{\θ}}_s}{2}\right)---\left(3\right)$

其中θ_s为分别代表极化状态s_r和s_s的在极化球上的两个点的圆弧的对向角。当θ_s＝0时，m_p(s_r,s_s)＝1，表示次用户接收机和次用户接收的信号为共轭状态，信号极化匹配；θ_s＝π时，m_p(s_r,s_s)＝0，表示次用户接收机和次用户接收的信号为正交状态，信号极化不匹配；

同样地，主用户接收到的信号的接收功率：

$P=P_p\·m_p(s_r,s_p)=P_p\·{\cos }^2\left(\frac{{\mathrm{\θ}}_p}{2}\right)---\left(4\right)$

$m_p(s_r,s_p)={\cos }^2\left(\frac{{\mathrm{\θ}}_p}{2}\right)---\left(5\right)$

其中P_p是主用户信号的功率，θ_p为分别代表极化状态s_r和s_p的在极化球上的两个点的圆弧的对向角。当θ_p＝0时，m_p(s_r,s_p)＝1，表示次用户接收机和主用户接收的信号为共轭状态，信号极化匹配；θ_p＝π时，m_p(s_r,s_s)＝0，表示次用户接收机和主用户接收的信号为正交状态，信号极化不匹配；

次用户的信号与干扰加噪声比SINR：

$SINR=\frac{S}{P+N_0}$

$\begin{array}{c}=SNR\frac{{\cos }^2\left(\frac{{\mathrm{\θ}}_s}{2}\right)}{INR\·{\cos }^2\left(\frac{{\mathrm{\θ}}_p}{2}\right)+1}\\ =SNR\·D\end{array}---\left(6\right)$

其中N₀是噪声功率谱密度，SNR为信噪比P_s/N₀，INR为干扰噪声比P_p/N₀。

为了得到最优的SINR，问题可表示为：

$\underset{D_{max}}{max}SINR=SNR*D---\left(7\right)$

次用户通过调整接收天线与接收的信号角度α_s，获得最大的信号对干扰加噪声功率比(SINR)。也就是当θ_p＝θ_s+α_s，此时：

$D_{max}=\frac{1+INRsin2\left({\mathrm{\α}}_s\right)}{INR+1}---\left(8\right)$

上面公式(8)转换，推导过程如下：

$tan\left(\frac{{\mathrm{\θ}}_s}{2}\right)=\frac{cos\left(\frac{{\mathrm{\θ}}_s+{\mathrm{\α}}_s}{2}\right)sin\left(\frac{{\mathrm{\θ}}_s+{\mathrm{\α}}_s}{2}\right)}{{\cos }^2\left(\frac{{\mathrm{\θ}}_s+{\mathrm{\α}}_s}{2}\right)+\frac{1}{INR}}---\left(9\right)$

${\mathrm{\θ}}_s=2arctan\left(\frac{{\mathrm{sin\α}}_s}{1+2/INR-{\mathrm{cos\α}}_s}\right)---\left(10\right)$

$\begin{array}{c}tan\left(\frac{{\mathrm{\θ}}_s}{2}\right)=\frac{{\mathrm{sin\α}}_s}{1+2/INR-{\mathrm{cos\α}}_s}\\ =\frac{A\left({\mathrm{\α}}_s\right)}{B\left({\mathrm{\α}}_s\right)}\end{array}---\left(11\right)$

定义：

$\begin{array}{c}C\left({\mathrm{\α}}_s\right)=\sqrt{A{\left({\mathrm{\α}}_s\right)}^2+B{\left({\mathrm{\α}}_s\right)}^2}\\ =\sqrt{(4+8/INR){\sin }^2({\mathrm{\α}}_s/2)+{(2/INR)}^2}\end{array}---\left(12\right)$

$sin\left(\frac{{\mathrm{\θ}}_s}{2}\right)=\frac{2sin({\mathrm{\α}}_s/2)cos({\mathrm{\α}}_s/2)}{C\left({\mathrm{\α}}_s\right)}---\left(13\right)$

$cos\left(\frac{{\mathrm{\θ}}_s}{2}\right)=\frac{2(1/INR+{\sin }^2({\mathrm{\α}}_s/2\left)\right)}{C\left({\mathrm{\α}}_s\right)}---\left(14\right)$

$cos\left(\frac{{\mathrm{\θ}}_p}{2}\right)=cos\left(\frac{{\mathrm{\θ}}_s+{\mathrm{\α}}_s}{2}\right)=\frac{(2/INR)cos\left({\mathrm{\α}}_s/2\right))}{C\left({\mathrm{\α}}_s\right)}---\left(15\right)$

$\begin{array}{c}{D}_{\max }=\frac{{(2/INR+2{\sin }^2({\mathrm{\α}}_s/2\left)\right)}^2}{INR(2/INR)\cos ({\mathrm{\α}}_s/2){)}^2+C{\left({\mathrm{\α}}_s\right)}^2}\\ =\frac{1+INR\sin 2\left({\mathrm{\α}}_s\right)}{INR+1}\end{array}---\left(16\right)$

步骤2：通过信号能量检测法，得到能量检测决策统计数值Y₁服从自由度为2TB的卡方分布χ²

双极化天线下极化接收机制的能量检测方法如图3所示：

通过双极化天线接收的水平和垂直极化信号，首先进入带通滤波器BPF，调整接收天线的极化状态，得到接收信号的最优信噪比，然后通过频谱能量法，得到决策统计数值的分布，通过事先规定的决策值，确定主用户信号是否存在。

s_p为主用户发射信号向量；

y_ps为次用户在双极化天线下接收的信号向量；

发射的信号向量s_p由正交的水平极化部分s_pH和垂直极化部分s_pV组成，以Jones向量表示如下：

$s_p=\left[\begin{array}{c}{s}_{pH}\\ s_{pV}\end{array}\right]=\left[\begin{array}{c}{s}_{pH}^i+s_{pV}^q\\ s_{pH}^q+s_{pV}^i\end{array}\right]---\left(17\right)$

其中i和q是同相和正交相位部分。

定义如下符号：X_t为主用户发射天线的极化状态；X_r为主用户接收天线的极化状态；C_t为次用户发射天线的极化状态；C_r为次用户接收天线的极化状态，则发射天线的极化状态向量表示为：

$X_t=\left[\begin{array}{c}cos{\mathrm{\α}}_p\\ {\mathrm{sin\α}}_p{e}^{{\mathrm{j\φ}}_p}\end{array}\right]---\left(18\right)$

其中，φ_p∈(0,2π)，并可通过下列公式计算：

${\mathrm{\α}}_p=arctan\frac{\left|s_{pV}\right|}{\left|s_{pH}\right|}---\left(19\right)$

${\mathrm{\φ}}_p={\mathrm{\φ}}_{pV}-{\mathrm{\φ}}_{pH}=arctan\frac{\left|s_{pV}^q\right|}{\left|s_{pV}^i\right|}-arctan\frac{\left|s_{pH}^q\right|}{\left|s_{pH}^i\right|}---\left(20\right)$

在观测时间间隔t的情况下，检测问题可归为有2种假设如下：

H₀:y₀(t)＝n(t)(21)

H₁:y₁(t)＝s_p(t)+n(t)

其中，H₀为假设主用户不存在，H₁为假设主用户存在。n(t)为带通下高斯噪声向量信号，可以表示为：

n(t)＝n_i(t)cos2πf_ct-n_q(t)sin2πf_ct(22)

其中，f_c为参考频率，n_i(t)为信号同相部分，n_q(t)为信号正交相位部分。

假设通频带宽为B，那么n_i(t)和n_q(t)均占B/2，在一定时间T内，噪声信号能量可以表示如下：

$\stackrel{\‾}{Y_0}={\∫}_0^T{n}^2\left(t\right)dt\≈\frac{1}{2B}\underset{m=1}{\overset{TB}{\Σ}}({a_{im}}^2+{a_{qm}}^2)=N_o\underset{i=1}{\overset{TB}{\Σ}}({b_{im}}^2+{b_{qm}}^2)---\left(23\right)$

其中，a_im和a_qm分别代表1/B下的n_i(t)和n_q(t)的m次采样，b_qm＝a_qm/σ_i ²，σ_i ²＝var(a_im)＝var(a_qm)＝2N₀B，N₀是噪声功率谱密度。能量检测决策统计数值Y₀表示如下：

$Y_0\≈\underset{i=1}{\overset{TB}{\Σ}}{b_{im}}^2+\underset{i=1}{\overset{TB}{\Σ}}{b_{qm}}^2---\left(24\right)$

从上面公式(24)，可得到噪声下，信号决策统计值Y₀服从自由度为2TB的中心卡方分布χ²。同理经公式(23)和(24)计算可得

$\stackrel{\‾}{Y_1}={{\∫}_0^T\[n\left(t\right)+s_p\left(t\right)\]}^{2}dt$

在主信号和噪声都存在下，决策统计值Y₁服从自由度为2TB的非中心参数为2γ的卡方分布χ²。

Y₀和Y₁服从卡方分布χ²如下所示

$\left\{\begin{array}{c}{Y}_0~{{\mathrm{\χ}}_{2u}}^2,H_0\\ Y_1~{{\mathrm{\χ}}_{2u}}^2\left(2\mathrm{\γ}\right),H_1\end{array}\right.---\left(25\right)$

其中，2u为自由度，γ为次用户接收信号信噪比。

进一步，Y₀和Y₁的概率密度函数如下：

$\left\{\begin{array}{c}{Y}_0=\frac{1}{2^u\mathrm{\Γ}\left(u\right)}{y_0}^{u-1}{e}^{-y/2},H_0\\ Y_1=\frac{1}{2}{\left(\frac{y_1}{2\mathrm{\γ}}\right)}^{\frac{u-1}{2}}{e}^{-\frac{2\mathrm{\γ}+y_1}{2}}{I}_{u-1}\left(\sqrt{2{\mathrm{\γy}}_1}\right),H_1\end{array}\right.---\left(26\right)$

其中，Γ(.)为伽马函数，I_n(.)为n阶第一形式修改贝叶斯函数(ModifiedBesselFunctionoftheFirstKind)，其中y₀表示y₀(t)，其中y₁表示y₁(t)，

H₀假设下，主用户没有发射信号时，次用户探测到发射信号的概率(误警率)P_f表示如下：

$P_f=P_r(Y_0>{\mathrm{\λ}}_{th}|H_0)=\frac{\mathrm{\Γ}(u,\frac{{\mathrm{\λ}}_{th}}{2})}{\mathrm{\Γ}\left(u\right)}---\left(27\right)$

其中 $\mathrm{\Γ}(u,\frac{{\mathrm{\λ}}_{th}}{2})=\frac{1}{\mathrm{\Γ}\left(\frac{{\mathrm{\λ}}_{th}}{2}\right)}{\∫}_0^u{e}^{-t}{t}^{\frac{{\mathrm{\λ}}_{th}}{2}-1}dt$

H₁假设下，主用户发射信号时，次用户检测概率P_d表示如下：

$P_d=P_r(Y_1>{\mathrm{\λ}}_{th}|H_1)=Q_u(\sqrt{2\mathrm{\γ}},\sqrt{{\mathrm{\λ}}_{th}})---\left(28\right)$

Q_u(.,.)为广义MarcumQ函数，λ_th为决策阈值，为0到1之间的常数。

对小规模衰减，接收的信号服从瑞利分布，信噪比γ的指数概率密度函数表示如下：

$f\left(\mathrm{\γ}\right)=\frac{1}{\stackrel{\‾}{\mathrm{\γ}}}\exp \left(-\frac{\mathrm{\γ}}{\stackrel{\‾}{\mathrm{\γ}}}\right),\mathrm{\γ}\≥0---\left(29\right)$

其中，为次用户接收信号的平均信噪比。

H₁假设下，主用户发射信号的平均检测概率表示如下：

${\stackrel{\‾}{P}}_{dr}={\∫}_0^{\mathrm{\∞}}{Q}_u(\sqrt{2\mathrm{\γ}},\sqrt{\mathrm{\λ}})f\left(\mathrm{\γ}\right)d\mathrm{\γ}---\left(30\right)$

设上面公式可换为：

${\stackrel{\‾}{P}}_{dr}=\frac{1}{\stackrel{\‾}{\mathrm{\γ}}}{\∫}_0^{\mathrm{\∞}}{Q}_u(x,\sqrt{\mathrm{\λ}}){\mathrm{xe}}^{-\frac{x^2}{2\stackrel{\‾}{\mathrm{\γ}}}}dx---\left(31\right)$

则

${\stackrel{\‾}{P}}_{dr}=e^{\frac{-\mathrm{\λ}}{2}}\underset{n=0}{\overset{u-2}{\Σ}}\frac{1}{n!}{\left(\frac{\mathrm{\λ}}{2}\right)}^n+{\left(\frac{1+\stackrel{\‾}{\mathrm{\γ}}}{\stackrel{\‾}{\mathrm{\γ}}}\right)}^{u-1}\×\[e^{-\frac{\mathrm{\λ}}{2(1+\stackrel{\‾}{\mathrm{\γ}})}}-e^{\frac{\mathrm{\λ}}{2}}\underset{n=0}{\overset{u-2}{\Σ}}\frac{1}{n!}{\left(\frac{\mathrm{\λ}\stackrel{\‾}{\mathrm{\γ}}}{2(1+\stackrel{\‾}{\mathrm{\γ}})}\right)}^n\]---\left(32\right)$

步骤3：结合接收信号对干扰加噪声功率比的最大值，输入到接收信号的决策统计数值分布参数中，在一定时间T内和带宽B，通过决策统计数值大于事先确定的某个阈值概率(P_d≥90％)，即可确定授权用户(主用户)信号是否存在。

次用户接收天线极化状态C_r表示如下：

$C_r=PA=\left[\begin{array}{cc}1& 0\\ 0& e^{{\mathrm{j\φ}}_s}\end{array}\right]\left[\begin{array}{c}cos{\mathrm{\α}}_s\\ {\mathrm{sin\α}}_s\end{array}\right]=\left[\begin{array}{c}cos{\mathrm{\α}}_s\\ {\mathrm{sin\α}}_s{e}^{{\mathrm{j\φ}}_s}\end{array}\right]---\left(33\right)$

在最优极化接收状态下，能量检测决策统计数值Y₁:

$Y=\underset{i=1}{\overset{M}{\Σ}}|C_r^{\±}{y}_1{|}^2>{\mathrm{\λ}}_{th}\left(34\right)$

其中，M为采样数，为复杂共轭转置。λ_th为决策阈值。

次用户采用最优极化接收信号最优信噪比为maxSINR＝γ，γ即为信噪比SINR，将其带入到步骤2的决策统计数值分布参数，公式(32)中得到主用户发射信号的检测概率：

$\begin{array}{c}{\stackrel{\‾}{P}}_{dr}={\∫}_0^{\mathrm{\∞}}{Q}_u(\sqrt{2\mathrm{\γ}},\sqrt{{\mathrm{\λ}}_{th}})f\left(\mathrm{\γ}\right)d\mathrm{\γ}\\ =e^{\frac{-{\mathrm{\λ}}_{th}}{2}}\underset{n=0}{\overset{u-2}{\Σ}}\frac{1}{n!}{\left(\frac{{\mathrm{\λ}}_{th}}{2}\right)}^n+{\left(\frac{1+\stackrel{\‾}{\mathrm{\γ}}}{\stackrel{\‾}{\mathrm{\γ}}}\right)}^{u-1}\×\[e^{-\frac{{\mathrm{\λ}}_{th}}{2(1+\stackrel{\‾}{\mathrm{\γ}})}}-e^{\frac{{\mathrm{\λ}}_{th}}{2}}\underset{n=0}{\overset{u-2}{\Σ}}\frac{1}{n!}{\left(\frac{{\mathrm{\λ}}_{th}\stackrel{\‾}{\mathrm{\γ}}}{2(1+\stackrel{\‾}{\mathrm{\γ}})}\right)}^n\]\end{array}---\left(35\right)$

其中，u为每个T秒和带宽B的乘积，u＝TB。

如果P_f≥40％，则表示的可信度较低，卫星不使用该频段，停止卫星与卫星终端的通讯；如果表示主用户存在，停止卫星与卫星终端的通讯；如果表示主用户不存在，卫星与卫星终端的可以使用该频段通讯。

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

本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。

Claims (7)

1.基于最优化极化接收机制的频谱感知能量检测方法，其特征在于包括下列步骤：

(1)地面基站为某波段的授权主用户，卫星为所述波段的非授权用户，为次用户，卫星、卫星用户终端、地面基站的发送和接收均采用极化天线；次用户根据接收天线的极化状态，计算次用户的信噪比SINR，其中，S为次用户接收功率，P为主用户接收功率，N₀为噪声平均功率，得出SINR的最大值maxSINR；

(2)在一定时间T内和带宽B下，推导得出主用户信号存在时，主用户发射信号的平均检测概率关于信噪比SINR的公式；

(3)使用步骤(1)计算的maxSINR代替公式中的SINR，通过主用户发射信号的平均检测概率与第一阈值比较，确定主用户信号是否存在。

2.根据权利要求1所述的基于最优化极化接收机制的频谱感知能量检测技术，其特征在于：得出SINR的最大值maxSINR的具体方法为：次用户调整接收极化天线与接收的信号角度α_s，使得θ_p＝θ_s+α_s，其中θ_s为分别代表极化状态s_r和s_s的在极化球上的两个点的圆弧的对向角，θ_p为分别代表极化状态s_r和s_p的在极化球上的两个点的圆弧的对向角，s_r为次用户接收机的极化状态，s_s表示次用户发射信号的极化状态，s_p表示主用户极化状态；

在所述波段干扰对噪声比INR确定情况下，利用如下公式得到信噪比SINR的最大值maxSINR；

$\max SINR=SNR*\frac{1+INRsin2\left({\mathrm{\α}}_s\right)}{INR+1}$

其中，SNR为信噪比P_s/N₀，P_s是次用户信号的功率，INR为干扰噪声比P_p/N₀，P_p是主用户信号的功率。

3.根据权利要求1所述的基于最优化极化接收机制的频谱感知能量检测技术，其特征在于，步骤(2)中推导得出主用户发射信号的平均检测概率关于信噪比SINR的公式如下：

${\stackrel{\‾}{P}}_{dr}=e^{\frac{-{\mathrm{\λ}}_{th}}{2}}\underset{n=0}{\overset{u-2}{\Σ}}\frac{1}{n!}{\left(\frac{{\mathrm{\λ}}_{th}}{2}\right)}^n+{\left(\frac{1+\stackrel{\‾}{\mathrm{\γ}}}{\stackrel{\‾}{\mathrm{\γ}}}\right)}^{u-1}\×\[e^{-\frac{{\mathrm{\λ}}_{th}}{2\left(1+\stackrel{\‾}{\mathrm{\γ}}\right)}}-e^{\frac{{\mathrm{\λ}}_{th}}{2}}\underset{n=0}{\overset{u-2}{\Σ}}\frac{1}{n!}{\left(\frac{{\mathrm{\λ}}_{th}\stackrel{\‾}{\mathrm{\γ}}}{2\left(1+\stackrel{\‾}{\mathrm{\γ}}\right)}\right)}^n\]$

其中λ_th为决策阈值，γ表示信噪比SINR，为平均信噪比，u为时间T秒和带宽B的乘积，u＝TB，n表示第一形式修改贝叶斯函数的阶数。

4.根据权利要求1所述的基于最优化极化接收机制的频谱感知能量检测技术，其特征在于，第一阈值选择90％，当时，表示主用户存在，停止卫星与卫星用户终端通讯；如果表示主用户不存在，卫星与卫星用户终端能够使用该频段通讯。

5.根据权利要求3所述的基于最优化极化接收机制的频谱感知能量检测技术，其特征在于，步骤(2)还包括计算，主用户没有发射信号时，次用户探测到发射信号的概率P_f，公式如下：

$P_f=P_r(Y_0>{\mathrm{\λ}}_{th}|H_0)=\frac{\mathrm{\Γ}(u,\frac{{\mathrm{\λ}}_{th}}{2})}{\mathrm{\Γ}\left(u\right)}$

其中 $\mathrm{\Γ}\left(u\right)={\∫}_0^{\mathrm{\∞}}{e}^{-t}{t}^{u-1}dt,\mathrm{\Γ}(u,\frac{{\mathrm{\λ}}_{th}}{2})=\frac{1}{\mathrm{\Γ}\left(\frac{{\mathrm{\λ}}_{th}}{2}\right)}{\∫}_0^u{e}^{-t}{t}^{\frac{{\mathrm{\λ}}_{th}}{2}-1}dt,$ t为积分参数；

步骤(3)中通过主用户发射信号的平均检测概率与第一阈值比较的同时，还比较次用户探测到发射信号的概率P_f与第二阈值，比较方法如下：

如果P_f≥40％，则表示的可信度较低，卫星不使用该频段，停止卫星与卫星用户终端通讯；如果P_f<40％,表示主用户存在，停止卫星与卫星用户终端通讯；如果P_f<40％,表示主用户不存在，卫星与卫星用户终端能够使用该频段通讯。

6.根据权利要求1所述的基于最优化极化接收机制的频谱感知能量检测技术，其特征在于，步骤(2)中推导得出主用户信号存在时，主用户发射信号的平均检测概率关于信噪比SINR的公式的方法为，采用信号能量检测法。

7.根据权利要求1所述的基于最优化极化接收机制的频谱感知能量检测技术，其特征在于，所述波段为C波段、L波段或K波段。