CN112261684A - 一种基于mub的拟合优度检验的频谱感知方法 - Google Patents
一种基于mub的拟合优度检验的频谱感知方法 Download PDFInfo
- Publication number
- CN112261684A CN112261684A CN202011314739.1A CN202011314739A CN112261684A CN 112261684 A CN112261684 A CN 112261684A CN 202011314739 A CN202011314739 A CN 202011314739A CN 112261684 A CN112261684 A CN 112261684A
- Authority
- CN
- China
- Prior art keywords
- mub
- goodness
- matrix
- fit test
- sample set
- 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
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W24/00—Supervisory, monitoring or testing arrangements
- H04W24/06—Testing, supervising or monitoring using simulated traffic
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
本发明公开了一种基于MUB的拟合优度检验的频谱感知方法,包括步骤如下:S1:基于MUB构造MUB矩阵,将采集到的样本集经过MUB矩阵转换,得到M2个新的样本集Z,M表示MUB矩阵的大小;S2:接收端对M2个新的样本集Z进行拟合优度检验,计算拟合优度检验的判决统计量;S3:根据给定的虚警概率Pfa要求和拟采用的数据样本长度n,确定检测门限γ,结合判决统计量,并依据判决规则来判断频谱使用状态。本发明通过在拟合优度检验中引入MUB方法,充分利用拟合优度检验在频谱检测中具有优秀检测性能的优势,同时利用MUB减少了需要的样本数量,提高了检测效率。
Description
技术领域
本发明涉及无线通信技术领域,更具体的,涉及一种基于MUB的拟合优度检验的频谱感知方法。
背景技术
随着无线通信技术飞速发展,由于大量的无线通信需求需要大量的无线频谱,频谱资源分配紧张的问题日益凸显。而另一方面,频段的平均利用率却很低。造成这一现象的很大一部分原因在于目前国际上所普遍采用的静态频谱分配机制,即为特定的无线通信系统或业务划定一段固定的频率范围。此举虽可有效避免不同系统或业务间的相互干扰,却也使得本就稀缺的无线频谱资源浪费严重。为了缓解上述频谱资源供需间的矛盾,近年来,认知无线电(Cognitive Radio,CR) 技术,作为一种智能的频谱共享技术,在此契机下应运而生。各种频谱感知技术已被提出。
基于拟合优度检验的算法直接从统计分布函数着手,能够充分挖掘数据所携带的信息,做出正确的判决。2009年,文献[H.Wang,E.H.Yang,Z.Zhao,W.Zhang, Spectrumsensing in cognitive radio using goodness of fit testing,IEEE Trans. WirelessCommun,2009,8(11):5427-5430]首次将频谱检测问题表述为非参数的拟合优度检验问题,而后,此方法得到了广泛的研究。其方法的基本思路是:当没有信号传输时,观测样本应该是来自独立的噪声分布,而当信号出现时,肯定会改变观测样本的分布,因此检测信号的出现与否等价于检测观测样本是否服从噪声分布,这便是一个典型的拟合优度检验问题。由此可知,基于拟合优度检验的频谱检测方法是考虑信号的整体特征,而非信号的某个参数特征,并且,它不需要任何信号的先验信息。
发明内容
本发明为了解决现有技术频谱感知的性能低下,且拟合优度检验需要较多样本的问题,提供了一种基于MUB的拟合优度检验的频谱感知方法,其能提高频谱感知的性能,同时提出基于MUB的样本扩充方法,克服了拟合优度检验需要较多样本的不足,提高检测效率,从而获得满意的频谱检测效果。
为实现上述本发明目的,采用的技术方案如下:一种基于MUB的拟合优度检验的频谱感知方法,所述的方法包括步骤如下:
S1:基于MUB构造MUB矩阵,将采集到的样本集经过MUB矩阵转换,得到M2个新的样本集Z;其中M表示MUB矩阵的大小;
S2:接收端对M2个新的样本集Z进行拟合优度检验,计算拟合优度检验的判决统计量;
S3:确定检测门限γ,结合判决统计量,并依据判决规则来判断频谱使用状态。
优选地,在步骤S1之前,先对MUB进行定义,具体如下:
若B1,B2,…,Bd是M×M阶d个复酉阵,即
记Bj=(bj1bj2…bjM),即bjk是矩阵Bj中的第k列,k=1、2、…、M,如果对于任何的i≠j,有
进一步地,根据M的不同值构造MUB矩阵,最终生成一个维度为M×M2的矩阵。
再进一步,当M=p1,p1为素数,MUB的构造如下:
当a,b都取定而x遍历数域即x依次取从0到p1-1时,得到一个p1维的向量va,b;当a取定而b,x遍历数域便得到p1个向量,将这p1个向量称作一组基Ba;当a,b,x遍历数域便得到p1个基,最终生成一个维度为M×M2的矩阵。
再进一步,当M=q=p2 s,p2为大于2的奇素数,s大于1的整数,MUB 的构造如下:令
再进一步,当M=q=2n,n>1时,MUB的构造如下:
再进一步,设{y1,y2,…,yM}为样本集,设矩阵Bi为MUB矩阵中的任一个矩阵,令
(zi1,zi2,…,ziM)=(y1,y2,…,yM)Bi
其中,i=1,2,…,M,从而得到M2个新的样本集Z
Z=(z11,z12,…,z1M,…,zM1,zM2,…,zMM)
基于这M2个新的样本集,进行拟合优度检验。
再进一步,步骤S2,所述的判决统计量计算公式如下:
式中,F*()分布函数如下所示:
式中,δ2为噪声方差。
再进一步地,在得到新的样本集Z之后,进行拟合优度检验之前,对新的样本集Z进行升序处理,再进行拟合优度检验。
再进一步,步骤S3中,判断频谱使用状态,具体如下:
其中,H0表示信道没有被占用,而H1表示信道被占用。
本发明的有益效果如下:
本发明通过在拟合优度检验中引入MUB方法,充分利用拟合优度检验在频谱检测中具有优秀检测性能的优势,同时利用MUB减少了需要的样本数量,提高了检测效率。具体而言,相比于非拟合优度检验方法,本发明所提方法具有更佳的检测性能;相比于传统的基于拟合优度的检验方法,本发明所提方法可以显著减少样本数量(或者在样本数量一致下,进一步显著提高检测性能)。
附图说明
图1是本实施例所述的方法的流程图。
图2是本实施例M=16时,采用KS与本实施例所述的方法检验仿真对比结果图。
图3是本实施例M=32时,采用KS与本实施例所述的方法检验检验仿真对比结果。
具体实施方式
下面结合附图和具体实施方式对本发明做详细描述。
实施例1
如图1所示,一种基于MUB的拟合优度检验的频谱感知方法,所述的方法包括步骤如下:
S1:基于MUB构造MUB矩阵,将采集到的样本集经过MUB矩阵转换,得到M2个新的样本集Z;
S2:接收端对M2个新的样本集Z进行拟合优度检验,计算拟合优度检验的判决统计量;
S3:根据虚警概率Pfa数据样本长度n和密度函数确定检测门限γ,结合判决统计量,并依据判决规则来判断频谱使用状态。
在一个具体的实施例中,在步骤S1之前,先对MUB进行定义,具体如下:
若B1,B2,…,Bd是M×M阶d个复酉阵,即
记Bj=(bj1bj2…bjM),即bjk是矩阵Bj中的第k列,如果对于任何的i≠j,有
在一个具体的实施例中,MUB矩阵的构造具体如下:
a)当M=p1,p1为素数,MUB的构造如下:
当a,b都取定而x遍历数域即x依次取从0到p1-1时,得到一个p1维的向量va,b;当a取定而b,x遍历数域便得到p1个向量,将这p1个向量称作一组基Ba;当a,b,x遍历数域便得到p1个基,最终生成一个维度为M×M2的矩阵。
b)当M=q=ps(p为大于2的奇素数,s为大于1的整数)时,MUB的构造具体如下:令
c)当M=q=2n(n>1)时MUB的构造具体如下:
在一个具体的实施例中,基于MUB的拟合优度(KS)检验,MUB矩阵中的每个酉矩阵都表示了空间中的一个座标架,并且不同座标架之间的相关系数都很小,仅为这样,样本经过MUB矩阵变换后,就可以产生几乎独立的新的样本。基于这些新样本,可以做出更加准确的判决。具体的数学过程如下所示,应用场景为加性高斯白噪声信道:
设{y1,y2,…,yM}为样本集,传统的拟合优度(KS)检验的判决统计量的计算如下所示:
本实施例基于MUB的拟合优度(KS)检验的判决统计量如下
设矩阵Bi为MUB矩阵中的任一个矩阵,令
(zi1,zi2,…,ziM)=(y1,y2,…,yM)Bi
其中i=1,2,…,M;从而得到M2个新的样本集Z
Z=(z11,z12,…,z1M,…,zM1,zM2,…,zMM)
对M2个新的样本集Z先进行升序处理之后,再进行拟合优度检验。
其判决统计量计算公式如下:
其中,分布函数和经验函数如下所示:
式中,δ2为噪声方差。
在一个具体的实施例中,步骤S3中,判断频谱使用状态,具体如下:
其中,H0表示信道没有被占用,而H1表示信道被占用。
在本实施例中,在仿真前必须确定好相应的门限值。本次仿真的虚警概率统一定为5%,通过10000次蒙特卡洛确定检验门限值。
表1两种检验方式的门限值
如图2所示,当M=16时,KS与本实施例所述的方法检验仿真对比结果图;如图3所示,当M=32时KS检验仿真对比结果。
从图2、图3中可以得到,当M分别为16和32时,在任意信噪比的情况下,基于MUB的拟合优度检验方法的检测概率总是高于传统的拟合优度检验,并且随着信噪比的提高,其检测概率更快的逼近1。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明权利要求的保护范围之内。
Claims (10)
1.一种基于MUB的拟合优度检验的频谱感知方法,其特征在于:所述的方法包括步骤如下:
S1:基于MUB构造MUB矩阵,将采集到的样本集经过MUB矩阵转换,得到M2个新的样本集Z;其中,M表示MUB矩阵的大小;
S2:接收端对M2个新的样本集Z进行拟合优度检验,计算拟合优度检验的判决统计量;
S3:确定检测门限γ,结合判决统计量,并依据判决规则来判断频谱使用状态。
3.根据权利要求2所述的基于MUB的拟合优度检验的频谱感知方法,其特征在于:根据M的不同值构造MUB矩阵,最终生成一个维度为M×M2的矩阵。
7.根据权利要求2~6任一项所述的基于MUB的拟合优度检验的频谱感知方法,其特征在于:设{y1,y2,…,yM}为样本集,设矩阵Bi为MUB矩阵中的任一个矩阵,令
(zi1,zi2,…,ziM)=(y1,y2,…,yM)Bi
其中,i=1,2,…,M,从而得到M2个新的样本集Z
Z=(z11,z12,…,z1M,…,zM1,zM2,…,zMM)
基于M2个新的样本集Z,进行拟合优度检验。
8.根据权利要求7所述的基于MUB的拟合优度检验的频谱感知方法,其特征在于:在得到新的样本集Z之后,进行拟合优度检验之前,对新的样本集Z进行升序处理,再进行拟合优度检验。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011314739.1A CN112261684B (zh) | 2020-11-20 | 2020-11-20 | 一种基于mub的拟合优度检验的频谱感知方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011314739.1A CN112261684B (zh) | 2020-11-20 | 2020-11-20 | 一种基于mub的拟合优度检验的频谱感知方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112261684A true CN112261684A (zh) | 2021-01-22 |
CN112261684B CN112261684B (zh) | 2023-06-16 |
Family
ID=74225029
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011314739.1A Active CN112261684B (zh) | 2020-11-20 | 2020-11-20 | 一种基于mub的拟合优度检验的频谱感知方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112261684B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101329400A (zh) * | 2008-07-30 | 2008-12-24 | 电子科技大学 | 一种基于拟合优度检验的雷达目标恒虚警检测方法 |
EP2086255A1 (en) * | 2008-02-01 | 2009-08-05 | Institut Eurecom | Process for sensing vacant sub-space over the spectrum bandwidth and apparatus for performing the same |
CN104320209A (zh) * | 2014-10-14 | 2015-01-28 | 宁波大学 | 一种基于拟合优度检验的频谱感知方法 |
US20150057041A1 (en) * | 2013-08-22 | 2015-02-26 | Nokia Corporation | Blind Spectrum Sensing Based on Maximum Correlation Coefficients and use Thereof |
-
2020
- 2020-11-20 CN CN202011314739.1A patent/CN112261684B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2086255A1 (en) * | 2008-02-01 | 2009-08-05 | Institut Eurecom | Process for sensing vacant sub-space over the spectrum bandwidth and apparatus for performing the same |
CN101329400A (zh) * | 2008-07-30 | 2008-12-24 | 电子科技大学 | 一种基于拟合优度检验的雷达目标恒虚警检测方法 |
US20150057041A1 (en) * | 2013-08-22 | 2015-02-26 | Nokia Corporation | Blind Spectrum Sensing Based on Maximum Correlation Coefficients and use Thereof |
CN104320209A (zh) * | 2014-10-14 | 2015-01-28 | 宁波大学 | 一种基于拟合优度检验的频谱感知方法 |
Non-Patent Citations (1)
Title |
---|
包志强;韩彦妮;崔妍蕊;: "采用T~2统计量和拟合优度的盲频谱感知算法", 信号处理, no. 01 * |
Also Published As
Publication number | Publication date |
---|---|
CN112261684B (zh) | 2023-06-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102013928B (zh) | 一种认知无线电系统中的快速频谱感知方法 | |
CN110932807B (zh) | 一种非高斯噪声下mimo系统的频谱感知方法及系统 | |
CN103297160A (zh) | 基于归一化特征值的拟合优度检验的频谱感知方法及装置 | |
CN107276696A (zh) | 一种低信噪比条件下信号检测算法 | |
CN107886113B (zh) | 一种基于卡方检验的电磁频谱噪声提取和滤波方法 | |
CN103346849B (zh) | 一种抵抗模仿主用户信号的恶意攻击的频谱感知方法 | |
CN112261684A (zh) | 一种基于mub的拟合优度检验的频谱感知方法 | |
CN114465681B (zh) | 一种用于电力物联网的多节点协作频谱感知方法及装置 | |
Rauniyar et al. | Cooperative adaptive threshold based energy and matched filter detector in cognitive radio networks | |
Rabelo et al. | On the Multivariate $\alpha $–$\mu $ Distribution: New Exact Analytical Formulations | |
CN113179143A (zh) | 一种基于最大最小值的融合频谱感知方法及系统 | |
Achtzehn et al. | IEEE 802.11 Wi-Fi access point density estimation with capture-recapture models | |
CN114759996B (zh) | 基于循环相关熵似然比检验的频谱感知方法 | |
Guedes et al. | Direct-conversion spectrum sensor impaired by symmetric α-stable and α-sub-gaussian noises | |
CN115021840B (zh) | 信息处理方法及装置、设备、存储介质 | |
Ali et al. | LTE signal detection using two-stage cooperative compressive sensing system | |
Huang et al. | Space-time correlation based fast regional spectrum sensing in cognitive radio | |
CN103560838B (zh) | 一种抑制直流偏置的能量检测方法 | |
CN116318444B (zh) | 二维频谱感知方法、装置、电子设备和存储介质 | |
CN107181548A (zh) | 一种压缩频谱感知性能提升方法 | |
Ismaeel et al. | Enhancement of Cognitive Radio Spectrum with Proposed Cascade Filter Model | |
CN106998235A (zh) | 主用户信号随机出现时基于循环平稳的频谱感知方法 | |
CN109672491B (zh) | 一种信号处理方法、装置、设备及计算机可读存储介质 | |
CN114448536B (zh) | 一种全双工频谱感知方法 | |
Chen et al. | Noise threshold estimation in spectrum monitoring data analysis application |
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 |