CN102694611B - 认知无线电系统中循环谱的自适应快速感知方法 - Google Patents
认知无线电系统中循环谱的自适应快速感知方法 Download PDFInfo
- Publication number
- CN102694611B CN102694611B CN201210180859.6A CN201210180859A CN102694611B CN 102694611 B CN102694611 B CN 102694611B CN 201210180859 A CN201210180859 A CN 201210180859A CN 102694611 B CN102694611 B CN 102694611B
- Authority
- CN
- China
- Prior art keywords
- data
- cyclic
- cognitive radio
- signal
- spectral density
- 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
Links
Landscapes
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
- Circuits Of Receivers In General (AREA)
Abstract
认知无线电系统中循环谱的自适应快速感知方法,涉及一种循环谱的自适应快速感知方法。它是为了解决在认知无线电系统中循环谱感知的感知时延长的问题。本发明通过数据叠加,自适应数据长度控制的方法,在保证其感知精度的情况下,缩短其感知时延。本发明由于采用数据重叠的方式进行循环谱估计,在减小对数据量要求的同时,也减小了估计的方差,同时,本发明采用了数据自适应长度选择,通过对信噪比的估计,选择合适的数据长度,即当信道条件好时,信噪比高,选择较短的数据计算就能准确估计出循环谱图,而反之如果信道条件不好,则选择较多的数据进行循环谱图的估计,减少了数据的计算量。本发明适用于认知无线电系统中循环谱的自适应快速感知。
Description
技术领域
本发明涉及一种循环谱的自适应快速感知方法。
背景技术
在认知无线电系统中,频谱感知通常采用能量检测、匹配检测、干扰温度检测、以及循环谱检测方法。其中,能量检测法虽然感知速度快,但是感知精度不高,且感知门限难以选择;匹配滤波检测法由于在通常情况下缺少先验信息,所以在认知感知中不采用该方法;干扰温度检测法由于需要对授权用户精确定位,否则容易对授权用户造成干扰。循环谱检测技术检测精度高,能够有效避免噪声的干扰,但是计算量大,检测时间长。
在通信系统中,由于信号经过编码,调制,加入循环前缀等加工,所以信号会呈现出循环平稳特性,而噪声由于其随机性,所以没有循环平稳性。根据这个性质,通过循环谱的感知,就能很好地将其区分开来。但是循环谱感知的感知时延较长。
发明内容
本发明是为了解决在认知无线电系统中循环谱感知的感知时延长的问题,从而提供一种认知无线电系统中循环谱的自适应快速感知方法。
认知无线电系统中循环谱的自适应快速感知方法,它由以下步骤实现:
步骤一、接收机接收一段认知无线电信号,将所述认知无线电信号以采样率fs进行采样,并量化,量化后截取一段长度为N1的接收信号序列x(n);
接收机继续接收一段认知无线电信号,将接收到的认知无线电信号进行采样和量化,并存入一号寄存器中;
步骤二、记录步骤一中的接收信号序列x(n)中的位于后面的长度为P的重叠部分数据,将所述长度为P的重叠部分数据存入二号寄存器中;
步骤三、对步骤一中的接收信号序列x(n)进行快速傅里叶变换,获得变换后序列X(k);
步骤四、将步骤三获得的变换后序列X(k)等分成M段,M为正整数,并分别计算M段序列的和值;计算和值最小的一段序列的平均功率,并将该平均功率作为高斯白噪声的功率密度,记为:N0/2;
根据公式:
获得平均高斯白噪声功率;
步骤五、根据公式:
计算步骤一中的接收信号序列x(n)的循环谱密度函数Sα(k);
其中:
为单次计算的循环谱密度函数,参数I为重叠计算循环功率谱密度的次数,L为数据平滑次数,k为频率点,α为循环频率;l为整数;
步骤六、找到步骤五获得的循环谱密度函数Sα(k)在除α=0,k=0以外的幅值最大值的点,并以该点为中心、3dB带宽范围内的点作为估计信号平均功率的估计空间;并估计出该空间的方差σs;
步骤七、根据公式:
估计出步骤六获得的方差σs的信噪比SNR;
步骤八、将步骤七获得的信噪比SNR与预设的判决门限γ进行比较,若SNR≥γ,或者数据重叠的总次数大于预设值Z,则将此循环谱密度函数Sα(k)作为估计信号有无的谱密度函数,执行步骤十;否则,执行步骤九;
步骤九、从一号寄存器中读取一段长度为D的数据,其中D=N0-P;将该段数据与二号寄存器中的长度为P的重叠部分数据组成一段长度为N0的数据x′(n),并返回执行步骤三;
步骤十、对步骤八获得的循环谱密度函数Sα(k)进行通过公式:
获得循环频率点的能量,并与门限值ξ进行比较,若ε≥η,则有主用户存在,完成认知无线电系统中循环谱的自适应快速感知;否则,则没有主用户存在,完成认知无线电系统中循环谱的自适应快速感知。
步骤十中门限值ξ是根据公式:
ξ=η×SNR
获得的;式中:η为信噪比系数,η为实数。
步骤六、步骤七、步骤八和步骤九中处理的数据长度根据通过估计当前信噪比、设定的门限值以及重叠数据的计算次数自适应地设定。
本发明通过数据叠加,自适应数据长度控制的方法,在保证其感知精度的情况下,缩短其感知时延。
具体实施方式
具体实施方式一、认知无线电系统中循环谱的自适应快速感知方法,它由以下步骤实现:
步骤一、接收机接收一段认知无线电信号,将所述认知无线电信号以采样率fs进行采样,并量化,量化后截取一段长度为N1的接收信号序列x(n);
接收机继续接收一段认知无线电信号,将接收到的认知无线电信号进行采样和量化,并存入一号寄存器中;
步骤二、记录步骤一中的接收信号序列x(n)中的位于后面的长度为P的重叠部分数据,将所述长度为P的重叠部分数据存入二号寄存器中;
步骤三、对步骤一中的接收信号序列x(n)进行快速傅里叶变换,获得变换后序列X(k);
步骤四、将步骤三获得的变换后序列X(k)等分成M段,M为正整数,并分别计算M段序列的和值;计算和值最小的一段序列的平均功率,并将该平均功率作为高斯白噪声的功率密度,记为:N0/2;
根据公式:
获得平均高斯白噪声功率;
步骤五、根据公式:
计算步骤一中的接收信号序列x(n)的循环谱密度函数Sα(k);
其中:
为单次计算的循环谱密度函数,参数I为重叠计算循环功率谱密度的次数,L为数据平滑次数,k为频率点,α为循环频率;l为整数;
步骤六、找到步骤五获得的循环谱密度函数Sα(k)在除α=0,k=0以外的幅值最大值的点,并以该点为中心、3dB带宽范围内的点作为估计信号平均功率的估计空间;并估计出该空间的方差σs;
步骤七、根据公式:
估计出步骤六获得的方差σs的信噪比SNR;
步骤八、将步骤七获得的信噪比SNR与预设的判决门限γ进行比较,若SNR≥γ,或者数据重叠的总次数大于预设值Z,则将此循环谱密度函数Sα(k)作为估计信号有无的谱密度函数,执行步骤十;否则,执行步骤九;
步骤九、从一号寄存器中读取一段长度为D的数据,其中D=N0-P;将该段数据与二号寄存器中的长度为P的重叠部分数据组成一段长度为N0的数据x′(n),并返回执行步骤三;
步骤十、对步骤八获得的循环谱密度函数Sα(k)进行通过公式:
获得循环频率点的能量,并与门限值ξ进行比较,若ε≥η,则有主用户存在,完成认知无线电系统中循环谱的自适应快速感知;否则,则没有主用户存在,完成认知无线电系统中循环谱的自适应快速感知。
步骤十中门限值ξ是根据公式:
ξ=η×SNR
获得的;式中:η为信噪比系数,η为实数。
步骤六、步骤七、步骤八和步骤九中处理的数据长度根据通过估计当前信噪比、设定的门限值以及重叠数据的计算次数自适应地设定。
步骤二和步骤五中通过对数据的叠加处理,分别计算出各自片段的循环谱密度函数,并将其相加求和,得到最终的循环谱图。
本发明由于采用数据重叠的方式进行循环谱估计,在减小对数据量要求的同时,也减小了估计的方差,并且就这一点而言,它的计算量是与传统的循环谱感知是相同的。同时,本发明采用了数据自适应长度选择,通过对信噪比的估计,选择合适的数据长度,即当信道条件好时,信噪比高,选择较短的数据计算就能准确估计出循环谱图,而反之如果信道条件不好,则选择较多的数据进行循环谱图的估计,从这个角度说,减少了数据的计算。
Claims (1)
1.认知无线电系统中循环谱的自适应快速感知方法,其特征是:它由以下步骤实现:
步骤一、接收机接收一段认知无线电信号,将所述认知无线电信号以采样率fs进行采样,并量化,量化后截取一段长度为N1的接收信号序列x(n);
接收机继续接收一段认知无线电信号,将接收到的认知无线电信号进行采样和量化,并存入一号寄存器中;
步骤二、记录步骤一中的接收信号序列x(n)中的位于后面的长度为P的重叠部分数据,将所述长度为P的重叠部分数据存入二号寄存器中;
步骤三、对步骤一中的接收信号序列x(n)进行快速傅里叶变换,获得变换后序列X(k);
步骤四、将步骤三获得的变换后序列X(k)等分成M段,M为正整数,并分别计算M段序列的和值;计算和值最小的一段序列的平均功率,并将该平均功率作为高斯白噪声的功率密度,记为:N0/2;
根据公式:
获得平均高斯白噪声功率;
步骤五、根据公式:
计算步骤一中的接收信号序列x(n)的循环谱密度函数Sα(k);
其中:
为单次计算的循环谱密度函数,参数I为重叠计算循环功率谱密度的次数,L为数据平滑次数,k为频率点,α为循环频率;l为整数;
步骤六、找到步骤五获得的循环谱密度函数Sα(k)在除α=0,k=0以外的幅值最大值的点,并以该点为中心、3dB带宽范围内的点作为估计信号平均功率的估计空间;并估计出该空间的方差σs;
步骤七、根据公式:
估计出步骤六获得的方差σs的信噪比SNR;
步骤八、将步骤七获得的信噪比SNR与预设的判决门限γ进行比较,若SNR≥γ,或者重叠计算循环功率谱密度的次数大于预设值Z,则将此循环谱密度函数Sα(k)作为估计信号有无的谱密度函数,执行步骤十;否则,执行步骤九;
步骤九、从一号寄存器中读取一段长度为D的数据,其中D=N0-P;将该段数据与二号寄存器中的长度为P的重叠部分数据组成一段长度为N0的数据x′(n),并返回执行步骤三;
步骤十、对步骤八获得的循环谱密度函数Sα(k)进行通过公式:
获得循环频率点的能量,并与门限值ξ进行比较,若ε≥η,则有主用户存在,完成认知无线电系统中循环谱的自适应快速感知;否则,则没有主用户存在,完成认知无线电系统中循环谱的自适应快速感知;
其中:门限值ξ是根据公式:
ξ=η×SNR
获得的;式中:η为信噪比系数,η为实数。。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210180859.6A CN102694611B (zh) | 2012-06-04 | 2012-06-04 | 认知无线电系统中循环谱的自适应快速感知方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210180859.6A CN102694611B (zh) | 2012-06-04 | 2012-06-04 | 认知无线电系统中循环谱的自适应快速感知方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102694611A CN102694611A (zh) | 2012-09-26 |
| CN102694611B true CN102694611B (zh) | 2014-06-11 |
Family
ID=46859906
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210180859.6A Expired - Fee Related CN102694611B (zh) | 2012-06-04 | 2012-06-04 | 认知无线电系统中循环谱的自适应快速感知方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102694611B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116054983A (zh) * | 2022-10-21 | 2023-05-02 | 上海微波技术研究所(中国电子科技集团公司第五十研究所) | 基于循环特征的认知网络协作频谱感知方法及系统 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008108795A1 (en) * | 2007-03-08 | 2008-09-12 | Thomson Licensing | Apparatus and method for sensing a multi-carrier signal using cyclostationarity |
| CN101630983A (zh) * | 2009-07-30 | 2010-01-20 | 哈尔滨工业大学 | 认知无线电中利用循环谱统计量的空闲频谱检测方法 |
| CN101834630A (zh) * | 2010-05-11 | 2010-09-15 | 南京邮电大学 | 一种基于能量-循环平稳特征的联合频谱检测方法 |
| CN101977091A (zh) * | 2010-10-19 | 2011-02-16 | 中国人民解放军信息工程大学 | 电磁频谱监测方法及系统 |
-
2012
- 2012-06-04 CN CN201210180859.6A patent/CN102694611B/zh not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008108795A1 (en) * | 2007-03-08 | 2008-09-12 | Thomson Licensing | Apparatus and method for sensing a multi-carrier signal using cyclostationarity |
| CN101630983A (zh) * | 2009-07-30 | 2010-01-20 | 哈尔滨工业大学 | 认知无线电中利用循环谱统计量的空闲频谱检测方法 |
| CN101834630A (zh) * | 2010-05-11 | 2010-09-15 | 南京邮电大学 | 一种基于能量-循环平稳特征的联合频谱检测方法 |
| CN101977091A (zh) * | 2010-10-19 | 2011-02-16 | 中国人民解放军信息工程大学 | 电磁频谱监测方法及系统 |
Non-Patent Citations (2)
| Title |
|---|
| 申达,甘小莺,许苗,钱良.认知无线电能量检测与循环谱检测算法研究.《信息技术》.2010,(第5期),全文. * |
| 谭学治,姜靖,孙洪剑.认知无线电的频谱感知技术研究.《信息安全与通信保密》.2007,(第3期),全文. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102694611A (zh) | 2012-09-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103346845B (zh) | 基于快速傅里叶变换的盲频谱感知方法及装置 | |
| CN102185666B (zh) | 复杂电磁环境下的多模式组合干扰预测系统及方法 | |
| CN104301280B (zh) | 一种基于ofdm的电力线载波周期脉冲噪声检测和抑制方法 | |
| CN103297159A (zh) | 一种频谱感知方法、频谱感知装置 | |
| CN102594472A (zh) | 一种基于小波分解阈值去噪的无线信道测量的方法及其装置 | |
| CN101577564A (zh) | 基于判决门限自适应的信号频谱感知与检测方法 | |
| CN105553896A (zh) | 宽带多频稀疏信号的非均匀采样与重建方法 | |
| CN102291186A (zh) | 一种基于信号到达方向估计的频谱感知方法 | |
| CN105025583A (zh) | 基于能量与协方差检测的分步频谱感知方法 | |
| CN104218973A (zh) | 基于Myriad滤波的跳频信号参数估计方法 | |
| CN103297160A (zh) | 基于归一化特征值的拟合优度检验的频谱感知方法及装置 | |
| CN102594471A (zh) | 一种基于分形盒维数和三阶循环累积量的频谱感知方法 | |
| CN102710349B (zh) | 一种脉冲干扰环境下基于数据挑选的频谱感知方法 | |
| CN103281141A (zh) | 一种盲频谱感知方法及装置 | |
| CN102006609B (zh) | 一种自适应序贯的协同频谱检测方法 | |
| CN103905129B (zh) | 基于谱型分析的信号检测及信号信息判读方法 | |
| CN103780324B (zh) | 动态频谱接入方法 | |
| CN102694611B (zh) | 认知无线电系统中循环谱的自适应快速感知方法 | |
| CN102111228B (zh) | 一种基于循环对称性的认知无线电频谱感知方法 | |
| CN110233687B (zh) | 一种多路数据信号联合判决检测方法 | |
| CN102104560B (zh) | 信道估计方法及装置 | |
| CN107040269A (zh) | 基于方差中值滤波的极/超低频信道大气噪声抑制方法 | |
| CN104363065B (zh) | 基于非高斯性测度的无线通信系统频谱感知方法 | |
| CN103078661A (zh) | 一种基于迭代门限的扩频系统干扰抑制方法 | |
| CN105813089A (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 |
Granted publication date: 20140611 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |