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

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

Info

Publication number
CN107733541A
CN107733541A CN201711226944.0A CN201711226944A CN107733541A CN 107733541 A CN107733541 A CN 107733541A CN 201711226944 A CN201711226944 A CN 201711226944A CN 107733541 A CN107733541 A CN 107733541A
Authority
CN
China
Prior art keywords
principal component
matrix
predetermined number
eigenvalue
covariance matrix
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.)
Granted
Application number
CN201711226944.0A
Other languages
English (en)
Other versions
CN107733541B (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

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/30Monitoring; Testing of propagation channels
    • H04B17/382Monitoring; Testing of propagation channels for resource allocation, admission control or handover

Abstract

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

Description

频谱感知的方法、装置、设备以及计算机可读存储介质
技术领域
本发明涉及认知无线电技术领域,特别是涉及一种基于主成分分 析的最大最小特征值频谱感知的方法、装置、设备以及计算机可读存 储介质。
背景技术
随着无线电技术的进步和社会的发展,各行各业对无线电频率的 需求越来越多,频谱资源匮乏的问题日益严重。目前,频谱资源采用 国家统一分配授权的管理模式,将频谱分为两种类型:授权频段和非 授权频段。其中,授权频段占据着大部分频谱资源,如电视广播频段, 但不少授权频段处于空闲状态;开放使用的非授权频段占整个频谱资 源的很少一部分,如无线局域网、无线城域网等无线网络大多使用非 授权频段在工作,该频段上的用户很多,业务量也很大,无线电频段 已基本趋于饱和。据美国联邦通信委员会的研究表明,已分配的频谱 利用率为15%~85%。所以说,频谱资源的匮乏并不是真正意义上的频谱资源不足,更多是由频谱利用率过低造成的。
比较经典的频谱感知的方法有能量检测(Energy Detection,ED) 算法,该算法是一种最基本、最简单的信号检测方法,它无需知道信 号的先验信息,复杂度低,容易实现,在现阶段的研究中广泛应用于 信号检测中。能量检测方法是一种非连贯性、非相干的检测算法,检 测的对象是特定频带内的频谱能量,而与具体信号的波形无关,特别 适用于未知信号的检测。
经典能量检测的缺点十分明显,如检测性能较差,不能区分信号 和噪声,在低信噪比情况下容易淹没在噪声当中,并且在噪声不确定 情况下很难用于信号检测。因此,能量检测只适用于对检测性能要求 不太高的认知无线电系统中,比如IEEE802.22系统。
综上所述可以看出,如何提高频谱检测器的检测性能是目前有待 解决的问题。
发明内容
本发明的目的是提供一种基于主成分分析的最大最小特征值频 谱感知的方法、装置、设备以及计算机可读存储介质,已解决现有技 术中频谱检测器检测性能较低的问题。
为解决上述技术问题,本发明提供一种基于主成分分析的最大最 小特征值频谱感知的方法,包括:对次用户的接收端接收到的信号进 行采样,获取预设个数的采样信号矩阵;利用主成分分析算法,分别 得到与预设个数的采样信号矩阵对应的样本协方差矩阵,以便分别计 算预设个数的样本协方差矩阵的最大特征值;分别计算与预设个数的 最大特征值对应的特征向量,以便分别提取与预设个数的特征向量对 应的主成分;将预设个数的主成分组建为主成分矩阵,得到与所述主 成分矩阵对应的主成分协方差矩阵;将所述主成分协方差矩阵的最大 特征值与最小特征值的比值作为最大最小特征值算法的检验统计量, 通过所述检验统计量与自适应阈值的比较结果,判断所述次用户的接 收端信号中是否存在主用户。
优选地,所述对次用户的接收端接收到的信号进行采样,获取预 设个数的采样信号矩阵包括:
对次用户的m个接收天线的观测值进行采样,获取采样信号矩阵 其中xi=[xi(1),xi(2), … ,xi(N)]表示第i个接收天线观测值的 1×N向量,N为接收样本的总数。
优选地,所述利用主成分分析算法,分别得到与预设个数的采样 信号矩阵对应的样本协方差矩阵,以便分别计算预设个数的样本协方 差矩阵的最大特征值包括:
利用分别估算与预设个数的采样信号矩阵对应的样 本协方差矩阵C;分别计算预设个数的样本协方差矩阵的最大特征 值
其中,矩阵C为m个天线上接收的PU信号的协方差矩阵Cs和噪 声分量的协方差矩阵Cw的和;C的特征值为Cs和Cw的特征值的和 λ(C)=λ(Cs)+λ(Cw);Cw的特征值等于Cs的特征值为零,除一个等于Cs的轨迹C的特征值为C的最大特征值为
优选地,所述分别计算与预设个数的最大特征值对应的特征向 量,以便分别提取与预设个数的特征向量对应的主成分包括:
通过(C-λmI2)vm=0分别计算预设个数的矩阵C的最大特征值λm对 应的特征向量vm,I2为单位矩阵;根据分别计算预设个数的最 大特征向量λm对应的主成分Pm
优选地,所述将预设个数的主成分组建为主成分矩阵,得到与所 述主成分矩阵对应的主成分协方差矩阵包括:将预设个数的主成分Pm组建为主成分矩阵Pval,通过计算得到与所述主成分矩阵对应的主成分 协方差矩阵。
本发明还提供了一种基于主成分分析的最大最小特征值频谱感 知的装置,包括:
采样模块,用于对次用户的接收端接收到的信号进行采样,获取 预设个数的采样信号矩阵;
计算模块,用于利用主成分分析算法,分别得到与预设个数的采 样信号矩阵对应的样本协方差矩阵,以便分别计算预设个数的样本协 方差矩阵的最大特征值;
提取模块,用于分别计算与预设个数的最大特征值对应的特征向 量,以便分别提取与预设个数的特征向量对应的主成分;
组建模块,用于将预设个数的主成分组建为主成分矩阵,得到与 所述主成分矩阵对应的主成分协方差矩阵;
判断模块,用于将所述主成分协方差矩阵的最大特征值与最小特 征值的比值作为最大最小特征值算法的检验统计量,通过所述检验统 计量与自适应阈值的比较结果,判断所述次用户的接收端信号中是否 存在主用户。
优选地,所述采样模块具体用于:
对次用户的m个接收天线的观测值进行采样,获取采样信号矩阵 其中xi=[xi(1),xi(2), … ,xi(N)]表示第i个接收天线观测值的 1×N向量,N为接收样本的总数。
优选地,所述计算模块具体用于:
利用分别估算与预设个数的采样信号矩阵对应的样 本协方差矩阵C;分别计算预设个数的样本协方差矩阵的最大特征 值
其中,矩阵C为m个天线上接收的PU信号的协方差矩阵Cs和噪 声分量的协方差矩阵Cw的和;C的特征值为Cs和Cw的特征值的和 λ(C)=λ(Cs)+λ(Cw);Cw的特征值等于Cs的特特征值为零,除一个等于Cs的轨迹C的特征值为C的最大特征值为
本发明还提供了一种基于主成分分析的最大最小特征值频谱感 知的设备,包括:
存储器,用于存储计算机程序;处理器,用于执行所述计算机程 序时实现上述一种基于主成分分析的最大最小特征值频谱感知的方法 的步骤。
本发明还提供了一种计算机可读存储介质,所述计算机可读存储 介质上存储有计算机程序,所述计算机程序被处理器执行时实现上述 一种基于主成分分析的最大最小特征值频谱感知的方法的步骤。
本发明所提供的一种基于主成分分析的最大最小特征值频谱感 知的方法、装置、设备以及计算机可读存储介质,对次用户的接收端 接收到的信号进行采样,得到预设个数的采样信号矩阵,分别对采集 到的采样信号矩阵运用主成分分析算法,得到所述采样信号矩阵的样 本协方差矩阵,计算所述样本协方差矩阵的最大特征值,求取所述最 大特征值对应的特征向量,提取所述特征向量对应的主成分。在存在 授权用户的环境中采集到的采样信号矩阵,经过主成分分析生产的主 成分中,只有一个主成分分量包含滤波的授权用户,该主成分分量对 应于最大特征值;因此,对次用户接收端接收到的预设个数的采样信号矩阵进行主成分分析后,提取预设个数的样本协方差矩阵最大特征 值对应的特征向量所对应的主成分,并且将提取到的预设个数的主成 分组件为主成分矩阵作为最大最小特征值算法的输入值,生成主成分 协方差矩阵,计算所述主成分协方差矩阵的最大特征值和最小特征值 的比值,将所述比值作为检验统计量,比较所述检验统计量与自适应 阈值的大小,从而判断次用户的接收端接收到的信号是否存在频谱信 号。利用本发明所提供的方法进行频谱感知时,先对采集到的信号运 用主成分分析算法,提高采样信号的信噪比,在对经过处理提高信噪 比后的采样信号进行最大最小特征值分析,大大的提高了最大最小特 征值算法的检测性能。因此本发明所提供的基于主成分分析的最大最 小特征值频谱感知的方法提高了频谱感知的检测性能。
附图说明
为了更清楚的说明本发明实施例或现有技术的技术方案,下面将 对实施例或现有技术描述中所需要使用的附图作简单的介绍,显而易 见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普 通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附 图获得其他的附图。
图1为本发明所提供的基于主成分分析的最大最小特征值频谱感 知的方法的第一种具体实施例的流程图;
图2为本发明所提供的基于主成分分析的最大最小特征值频谱感 知的方法的第二种具体实施例的流程图;
图3为不同的频谱感知算法检测性能的实验结果图;
图4为本发明实施例提供的一种基于主成分分析的最大最小特征 值频谱感知的装置的结构框图。
具体实施方式
本发明的核心是提供一种基于主成分分析的最大最小特征值频 谱感知的方法、装置、设备以及计算机可读存储介质,提高了频谱检 测器的检测性能。
为了使本技术领域的人员更好地理解本发明方案,下面结合附图
请参考图1,图1为本发明所提供的基于主成分分析的最大最小 特征值频谱感知的方法的第一种具体实施例的流程图;具体操作步骤 如下:
步骤S101:对次用户的接收端接收到的信号进行采样,获取预设 个数的采样信号矩阵;
在本实施例中,可以选择获取两个采样信号矩阵,也可以选择获 取多个的采样信号矩阵。
假设在认知无线电网络中有1个次用户(Secondary user非授权用户或认知 用户)和m个接收天线;对次用户的m个接收天线的观测值进行采样,获取采 样信号矩阵其中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, 方差为的高斯白噪声信号(AWGN)。hi是PU和第i个SU接收天 线之间的信道增益。
采样信号矩阵为m×N矩阵,协方差矩阵C为m×m矩阵:
由于η2=η,为了简单起见,将η2代入η;协方差矩阵C可以估计 为:
步骤S103:分别计算与预设个数的最大特征值对应的特征向量, 以便分别提取与预设个数的特征向量对应的主成分;
通过使用PU信号和噪声之间的独立性假设,矩阵C可以被写为 m个天线上接收的PU信号的协方差矩阵Cs和噪声分量的协方差矩阵 Cw的和C=Cs+Cw矩阵C的特征值为Cs和Cw的特征值的和 λ(C)=λ(Cs)+λ(Cw);
Cw的特征值等于Cs的特征值为零,除一个等于Cs的轨迹这是因为Cs是一阶矩阵。因此,矩阵C的特征值为:
特征值对应的特征向量可以通过求解(C-λiI2)vi=0等式得到,其中 是对应于第i个特征值的第i个特征向量,I2是单位矩阵。一旦找到特 征向量,主成分分量PC可以通过获得。
将样本协方差的最大特征值λm代入(C-λiI2)vi=0,通过 (C-λmI2)vm=0计算矩阵C的最大特征值λm对应的特征向量vm,根据 分别计算预设个数的最大特征向量λm对应的主成分Pm
在H1下,PCA产生m个主成分(PCs),其中只有一个分量包含滤 波的PU信号,这个分量Pm对应于最大特征值其 中代表PU的功率,表示存在于Pm中的噪声分量信号的功 率。其他m-1分量是在m个天线观察到的噪声的混合。
在H0下,矩阵C产生一个对角矩阵:C0=CW;由于则C0的特征值为:由于C0是一个对角矩阵,所以m×m 个单位矩阵Im可以是收集特征向量的矩阵:
通过可知,H0下的PC只不过是噪声成分。然而,由于m 个SU接收天线的m个噪声分量是白高斯和独立的,所以该组特征向 量的任何旋转不会影响PC在H0下的统计特性。因此,Pi成为wi的线 性组合,并且保持为高斯白噪声。
步骤S104:将预设个数的主成分组建为主成分矩阵,得到与所述 主成分矩阵对应的主成分协方差矩阵;
将预设个数的主成分Pm组建为主成分矩阵Pval,通过计算得到与所 述主成分矩阵对应的主成分协方差矩阵。
步骤S105:将所述主成分协方差矩阵的最大特征值与最小特征值 的比值作为最大最小特征值算法(MME)的检验统计量,通过所述检 验统计量与自适应阈值的比较结果,判断所述次用户的接收端信号中 是否存在主用户。
在本实施例中,采样信号矩阵通过PCA产生多达m个分量(相 同数目的观察值),检测器必须选择验证的分量来执行频谱感知。如上 所述,在H0下的主成分分量是等效的,因为wi是具有相同方差的 AWGN。与H0不同,H1导致非等效PCs,仅仅是噪声的分 量,而是PU信号和噪声的组合。因此,通过对应用测试统计量, SU能够诊断通道状态。所以,SU应该能够选择良好的PCA输出,导 致对PU状态的有效判定。
请参考图2,图2为本发明所提供的基于主成分分析的最大最小 特征值频谱感知的方法的第一种具体实施例的流程图。
在上述实施例的基础上,本实施例对次用户的m个接收天线接收 的信号进行两次采集,得到第一采样信号矩阵和第二采样信号矩阵; 当然在本申请的其他实施例中,也可以对次用户的接收端进行多次采 样,如三次,四次。本实施例的具体操作步骤如下:
步骤S201:对次用户的M个接收天线接收到的信号进行两次采 样,得到第一采样信号矩阵和第二采样信号矩阵;
步骤S202:利用PCA算法估算所述第一采样信号矩阵的样本协 方差矩阵C1,利用PCA算法估算所述第二采样信号矩阵的样本协方差 矩阵C2
步骤S203:利用计算C1的最大特征值计算C2的最大特征值
步骤S204:利用(C-λmI2)vm=0计算对应的特征向量计算对应的特征向量
步骤S205:利用计算对应的主成分分量计算对 应的主成分分量
步骤S206:将组建为主成分矩阵Pval,获得主成分协方差 矩阵;计算所述主成分协方差矩阵的最大特征值λmax和最小特征值λmin, 将作为最大最小特征值算法的检验统计量;
步骤S207:比较与自适应阈值的大小,判断次用户m个接收天 线接收的信号中是否存在主用户信号。
一旦SU选择适当的PC,则通过MME计算测试统计量T,并与 自适应阈值ξ进行比较以对PU状态做出决定。因此,新通道假设可 以表示为:
其中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 (10)

1.一种基于主成分分析的最大最小特征值频谱感知的方法,其特征在于,包括:
对次用户的接收端接收到的信号进行采样,获取预设个数的采样信号矩阵;
利用主成分分析算法,分别得到与预设个数的采样信号矩阵对应的样本协方差矩阵,以便分别计算预设个数的样本协方差矩阵的最大特征值;
分别计算与预设个数的最大特征值对应的特征向量,以便分别提取与预设个数的特征向量对应的主成分;
将预设个数的主成分组建为主成分矩阵,得到与所述主成分矩阵对应的主成分协方差矩阵;
将所述主成分协方差矩阵的最大特征值与最小特征值的比值作为最大最小特征值算法的检验统计量,通过所述检验统计量与自适应阈值的比较结果,判断所述次用户的接收端信号中是否存在主用户。
2.如权利要求1所述的方法,其特征在于,所述对次用户的接收端接收到的信号进行采样,获取预设个数的采样信号矩阵包括:
对次用户的m个接收天线的观测值进行采样,获取采样信号矩阵其中xi=[xi(1),xi(2),…,xi(N)]表示第i个接收天线观测值的1×N向量,N为接收样本的总数。
3.如权利要求2所述的方法,其特征在于,所述利用主成分分析算法,分别得到与预设个数的采样信号矩阵对应的样本协方差矩阵,以便分别计算预设个数的样本协方差矩阵的最大特征值包括:
利用分别估算与预设个数的采样信号矩阵对应的样本协方差矩阵C;
分别计算预设个数的样本协方差矩阵的最大特征值
其中,矩阵C为m个天线上接收的PU信号的协方差矩阵Cs和噪声分量的协方差矩阵Cw的和;C的特征值为Cs和Cw的特征值的和λ(C)=λ(Cs)+λ(Cw);Cw的特征值等于Cs的特征值为零,除一个等于Cs的轨迹C的特征值为C的最大特征值为
4.如权利要求3所述的方法,其特征在于,所述分别计算与预设个数的最大特征值对应的特征向量,以便分别提取与预设个数的特征向量对应的主成分包括:
通过(C-λmI2)vm=0分别计算预设个数的矩阵C的最大特征值λm对应的特征向量vm,I2为单位矩阵;
根据分别计算预设个数的最大特征向量λm对应的主成分Pm
5.如权利要求4所述的方法,其特征在于,所述将预设个数的主成分组建为主成分矩阵,得到与所述主成分矩阵对应的主成分协方差矩阵包括:
将预设个数的主成分Pm组建为主成分矩阵Pval,通过计算得到与所述主成分矩阵对应的主成分协方差矩阵。
6.一种基于主成分分析的最大最小特征值频谱感知的装置,其特征在于,包括:
采样模块,用于对次用户的接收端接收到的信号进行采样,获取预设个数的采样信号矩阵;
计算模块,用于利用主成分分析算法,分别得到与预设个数的采样信号矩阵对应的样本协方差矩阵,以便分别计算预设个数的样本协方差矩阵的最大特征值;
提取模块,用于分别计算与预设个数的最大特征值对应的特征向量,以便分别提取与预设个数的特征向量对应的主成分;
组建模块,用于将预设个数的主成分组建为主成分矩阵,得到与所述主成分矩阵对应的主成分协方差矩阵;
判断模块,用于将所述主成分协方差矩阵的最大特征值与最小特征值的比值作为最大最小特征值算法的检验统计量,通过所述检验统计量与自适应阈值的比较结果,判断所述次用户的接收端信号中是否存在主用户。
7.如权利要求6所述的装置,其特征在于,所述采样模块具体用于:
对次用户的m个接收天线的观测值进行采样,获取采样信号矩阵其中xi=[xi(1),xi(2),…,xi(N)]表示第i个接收天线观测值的1×N向量,N为接收样本的总数。
8.如权利要求7所述的装置,其特征在于,所述计算模块具体用于:
利用分别估算与预设个数的采样信号矩阵对应的样本协方差矩阵C;
分别计算预设个数的样本协方差矩阵的最大特征值
其中,矩阵C为m个天线上接收的PU信号的协方差矩阵Cs和噪声分量的协方差矩阵Cw的和;C的特征值为Cs和Cw的特征值的和λ(C)=λ(Cs)+λ(Cw);Cw的特征值等于Cs的特征值为零,除一个等于Cs的轨迹C的特征值为C的最大特征值为
9.一种基于主成分分析的最大最小特征值频谱感知的设备,其特征在于,包括:
存储器,用于存储计算机程序;
处理器,用于执行所述计算机程序时实现如权利要求1至5任一项所述一种基于主成分分析的最大最小特征值频谱感知的方法的步骤。
10.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现如权利要求1至5任一项所述一种基于主成分分析的最大最小特征值频谱感知的方法的步骤。
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 true CN107733541A (zh) 2018-02-23
CN107733541B 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)

Cited By (10)

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

Citations (5)

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

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140133848A1 (en) * 2012-11-15 2014-05-15 Mitsubishi Electric Research Laboratories, Inc. Adaptively Coding and Modulating Signals Transmitted Via Nonlinear Channels
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
韩新胜: "基于随机矩阵理论的频谱感知技术研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 *

Cited By (12)

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

Also Published As

Publication number Publication date
CN107733541B (zh) 2021-07-06

Similar Documents

Publication Publication Date Title
CN107733541A (zh) 频谱感知的方法、装置、设备以及计算机可读存储介质
Xu et al. Estimating wind speed probability distribution by diffusion-based kernel density method
CN104135327B (zh) 基于支持向量机的频谱感知方法
CN106549722B (zh) 一种基于历史感知信息的双门限能量检测方法
CN107395301A (zh) 一种基于k均值算法的频谱感知方法
US20130225105A1 (en) Method, Apparatus and Computer Program Product for Identifying Frequency Bands, and a Method, Apparatus and Computer Program Product for Evaluating Performance
CN108322277B (zh) 一种基于协方差矩阵反特征值的频谱感知方法
De Vito A review of wideband spectrum sensing methods for cognitive radios
CN106533577A (zh) 一种基于能量检测的非高斯噪声抑制方法
CN107370521A (zh) 一种认知无线电多用户协作频谱感知方法
CN107682103B (zh) 一种基于最大特征值和主特征向量的双特征频谱感知方法
CN103491551B (zh) 一种基于特征向量的加权协作频谱感知方法
WO2016060620A1 (en) Frame based spike detection module
CN105429719A (zh) 基于功率谱和多尺度小波变换分析强干扰信号检测方法
CN110932806A (zh) 一种alpha稳定噪声衰落信道下多天线频谱感知方法
CN107808100B (zh) 针对特定测试样本的隐写分析方法
CN106788817A (zh) 一种基于贝叶斯准则和能量检测法的频谱感知方法
CN107171752A (zh) 一种认知无线电的频谱感知方法及系统
CN108768563A (zh) 一种协作频谱感知方法及相关装置
CN106972900B (zh) 基于广义t2统计量的盲频谱感知方法
CN109309538A (zh) 一种频谱感知方法、装置、设备、系统及存储介质
CN108400826A (zh) 一种基于循环矩阵特征值的频谱感知方法
CN110190917B (zh) 一种LTE230MHz电力无线专网的频谱空洞感知方法、装置及设备
Alammar et al. A minesweeper algorithm for improved signal area estimation in spectrum aware systems
CN105119668B (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