CN107517089B - 一种基于频谱特征的宽带频谱感知方法 - Google Patents
一种基于频谱特征的宽带频谱感知方法 Download PDFInfo
- Publication number
- CN107517089B CN107517089B CN201710795476.2A CN201710795476A CN107517089B CN 107517089 B CN107517089 B CN 107517089B CN 201710795476 A CN201710795476 A CN 201710795476A CN 107517089 B CN107517089 B CN 107517089B
- Authority
- CN
- China
- Prior art keywords
- spectrum
- signal
- sampling
- matrix
- samples
- 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
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
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/391—Modelling the propagation channel
- H04B17/3911—Fading models or fading generators
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
本发明公开了一种基于频谱特征的宽带频谱感知方法,属于认知无线电的频谱感知领域。该方法利用信号频谱特征检测来检测主用户占用频段并估计功率水平的一种方法。根据实际中频谱资源利用率较低的特点,对可能存在的频率位置和相应的功率值水平进行估计,该方法在原先无衰落信道模型基础上加入瑞利衰落信道的信道模型,以求出本方法在瑞利衰落信道,利用基于频谱特征的技术,以加权正交匹配追踪(Weighted Orthogonal Matching Pursuit,WOMP)算法为信号重建算法估计信号功率值水平和频谱占用情况。
Description
技术领域
本发明涉及一种基于频谱特征的宽带频谱感知方法,属于认知无线电(CognitiveRadio,CR)的频谱感知(Spectrum Sensing,SS)技术领域。
背景技术
由于被分配的频谱没有得到充分利用,所以现行的基于特定地理区域内独占特定授权频带的规定是低效的。为了提高当前稀缺的频谱资源的利用效率,频谱管理制度向更加智能和灵活的方向迈进。在这一背景下,认知无线电技术得以推广。通过使未授权系统有机会利用未使用的许可频带,CR解决了频谱资源有限的问题。在CR的各项技术中,频谱感知技术可用于检测某频段上许可用户是否存在。频谱感知技术包括传统感知技术和压缩感知技术(Compressive sensing,CS)。传统的信号获取与处理过程包括:采样、压缩、传输、解压缩四个部分,其采样过程必须遵循奈奎斯特采样定理。但是这种方式采样数据量大,先采样后压缩,浪费了大量的传感元、时间和存储空间。而压缩感知技术则是利用信号在某一域上的稀疏性,将数据采集和压缩合二为一,从而允许采样率明显低于奈奎斯特速率,减少传感元数量,降低采样速率,减少传输时延。
目前的频谱感知技术有:基于匹配滤波器的检测、能量检测、特征检测、基于自相关的检测等技术。检测性能方面比较优秀的技术是基于信号自相关特性的检测技术,自相关特征中包含可用的信号的频谱特征,能够实现显著超越盲目技术的检测性能。然而在实际的CS应用中,衰落信道带来的影响使得干扰具有非规律性,不能被考虑进去。而本发明能够很好地解决上面的问题。
发明内容
本发明目的在于针对上述现有技术的不足,提出了一种基于频谱特征的宽带频谱感知方法,该方法引入瑞利衰落信道,考虑它带来的小尺度衰落的影响。该方法是将经过衰落信道后的接收信号的自相关函数,与主用户的自相关函数比较,通过选择合适的过完备字典,保证残差的半正定特性,最终直接估计主用户的频谱占用情况,同时检测出对应的功率水平。
本发明解决其技术问题所采取的技术方案是:一种基于频谱特征的宽带频谱感知方法,该方法采用瑞利衰落信道模型,将主用户的频谱特征当做唯一的先验条件,通过比较接收信号的功率谱密度与主用户的频谱信息,求出信号功率值和频谱占用情况。包括如下:
首先建立基于瑞利衰落信道的信号传输模型;其次根据主用户和接收信号的自相关函数建立函数关系;然后对信号进行压缩采样;再利用基于频谱特征的压缩感知技术,采用相关匹配策略,通过保证残差自相关函数的半正定特性,利用l1-最小化公式推导出算法中的权值表示;最后利用加权正交匹配追踪(Weighted Orthogonal Matching Pursuit,WOMP)算法识别主用户占用的频谱,估计对应的功率值水平。具体包括:
1)建立基于瑞利衰落信道的信号传输模型:
无线信号经过单径瑞利衰落及加性高斯白噪声(additive white gaussiannoise,AWGN)的干扰,到达接收机,设置如下参数:x(t)是发送信号;α(t)是无线信道的衰落因子;为导致的传输损耗;ζ(t)为瑞利衰落随机过程;n(t)为加性高斯白噪声;i(t)为干扰。
2)接收信号经过多陪集采样后自相关函数模型:
参数为:设表示主用户可能的频率位置的集合,并且令表示能包含的网格。γi(τ)和γn(τ)分别表示干扰和噪声的自相关函数,p(ωm)是对应于主要用户的的发射功率。γcm(τ,ωm)表示参考信号xc(t)在载波频率ωm处的参考自相关函数。实际情况下,因为许多频带没有被占用,所以频谱具有稀疏性。本发明定义稀疏向量p为检测函数的输出,零值直接表明在ω处不存在主用户,非零值直接表明在ω处主用户存在。
本发明对通过瑞利信道后的信号y(t)进行压缩采样,采用的是多陪集(Multi-coset,MC)采样:给定接收的多频带信号y(t),在一定时刻ti(n)得到MC采样样本,其中L>0是适当的整数,i=1,2...,k,n∈Z。集合{ci}包含从{0,1...,L-1}中选择的k个不同的整数。注意到MC采样过程可以被等效为经典的以块的形式进行奈奎斯特采样,从每块的L个样本中周期性保留k个样本。未丢弃的样本由集合{ci}指定。每个序列的周期等于LT。
本发明将zf看做y(t)中均匀奈奎斯特采样样本的第f块,每块长度为L。多陪集采样后的总的样本数据可用Y表示,总块数为Nf,为了将采集的样本与原始奈奎斯特采样信号相关联,令yf表示每块的采样样本,每块有k个非均匀样本,Φ∈Ck×L是一个矩阵,它的作用是从zf中非均匀选择k个样本。本发明中Φ可以通过对单位矩阵IL非均匀抽取k行实现。
3)基于相关匹配的压缩频谱感知:
令{Ry(k)}表示奈奎斯特采样信号{y(n)}的相关性,Ry是维度为L×L的相关矩阵。本发明将{xc(n)}定义为奈奎斯特采样候选信号,{Rc(k)}定义为{xc(n)}的相关函数。相关矩阵表示采样样本相关性的平均值。基带参考自相关是从主用户的频谱特征中抽取出来的。将纯音的自相关函数用e(ωm)eH(ωm)给出,e(ωm)∈CL×1是频率转向向量,利用e(ωm)在感知频率ωm处形成的秩为1的矩阵对Rc调制。调制参考相关性Rcm(ωm)必须也以对信号同样的方法进行压缩,比如利用相同的采Ry样矩阵Φ。
本发明将自相关函数模型重写成稀疏表示,是k2×1维的向量,由的列级联形成。列的级联由运算符vec()表示。B包含主信号的频谱信息,并被定义为diag(rc),其中rc是vec(Rc)。矩阵S定义了扫频网格。过完备字典的稀疏表示由矩阵BS给出。
4)重构算法中权值的选择包括:
在根据残差矩阵相关匹配算法选择的方法中,必须满足以下条件:
1)p必须是稀疏的;
3)对稀疏向量施加稀疏约束。
有益效果:
1、本发明不需完整重建接收信号的功率谱密度,可以直接检测主用户占用哪段频率,还可以识别主要用户的发射功率,具有较强的抗干扰能力。
2、本发明可以用于在瑞利衰落信道中检测主用户所占频段,具有较好的检测性能。
3、本发明经过瑞利衰落信道后,可以较准确地恢复信号特征,对干扰具有很好地鲁棒性。
4、本发明引入瑞利衰落信道模型,将基于频谱特征的频谱感知技术与小尺度衰落联系起来。
5、本发明将衰落信道对主用户频谱特征的影响加到相关匹配中,估计频谱位置,检测信号功率水平。
附图说明
图1为系统模型示意图。
图2为本发明实施例中的公式15形象表示示意图。
具体实施方式
下面结合说明书附图对本发明创造作进一步的详细说明。
如图1所示,本发明提供了一种宽带频谱感知中基于频谱特征检测的方法,该方法包括如下:
1)系统模型
无线信号经过单径瑞利衰落信道后,加入了加性高斯白噪声和干扰后信号为x(t)。
y(t)=α(t)x(t)+i(t)+n(t) 式(1)
x(t)是发送信号;n(t)为加性高斯白噪声;i(t)为干扰;α(t)是无线信道的衰落因子,可用式(2)表示:
其中γε(τ)表示噪声和干扰和的自相关函数,p(ωm)是对应于主用户的的发射功率。γcm(τ,ωm)表示参考信号xc(t)在载波频率ωm处的参考自相关函数。本发明定义稀疏向量p,非零值表明在ω处主用户存在。p由式(4)表示:
p=[p(ω0)p(ω1)…p(ωM-1)]T 式(4)
2)多陪集采样
给定接收的多频带信号y(t),在特定时刻ti(n)得到MC采样样本,ti(n)用式(5)表示:
ti(n)=(nL+ci)T 式(5)
其中L>0是适当的整数,i=1,2...,k,n∈Z。集合{ci}包含从{0,1...,L-1}中选择的k个不同的整数。每个序列的周期等于LT。注意到MC采样过程可被视为经典的以块进行奈奎斯特采样,从每块的L个样本中周期性舍弃(L-k)个样本。本发明将zf看做y(t)中均匀奈奎斯特采样样本的第f块,每块长度为L,其中由式(6)表示。
多陪集采样后的总的样本数据可用Y表示,式(8)中yf表示每块的采样样本:
根据式(6)和式(8)可以将yf和zf联系起来,用式(9)表示,其中Φ∈Ck×L是一个矩阵,它的作用是从zf中非均匀选择k个样本。Φ可以通过对单位矩阵IL非均匀抽取k行实现。
yf=Φzf 式(9)
3)基于相关匹配的频谱感知
将纯音的自相关函数用e(ωm)eH(ωm)给出,其中e(ωm)∈CL×1是频率转向向量,定义为
本发明的目的在于估计在中包含某些频率的功率有多大。因此本发明设计一种频谱分析仪,它的脉冲响应不是单独频率,而是预先设定的频谱形状。为了获得特定频谱形状的频率位置,用e(ωm)eH(ωm)形成的秩为1的矩阵对基带参考自相关调制,得到式(12)。
Rcm(ωm)=[Rc⊙e(ωm)eH(ωm)] 式(12)
⊙表示两个矩阵的对应元素乘积。调制参考(自)相关性Rcm(ωm)也必须以对信号同样的方法进行压缩,进行如下操作:
根据式(3)和式(13)中定义的数据自相关矩阵,可以得到相应的模型为:
式(14)中的模型可以改写成如下的稀疏表示,
式(15)中kron(Φ,Φ)表示Kronecker积。定义B是diag(rc),它包含主信号的频谱信息变量,其中rc可用式(17)表示。
rc=vec(Rc) 式(17)
同样的,ε可用式(18)来表示。
此外,式(15)中矩阵S定义了扫频网格,用式(19)表示。
S=[s(ω0)s(ω1)…s(ωM-1)] 式(19)
s(ωm)=vec(e(ωm)eH(ωm)) 式(20)
式(20)为S中每个元素表示形式。
式(15)给出的过完备字典的稀疏表示由矩阵BS给出,它的列包含式(3)中引入的不同候选自相关函数的信息。由于k<L,可将式(15)的问题转换为以下问题:
接下来推导这一限制的重要性。
①残差自相关的半正定特性
I-p(ωm)U(ωm)Λ(ωm)UH(ωm)≥0
→I-p(ωm)Λ(ωm)≥0 式(24)
上式定义了k个不同的约束条件,其中λ1(ωm)≥λ2(ωm)≥...≥λk(ωm)是Λ(ωm)的特征值。注意,如果对于q=1满足式(24),那么对任意q均满足式(25),进而可以推导出式(26)。
I-p(ωm)λq(ωm)≥0 q=1,2,...,k 式(25)
②施加稀疏约束
为了更准确的得到p的值,本发明选择系数序列具有最小l1范数(的p的值)
其中w0,...,wM-1是正权值,因为p(ωm)的值大于等于零,为了简单起见,去除绝对值。β的值取决于所选择的权重。纠正l1范数大小的依赖性的一种方法是强制每个乘积wmp(ωm)等于1。理想的情况下,
式(28)中的wm的值可以由式(30)得出,那么式(27)中β的值近似等于被检查频带中存在的主要用户的数量K,决定了向量p的稀疏度水平。
wm=λ1(ωm) 式(30)
4)稀疏信号重构算法
在正交匹配追踪(Orthogonal Matching Pursuit,OMP)算法的基础上,通过增加权值可以更好地检测出信号的频谱位置和功率值水平。该算法的具体步骤包括如下:
步骤4:Λt=Λt-1∪{λt}
步骤5:At=[At-1aλt],A0是空矩阵。
步骤9:t=t+1
步骤10:end while
步骤11:返回pt
Claims (6)
1.一种基于频谱特征的宽带频谱感知方法,其特征在于,所述方法采用瑞利衰落信道模型,将主用户的频谱特征当做唯一的先验条件,通过比较接收信号的功率谱密度与主用户的频谱信息,求出信号功率值和频谱占用情况;
所述方法首先建立基于瑞利衰落信道的信号传输模型,其次根据主用户和接收信号的自相关函数建立函数关系,然后对信号进行压缩采样,再利用基于频谱特征的压缩感知技术,采用相关匹配策略,通过保证残差自相关函数的半正定特性,利用l1-最小化公式推导出算法中的权值表示,最后利用加权正交匹配追踪(Weighted Orthogonal MatchingPursuit,WOMP)算法识别主用户占用的频谱,估计对应的功率值水平;
所述方法接收信号经过多陪集采样后自相关函数模型,包括:设表示主用户可能的频率位置的集合,并且令表示能包含的网格,γi(τ)和γn(τ)分别表示干扰和噪声的自相关函数,p(ωm)是对应于主要用户的的发射功率,γcm(τ,ωm)表示参考信号xc(t)在载波频率ωm处的参考自相关函数,定义稀疏向量p为检测函数的输出,零值直接表明在ω处不存在主用户,非零值直接表明在ω处主用户存在;
所述方法对通过瑞利信道后的信号y(t)进行压缩采样,采用的是多陪集(Multi-coset,MC)采样:给定接收的多频带信号y(t),在一定时刻ti(n)得到MC采样样本,其中L>0是适当的整数,i=1,2...,k,n∈Z,集合{ci}包含从{0,1...,L-1}中选择的k个不同的整数,注意到MC采样过程被等效为经典的以块的形式进行奈奎斯特采样,从每块的L个样本中周期性保留k个样本,未丢弃的样本由集合{ci}指定,每个序列的周期等于LT;
所述方法将zf看做y(t)中均匀奈奎斯特采样样本的第f块,每块长度为L,多陪集采样后的总的样本数据可用Y表示,总块数为Nf,为了将采集的样本与原始奈奎斯特采样信号相关联,令yf表示每块的采样样本,每块有k个非均匀样本,Φ∈Ck×L是一个矩阵,它的作用是从zf中非均匀选择k个样本,Φ通过对单位矩阵IL非均匀抽取k行实现。
3.根据权利要求1所述的一种基于频谱特征的宽带频谱感知方法,其特征在于,所述方法基于相关匹配的压缩频谱感知,包括:令{Ry(k)}表示奈奎斯特采样信号{y(n)}的相关性,Ry是维度为L×L的相关矩阵,将{xc(n)}定义为奈奎斯特采样候选信号,{Rc(k)}定义为{xc(n)}的相关函数,相关矩阵表示采样样本相关性的平均值,基带参考自相关是从主用户的频谱特征中抽取出来的,将纯音的自相关函数用e(ωm)eH(ωm)给出,e(ωm)∈CL×1是频率转向向量,利用e(ωm)在感知频率ωm处形成的秩为1的矩阵对Rc调制,调制参考相关性Rcm(ωm)必须也以对信号同样的方法进行压缩,利用相同的采Ry样矩阵Φ。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710795476.2A CN107517089B (zh) | 2017-09-06 | 2017-09-06 | 一种基于频谱特征的宽带频谱感知方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710795476.2A CN107517089B (zh) | 2017-09-06 | 2017-09-06 | 一种基于频谱特征的宽带频谱感知方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107517089A CN107517089A (zh) | 2017-12-26 |
CN107517089B true CN107517089B (zh) | 2020-09-01 |
Family
ID=60724682
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710795476.2A Active CN107517089B (zh) | 2017-09-06 | 2017-09-06 | 一种基于频谱特征的宽带频谱感知方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107517089B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109547134B (zh) * | 2018-12-30 | 2021-10-26 | 清华大学 | 基于多发射功率等级的协作频谱感知方法和装置 |
CN110311744A (zh) * | 2019-07-10 | 2019-10-08 | 西南交通大学 | 一种基于Catboost算法的信道环境自适应频谱感知方法 |
CN110611542B (zh) * | 2019-08-28 | 2021-06-22 | 南京邮电大学 | 一种基于改进融合算法的三阈值协作频谱感知方法 |
CN110677204B (zh) * | 2019-11-20 | 2021-03-02 | 南京邮电大学 | 一种存在智能干扰的频谱感知时间优化方法 |
CN110932753B (zh) * | 2019-12-11 | 2021-04-09 | 电子科技大学 | 一种基于智能决策的变换域自适应通信传输方法 |
CN112543073B (zh) * | 2020-11-27 | 2022-03-08 | 电子科技大学 | 基于次奈奎斯特采样的联合宽带频谱感知和载频估计方法 |
CN113037406B (zh) * | 2020-12-29 | 2022-07-05 | 杭州电子科技大学 | 一种时频特性提取及压缩感知融合的高效协作频谱感知方法 |
CN116318476B (zh) * | 2023-05-23 | 2023-09-19 | 北京星河亮点技术股份有限公司 | 基于功率比较的信道占用状态预测方法及装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101420758A (zh) * | 2008-11-26 | 2009-04-29 | 北京科技大学 | 一种认知无线电中对抗模仿主用户攻击的方法 |
CN103001684A (zh) * | 2012-12-18 | 2013-03-27 | 清华大学 | 多天线认知无线电通信系统中的鲁棒干扰抑制方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104682963A (zh) * | 2015-03-03 | 2015-06-03 | 北京邮电大学 | 一种信号循环平稳特性的重构方法 |
CN105915473B (zh) * | 2016-05-26 | 2019-07-12 | 中南大学 | 一种基于压缩感知技术的ofdm系统参数化信道估计及均衡方法 |
-
2017
- 2017-09-06 CN CN201710795476.2A patent/CN107517089B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101420758A (zh) * | 2008-11-26 | 2009-04-29 | 北京科技大学 | 一种认知无线电中对抗模仿主用户攻击的方法 |
CN103001684A (zh) * | 2012-12-18 | 2013-03-27 | 清华大学 | 多天线认知无线电通信系统中的鲁棒干扰抑制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107517089A (zh) | 2017-12-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107517089B (zh) | 一种基于频谱特征的宽带频谱感知方法 | |
CN110061792B (zh) | 一种基于变分模态分解的频谱感知方法 | |
WO2005120000A1 (en) | Fourier-transform based linear equalization for mimo cdma downlink | |
CN103051403A (zh) | 基于多个mwc分布式亚奈奎斯特采样联合重构的频谱感知方法 | |
Leinonen et al. | Compressed sensing with applications in wireless networks | |
CN110768754B (zh) | 一种信号检测方法及装置 | |
CN103532645A (zh) | 一种观测矩阵优化的压缩频谱感知方法 | |
CN105024957B (zh) | 一种直扩码分多址上行链路信道估计方法及装置 | |
CN109818645B (zh) | 基于信号检测与支撑集辅助的叠加csi反馈方法 | |
CN103139828B (zh) | 一种宽带频谱感知装置及方法 | |
CN103346984A (zh) | 基于bsl0的局部聚类稀疏信道估计方法 | |
WO2009032656A1 (en) | Receiver and method for blind adaptive thresholding using noise estimator | |
CN104682963A (zh) | 一种信号循环平稳特性的重构方法 | |
CN114641972A (zh) | 估计过载通信信道中的发射符号向量的方法 | |
CN106878997A (zh) | 一种分布式网络中的贝叶斯压缩频谱感知方法 | |
CN113037406B (zh) | 一种时频特性提取及压缩感知融合的高效协作频谱感知方法 | |
CN115201759B (zh) | 一种基于奇异值分解的雷达嵌入式通信波形设计方法 | |
CN114978818B (zh) | 一种基于压缩感知的自适应信道估计方法及估计系统 | |
CN104980238B (zh) | 一种基于群稀疏的协作频谱感知方法 | |
Hu et al. | Adaptive and blind wideband spectrum sensing scheme using singular value decomposition | |
El-Khamy et al. | C9. a new wideband spectrum sensing technique based on compressive sensing, wavelet edge detection and parallel processing | |
CN114244458A (zh) | 一种亚奈奎斯特采样前端的全盲频谱感知方法 | |
CN108199755B (zh) | 一种无线通信的码本和数据导频功率比的检测方法及系统 | |
Astaiza et al. | Compressive local wideband spectrum sensing algorithm for multiantenna cognitive radios | |
Fu et al. | Simplified cyclostationary detector using compressed sensing |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
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 |