CN102223195B - 一种频谱感知方法 - Google Patents

一种频谱感知方法 Download PDF

Info

Publication number
CN102223195B
CN102223195B CN2011102152664A CN201110215266A CN102223195B CN 102223195 B CN102223195 B CN 102223195B CN 2011102152664 A CN2011102152664 A CN 2011102152664A CN 201110215266 A CN201110215266 A CN 201110215266A CN 102223195 B CN102223195 B CN 102223195B
Authority
CN
China
Prior art keywords
frequency spectrum
parameter
sensing method
spectrum sensing
signal
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.)
Expired - Fee Related
Application number
CN2011102152664A
Other languages
English (en)
Other versions
CN102223195A (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN2011102152664A priority Critical patent/CN102223195B/zh
Publication of CN102223195A publication Critical patent/CN102223195A/zh
Application granted granted Critical
Publication of CN102223195B publication Critical patent/CN102223195B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明属于认知无线电技术领域,公开了一种频谱感知方法。针对现有的固定参数SR系统在频谱感知过程中存在的问题,本发明的频谱感知方法首先固定SR系统的参数a,利用接收信号r(t)获得噪声方差
Figure DDA0000079628300000011
最后由a和
Figure DDA0000079628300000012
的值确定SR系统的另一个参数b。本发明的感知方法通过调整SR系统参数从而使得随机共振系统与接收信号特性相匹配,具有自适应的特点,更容易产生共振,明显提高了输出信号的SNR,提高了检测概率,更容易感知授权用户是否存在,由此提高了频谱感知的性能。

Description

一种频谱感知方法
技术领域
本发明属于认知无线电技术领域,具体涉及一种频谱感知方法。
背景技术
认知无线电(Cognitive Radio,CR)技术是为了解决目前面临的频谱资源稀缺提出来的,而频谱感知算法是认知无线电的关键技术之一。为了避免认知无线电系统对授权用户产生有害的干扰,要求频谱感知方法能够在低信噪比下可靠地检测出授权用户信号来。能量检测是常用的一种频谱感知方法,但是在低信噪比(Signal-to-Noise Ratio,SNR)下检测概率较低,不满足认知无线电技术的要求。
研究发现把要感知的接收信号通过参数匹配的随机共振(Stochastic Resonance,SR)系统能够明显的提高其输出的SNR,于是把SR理论应用于能量检测(Energy Detection,ED)方法中能够提高检测概率,从而提高了认知无线电频谱感知的性能,具体如下:
频谱感知以二元假设检验为模型: H 1 : r ( t ) = s ( t ) + n ( t ) H 0 : r ( t ) = n ( t ) , 其中,假设H1表示存在授权用户信号,H0则表示不存在授权用户信号;r(t)表示要感知的接收信号,s(t)表示授权用户信号,n(t)表示零均值、方差为
Figure GDA00003110374800015
的加性白高斯噪声(Additive White Gaussian Noise,AWGN)。
在认知无线电的频谱感知中通常采用如下的非线性langevin方程描述双稳态SR系统:
Figure GDA00003110374800012
其中,V(x)是双稳态SR系统的势函数,其表达式为:
Figure GDA00003110374800013
其中,a,b为SR系统参数,把要感知的接收信号r(t)=s(t)+n(t)代入langevin方程中求解得到x(t),x(t)通过能量检测,判断授权用户是否信号存在,即通过
Figure GDA00003110374800014
计算得到检验统计量T(r),把T(r)与门限r进行比较,判断授权用户信号是否存在。这里,N表示每次频谱感知累积样本数,也即是对r(t)采样的采样点数,从采样序列中取N个采样点用来计算统计量T(r),l是计算T(r)的开始的采样序列序号。
随机共振系统的参数要求与接受信号匹配才能较好的提高输出SNR,但是由于外界的干扰使得接收信号r(t)不断变化,所以固定参数的SR系统很难与接收信号良好匹配,从而使得SR系统输出SNR的提高不能达到期望的效果,并且很多时候检测性能不能满足频谱感知的实际需要。
发明内容
本发明的目的是为了解决现有的固定参数SR系统在频谱感知过程中存在的问题,提出了一种频谱感知方法。
本发明的技术方案是:一种频谱感知方法,包括如下步骤:
S1.通过接收设备接受要感知的接收信号r(t),其中,r(t)=s(t)+n(t),s(t)为授权用户信号;n(t)是均值为零方差为
Figure GDA00003110374800021
的高斯白噪声;
S2.确定SR系统参数b:所述SR系统通过langevin方程进行描述,其中,
Figure GDA00003110374800023
固定参数a,使参数a满足绝热近似条件,然后根据接收信号r(t)获得噪声方差
Figure GDA00003110374800024
由a和
Figure GDA00003110374800025
的值确定参数b;
S3.接收信号r(t)通过langevin方程求得x(t),x(t)通过能量检测,判断授权用户信号是否存在。
本发明的有益效果:本发明的频谱感知方法通过把随机共振理论应用到能量检测中,首先固定SR系统的参数a,利用接收信号r(t)获得噪声方差
Figure GDA00003110374800026
最后由a和
Figure GDA00003110374800027
的值确定SR系统的另一个参数b。本发明的感知方法通过调整SR系统参数从而使得随机共振系统与接收信号特性相匹配,具有自适应的特点,所以更容易产生共振,可以明显提高输出信号的SNR,提高了检测概率,更容易感知授权用户是否存在,由此提高了频谱感知的性能。
附图说明
图1为本发明的频谱感知方法示意图。
图2为本发明的频谱感知方法流程示意图。
图3为三种基于能量检测的频谱感知方法性能仿真比较示意图。
具体实施方式
下面结合图1和图2对本发明的频谱感知方法进行阐述,具体包括以下步骤:
S1.通过接收设备接受要感知的接收信号r(t),其中,r(t)=s(t)+n(t),s(t)为授权用户信号;n(t)是均值为零方差为
Figure GDA00003110374800028
的高斯白噪声;
S2.确定SR系统参数b:所述SR系统通过langevin方程
Figure GDA00003110374800029
进行描述,其中,
Figure GDA00003110374800031
固定参数a,使参数a满足绝热近似条件,然后根据接收信号r(t)获得噪声方差
Figure GDA00003110374800032
由a和
Figure GDA00003110374800033
的值确定参数b。参数b的具体确定过程如下:
利用绝热近似(Adiabatic Approximation)理论,当信号r(t)=s(t)+n(t)通过langevin方程定义的双稳态SR系统时,输出信号x(t)的SNR为:
SNR o = ( 2 a A m 2 c 2 σ n 4 e - 2 U 0 / σ n 2 ) ( 1 - 4 a 2 A m 2 c 2 π 2 σ n 2 e - 4 U 0 / σ n 2 2 a 2 π 2 e - 4 U 0 / σ n 2 + ( 2 π f s ) 2 ) - 1 ≈ 2 a A m 2 c 2 σ n 4 e - 2 U 0 / σ n 2
其中,a是SR系统参数,Am是授权用户信号s(t)的幅度,c是双稳态SR系统的势阱点,
Figure GDA00003110374800035
是高噪声的方差,U0=a2/(4b)是当Am=0时的双稳态SR系统的势垒高度。具体可参考文献:McNamara B,Wiesenfeld K.Theory of stochastic resonance,Physical Review A,1989,39(9):4854-4869。
为了满足绝热近似理论,要求a?πfs,其中fs为授权用户信号s(t)的频率。
由于输入信号的平均信噪比为:
Figure GDA00003110374800036
因此,当发生随机共振时,接收信号r(t)经过双稳态SR系统后的输出信噪比增益为:
Figure GDA00003110374800037
令k=a2/b,显然有k>0,则即给定噪声方差,输出信号SNR增益ηSNR是系统参数k的非线性函数。
ηSNR对k的二阶导数为:因此,为了使ηSNR是关于k的下凹函数,以便取得唯一的极大值,要求:
Figure GDA000031103748000310
于是,最大化SNR增益的最优k的取值满足: k op = arg max k &eta; SNR s . t . , 0 < k < 4 &sigma; n 2 , 求解上式可得:
Figure GDA000031103748000312
于是最大化SNR增益的双稳态SR系统的参数需满足 b = a 2 / ( 2 &sigma; n 2 ) a ? &pi; f s .
在本发明中,对于上式得到的a和b的关系可以通过一个调节系数h进行调整,即 b = ha 2 / ( 2 &sigma; n 2 ) .
这里只是给出了一种确立a和b关系的技术手段,本领域的普通技术人员应该意识到还可以通过其它方法来确定a和b关系。
S3.接收信号r(t)通过langevin方程求得x(t),x(t)通过能量检测,判断授权用户信号是否存在:具体如下:计算得到检验统计量T(r),把T(r)与门限r比较,如果T(r)>r则认为授权用户信号存在,即授权用户正在使用这一频段;如果T(r)<r则认为授权用户信号不存在,即授权用户没有使用这一频段。
下面对本发明方法进行仿真测试。仿真的参数为:输入正弦信号幅度Am=1,频率设为fs=0.01Hz,采样周期设为Δt=0.02s进行离散化,每次频谱感知累积样本数为N=1000;双稳态系统参数
Figure GDA00003110374800043
在仿真中假设噪声方差
Figure GDA00003110374800044
已知。
仿真表明:当调节系数h在区间[0.8,1.2]进行变动时,检测性能有很好的提高,在h=1时,性能得到最大的提高。
关于绝热近似条件,a?πfs,其中,在
Figure GDA00003110374800045
Figure GDA00003110374800046
时检测性能得到很好的提高,
Figure GDA00003110374800047
Figure GDA00003110374800048
时提高的性能略次于
Figure GDA000031103748000410
时的性能的提高,但也得到了明显的检测性能的提高,但是当
Figure GDA000031103748000411
是检测性能会明显下降。本领域的普通技术人员可以根据实际需要进行具体选择。
与固定参数的SR ED频谱感知方法以及传统ED频谱感知方法比较的结果如图3所示,这里固定参数的SR ED频谱感知方法选取a=2,b=1。从图3可以看出,当取虚警概率Pf=0.05时,本发明的频谱感知方法(选取h=1)、固定参数SR ED频谱感知方法和传统ED频谱感知方法的检测概率进行比较得出,当检测概率Pd=0.9时,本发明的方法相对于传统ED频谱感知方法有约10dB的性能增益,而基于固定参数SR的ED方法检测概率没有达到0.9。
本领域普通技术人员可以理解,实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件完成,所述的程序可以存储于可读存储介质中,例如只读存储器、随机存取存储器、磁盘、光盘等。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。

Claims (3)

1.一种频谱感知方法,包括如下步骤:
S1.通过接收设备接受要感知的接收信号r(t),其中,r(t)=s(t)+n(t),s(t)为授权用户信号;n(t)是均值为零方差为
Figure FDA00002933236100011
的高斯白噪声;
S2.确定SR系统参数b:所述SR系统通过langevin方程
Figure FDA00002933236100012
进行描述,其中,
Figure FDA00002933236100013
固定参数a,使参数a满足绝热近似条件,然后根据接收信号r(t)获得噪声方差
Figure FDA00002933236100014
由a和的值确定参数b;
S3.接收信号r(t)通过langevin方程求得x(t),x(t)通过能量检测,判断授权用户信号是否存在。
2.根据权利要求1所述的频谱感知方法,其特征在于,步骤S2所述的由a和
Figure FDA00002933236100016
的值确定参数b具体为:
Figure FDA00002933236100017
其中,h为调节系数,范围为:[0.8,1.2]。
3.根据权利要求2所述的频谱感知方法,其特征在于,所述的调节系数h为1。
CN2011102152664A 2011-07-29 2011-07-29 一种频谱感知方法 Expired - Fee Related CN102223195B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2011102152664A CN102223195B (zh) 2011-07-29 2011-07-29 一种频谱感知方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2011102152664A CN102223195B (zh) 2011-07-29 2011-07-29 一种频谱感知方法

Publications (2)

Publication Number Publication Date
CN102223195A CN102223195A (zh) 2011-10-19
CN102223195B true CN102223195B (zh) 2013-09-25

Family

ID=44779632

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2011102152664A Expired - Fee Related CN102223195B (zh) 2011-07-29 2011-07-29 一种频谱感知方法

Country Status (1)

Country Link
CN (1) CN102223195B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102710348B (zh) * 2012-05-08 2014-01-22 电子科技大学 一种基于周期图的能量检测方法
CN102710347B (zh) * 2012-05-08 2014-01-22 电子科技大学 一种基于偏差修正的能量检测方法
CN104850752A (zh) * 2015-05-29 2015-08-19 桂林电子科技大学 一种基于自适应随机共振的参数搜索范围确定方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101521896A (zh) * 2009-04-17 2009-09-02 东南大学 认知无线电中的基于似然比的协同频谱感知方法
CN102083101A (zh) * 2011-01-25 2011-06-01 东南大学 一种认知无线电传感器网络信息传输方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101408929B1 (ko) * 2008-02-27 2014-06-18 삼성전자주식회사 인지 무선 기반의 무선통신 시스템에서 제어 정보 송수신장치 및 방법

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101521896A (zh) * 2009-04-17 2009-09-02 东南大学 认知无线电中的基于似然比的协同频谱感知方法
CN102083101A (zh) * 2011-01-25 2011-06-01 东南大学 一种认知无线电传感器网络信息传输方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陈伟.认知无线电中的频谱感知与干扰评估技术研究.《认知无线电中的频谱感知与干扰评估技术研究》.2011,第9页,第23-25页,第30-34页. *

Also Published As

Publication number Publication date
CN102223195A (zh) 2011-10-19

Similar Documents

Publication Publication Date Title
CN102608553B (zh) 一种基于自适应随机共振的微弱信号提取方法
CN106130942B (zh) 一种基于循环谱的无线通信信号调制识别及参数估计方法
US20160197748A1 (en) Blind spectrum sensing method and device based on fast fourier transform
CN107612587A (zh) 一种用于跳频非合作通信中跳频信号的参数估计方法
CN101291055A (zh) 一种输电线路故障行波初始波头到达时刻精确标定方法
CN102253421B (zh) 一种红外探测系统
CN102223195B (zh) 一种频谱感知方法
CN101848177B (zh) 基于变频的双稳态最优随机共振单频微弱信号检测方法
CN102645589B (zh) 脉冲检测方法及系统
JP2014180900A (ja) 車両システム、車載装置、及び携帯機
CN101854219B (zh) 基于随机共振的能量检测装置及其检测方法
CN106170139B (zh) 一种频谱检测方法及系统
CN102710347B (zh) 一种基于偏差修正的能量检测方法
CN103281141A (zh) 一种盲频谱感知方法及装置
CN102236944B (zh) 一种适用于报警领域的红外探测系统
CN103888201B (zh) 一种利用空间分集的协作频谱感知方法
Lei et al. Spectrum sensing based on goodness of fit tests
CN104734793A (zh) 基于p次方的无线协作频谱感知的能量检测方法
CN102970083B (zh) 基于噪声自适应随机共振系统的频谱感知方法
CN105680964B (zh) 一种频谱感知方法和频谱感知系统、客户端和服务端
CN104320360B (zh) 一种基于分数阶傅里叶变换的线性调频信号时延估计方法
CN106254007A (zh) 一种应用于lte230系统的频谱感知方法和装置
CN102710348B (zh) 一种基于周期图的能量检测方法
CN102236945B (zh) 适用于报警领域的红外探测系统
CN102148650B (zh) 基于检测率与虚警率加权合并的能量检测器的检测方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20130925

Termination date: 20160729