首页 > 专利 > 杭州电子科技大学 > 一种基于压缩感知和支持向量机的高效频谱感知方法专利详情

一种基于压缩感知和支持向量机的高效频谱感知方法   0    0

有效专利 查看PDF
专利申请流程有哪些步骤?
专利申请流程图
申请
申请号:指国家知识产权局受理一件专利申请时给予该专利申请的一个标示号码。唯一性原则。
申请日:提出专利申请之日。
2020-12-30
申请公布
申请公布指发明专利申请经初步审查合格后,自申请日(或优先权日)起18个月期满时的公布或根据申请人的请求提前进行的公布。
申请公布号:专利申请过程中,在尚未取得专利授权之前,国家专利局《专利公报》公开专利时的编号。
申请公布日:申请公开的日期,即在专利公报上予以公开的日期。
2021-06-04
授权
授权指对发明专利申请经实质审查没有发现驳回理由,授予发明专利权;或对实用新型或外观设计专利申请经初步审查没有发现驳回理由,授予实用新型专利权或外观设计专利权。
2022-07-22
预估到期
发明专利权的期限为二十年,实用新型专利权期限为十年,外观设计专利权期限为十五年,均自申请日起计算。专利届满后法律终止保护。
2040-12-30
基本信息
有效性 有效专利 专利类型 发明专利
申请号 CN202011618079.6 申请日 2020-12-30
公开/公告号 CN112821968B 公开/公告日 2022-07-22
授权日 2022-07-22 预估到期日 2040-12-30
申请年 2020年 公开/公告年 2022年
缴费截止日
分类号 H04B17/382G06N20/10G06K9/62 主分类号 H04B17/382
是否联合申请 独立申请 文献类型号 B
独权数量 1 从权数量 5
权利要求数量 6 非专利引证数量 1
引用专利数量 1 被引证专利数量 0
非专利引证 1、2011.02.17林畅.认知无线电系统中基于压缩感知的宽带频谱感知技术研究《.中国优秀硕士学位论文全文数据库(信息科技辑)》.2016,吕斌等.采用支持向量机的宽带频谱感知算法《.信号处理》.2014,(第12期),;
引用专利 US2011040711A 被引证专利
专利权维持 1 专利申请国编码 CN
专利事件 事务标签 公开、实质审查、授权
申请人信息
申请人 第一申请人
专利权人 杭州电子科技大学 当前专利权人 杭州电子科技大学
发明人 包建荣、鲁彪、刘超、姜斌、吴俊、唐向宏、邱雨、朱芳 第一发明人 包建荣
地址 浙江省杭州市经济技术开发区白杨街道2号大街1158号 邮编 310018
申请人数量 1 发明人数量 8
申请人所在省 浙江省 申请人所在市 浙江省杭州市
代理人信息
代理机构
专利代理机构是经省专利管理局审核,国家知识产权局批准设立,可以接受委托人的委托,在委托权限范围内以委托人的名义办理专利申请或其他专利事务的服务机构。
浙江千克知识产权代理有限公司 代理人
专利代理师是代理他人进行专利申请和办理其他专利事务,取得一定资格的人。
周希良
摘要
本发明公开了一种基于压缩感知和支持向量机的高效频谱感知方法,包括以下步骤:S1:输入待感知的接收信号;S2:对待感知接收信号通过压缩感知采样处理得到测量矩阵Y,并对测量矩阵Y的采样协方差矩阵对角线元素判决,并取均值得到特征统计量Tk;S3:通过特征统计量与门限对比获取待感知接收信号的标签,并将获取到的标签和特征统计量组成的样本训练集;S4:将样本训练集输入到支持向量机SVM分类器中训练,得到频谱分类器;S5:将收集的数据输入至频谱分类器中处理,得到分类结果。本发明无需重构原信号,只需通过协方差对角线元素作为判决统计量,极大降低复杂度,且提高检测概率。
  • 摘要附图
    一种基于压缩感知和支持向量机的高效频谱感知方法
  • 说明书附图:图1
    一种基于压缩感知和支持向量机的高效频谱感知方法
  • 说明书附图:图2
    一种基于压缩感知和支持向量机的高效频谱感知方法
  • 说明书附图:图3
    一种基于压缩感知和支持向量机的高效频谱感知方法
  • 说明书附图:图4
    一种基于压缩感知和支持向量机的高效频谱感知方法
  • 说明书附图:图5
    一种基于压缩感知和支持向量机的高效频谱感知方法
  • 说明书附图:图6
    一种基于压缩感知和支持向量机的高效频谱感知方法
  • 说明书附图:图7
    一种基于压缩感知和支持向量机的高效频谱感知方法
法律状态
序号 法律状态公告日 法律状态 法律状态信息
1 2022-07-22 授权
2 2021-06-04 实质审查的生效 IPC(主分类): H04B 17/382 专利申请号: 202011618079.6 申请日: 2020.12.30
3 2021-05-18 公开
权利要求
权利要求书是申请文件最核心的部分,是申请人向国家申请保护他的发明创造及划定保护范围的文件。
1.一种基于压缩感知和支持向量机的高效频谱感知方法,其特征在于,包括以下步骤:
S1:输入待感知的接收信号;
S2:对待感知接收信号通过压缩感知采样处理得到测量矩阵Y,并对测量矩阵Y的采样协方差矩阵对角线元素判决,并取均值得到特征统计量Tk;
S3:通过特征统计量与门限对比获取待感知接收信号的标签,并将获取到的标签和特征统计量组成的样本训练集;
S4:将样本训练集输入到支持向量机SVM分类器中训练,得到频谱分类器;
S5:将收集的数据输入至频谱分类器中处理,得到分类结果;
步骤S2中,测量矩阵Y的采样协方差矩阵对角线元素判决并取均值得到特征统计量为通过对采样协方差矩阵对角线元素进行去优化判决去噪处理并对其取均值;其中,采样协方差矩阵R表示为:
T
其中,R为N×N维矩阵;Y的维度为1×N,为经过压缩测量得到的测量矩阵;Y 为Y的转置矩阵,N为测量矩阵的长度;
假设存在K个次用户SU协作频谱感知,得到K个采样协方差矩阵;
H0和H1分别表示主用户PU是否存在的假设条件,表示为:
其中,s(k)和n(k)分别表示第k个接收的主用户PU信号序列和均值为零方差为 的加性高斯白噪声,k=1,2,...N;N表示接收信号总样本数;h(k)表示第k个主用户PU信号序列的信道增益;x(k)表示认知用户SU的接收信号;
取采样协方差矩阵对角线元素rii,与其对应的噪声水平σi比较,去噪处理,如果rii≥σi,表示当前位置信号为正常信号,反之,rii<σi,判定此处为噪声信号,将其对角线元素置为零;其过程表示为:
去噪后的采样协方差矩阵对角线元素均值作为特征统计量,具体表示为:
k
其中,Tk表示第k个感知用户的到的特征统计量,rii为第k个感知用户压缩检测采样协方差矩阵的对角线元素,且k=1,2,…N,N为感知用户的数量。

2.根据权利要求1所述的基于压缩感知和支持向量机的高效频谱感知方法,其特征在于,步骤S2中,对待感知接收信号通过压缩感知采样处理为对待感知接收信号降维处理和主成分提取,具体为:
对信号x=[x1,x2,…xN],其中,N为整数,表示信号长度;xi实数,表示信号x的第i个元素;在N维实数集RN中,存在正交矩阵,为N×N维基矩阵,使得任何信号均由基矩阵线性表示;其中,ψ=[ψ1,ψ2,…,ψN],ψi为N×1维基向量;
同时,信号x还表示为:
x=ψθ                           (1)
或者
其中,θi为信号x和ψi的内积结果,即
θi=<x,ψi>=x1ψi1+x2ψi2+…+xψN iN,θ为信号x在正交基变换矩阵ψ下的投影;若经过变换后的θ中的元素θi仅有K个非零元素时,i=1,2,…N,称信号x为在此正交基变换下K稀疏的稀疏信号,K<若信号x通过变换后满足K稀疏特性,则称信号x在正交变换基矩阵ψ下是能压缩的;信号压缩在满足稀疏特征的条件下,通过M×N维观测矩阵Φ对信号x线性测量,观测矩阵Φ与正交变换基矩阵ψ满足不相干特性,线性测量过程如下:
Y=Φx                           (3)
具体测量过程写为:
Y=Φψx                          (4)
其中,测量矩阵Y为M×1维信号,且M<<N,测量矩阵Y中具有能够表征信号x中的全部信息量;
使得测量矩阵Y能够重构出信号x的充要条件为感知矩阵A=Φψ满足有限等距条件;观测矩阵Φ和正交变换基矩阵Ψ满足不相干性质,不相干性质表示为:
其中,μ(Φ,ψ)为矩阵Φ和Ψ的相关度,且 Φi为观测矩阵的第i行,
1≤i≤M,ψj为正交变换基矩阵的第j列,1≤j≤N;
在认知无线电场景满足稀疏条件下,将接收到信号进行压缩测量,得到测量矩阵Y。

3.根据权利要求1所述的基于压缩感知和支持向量机的高效频谱感知方法,其特征在于,步骤S3中得到特征统计量与门限对比获取待感知接收信号的标签,具体为:
将特征统计量Tk与门限γk比较;若Tk>γk,则将标签f置为+1;若否,则将标签f置为‑1。

4.根据权利要求3所述的基于压缩感知和支持向量机的高效频谱感知方法,其特征在于,步骤S3中将获取得到的标签和得到的特征统计量组成样本训练集,具体为:
比较所有Tk与门限的大小,并计算+1和‑1的数量,如果+1的数量大于‑1的数量,则将特征统计量Tk的标签f置为+1;否则,将特征统计量Tk的标签置f为‑1;将特征统计量Tk和特征统计量Tk对应的标签组成样本训练集G={Tk,f}。

5.根据权利要求4所述的基于压缩感知和支持向量机的高效频谱感知方法,其特征在于,步骤S4中将样本训练集输入到支持向量机SVM分类器中训练得到频谱分类器,具体为:
N
设(T1,f1),(T2,f2)...,(TN,fN),Ti∈R为训练样本数据,fi∈{+1,‑1}为Ti对应的标签;
其中,(Ti,fi)为关联数据的组合,表示接收数据和其对应的标签;N表示样本维度,则最大间隔超平面表示为:
w·T+d=0                       (10)
其中,w为与超平面垂直的向量,即法向量;d表示相对于原点的偏移量;T为特征统计量;
超平面两侧的样本满足以下约束条件:
w·Ti+d≥0(fi=+1)                     (11)
w·Ti+d≤0(fi=‑1)                     (12)
在超平面中增加映射函数φ(T),将T映射到高维空间,得到样本点的决策函数,表示为:
f(T)=sign(w·φ(T)+b)                    (13)
2 2
由超平面决定的分类间隔为2/||w||;||g||表示L2范数,优化目标为分类间隔;
2
当||w||最小时,分类间隔为最大,则优化目标函数为:
分类间隔上的点称为支持向量,为缓解在高维空间的过度拟合,增加一个松弛变量ξ调整间隔,则优化的超平面间隔表示为:
其中,ξi为实数,表示第i个松弛变量;C表示惩罚参数,用于限制ξi;
对优化的超平面间隔方程采用拉格朗日优化处理,表示为:
其中,拉格朗日因子取αi≥0,βi≥0,i=1,2,...,L;
对La(w,b,α,β)中w和d两个变量求偏导数,并令偏导数为零,表示为:
超平面优化方程表示为:
其中,K(Ti,Tj)为核函数,是φ(x)在特征空间的内积,表示为K(Ti,Tj)=<φ(Ti),φ(Tj)>;
采用序列最小最优化方法解决支持向量机SVM分类器中核心优化问题得到拉格朗日因子α和偏移量b,则分类函数即频谱分类器表示为:
其中,sign(﹒)为符号函数。

6.根据权利要求1所述的基于压缩感知和支持向量机的高效频谱感知方法,其特征在于,步骤S5中得到的分类结果为次用户SU信号感知到的频谱占用情况;频谱占用情况为:输入测试样本,若频谱分类器输出的结果为+1,则主用户PU信号频谱被占用;若频谱分类器输出的结果为‑1,则主用户PU信号频谱未被占用。
说明书

技术领域

[0001] 本发明涉及数字通信技术领域,尤其涉及一种基于压缩感知和支持向量机的高效频谱感知方法。

背景技术

[0002] 近年来,随着无线电通信技术的飞速发展,无线电频谱资源的需求迅速攀升,另外,由于无线电频谱资源主要采用固定分配模式,使得频谱利用率低、频谱资源浪费的问题十分严重。图1为典型的认知无线网络(CRN)系统架构示意图。典型的认知无线电网络(CNR)由主用户(PU)和次用户(SU)组成,通常假设PU和SU的无线网络通信在物理上是分隔的,SU无法直接获得PU的信道状态。在此系统中,PU对所占信道具有优先使用权,认知基站(CBS)首先通过检测信道中的PU信号,判断频谱中的空闲信道,然后,它发送PU接收机(PU‑R)的状态并确定空闲频谱。直到PU不再占用频谱,SU可以重新使用该频谱。若SU正使用的频谱被PU访问,则SU要退出该频谱转入缓存中,认知设备同时检测别的空闲频谱。
[0003] 关于认知无线电中的频谱感知技术的研究,已经有大量学者做了许多工作,针对不同的应用场景,提出了多种类型的频谱感知算法。传统的频谱感知算法主要针对单用户的窄带信号检测,包括能量检测算法,匹配滤波算法,循环平稳特征检测算法。能量检测算法是最常用的频谱感知方法。能量检测算法简单,所需时间短,但是会受到噪声峰值影响,不适合低信噪比的环境。匹配滤波检测算法精度高,所需时间短,然而需要授权用户的先验信息。循环平稳特征检查测算法检测精度高,支持在低信噪比环境下检测;但其算法复杂度高,检测时间最长。一般来说,频谱感知方法需要信号和噪声的先验信息,这限制了实际应用。近年来,随机矩阵理论的兴起,越来越多的基于协方差矩阵的盲频谱检测算法应用于频谱感知。传统盲检测方法包括:协方差绝对值(CAV)检测、最大最小特征值法(MME)检测等。这些算法可以消除一定的噪声不确定度影响,提高检测概率。但是MME和CAV的门限均是渐进的,门限的不精确性会影响算法检测性能
[0004] 认知无线电(CR)技术的出现,使次用户(SU)能智能接入未被占用的空闲频谱,大大的提高的频谱利用率。其中,频谱感知作为CR关键,能够精确智能的识别出空闲频谱,充分利用频谱资源,有效提高频谱利用率。其中,传统单用户频谱感知技术因其局限性不适合实际复杂环境。为此,多用户协作频谱感知技术应运而生。它通过对多用户感知结果融合,能有效提高频谱感的性能。其中,基于压缩感知技术的频谱感知在检测性能方面有着显著的效果,但传统压缩频谱感知技术需要对原始信号重构,大大提高了频谱感知的复杂度。而利用支持向量机可以省去信号重构步骤,直接对压缩测量后的矩阵处理,有效降低了算法复杂度。
[0005] 因此,本发明基于现有技术当中的受限等距特性(RIP)、序列最小最优化方法以及支持向量机,提出了一种基于压缩感知和支持向量机的高效频谱感知方法,其中,[0006] 受限等距特性(RIP):
[0007] 受限等距特性是观测矩阵Φ能够精确重构原始信号的充要条件,指对任意信号X,观测矩阵Φ均可满足下式:
[0008]
[0009] 其中,K为信号稀疏度,δ为足够小的一个数值,观测矩阵Φ的维度为M×N。
[0010] 序列最小最优化方法:
[0011] 序列最小优化(SMO)是一种用于解决支持向量机训练过程中系数优化问题的迭代算法,SMO算法的中心思想就是每次选出两个α优化,然后固定其他的α值。重复此过程,直到达到某个终止条件程序退出并得到我们需要的优化结果。
[0012] 支持向量机:
[0013] 支持向量机SVM是一类按监督学习方式对数据二分类的广义线性分类器。其通过二次规划来解决模式识别及分类问题。SVM分类器的最主要思想是找到线性可分离特征空间的最优超平面,也就是要使正确分类的概率达到最大,并达到最大超平面间隔。最大超平面方程表达式为(w·x)+b=0,其中,w为超平面的法向量,b是相对于原点的偏移量。在超平面两侧的点表示被分类的两种情况,若以y表示分类类别,若yi=+1,则有(w·xi)+b≥0;若2
yi=‑1,则有(w·xi)+b≤0。由超平面决定的分类间隔为2/||w|| ,分类间隔上的点成为支
2
持向量,最大超平面即分类间隔最大,也就是满足上述条件的情况下使得2/||w||最大,即发明内容
[0014] 本发明的目的是针对现有压缩频谱感知技术的缺陷,提供了一种基于压缩感知和支持向量机的高效频谱感知方法,在低信噪比的条件下依然能有较高的频谱识别率,无需信号重构,有效降低运算复杂度,同时,对支持向量机判决门限的优化,使得频谱感知更加准确。
[0015] 为了实现以上目的,本发明采用以下技术方案:
[0016] 一种基于压缩感知和支持向量机的高效频谱感知方法,包括以下步骤:
[0017] S1:输入待感知的接收信号;
[0018] S2:对待感知接收信号通过压缩感知采样处理得到测量矩阵Y,并对测量矩阵Y的采样协方差矩阵对角线元素判决并取均值得到特征统计量;
[0019] S3:通过特征统计量与门限对比获取待感知接收信号的标签,并将获取到的标签和得 到的 特 征统 计 量组 成样 本 训练 集 ;门限 为已 知 给 定的 ;门 限 为‑1M为采样数量, 为噪声方差。其中,Q (·)为Q函数的
互逆函数并且有 Pf为虚警概率;
[0020] S4:将组成的样本训练集输入到支持向量机SVM分类器中训练,得到频谱分类器;
[0021] S5:将收集的数据输入至频谱分类器中处理,得到分类结果。
[0022] 进一步的,所述步骤S2中对待感知接收信号通过压缩感知采样处理是对待感知接收信号降维处理和主成分提取。所述对待感知接收信号降维处理和特征提取具体为:
[0023] 对信号x=[x1,x2,…xN]。其中,N为整数,表示离散信号长度。xi实数,表示信号x的第i个元素。在N维实数集RN中,存在正交矩阵,为N×N维基矩阵,使得任何信号均可由基矩阵线性表示。其中,ψ=[ψ1,ψ2,…,ψN],ψi为N×1维基向量,同时,信号x还可表示为:
[0024] x=ψθ                 (1)
[0025] 或者
[0026]
[0027] 其中,θi为信号x和ψi的内积结果,即θi=<x,ψi>=x1ψi1+x2ψi2+…+xNψiN,θ为信号x在正交基变换矩阵ψ下的投影。若经过变换后的θ中的元素θi(i=1,2,…N)仅有K(K<
[0028] 若信号x通过变换后满足K稀疏特性,则称信号x在正交变换基矩阵ψ下是可压缩的。在实际的信号处理中我们要将信号压缩,即可在满足稀疏特征的条件下,通过M×N维观测矩阵Φ对信号x线性测量,观测矩阵Φ与正交基变换矩阵ψ满足不相干特性,线性测量过程如下:
[0029] Y=Φx                  (3)
[0030] 具体测量过程可写为:
[0031] Y=Φψx                 (4)
[0032] 其中,测量矩阵Y为M×1维信号,且M<<N,测量矩阵Y中具有能够表征信号x中的全部信息量,故可以用上述压缩测量来实现由N维数据到到M维数据的压缩采样,且不丢失信号x的信息。
[0033] 另外,使得测量矩阵能够重构出信号x的充要条件为感知矩阵A=Φψ满足有限等距条件(RIP)。研究表明,观测矩阵Φ和正交基变换矩阵Ψ满足不相干性质,则感知矩阵A=Φψ大概率满足背景技术所述的RIP性质,故都可作为观测矩阵的衡量标准。不相干性质可表示为:
[0034]
[0035] 其中,μ(Φ,ψ)为矩阵Φ和Ψ的相关度,且 Φi(1≤i≤M)为观测矩阵的第i行,ψj(1≤j≤N)为正交基变换矩阵的第j列。
[0036] 实际上,很多随机矩阵都可与变换矩阵较好的满足RIP性质,如高斯矩阵、伯努利矩阵、傅里叶变换矩阵等。而在认知无线电场景下,频谱上存在频谱空洞,需要次用户准确的判断频谱空洞并智能接入,此场景满足稀疏条件,故可将接收到信号上述压缩观测,得到测量矩阵Y。
[0037] 进一步的,所述步骤S2中测量矩阵Y的采样协方差矩阵对角线元素判决并取均值得到特征统计量为通过对采样协方差矩阵对角线元素进行去优化判决去噪处理并对其取均值,具体为:
[0038] 采样协方差矩阵表示为:
[0039]
[0040] 其中,R为N×N维矩阵,Y的维度为1×N,为经过压缩测量得到的测量矩阵,YT为Y的转置矩阵,N为测量矩阵的长度。
[0041] 假设存在K个次用户(SU)协作频谱感知,故在此得到K个采样协方差矩阵。
[0042] H0和H1分别表示主用户PU是否存在的假设条件,表示为:
[0043]
[0044] 其中,s(k)和n(k)(k=1,2,...N)分别表示第k个接收的主用户PU信号序列和均值为零方差为 的加性高斯白噪声;N表示接收信号总样本数;h(k)表示第k个主用户PU信号序列的信道增益;x(k)表示认知用户SU的接收信号。
[0045] 对压缩观测信号的采样协方差矩阵而言,对角线元素为各个变量的方差,如果频谱检测过程中噪声过大,会对对角线元素方差值产生较大影响,从而对判决结果产生很大影响,因此在本步骤中首先对采样协方差矩阵对角线元素判决优化。
[0046] 取采样协方差矩阵对角线元素rii,与其对应的噪声水平σi比较,去噪处理,如果rii≥σi,表示当前位置信号为正常信号,反之,rii<σi,判定此处为噪声信号,将其对角线元素置为零。其过程可表示为:
[0047]
[0048] 去噪后的采样协方差矩阵对角线元素均值作为特征统计量,具体表示为:
[0049]
[0050] 其中,Tk表示第k个感知用户的到的特征统计量,riik为第k个感知用户压缩检测采样协方差矩阵的对角线元素,且k=1,2,…=N,N为感知用户的数量。
[0051] 进一步地,所述步骤S3中得到特征统计量与门限对比获取待感知接收信号的标‑1签,具体为将判决统计量Tk与门限 比较。其中,Q (·)
为Q函数的互逆函数并且有 Pf为虚警概率。若Tk>γk,则
将标签f置为+1;若否,则将标签f置为‑1。
[0052] 进一步的,所述步骤S3中将所得标签和特征统计量组成样本训练集,具体为:
[0053] 比较所有Tk与门限的大小,并计算+1和‑1的数量,如果+1的数量大于‑1的数量,则将特征统计量Tk的标签f置为+1;否则,将特征统计量Tk的标签置f为‑1。将特征统计量Tk和特征统计量Tk对应的标签组成样本训练集G={Tk,f}。
[0054] 进一步的,所述步骤S4中将样本训练集输入到支持向量机SVM分类器中训练得到频谱分类器,具体为:
[0055] 设(T1,f1),(T2,f2)...,(TN,fN),Ti∈RN为训练样本数据,fi∈{+1,‑1}为Ti对应的标签;其中,(Ti,fi)为关联数据的组合,表示接收数据和其对应的标签;N表示样本维度,则最大间隔超平面表示为:
[0056] w·T+d=0                       (10)
[0057] 其中,w为与超平面垂直的向量,即法向量;d表示相对于原点的偏移量;T为特征统计量。
[0058] 超平面两侧的样本满足以下约束条件:
[0059] w·Ti+d≥0(fi=+1)                     (11)
[0060] w·Ti+d≤0(fi=‑1)                     (12)
[0061] 在超平面中增加映射函数φ(T),将T映射到高维空间,得到样本点的决策函数,表示为:
[0062] f(T)=sign(w·φ(T)+b)                    (13)
[0063] 由超平面决定的分类间隔为2/||w||2;||·||2表示L2范数,优化目标为分类间隔;
[0064] 当||w||2最小时,分类间隔为最大,则优化目标函数为:
[0065]
[0066] 分类间隔上的点称为支持向量,为缓解在高维空间的过度拟合,实际中通常增加一个松弛变量ξ调整间隔,则优化的超平面间隔表示为:
[0067]
[0068] 其中,ξi为实数,表示第i个松弛变量;C表示惩罚参数,用于限制ξi;
[0069] 对优化的超平面间隔方程采用拉格朗日优化处理,表示为:
[0070]
[0071] 其中,拉格朗日因子取αi≥0,βi≥0,i=1,2,...,L;
[0072] 对La(w,b,α,β)中w和d两个变量求偏导数,并令偏导数为零,表示为:
[0073]
[0074] 超平面优化方程表示为:
[0075]
[0076] 其中,K(Ti,Tj)为核函数,是φ(x)在特征空间的内积,表示为K(Ti,Tj)=<φ(Ti),φ(Tj)>;
[0077] 采用背景技术所述的序列最小最优化方法解决支持向量机SVM分类器中核心优化问题得到拉格朗日因子α和偏移量b,则分类函数即频谱分类器表示为:
[0078]
[0079] 其中,sign(﹒)为符号函数。
[0080] 进一步的,所述步骤S5中得到的分类结果为次用户SU信号感知到的频谱占用情况;所述频谱占用情况为:输入测试样本,若频谱分类器输出的结果为+1,则主用户PU信号频谱被占用;若频谱分类器输出的结果为‑1,则主用户PU信号频谱未被占用。
[0081] 与现有技术相比,本发明利用压缩感知技术对感知信号压缩测量,有效地降低数据维数,且基于支持向量机技术提出一种非重构的压缩频谱感知算法,并且对SVM判决统计量优化,经SVM训练寻找最优超平面,使得PU和噪声间隔达到最大。在低信噪比下具有较高的检测率,有效的提高频谱利用率,具有较高应用价值。

实施方案

[0089] 以下通过特定的具体实例说明本发明的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本发明的精神下各种修饰或改变。需说明的是,在不冲突的情况下,以下实施例及实施例中的特征可以相互组合。
[0090] 本实施例提供了一种基于支持向量机的高效频谱感知方法,如图2所示,包括以下步骤:
[0091] S1:输入待感知的接收信号;
[0092] S2:对待感知接收信号通过压缩感知采样处理得到测量矩阵Y,并对测量矩阵Y的采样协方差矩阵对角线元素判决并取均值得到特征统计量;
[0093] S3:通过特征统计量与门限对比获取待感知接收信号的标签,并将获取到的标签和得到的特征统计量组成样本训练集;
[0094] S4:将组成的样本训练集输入到支持向量机SVM分类器中训练,得到频谱分类器;
[0095] S5:将收集的数据输入至频谱分类器中处理,得到分类结果。
[0096] 本实施例具体采用以下技术方案:提取采样协方差矩阵对角线元素平均值作为特征,将其与已知门限比较生成标签。压缩观测对感知信号降维处理,对采样协方差矩阵对角线元素优化判决得到特征统计量,最后将统计量和标签组成训练样本输入到SVM分类器中训练得到频谱分类器。测试数据输入到频谱分类器中,输出“+1”表示PU信号存在,输出“‑1”表示PU信号不存在。首先,传统压缩频谱感知算法需要对原信号重构,算法整体复杂度很高,本实施例将支持向量机技术引入压缩频谱感知算法中,对压缩观测后的信号直接SVM分解,有效降低运算复杂度。另外,本实施例频谱感知的结果只有频谱是否被占用两种情况,满足SVM分类的性质。其次,本实施例通过对采样协方差矩阵对角线元素优化判决,有效提高频谱检测效率。因此,本实施例相比于传统SVM频谱感知算法,极大的减少样本维度,有效减少复杂度,提高频谱感知效率。
[0097] 步骤S2中,对待感知接收信号通过压缩感知采样处理是对待感知接收信号降维处理和主成分提取,以及测量矩阵Y的采样协方差矩阵对角线元素判决并取均值得到特征统计量的流程如图3所示。
[0098] 对待感知接收信号通过压缩感知采样处理是对待感知接收信号降维处理和主成分提取,具体为:
[0099] 对信号x=[x1,x2,…xN]。其中,N为整数,表示离散信号长度。xi实数,表示信号x的第i个元素。在N维实数集RN中,存在正交矩阵,为N×N维基矩阵,使得任何信号均可由基矩阵线性表示。其中,ψ=[ψ1,ψ2,…,ψN],ψi为N×1维基向量,同时,信号x还可表示为:
[0100] x=ψθ                 (20)
[0101] 或者
[0102]
[0103] 其中,θi为信号x和ψi的内积结果,即θi=<x,ψi>=x1ψi1+x2ψi2+…+xNψiN,θ为时域信号x在正交基变换矩阵ψ下的投影。若经过变换后的θ中的元素θi(i=1,2,…N)仅有K(K<
[0104] 若信号x通过变换后满足K稀疏特性,则称信号x在正交基变换矩阵ψ下是可压缩的。在实际的信号处理中我们要将信号压缩,即可在满足稀疏特征的条件下,通过M×N维观测矩阵Φ对信号x线性测量,观测矩阵Φ与正交基变换矩阵ψ满足不相干特性,线性测量过程如下:
[0105] Y=Φx                  (22)
[0106] 具体测量过程可写为:
[0107] Y=Φψx               (23)
[0108] 其中,测量矩阵Y为M×1维信号,且M<<N,测量矩阵Y中具有能够表征信号x中的全部信息量,故可以用上述压缩观测来实现由N维数据到到M维数据的压缩采样,且不丢失信号x的信息。
[0109] 另外,使得测量矩阵Y能够重构出信号x的充要条件为感知矩阵A=Φψ满足有限等距条件(RIP)。研究表明,观测矩阵Φ和正交基变换矩阵Ψ满足不相干性质,则感知矩阵A=Φψ大概率满足背景技术所述的RIP性质,故都可作为观测矩阵的衡量标准。不相干性质可表示为:
[0110]
[0111] 其中,μ(Φ,ψ)为矩阵Φ和Ψ的相关度,且 Φi(1≤i≤M)为观测矩阵的第i行,ψj(1≤j≤N)为变换基矩阵的第j列。
[0112] 实际上,很多随机矩阵都可与变换矩阵较好的满足RIP性质,如高斯矩阵、伯努利矩阵、傅里叶变换矩阵等。而在认知无线电场景下,频谱上存在频谱空洞,需要次用户准确的判断频谱空洞并智能接入,此场景满足稀疏条件,故可将接收到信号上述压缩观测,得到测量矩阵Y。
[0113] 测量矩阵Y的采样协方差矩阵对角线元素判决并取均值得到特征统计量为通过对采样协方差矩阵对角线元素进行去优化判决去噪处理并对其取均值;
[0114] 其中,采样协方差矩阵R表示为:
[0115]
[0116] 其中,R为N×N维矩阵,Y的维度为1×N,为经过压缩测量得到的测量矩阵,YT为Y的转置矩阵,N为测量矩阵的长度。
[0117] 假设存在K个次用户(SU)协作频谱感知,故在此得到K个采样协方差矩阵。
[0118] H0和H1分别表示主用户PU是否存在的假设条件,表示为:
[0119]
[0120] 其中,s(k)和n(k)(k=1,2,...N)分别表示第k个接收的主用户PU信号序列和均值为零方差为 的加性高斯白噪声;N表示接收信号总样本数;h(k)表示第k个主用户PU信号序列的信道增益;x(k)表示认知用户SU的接收信号。
[0121] 对压缩观测信号的采样协方差矩阵而言,对角线元素为各个变量的方差,如果频谱检测过程中噪声过大,会对对角线元素方差值产生较大影响,从而对判决结果产生很大影响,因此在本步骤中首先对采样协方差矩阵对角线元素判决优化。
[0122] 取采样协方差矩阵对角线元素rii,与其对应的噪声水平σi比较,去噪处理,如果rii≥σi,表示当前位置信号为正常信号,反之,rii<σi,判定此处为噪声信号,将其对角线元素置为零。其过程可表示为:
[0123]
[0124] 去噪后的采样协方差矩阵对角线元素均值作为特征统计量,具体表示为:
[0125]
[0126] 其中,Tk表示第k个感知用户的到的特征统计量,riik为第k个感知用户压缩检测采样协方差矩阵的对角线元素,且k=1,2,…=N,N为感知用户的数量。
[0127] 步骤S3中得到特征统计量与门限对比获取待感知接收信号的标签的流程图如图4所示,具体为:
[0128] 判决统计量Tk与门限 比较。其中,M为采样数‑1
量, 为噪声方差。Q (·)为Q函数的互逆函数并且有 Pf
为虚警概率。若Tk>γk,则将标签f置为+1;若否,则将标签f置为‑1。比较所有Tk与门限的大小,并计算+1和‑1的数量,如果+1的数量大于‑1的数量,则将特征统计量Tk的标签f置为+1;
否则,将特征统计量Tk的标签置f为‑1。将特征统计量Tk和特征统计量Tk对应的标签组成样本训练集G={Tk,f}。
[0129] 骤S4中,将样本训练集输入到支持向量机SVM分类器中训练得到频谱分类器的流N程图如图5所示。设(T1,f1),(T2,f2)...,(TN,fN),Ti∈R为训练样本数据,fi∈{+1,‑1}为Ti对应的标签;其中,(Ti,fi)为关联数据的组合,表示接收数据和其对应的标签;N表示样本维度,则最大间隔超平面表示为:
[0130] w·T+d=0                 (29)
[0131] 其中,w为与超平面垂直的向量,即法向量,即;d表示相对于原点的偏移量;T为特征统计量。
[0132] 超平面两侧的样本满足以下约束条件:
[0133] w·Ti+d≥0(fi=+1)              (30)
[0134] w·Ti+d≤0(fi=‑1)             (31)
[0135] 在超平面中增加映射函数φ(T),将T映射到高维空间,得到样本点的决策函数,表示为:
[0136] f(T)=sign(w·φ(T)+b)            (32)
[0137] 由超平面决定的分类间隔为2/||w||2;||·||2表示L2范数,优化目标为分类间隔;
[0138] 当||w||2最小时,分类间隔为最大,则优化目标函数为:
[0139]
[0140] 分类间隔上的点称为支持向量,为缓解在高维空间的过度拟合,实际中通常增加一个松弛变量ξ调整间隔,则优化的超平面间隔表示为:
[0141]
[0142] 其中,ξi为实数,表示第i个松弛变量;C表示惩罚参数,用于限制ξi;
[0143] 对优化的超平面间隔方程采用拉格朗日优化处理,表示为:
[0144]
[0145] 其中,拉格朗日因子取αi≥0,βi≥0,i=1,2,...,L;
[0146] 对La(w,b,α,β)中w和d两个变量求偏导数,并令偏导数为零,表示为:
[0147]
[0148] 超平面优化方程表示为:
[0149]
[0150] 其中,K(Ti,Tj)为核函数,是φ(x)在特征空间的内积,表示为K(Ti,Tj)=<φ(Ti),φ(Tj)>;
[0151] 采用背景技术所述的序列最小最优化方法解决支持向量机SVM分类器中核心优化问题得到拉格朗日因子α和偏移量b,则分类函数即频谱分类器表示为:
[0152]
[0153] 其中,sign(﹒)为符号函数。
[0154] 步骤S5中,得到的分类结果为次用户SU信号感知到的频谱占用情况流程如图6所示。将认知用户收集到的数据先压缩观测处理,数据压缩降低数据维度后对压缩感测数据采样协方差矩阵对角线元素优化判决得到特征统计量T,将统计量输入到SVM分类器中,得到频谱分类器,将测试数据输入到训练好的频谱分类器中,得到PU信号的频谱占用情况。若频谱分类器输出为+1,则频谱被占用;若输出为‑1,则表示频谱未被占用。
[0155] 如图7所示为SVM线性最大间隔超平面示意图。SVM的关键是通过最大化分类器间隔余量,使误差之和最小化,从而使泛化能力最大化。但在低信噪比环境,其无法通过线性超平面去实现频谱感知。通常使用核函数(即 表示的非线性映射函数)将输入低维空间映射到特征空间的高维点积空间,则可改善低维空间线性不可分离性。
[0156] 综上所述,本实施例提出一种低信噪比下基于压缩感知和支持向量机高(SVM)的高效频谱感知,其包括了压缩观测、采样协方差矩阵对角线元素平均值优化统计量和SVM训练和测试过程。针对传统压缩频谱感知需要重构原始信号从而计算量大的问题,本方法无需重构原始信号,对压缩观后的信号直接SVM分类处理。在低信噪比的情况下,实现高效快速的感知频谱,提高频谱利用率,具有较高应用价值。本实施例利用压感知对输入信号压缩观测处理,有效的降低了信号维度,减小后续训练测试复杂度,且利用采样协方差矩阵对角线元素来构造特征统计量,减提高检测效率,最后通过SVM算法将噪声和PU信号分离,有效提高频谱检测概率,具有较高的应用价值。
[0157] 尽管已清晰描述了本发明的实施例,但对本领域的技术人员而言,可在不脱离本发明方法原理和思想的情况下,对这些实施例开展多种变化、修改、替换和变型,则本发明的范围由所附权利要求及其等同限定。即可根据实际情况改变本发明方法所述方法中认知无线电的主信号PU带宽、次用户PU数量、采样频率、采样时间等模型参数,观测矩阵的选取、训练信号和测试信号的产生方式,测试样本的维数,采样协方差矩阵的构造方式,核函数的参数以及SVM中的惩罚因子等参数。仍属于本发明所述方法的范畴,仍受本发明保护。
[0158] 注意,上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解,本发明不限于这里所述的特定实施例,对本领域技术人员来说能够各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此,虽然通过以上实施例对本发明了较为详细的说明,但是本发明不仅仅限于以上实施例,在不脱离本发明构思的情况下,还可以包括更多其他等效实施例,而本发明的范围由所附的权利要求范围决定。

附图说明

[0082] 图1是典型的认知无线网络系统架构示意图;
[0083] 图2是一种基于压缩感知和支持向量机的高效频谱感知方法流程图;
[0084] 图3是压缩观测和采样协方差矩阵对角线元素优化判决构造统计量流程图;
[0085] 图4是标签的生成与统计量构成训练集流程图;
[0086] 图5是训练频谱分类器流程图;
[0087] 图6是测试频谱分类器流程图;
[0088] 图7是SVM最大线性间隔超平面示意图。
版权所有:盲专网 ©2023 zlpt.xyz  蜀ICP备2023003576号