CN107733541B - 频谱感知的方法、装置、设备以及计算机可读存储介质 - Google Patents

频谱感知的方法、装置、设备以及计算机可读存储介质 Download PDF

Info

Publication number
CN107733541B
CN107733541B CN201711226944.0A CN201711226944A CN107733541B CN 107733541 B CN107733541 B CN 107733541B CN 201711226944 A CN201711226944 A CN 201711226944A CN 107733541 B CN107733541 B CN 107733541B
Authority
CN
China
Prior art keywords
principal component
preset number
matrix
maximum
matrixes
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201711226944.0A
Other languages
English (en)
Other versions
CN107733541A (zh
Inventor
孙晨皓
王永华
万频
杜艺期
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangdong University of Technology
Original Assignee
Guangdong University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Guangdong University of Technology filed Critical Guangdong University of Technology
Priority to CN201711226944.0A priority Critical patent/CN107733541B/zh
Publication of CN107733541A publication Critical patent/CN107733541A/zh
Application granted granted Critical
Publication of CN107733541B publication Critical patent/CN107733541B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/30—Monitoring; Testing of propagation channels
    • H04B17/382—Monitoring; Testing of propagation channels for resource allocation, admission control or handover

Landscapes

  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Monitoring And Testing Of Transmission In General (AREA)

Abstract

本发明公开了一种频谱感知的方法、装置、设备及存储介质，包括：对次用户的接收端接收到的信号进行采样，获取预设个数的采样信号矩阵；利用主成分分析法，分别得到与预设个数的采样信号矩阵对应的样本协方差矩阵，分别计算预设个数的样本协方差矩阵最大特征值对应的特征向量，以便分别提取与预设个数的特征向量对应的主成分；将预设个数的主成分组建为主成分矩阵，得到与所述主成分矩阵对应的主成分协方差矩阵；将所述主成分协方差矩阵的最大特征值与最小特征值的比值作为最大最小特征值算法的检验统计量，判断所述次用户的接收端信号中是否存在主用户。利用本发明所提供的方法、装置、设备及存储介质，提高了频谱检测器的检测性能。

Description

频谱感知的方法、装置、设备以及计算机可读存储介质
技术领域
本发明涉及认知无线电技术领域，特别是涉及一种基于主成分分析的最大最小特征值频谱感知的方法、装置、设备以及计算机可读存储介质。
背景技术
随着无线电技术的进步和社会的发展，各行各业对无线电频率的需求越来越多，频谱资源匮乏的问题日益严重。目前，频谱资源采用国家统一分配授权的管理模式，将频谱分为两种类型：授权频段和非授权频段。其中，授权频段占据着大部分频谱资源，如电视广播频段，但不少授权频段处于空闲状态；开放使用的非授权频段占整个频谱资源的很少一部分，如无线局域网、无线城域网等无线网络大多使用非授权频段在工作，该频段上的用户很多，业务量也很大，无线电频段已基本趋于饱和。据美国联邦通信委员会的研究表明，已分配的频谱利用率为15％～85％。所以说，频谱资源的匮乏并不是真正意义上的频谱资源不足，更多是由频谱利用率过低造成的。
比较经典的频谱感知的方法有能量检测(Energy Detection，ED)算法，该算法是一种最基本、最简单的信号检测方法，它无需知道信号的先验信息，复杂度低，容易实现，在现阶段的研究中广泛应用于信号检测中。能量检测方法是一种非连贯性、非相干的检测算法，检测的对象是特定频带内的频谱能量，而与具体信号的波形无关，特别适用于未知信号的检测。
经典能量检测的缺点十分明显，如检测性能较差，不能区分信号和噪声，在低信噪比情况下容易淹没在噪声当中，并且在噪声不确定情况下很难用于信号检测。因此，能量检测只适用于对检测性能要求不太高的认知无线电系统中，比如IEEE802.22系统。
综上所述可以看出，如何提高频谱检测器的检测性能是目前有待解决的问题。
发明内容
本发明的目的是提供一种基于主成分分析的最大最小特征值频谱感知的方法、装置、设备以及计算机可读存储介质，已解决现有技术中频谱检测器检测性能较低的问题。
为解决上述技术问题，本发明提供一种基于主成分分析的最大最小特征值频谱感知的方法，包括：对次用户的接收端接收到的信号进行采样，获取预设个数的采样信号矩阵；利用主成分分析算法，分别得到与预设个数的采样信号矩阵对应的样本协方差矩阵，以便分别计算预设个数的样本协方差矩阵的最大特征值；分别计算与预设个数的最大特征值对应的特征向量，以便分别提取与预设个数的特征向量对应的主成分；将预设个数的主成分组建为主成分矩阵，得到与所述主成分矩阵对应的主成分协方差矩阵；将所述主成分协方差矩阵的最大特征值与最小特征值的比值作为最大最小特征值算法的检验统计量，通过所述检验统计量与自适应阈值的比较结果，判断所述次用户的接收端信号中是否存在主用户。
优选地，所述对次用户的接收端接收到的信号进行采样，获取预设个数的采样信号矩阵包括：
对次用户的m个接收天线的观测值进行采样，获取采样信号矩阵
Figure GDA0001543786970000021
其中xi＝[xi(1),xi(2),…,xi(N)]表示第i个接收天线观测值的1×N向量，N为接收样本的总数。
优选地，所述利用主成分分析算法，分别得到与预设个数的采样信号矩阵对应的样本协方差矩阵，以便分别计算预设个数的样本协方差矩阵的最大特征值包括：
利用
Figure GDA0001543786970000022
分别估算与预设个数的采样信号矩阵对应的样本协方差矩阵C；分别计算预设个数的样本协方差矩阵的最大特征值
Figure GDA0001543786970000023
其中，矩阵C为m个天线上接收的PU信号的协方差矩阵Cs和噪声分量的协方差矩阵Cw的和；C的特征值为Cs和Cw的特征值的和λ(C)＝λ(Cs)+λ(Cw)；
Figure GDA0001543786970000024
Cw的特征值
Figure GDA0001543786970000025
等于
Figure GDA0001543786970000026
Cs的特征值
Figure GDA0001543786970000027
为零，除一个等于Cs的
Figure GDA0001543786970000028
优选地，所述分别计算与预设个数的最大特征值对应的特征向量，以便分别提取与预设个数的特征向量对应的主成分包括：
通过(C-λmI2)vm＝0分别计算预设个数的矩阵C的最大特征值λm对应的特征向量vm，I2为单位矩阵；根据
Figure GDA0001543786970000031
分别计算预设个数的最大特征向量λm对应的主成分Pm。
优选地，所述将预设个数的主成分组建为主成分矩阵，得到与所述主成分矩阵对应的主成分协方差矩阵包括：将预设个数的主成分Pm组建为主成分矩阵Pval，通过计算得到与所述主成分矩阵对应的主成分协方差矩阵。
本发明还提供了一种基于主成分分析的最大最小特征值频谱感知的装置，包括：
采样模块，用于对次用户的接收端接收到的信号进行采样，获取预设个数的采样信号矩阵；
计算模块，用于利用主成分分析算法，分别得到与预设个数的采样信号矩阵对应的样本协方差矩阵，以便分别计算预设个数的样本协方差矩阵的最大特征值；
提取模块，用于分别计算与预设个数的最大特征值对应的特征向量，以便分别提取与预设个数的特征向量对应的主成分；
组建模块，用于将预设个数的主成分组建为主成分矩阵，得到与所述主成分矩阵对应的主成分协方差矩阵；
判断模块，用于将所述主成分协方差矩阵的最大特征值与最小特征值的比值作为最大最小特征值算法的检验统计量，通过所述检验统计量与自适应阈值的比较结果，判断所述次用户的接收端信号中是否存在主用户。
优选地，所述采样模块具体用于：
对次用户的m个接收天线的观测值进行采样，获取采样信号矩阵
Figure GDA0001543786970000032
其中xi＝[xi(1),xi(2),…,xi(N)]表示第i个接收天线观测值的1×N向量，N为接收样本的总数。
优选地，所述计算模块具体用于：
利用
Figure GDA0001543786970000033
分别估算与预设个数的采样信号矩阵对应的样本协方差矩阵C；分别计算预设个数的样本协方差矩阵的最大特征值
Figure GDA0001543786970000041
其中，矩阵C为m个天线上接收的PU信号的协方差矩阵Cs和噪声分量的协方差矩阵Cw的和；C的特征值为Cs和Cw的特征值的和λ(C)＝λ(Cs)+λ(Cw)；
Figure GDA0001543786970000042
Cw的特征值
Figure GDA0001543786970000043
等于
Figure GDA0001543786970000044
Cs的特征值
Figure GDA0001543786970000045
为零，除一个等于Cs的
Figure GDA0001543786970000046
本发明还提供了一种基于主成分分析的最大最小特征值频谱感知的设备，包括：
存储器，用于存储计算机程序；处理器，用于执行所述计算机程序时实现上述一种基于主成分分析的最大最小特征值频谱感知的方法的步骤。
本发明还提供了一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现上述一种基于主成分分析的最大最小特征值频谱感知的方法的步骤。
本发明所提供的一种基于主成分分析的最大最小特征值频谱感知的方法、装置、设备以及计算机可读存储介质，对次用户的接收端接收到的信号进行采样，得到预设个数的采样信号矩阵，分别对采集到的采样信号矩阵运用主成分分析算法，得到所述采样信号矩阵的样本协方差矩阵，计算所述样本协方差矩阵的最大特征值，求取所述最大特征值对应的特征向量，提取所述特征向量对应的主成分。在存在授权用户的环境中采集到的采样信号矩阵，经过主成分分析生产的主成分中，只有一个主成分分量包含滤波的授权用户，该主成分分量对应于最大特征值；因此，对次用户接收端接收到的预设个数的采样信号矩阵进行主成分分析后，提取预设个数的样本协方差矩阵最大特征值对应的特征向量所对应的主成分，并且将提取到的预设个数的主成分组件为主成分矩阵作为最大最小特征值算法的输入值，生成主成分协方差矩阵，计算所述主成分协方差矩阵的最大特征值和最小特征值的比值，将所述比值作为检验统计量，比较所述检验统计量与自适应阈值的大小，从而判断次用户的接收端接收到的信号是否存在频谱信号。利用本发明所提供的方法进行频谱感知时，先对采集到的信号运用主成分分析算法，提高采样信号的信噪比，在对经过处理提高信噪比后的采样信号进行最大最小特征值分析，大大的提高了最大最小特征值算法的检测性能。因此本发明所提供的基于主成分分析的最大最小特征值频谱感知的方法提高了频谱感知的检测性能。
附图说明
为了更清楚的说明本发明实施例或现有技术的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单的介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明所提供的基于主成分分析的最大最小特征值频谱感知的方法的第一种具体实施例的流程图；
图2为本发明所提供的基于主成分分析的最大最小特征值频谱感知的方法的第二种具体实施例的流程图；
图3为不同的频谱感知算法检测性能的实验结果图；
图4为本发明实施例提供的一种基于主成分分析的最大最小特征值频谱感知的装置的结构框图。
具体实施方式
本发明的核心是提供一种基于主成分分析的最大最小特征值频谱感知的方法、装置、设备以及计算机可读存储介质，提高了频谱检测器的检测性能。
为了使本技术领域的人员更好地理解本发明方案，下面结合附图和具体实施方式对本发明作进一步的详细说明。显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
请参考图1，图1为本发明所提供的基于主成分分析的最大最小特征值频谱感知的方法的第一种具体实施例的流程图；具体操作步骤如下：
步骤S101：对次用户的接收端接收到的信号进行采样，获取预设个数的采样信号矩阵；
在本实施例中，可以选择获取两个采样信号矩阵，也可以选择获取多个的采样信号矩阵。
假设在认知无线电网络中有1个次用户(Secondary user非授权用户或认知用户)和m个接收天线；对次用户的m个接收天线的观测值进行采样，获取采样信号矩阵
Figure GDA0001543786970000061
其中xi＝[xi(1),xi(2),…,xi(N)]表示第i个接收天线观测值的1×N向量，N为接收样本的总数。
步骤S102：利用主成分分析算法，分别得到与预设个数的采样信号矩阵对应的样本协方差矩阵，以便分别计算预设个数的样本协方差矩阵的最大特征值；
通过m个天线观测的采样信号矩阵，可以应用主成分分析法(PCA)来生成m个主成分(PCs)。当主用户(Primary user PU授权用户)存在时，只有一个主成分分量PC包含具有增强信噪比(SNR)的有噪声的PU信号，而其它主成分分量是噪声的线性组合。因此，SU应该能够选择适当的主成分分量PC来执行频谱感知。
关于PU存在和不存在的问题公式可以由如下经典的贝叶斯检测问题表示：Hη:xi＝ηhis+wi；其中η∈{0；1}，H0代表PU不存在，而H1表示PU存在。xi表示第i个SU接收天线观测值的1×N向量，N是接收样本的总数；s是包含PU信号的1×N向量。1×N向量wi表示均值为0，方差为
Figure GDA0001543786970000062
的高斯白噪声信号(AWGN)。hi是PU和第i个SU接收天线之间的信道增益。
采样信号矩阵为m×N矩阵，协方差矩阵C为m×m矩阵：
Figure GDA0001543786970000063
由于η2＝η，为了简单起见，将η2代入η；协方差矩阵C可以估计为：
Figure GDA0001543786970000064
步骤S103：分别计算与预设个数的最大特征值对应的特征向量，以便分别提取与预设个数的特征向量对应的主成分；
通过使用PU信号和噪声之间的独立性假设，矩阵C可以被写为m个天线上接收的PU信号的协方差矩阵Cs和噪声分量的协方差矩阵Cw的和C＝Cs+Cw；
Figure GDA0001543786970000065
矩阵C的特征值为Cs和Cw的特征值的和λ(C)＝λ(Cs)+λ(Cw)；
Cw的特征值
Figure GDA0001543786970000066
等于
Figure GDA0001543786970000067
Cs的特征值
Figure GDA0001543786970000068
为零，除一个等于Cs的轨迹
Figure GDA0001543786970000069
这是因为Cs是一阶矩阵。因此，矩阵C的特征值为：
Figure GDA0001543786970000071
特征值对应的特征向量可以通过求解(C-λiI2)vi＝0等式得到，其中是对应于第i个特征值的第i个特征向量，I2是单位矩阵。一旦找到特征向量，主成分分量PC可以通过
Figure GDA0001543786970000072
获得。
将样本协方差的最大特征值λm代入(C-λiI2)vi＝0，通过(C-λmI2)vm＝0计算矩阵C的最大特征值λm对应的特征向量vm，根据
Figure GDA0001543786970000073
分别计算预设个数的最大特征向量λm对应的主成分Pm。
在H1下，PCA产生m个主成分(PCs)，其中只有一个分量包含滤波的PU信号，这个分量Pm对应于最大特征值
Figure GDA0001543786970000074
其中
Figure GDA0001543786970000075
代表PU的功率，
Figure GDA0001543786970000076
表示存在于Pm中的噪声分量信号的功率。其他m-1分量是在m个天线观察到的噪声的混合。
在H0下，矩阵C产生一个对角矩阵：C0＝CW；由于
Figure GDA0001543786970000077
则C0的特征值为:
Figure GDA0001543786970000078
由于C0是一个对角矩阵，所以m×m个单位矩阵Im可以是收集特征向量的矩阵：
Figure GDA0001543786970000079
通过
Figure GDA00015437869700000710
可知，H0下的PC只不过是噪声成分。然而，由于m个SU接收天线的m个噪声分量是白高斯和独立的，所以该组特征向量的任何旋转不会影响PC在H0下的统计特性。因此，Pi成为wi的线性组合，并且保持为高斯白噪声。
步骤S104：将预设个数的主成分组建为主成分矩阵，得到与所述主成分矩阵对应的主成分协方差矩阵；
将预设个数的主成分Pm组建为主成分矩阵Pval，通过计算得到与所述主成分矩阵对应的主成分协方差矩阵。
步骤S105：将所述主成分协方差矩阵的最大特征值与最小特征值的比值作为最大最小特征值算法(MME)的检验统计量，通过所述检验统计量与自适应阈值的比较结果，判断所述次用户的接收端信号中是否存在主用户。
在本实施例中，采样信号矩阵通过PCA产生多达m个分量(相同数目的观察值)，检测器必须选择验证的分量来执行频谱感知。如上所述，在H0下的主成分分量
Figure GDA0001543786970000081
是等效的，因为wi是具有相同方差的AWGN。与H0不同，H1导致非等效PCs，
Figure GDA0001543786970000082
仅仅是噪声的分量，而
Figure GDA0001543786970000083
是PU信号和噪声的组合。因此，通过对
Figure GDA0001543786970000084
应用测试统计量，SU能够诊断通道状态。所以，SU应该能够选择良好的PCA输出，导致对PU状态的有效判定。
请参考图2，图2为本发明所提供的基于主成分分析的最大最小特征值频谱感知的方法的第一种具体实施例的流程图。
在上述实施例的基础上，本实施例对次用户的m个接收天线接收的信号进行两次采集，得到第一采样信号矩阵和第二采样信号矩阵；当然在本申请的其他实施例中，也可以对次用户的接收端进行多次采样，如三次，四次。本实施例的具体操作步骤如下：
步骤S201：对次用户的M个接收天线接收到的信号进行两次采样，得到第一采样信号矩阵和第二采样信号矩阵；
步骤S202：利用PCA算法估算所述第一采样信号矩阵的样本协方差矩阵C1，利用PCA算法估算所述第二采样信号矩阵的样本协方差矩阵C2；
步骤S203：利用
Figure GDA0001543786970000085
计算C1的最大特征值
Figure GDA0001543786970000086
计算C2的最大特征值
Figure GDA0001543786970000087
步骤S204：利用(C-λmI2)vm＝0计算
Figure GDA0001543786970000088
对应的特征向量
Figure GDA0001543786970000089
计算
Figure GDA00015437869700000810
对应的特征向量
Figure GDA00015437869700000811
步骤S205：利用
Figure GDA00015437869700000812
计算
Figure GDA00015437869700000813
对应的主成分分量
Figure GDA00015437869700000814
计算
Figure GDA00015437869700000815
对应的主成分分量
Figure GDA00015437869700000816
步骤S206：将
Figure GDA00015437869700000817
和
Figure GDA00015437869700000818
组建为主成分矩阵Pval，获得主成分协方差矩阵；计算所述主成分协方差矩阵的最大特征值λmax和最小特征值λmin，将
Figure GDA00015437869700000819
作为最大最小特征值算法的检验统计量；
步骤S207：比较与自适应阈值的大小，判断次用户m个接收天线接收的信号中是否存在主用户信号。
一旦SU选择适当的PC，则通过MME计算测试统计量T，并与自适应阈值ξ进行比较以对PU状态做出决定。因此，新通道假设可以表示为：
Figure GDA0001543786970000091
其中w对应于H0下获得的噪声分量，y和r分别表示在H1下存在于Pval中的PU信号和噪声。
在本实施例中先对采集到的信号运用PCA算法，提高采样信号的信噪比，在对经过处理提高信噪比后的采样信号进行MME分析，大大的提高了MME算法的检测性能。如图3所示，通过不同的频谱感知算法对同样的采样信号进行检测，求得检测概率pd＝P(η*＝1|η＝1)；试验中频谱感知的算法包括ED、MME、ED与PCA结合以及本实施例中PCA与MME结合的算法，通过实验结果，可以看出本实施例所提供的基于主成分分析的最大最小特征值频谱感知的方法相对于用信号能量作为特征有较好的检测性能。
请参考图4，图4为本发明实施例提供的一种基于主成分分析的最大最小特征值频谱感知的装置的结构框图；具体装置可以包括：
采样模,100，用于对次用户的接收端接收到的信号进行采样，获取预设个数的采样信号矩阵；
计算模块200，用于利用主成分分析算法，分别得到与预设个数的采样信号矩阵对应的样本协方差矩阵，以便分别计算预设个数的样本协方差矩阵的最大特征值；
提取模块300，用于分别计算与预设个数的最大特征值对应的特征向量，以便分别提取与预设个数的特征向量对应的主成分；
组建模块400，用于将预设个数的主成分组建为主成分矩阵，得到与所述主成分矩阵对应的主成分协方差矩阵；
判断模块500，用于将所述主成分协方差矩阵的最大特征值与最小特征值的比值作为最大最小特征值算法的检验统计量，通过所述检验统计量与自适应阈值的比较结果，判断所述次用户的接收端信号中是否存在主用户。
本实施例的基于主成分分析的最大最小特征值频谱感知的装置用于实现前述的基于主成分分析的最大最小特征值频谱感知的方法，因此基于主成分分析的最大最小特征值频谱感知的装置中的具体实施方式可见前文中的基于主成分分析的最大最小特征值频谱感知的方法的实施例部分，例如，采样模块100，计算模块200，提取模块300，组建模块400，判断模块500，分别用于实现上述基于主成分分析的最大最小特征值频谱感知的方法中步骤S101，S102，S103，S104和S105，所以，其具体实施方式可以参照相应的各个部分实施例的描述，在此不再赘述。
本发明具体实施例还提供了一种基于主成分分析的最大最小特征值频谱感知的设备，包括：存储器，用于存储计算机程序；处理器，用于执行所述计算机程序时实现上述一种基于主成分分析的最大最小特征值频谱感知的方法的步骤。
本发明具体实施例还提供了一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现上述一种基于主成分分析的最大最小特征值频谱感知的方法的步骤。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其它实施例的不同之处，各个实施例之间相同或相似部分互相参见即可。对于实施例公开的装置而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
专业人员还可以进一步意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、计算机软件或者二者的结合来实现，为了清楚地说明硬件和软件的可互换性，在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。
结合本文中所公开的实施例描述的方法或算法的步骤可以直接用硬件、处理器执行的软件模块，或者二者的结合来实施。软件模块可以置于随机存储器(RAM)、内存、只读存储器(ROM)、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其它形式的存储介质中。
以上对本发明所提供的基于主成分分析的最大最小特征值频谱感知的方法、装置、设备以及计算机可读存储介质进行了详细介绍。本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以对本发明进行若干改进和修饰，这些改进和修饰也落入本发明权利要求的保护范围内。

Claims (5)

1.一种基于主成分分析的最大最小特征值频谱感知的方法，其特征在于，包括：
对次用户的接收端接收到的信号进行采样，获取预设个数的采样信号矩阵；
利用主成分分析算法，分别得到与预设个数的采样信号矩阵对应的样本协方差矩阵，以便分别计算预设个数的样本协方差矩阵的最大特征值；
分别计算与预设个数的最大特征值对应的特征向量，以便分别提取与预设个数的特征向量对应的主成分；
将预设个数的主成分组建为主成分矩阵，得到与所述主成分矩阵对应的主成分协方差矩阵；
将所述主成分协方差矩阵的最大特征值与最小特征值的比值作为最大最小特征值算法的检验统计量，通过所述检验统计量与自适应阈值的比较结果，判断所述次用户的接收端信号中是否存在主用户；
所述对次用户的接收端接收到的信号进行采样，获取预设个数的采样信号矩阵包括：
对次用户的m个接收天线的观测值进行采样，获取采样信号矩阵
Figure FDA0002969276040000011
其中xi＝[xi(1),xi(2),...,xi(N)]表示第i个接收天线观测值的1×N向量，N为接收样本的总数；
所述利用主成分分析算法，分别得到与预设个数的采样信号矩阵对应的样本协方差矩阵，以便分别计算预设个数的样本协方差矩阵的最大特征值包括：
利用
Figure FDA0002969276040000012
分别估算与预设个数的采样信号矩阵对应的样本协方差矩阵C；
分别计算预设个数的样本协方差矩阵的最大特征值
Figure FDA0002969276040000013
其中，矩阵C为m个天线上接收的PU信号的协方差矩阵Cs和噪声分量的协方差矩阵Cw的和；C的特征值为Cs和Cw的特征值的和λ(C)＝λ(Cs)+λ(Cw)；
Figure FDA0002969276040000014
Cw的特征值λi w等于
Figure FDA0002969276040000015
其中，1≤i≤m，Cs的特征值
Figure FDA0002969276040000016
为零，除一个等于Cs的轨迹
Figure FDA0002969276040000017
C的特征值为
Figure FDA0002969276040000021
C的最大特征值为
Figure FDA0002969276040000022
Figure FDA0002969276040000023
表示存在于Pm中的噪声分量信号的功率，hi表示PU和第i个SU接收天线之间的信道增益，s表示包含PU信号的1×N向量，
Figure FDA0002969276040000024
表示包含PU信号的1×N向量的方差；
所述分别计算与预设个数的最大特征值对应的特征向量，以便分别提取与预设个数的特征向量对应的主成分包括：
通过(C-λmI2)vm＝0分别计算预设个数的矩阵C的最大特征值λm对应的特征向量vm，I2为单位矩阵；
根据
Figure FDA0002969276040000025
分别计算预设个数的最大特征向量λm对应的主成分Pm。
2.如权利要求1所述的方法，其特征在于，所述将预设个数的主成分组建为主成分矩阵，得到与所述主成分矩阵对应的主成分协方差矩阵包括：
将预设个数的主成分Pm组建为主成分矩阵Pval，通过计算得到与所述主成分矩阵对应的主成分协方差矩阵。
3.一种基于主成分分析的最大最小特征值频谱感知的装置，其特征在于，包括：
采样模块，用于对次用户的接收端接收到的信号进行采样，获取预设个数的采样信号矩阵；
计算模块，用于利用主成分分析算法，分别得到与预设个数的采样信号矩阵对应的样本协方差矩阵，以便分别计算预设个数的样本协方差矩阵的最大特征值；
提取模块，用于分别计算与预设个数的最大特征值对应的特征向量，以便分别提取与预设个数的特征向量对应的主成分；
组建模块，用于将预设个数的主成分组建为主成分矩阵，得到与所述主成分矩阵对应的主成分协方差矩阵；
判断模块，用于将所述主成分协方差矩阵的最大特征值与最小特征值的比值作为最大最小特征值算法的检验统计量，通过所述检验统计量与自适应阈值的比较结果，判断所述次用户的接收端信号中是否存在主用户；
所述采样模块具体用于：
对次用户的m个接收天线的观测值进行采样，获取采样信号矩阵
Figure FDA0002969276040000026
其中xi＝[xi(1),xi(2),...,xi(N)]表示第i个接收天线观测值的1×N向量，N为接收样本的总数；
所述计算模块具体用于：
利用
Figure FDA0002969276040000031
分别估算与预设个数的采样信号矩阵对应的样本协方差矩阵C；
分别计算预设个数的样本协方差矩阵的最大特征值
Figure FDA0002969276040000032
其中，矩阵C为m个天线上接收的PU信号的协方差矩阵Cs和噪声分量的协方差矩阵Cw的和；C的特征值为Cs和Cw的特征值的和λ(C)＝λ(Cs)+λ(Cw)；
Figure FDA0002969276040000033
Cw的特征值λi w等于
Figure FDA0002969276040000034
其中，1≤i≤m，Cs
Figure FDA0002969276040000035
的噪声分量信号的功率，hi表示PU和第i个SU接收天线之间的信道增益，s表示包含PU信号的1×N向量，
Figure FDA0002969276040000036
表示包含PU信号的1×N向量的方差；
所述提取模块具体用于：
通过(C-λmI2)vm＝0分别计算预设个数的矩阵C的最大特征值λm对应的特征向量vm，I2为单位矩阵；
根据
Figure FDA0002969276040000037
分别计算预设个数的最大特征向量λm对应的主成分Pm。
4.一种基于主成分分析的最大最小特征值频谱感知的设备，其特征在于，包括：
存储器，用于存储计算机程序；
处理器，用于执行所述计算机程序时实现如权利要求1至2任一项所述一种基于主成分分析的最大最小特征值频谱感知的方法的步骤。
5.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现如权利要求1至2任一项所述一种基于主成分分析的最大最小特征值频谱感知的方法的步骤。
CN201711226944.0A 2017-11-29 2017-11-29 频谱感知的方法、装置、设备以及计算机可读存储介质 Active CN107733541B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711226944.0A CN107733541B (zh) 2017-11-29 2017-11-29 频谱感知的方法、装置、设备以及计算机可读存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711226944.0A CN107733541B (zh) 2017-11-29 2017-11-29 频谱感知的方法、装置、设备以及计算机可读存储介质

Publications (2)

Publication Number Publication Date
CN107733541A CN107733541A (zh) 2018-02-23
CN107733541B true CN107733541B (zh) 2021-07-06

Family

ID=61220088

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711226944.0A Active CN107733541B (zh) 2017-11-29 2017-11-29 频谱感知的方法、装置、设备以及计算机可读存储介质

Country Status (1)

Country Link
CN (1) CN107733541B (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108494509A (zh) * 2018-03-09 2018-09-04 广东工业大学 一种协作频谱感知方法及装置
CN108462544B (zh) * 2018-03-27 2021-09-17 广东工业大学 一种频谱感知方法及装置
CN108768563A (zh) * 2018-05-17 2018-11-06 广东工业大学 一种协作频谱感知方法及相关装置
CN108736992A (zh) * 2018-05-17 2018-11-02 广东工业大学 一种协作频谱感知方法及相关装置
CN109004997A (zh) * 2018-08-09 2018-12-14 广东工业大学 一种频谱感知方法、装置和计算机可读存储介质
CN109309538A (zh) * 2018-08-28 2019-02-05 广东工业大学 一种频谱感知方法、装置、设备、系统及存储介质
CN110048788B (zh) * 2019-03-15 2021-08-24 广东工业大学 一种基于聚类算法的联合频谱感知方法
CN111062003A (zh) * 2019-12-13 2020-04-24 武汉轻工大学 样本总体协方差判定方法、装置、设备及存储介质
CN112333120B (zh) * 2020-11-09 2021-08-24 电子科技大学 基于pca的信道增益矩阵的特征提取方法
CN115728383B (zh) * 2022-10-14 2023-12-05 港珠澳大桥管理局 桥梁结构损伤定位方法、装置、计算机设备和介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103795479A (zh) * 2014-01-28 2014-05-14 南京邮电大学 一种基于特征值的协作频谱感知方法
CN104135327A (zh) * 2014-07-10 2014-11-05 上海大学 基于支持向量机的频谱感知方法
CN105357676A (zh) * 2015-12-10 2016-02-24 上海交通大学 基于矩阵填充的室内空闲频谱检测方法
CN106950446A (zh) * 2017-03-22 2017-07-14 国网上海市电力公司 基于主成分分析法的电网异常快速检测与定位方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9077508B2 (en) * 2012-11-15 2015-07-07 Mitsubishi Electric Research Laboratories, Inc. Adaptively coding and modulating signals transmitted via nonlinear channels

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103795479A (zh) * 2014-01-28 2014-05-14 南京邮电大学 一种基于特征值的协作频谱感知方法
CN104135327A (zh) * 2014-07-10 2014-11-05 上海大学 基于支持向量机的频谱感知方法
CN105357676A (zh) * 2015-12-10 2016-02-24 上海交通大学 基于矩阵填充的室内空闲频谱检测方法
CN106950446A (zh) * 2017-03-22 2017-07-14 国网上海市电力公司 基于主成分分析法的电网异常快速检测与定位方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于随机矩阵理论的频谱感知技术研究;韩新胜;《中国优秀硕士学位论文全文数据库信息科技辑》;20170215(第2017年02期);第1章1.3节，第2章2.3节，第4章4.1-4.2节 *

Also Published As

Publication number Publication date
CN107733541A (zh) 2018-02-23

Similar Documents

Publication Publication Date Title
CN107733541B (zh) 频谱感知的方法、装置、设备以及计算机可读存储介质
Köse et al. RF fingerprinting of IoT devices based on transient energy spectrum
Nadler et al. Performance of eigenvalue-based signal detectors with known and unknown noise level
CN111600667B (zh) 一种基于cnn-lstm的频谱感知方法
CN108462544B (zh) 一种频谱感知方法及装置
CN107395301B (zh) 一种基于k均值算法的频谱感知方法与装置
CN107979431B (zh) 基于黎曼中值的频谱感知的方法、装置以及设备
CN107426736B (zh) 一种认知无线电的频谱感知方法及系统
CN108322277B (zh) 一种基于协方差矩阵反特征值的频谱感知方法
CN107360577B (zh) 一种基于机器学习的频谱感知方法及装置
Lin et al. Generalized FMD detection for spectrum sensing under low signal-to-noise ratio
Luís et al. RF-spectrum opportunities for cognitive radio networks operating over GSM channels
Ikuma et al. A comparison of three classes of spectrum sensing techniques
CN107171752B (zh) 一种认知无线电的频谱感知方法及系统
CN109525339B (zh) 认知无线电的频谱感知方法、装置、设备以及存储介质
Hamid et al. Spectrum sensing through spectrum discriminator and maximum minimum eigenvalue detector: A comparative study
EP2086255B1 (en) Process for sensing vacant sub-space over the spectrum bandwidth and apparatus for performing the same
CN112994813B (zh) 自适应采样的频谱感知方法及相关装置
CN114268393B (zh) 一种基于连通分量个数特征的认知无线电频谱感知方法
Sindhubargavi et al. Spectrum sensing using energy detection technique for cognitive radio networks using PCA technique
Subekti et al. A Jarque-Bera test based spectrum sensing for cognitive radio
Iwata et al. A study on Welch FFT segment size selection method for spectrum awareness
CN109286937B (zh) 利用小特征值估计噪声功率的协方差矩阵频谱感知方法
CN106100775B (zh) 基于邻接矩阵的ofdm频谱感知方法
CN113132033A (zh) 基于多项式插值处理的通信干扰检测方法和装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant