CN107576848B - 一种基于频谱分析的模板设定与模板检测方法 - Google Patents
一种基于频谱分析的模板设定与模板检测方法 Download PDFInfo
- Publication number
- CN107576848B CN107576848B CN201710888152.3A CN201710888152A CN107576848B CN 107576848 B CN107576848 B CN 107576848B CN 201710888152 A CN201710888152 A CN 201710888152A CN 107576848 B CN107576848 B CN 107576848B
- Authority
- CN
- China
- Prior art keywords
- template
- maximum
- point
- signal
- reference value
- 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
Landscapes
- Monitoring And Testing Of Transmission In General (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种模板设定、模板检测的实现方法,它涉及了一种在频谱监测领域用于实现模板检测功能的方法。它由FPGA进行下变频转换和快速傅里叶变化生成频谱数据,对频谱数据进行平滑处理生成频谱模板基准值,判断频谱模板基准值的信号类型,根据不同信号类型采用不同方法生成上下频谱模板。对待识别信号进行预处理,将预处理后的待识别信号与上下频谱模板进行比较和匹配,超出上下频谱模板范围判定为异常信号,并给出告警信号,同时对告警时异常频谱进行存储。本发明具有能够实时分析、实时显示等特点,提取信号特征运算量小,实时处理快,检测精度高等优点。特别适合用于卫星通信领域对瞬时异常信号进行快速检测,便于查找和定位问题。
Description
技术领域
本发明属于频谱监测领域,具体涉及一种模板设定、模板检测的实现方法。
背景技术
频谱监测技术是用来监测和管理频段信息的主要手段。随着无线电的发展,频谱资源的利用也在不断增加,致使频谱资源的使用变得日益紧张。大量设备工作于相邻频段,易产生相互干扰,提升频谱资源的利用率,监测异常信号已经成为了面临的重大问题。
对异常信号的识别是频谱监测主要内容之一。常见的信号识别是对信号特征进行提取,利用模式识别的方式进行信号识别。现有的识别信号多为调制信号,没有考虑实际复杂的电磁环境对其信号的影响。且当前频谱监测中信号识别是对信号特征进行提取,提取其信号包络,根据信号特征与包络进行匹配,设置相关系数的计算,根据计算结果来对信号进行识别,识别过程繁琐,且计算量大,在嵌入式设备应用较难,降低了识别速度。
发明内容
针对以上不足,本发明所要解决的技术问题是提供一种基于模板设定、频谱模板检测的功能实现方法。
本发明的目的是通过以下步骤来实现的:
步骤1:将模拟信号进行AD转换、下变频转换和快速傅立叶变换,形成频谱数据;
步骤2:将生成的频谱数据进行平均处理得到平均后的频谱数据,即为频谱数据的模板基准值为频谱数据的位置;
步骤3:对模板基准值进行极大值搜索,得到极大值所在位置j和极大值将极大值点周围(2m-4)个点的功率值与极大值进行比较并计算极大值周边(2m-4)个点的均值得到特征值;分别判断极大值与极大值点周围(2m-4)个点的功率值差值是否大于等于10且特征值是否大于设定阈值,若是则模板基准值为单载波信号,执行步骤4,否则模板基准值为宽带信号,执行步骤5;其中,m≥5且为整数;
步骤4:对极大值位置附近(j-m)点到(j+m)点位置的任一点k,选取(k-2)到(k+2)五个点的功率值进行极大值搜索,并将k点的功率值设为搜索到的极大值,生成新的上频谱模板基准值并对极大值位置附近(k-1)到(k+1)三个点的功率值进行极小值搜索,将k点的功率值设为搜索到的极小值,生成新的下频谱模板基准值设定上模板差值为a,下模板差值为b,则上频谱模板为下频谱模板为执行步骤6;其中,k为(j-m)点到(j+m)范围内任意一整数点;
步骤5:设定上模板差值为a,下模板差值为b,则上频谱模板为下频谱模板为
步骤6:根据预设带宽对待识别信号进行平滑处理,获得预处理的识别信号,对待识别信号和频谱模板的功率值逐点比对,如果待识别信号高于上频谱模板则给出功率超上限告警,如果待识别信号低于下频谱模板则给出功率超下限告警,并记录存储告警发生时的待识别信号。
其中,步骤3具体为:
对模板基准值进行极大值搜索,得到极大值的所在位置j和极大值将极大值点周围(j-m)、(j-(m-1))、…、(j-3)、(j+3)、…、(j+(m-1))、(j+m)即(2m-4)个点的功率值与极大值比较,并计算该(2m-4)个点均值并计算得到特征值;计算公式如下:
计算特征值θ
判断判断极大值与极大值点周围(2m-4)个点的功率值差值是否大于等于10,且特征值是否大于设定阈值,即
其中i∈[j-m,j-3]∪[j+3,j+m],i∈Z+,m≥5且为整数;
θ>θ1 (2)
其中,θ1为设定阈值;
当公式(1)(2)同时满足时,判定模板基准值为单载波信号,执行步骤4,否则为宽带信号,执行步骤5。
本发明与背景技术相比具有如下优点:
1.本发明频谱模板设定方法较之常见方法,其提取信号特征运算量小,提取包络快,对其应用的硬件平台配置要求低,有利于低成本控制。
2.本发明对待测信号实时分析,实时识别速度快,能够快速识别异常偶发信号,识别率高,并且对异常信号能够及时存储,便于分析。
3.本发明制造的设备具有集成化程度高、结构简单、可靠性高、人机交互界面友好、使用方便等特点。
附图说明
图1是本发明的模板设定检测方法流程框图。
图2是单载波信号用本发明的方法产生的上下频谱模板。
图3是宽带信号用本发明的方法产生的上下频谱模板。
图4是待识别信号与模板信号匹配检测正常的原理图。
图5是待识别信号与模板信号匹配检测超上限告警的原理图。
图6是待识别信号与模板信号匹配检测超下限告警的原理图。
具体实施方式
本发明的实施步骤:如图1所示
步骤1;频谱数据生成
系统前端通过FPGA来实现AD数据采集、数字下变频器和快速傅立叶变换。频谱数据生成单元通过DMA模式读取FPGA的频谱数据。在存储区开辟两个数据区A和B,区域大小为615点数据,每个数据为16比特。采用乒乓结构方式将频谱数据接收存储。设当前存储新的数据到数据区A,同时对数据区B的数据进行存储和计算处理,数据区B处理完成后,等待数据区A接收完成,一旦接收完成,数据区A转为存储、计算处理流程,同时数据区B转为接收新的数据流程,如此数据区A、B轮流循环接收数据,处理数据。
步骤2:对生成的频谱数据进行平滑处理。
对生成的频谱数据进行平滑处理,对其进行平均运算,设频谱数据为Pn(i),则平均后的频谱数据,即为频谱数据的模板基准值
其中N为频谱数据的平均次数,i为频谱数据的位置。
步骤3:设定上下模板的差值,生成上下频谱模板。
首先判定模板基准值是否为单载波信号,对模板基准值进行极大值搜索,得到极大值的所在位置j和极大值判定极大值点周围(j-m)、(j-(m-1))、…、(j-3)、(j+3)、…、(j+(m-1))、(j+m)即(2m-4)个点的功率值与极大值比较,其中,m≥5且为整数;
满足以下条件:
其中i∈[j-m,j-3]∪[j+3,j+m],i∈Z+,10≥m≥5;
且计算极大值周边(2m-4)各点的均值
则选取特征值θ
设置一个阈值θ1,判断
θ>θ1 (2)
当公式(1)(2)同时满足时,判定模板基准值为单载波信号。如果不满足上述条件为宽带信号。
若为单载波信号,对极大值位置附近(j-m)点到(j+m)点位置的模板基准值进行处理,得到新的模板基准值。当处理其中某个点k时,对其(k-2)到(k+2)五个点的功率值进行极大值搜索,得到的极大值为k点的功率值,依次类推,当所有点处理完成后,生成新的上频谱模板基准值同理,当处理其中某个点k时,对其(k-1)到(k+1)三个点的功率值进行极小值搜索,得到的极小值为k点的功率值,当所有点处理完成后,生成新的下频谱模板基准值设定上模板差值为a,下模板差值为b,则上频谱模板为 下频谱模板为其中,k为(j-m)点到(j+m)范围内任意一整数点。单载波信号生成的上下频谱模板如图2所示。
若为宽带信号,则设定上模板差值为a,下模板差值为b,则上频谱模板为下频谱模板为宽带信号生成的上下频谱模板如图3所示。
频谱模板可以根据频谱的特征(如频率、SPAN)和类型形成模板库。在模板检测匹配过程中,可以直接调用模板库的模板进行匹配。
步骤4:模板检测匹配。
根据预设频率带宽对待识别信号进行平滑处理,获得预处理的识别信号,将预处理后的待识别信号与频谱模板进行比较,判断待识别信号是否正常。通过对待识别信号和模板信号的功率值逐点比对,如果高于上频谱模板则给出超上限告警,如果低于下频谱模板则给出超下限告警。并记录存储告警发生时的待识别信号。如图4所示,正常信号的频谱位于上频谱模板与下频谱模板之间。待识别信号功率过高时,如图5所示,频谱信号超过上频率模板,监测到后发出超上限告警。待识别信号功率低于下频谱模板时时,如图6所示,监测到后发出超下限告警。
Claims (2)
1.一种基于频谱分析的模板设定与模板检测方法,其特征在于,包括以下步骤:
步骤1:将模拟信号进行AD转换、下变频转换和快速傅立叶变换,形成频谱数据;
步骤2:将生成的频谱数据进行平均处理得到平均后的频谱数据,即为频谱数据的模板基准值i为频谱数据的位置;
步骤3:对模板基准值进行极大值搜索,得到极大值所在位置j和极大值将极大值点周围2m-4个点的功率值与极大值进行比较并计算极大值周边2m-4个点的均值根据极大值和均值得到特征值;分别判断极大值与极大值点周围2m-4个点的功率值差值是否大于等于10且特征值是否大于设定阈值,若是则模板基准值为单载波信号,执行步骤4,否则模板基准值为宽带信号,执行步骤5;其中,m≥5且为整数;
步骤4:对极大值位置附近j-m点到j+m点位置的任一点k,选取k-2到k+2五个点的功率值进行极大值搜索,并将k点的功率值设为搜索到的极大值,生成新的上频谱模板基准值并对极大值位置附近k-1到k+1三个点的功率值进行极小值搜索,将k点的功率值设为搜索到的极小值,生成新的下频谱模板基准值设定上模板差值为a,下模板差值为b,则上频谱模板为下频谱模板为执行步骤6;其中,k为j-m点到j+m范围内任意一整数点;
步骤5:设定上模板差值为a,下模板差值为b,则上频谱模板为下频谱模板为
步骤6:根据预设带宽对待识别信号进行平滑处理,获得预处理的识别信号,对待识别信号和频谱模板的功率值逐点比对,如果待识别信号高于上频谱模板则给出功率超上限告警,如果待识别信号低于下频谱模板则给出功率超下限告警,并记录存储告警发生时的待识别信号。
2.根据权利要求1所述的一种基于频谱分析的模板设定与模板检测方法,其特征在于,步骤3具体为:
对模板基准值进行极大值搜索,得到极大值的所在位置j和极大值将极大值点周围j-m、j-(m-1)、...、j-3、j+3、...、j+(m-1)、j+m即2m-4个点的功率值与极大值比较,并计算该2m-4个点均值并计算得到特征值;计算公式如下:
计算特征值θ
判断极大值与极大值点周围2m-4个点的功率值差值是否大于等于10,且特征值是否大于设定阈值,即
其中i∈[j-m,j-3]∪[j+3,j+m],i∈Z+,m≥5且为整数;
θ>θ1 (2)
其中,θ1为设定阈值;
当公式(1)(2)同时满足时,判定模板基准值为单载波信号,执行步骤4,否则为宽带信号,执行步骤5。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710888152.3A CN107576848B (zh) | 2017-09-27 | 2017-09-27 | 一种基于频谱分析的模板设定与模板检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710888152.3A CN107576848B (zh) | 2017-09-27 | 2017-09-27 | 一种基于频谱分析的模板设定与模板检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107576848A CN107576848A (zh) | 2018-01-12 |
CN107576848B true CN107576848B (zh) | 2019-07-26 |
Family
ID=61038648
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710888152.3A Active CN107576848B (zh) | 2017-09-27 | 2017-09-27 | 一种基于频谱分析的模板设定与模板检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107576848B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108414832B (zh) * | 2018-02-09 | 2020-11-03 | 中国人民解放军31007部队 | 一种基于图像处理的异常信号频谱参数检测方法 |
CN116015497A (zh) * | 2022-12-27 | 2023-04-25 | 成都爱科特科技发展有限公司 | 一种信号频谱模板匹配方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1848240A (zh) * | 2005-04-12 | 2006-10-18 | 佳能株式会社 | 基于离散对数傅立叶变换的基音检测方法、设备和介质 |
CN1941762A (zh) * | 2005-09-27 | 2007-04-04 | 华为技术有限公司 | 一种分配正交多路频分复用系统中保护子载波的方法 |
CN101173870A (zh) * | 2007-09-27 | 2008-05-07 | 上海交通大学 | 复杂噪声基底下短时特征声信号时频域辨识与检测方法 |
CN102466764A (zh) * | 2010-11-03 | 2012-05-23 | 北京普源精电科技有限公司 | 一种频谱超限测量模板的生成方法和装置 |
JP2014016552A (ja) * | 2012-07-10 | 2014-01-30 | Pioneer Electronic Corp | 音声信号処理方法、音声信号処理装置およびプログラム |
CN104640139A (zh) * | 2015-02-12 | 2015-05-20 | 成都大公博创信息技术有限公司 | 一种高铁通信的干扰检测方法 |
JP2017032488A (ja) * | 2015-08-05 | 2017-02-09 | 株式会社リコー | 情報処理装置、診断方法、およびプログラム |
CN106646012A (zh) * | 2016-09-09 | 2017-05-10 | 中国舰船研究设计中心 | 一种用频装备发射谱包络提取方法 |
CN106792706A (zh) * | 2016-12-09 | 2017-05-31 | 中国人民解放军国防科学技术大学 | 一种基于无线网信号频谱非平稳特性的非法入侵检测方法 |
CN106771589A (zh) * | 2016-12-29 | 2017-05-31 | 中国电子科技集团公司第五十四研究所 | 一种频谱分析实时瀑布图的实现方法 |
-
2017
- 2017-09-27 CN CN201710888152.3A patent/CN107576848B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1848240A (zh) * | 2005-04-12 | 2006-10-18 | 佳能株式会社 | 基于离散对数傅立叶变换的基音检测方法、设备和介质 |
CN1941762A (zh) * | 2005-09-27 | 2007-04-04 | 华为技术有限公司 | 一种分配正交多路频分复用系统中保护子载波的方法 |
CN101173870A (zh) * | 2007-09-27 | 2008-05-07 | 上海交通大学 | 复杂噪声基底下短时特征声信号时频域辨识与检测方法 |
CN102466764A (zh) * | 2010-11-03 | 2012-05-23 | 北京普源精电科技有限公司 | 一种频谱超限测量模板的生成方法和装置 |
JP2014016552A (ja) * | 2012-07-10 | 2014-01-30 | Pioneer Electronic Corp | 音声信号処理方法、音声信号処理装置およびプログラム |
CN104640139A (zh) * | 2015-02-12 | 2015-05-20 | 成都大公博创信息技术有限公司 | 一种高铁通信的干扰检测方法 |
JP2017032488A (ja) * | 2015-08-05 | 2017-02-09 | 株式会社リコー | 情報処理装置、診断方法、およびプログラム |
CN106646012A (zh) * | 2016-09-09 | 2017-05-10 | 中国舰船研究设计中心 | 一种用频装备发射谱包络提取方法 |
CN106792706A (zh) * | 2016-12-09 | 2017-05-31 | 中国人民解放军国防科学技术大学 | 一种基于无线网信号频谱非平稳特性的非法入侵检测方法 |
CN106771589A (zh) * | 2016-12-29 | 2017-05-31 | 中国电子科技集团公司第五十四研究所 | 一种频谱分析实时瀑布图的实现方法 |
Non-Patent Citations (2)
Title |
---|
TD-LTE-Advanced终端测试仪频谱发射模板测试方法与实现;王先鹏 等;《安徽电子信息职业技术学院学报》;20161231;第37-40页 |
基于FPGA的数字下变频器的设计;张红涛 等;《现代电子技术》;20140801;第73-75页 |
Also Published As
Publication number | Publication date |
---|---|
CN107576848A (zh) | 2018-01-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105425123A (zh) | 一种紫外成像与红外成像协同检测电力设备故障的方法及系统 | |
CN108830162B (zh) | 无线电频谱监测数据中的时序模式序列提取方法及存储方法 | |
CN107576848B (zh) | 一种基于频谱分析的模板设定与模板检测方法 | |
CN109470957B (zh) | 一种基于弗雷歇距离的非侵入式电器识别方法 | |
CN111626360B (zh) | 用于检测锅炉故障类型的方法、装置、设备和存储介质 | |
CN108489996A (zh) | 一种绝缘子的缺陷检测方法、系统及终端设备 | |
CN109039503A (zh) | 一种频谱感知方法、装置、设备及计算机可读存储介质 | |
CN110298369A (zh) | 一种电力系统不良数据的辨识方法及系统 | |
Honghao et al. | Spectrum anomalies autonomous detection in cognitive radio using hidden markov models | |
CN110069961A (zh) | 一种物体检测方法和装置 | |
CN117556377B (zh) | 用于站所自动化终端的多源数据融合处理方法 | |
CN116090709B (zh) | 基于金属加工监控的偷排预警分析方法及系统 | |
CN105974223B (zh) | 一种用于在线检测用电设备工作状态的方法及系统 | |
CN110555619A (zh) | 一种基于智能配电网的供电能力评估方法 | |
CN114236234A (zh) | 一种基于基波和谐波混合判据的电器特征辨识方法 | |
CN115800552A (zh) | 一种超级电容运行电力调频用智能调控系统及方法 | |
CN109146265A (zh) | 基于动态规划时间序列相似性的三相不平衡源溯源方法 | |
CN117851907B (zh) | 一种基于物联网技术的水闸渗流监测方法 | |
CN110967585A (zh) | 恶性负载的识别方法及装置 | |
CN116861316B (zh) | 一种电器监测方法及装置 | |
CN102932077A (zh) | 一种频谱占用度的测量方法 | |
CN115912359B (zh) | 基于大数据的数字化安全隐患识别排查治理方法 | |
CN116502149A (zh) | 基于电流特征传导的低压配电网户变关系识别方法及系统 | |
CN110687379A (zh) | 一种拓扑可配置的非侵入式楼宇用电设备监测分析系统 | |
CN107135126B (zh) | 基于子流分形指数的流量在线识别方法 |
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 |