CN106771589B - 一种频谱分析实时瀑布图的实现方法 - Google Patents

一种频谱分析实时瀑布图的实现方法 Download PDF

Info

Publication number
CN106771589B
CN106771589B CN201611243290.8A CN201611243290A CN106771589B CN 106771589 B CN106771589 B CN 106771589B CN 201611243290 A CN201611243290 A CN 201611243290A CN 106771589 B CN106771589 B CN 106771589B
Authority
CN
China
Prior art keywords
frequency spectrum
data
spectrum data
real
time
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
Application number
CN201611243290.8A
Other languages
English (en)
Other versions
CN106771589A (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.)
CETC 54 Research Institute
Original Assignee
CETC 54 Research Institute
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 CETC 54 Research Institute filed Critical CETC 54 Research Institute
Priority to CN201611243290.8A priority Critical patent/CN106771589B/zh
Publication of CN106771589A publication Critical patent/CN106771589A/zh
Application granted granted Critical
Publication of CN106771589B publication Critical patent/CN106771589B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/16—Spectrum analysis; Fourier analysis

Landscapes

  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
  • Circuits Of Receivers In General (AREA)

Abstract

本发明公开了一种实时瀑布图的实现方法，它涉及了一种在频谱监测领域用于实现实时瀑布图的方法。它由频谱数据生成单元通过DMA模式读取频谱数据，采用乒乓操作方式将频谱数据循环存储在数据区。频谱数据处理单元采用加权因子的方法将每帧频谱数据的相同位置的数据进行加权计算。按照时间顺序将将经过加权因子权衡的频谱数据的幅度根据自定义的颜色‑幅度比例转换成颜色数据并显示在显示屏上。本发明具有能够实时采集、实时分析、实时显示等特点，能够不丢失频谱的重要信息特征，并且有效的避免了硬件之间处理速度的差异造成的丢帧问题。特别适合用于制作频谱监测领域中检测快速跳变信号等瞬时信号的实时频谱仪。

Description

一种频谱分析实时瀑布图的实现方法
技术领域
本发明属于频谱监测领域，具体涉及一种实时瀑布图的实现方法。
背景技术
在频谱分析中，为了方便用户从频谱数据中获取更多消息，一般都会提供多条轨迹，如最大保持、最小保持、0-SPAN模式等。这都是基于幅度和频率、或者幅度和时间的二维关系。时域或者频域分析只能分析一帧长度的信号，特别是在非平稳信号分析中，如果需要查看幅度、频率随时间变化的三维图，对动态信号显示显得非常有用，需要一种新的显示方法瀑布图。由于瀑布图方法数据量大，常见的瀑布图分析会采用抽取非实时处理方法，缺点是不能实时显示频谱信息，会遗漏频谱的重要信息特征。
发明内容
针对以上不足，本发明提供一种可在一个界面同时显示幅度、频率和时间三者关系的实时瀑布图的功能实现方法。本发明通过DMA模式采用乒乓操作方式读取频谱数据，采用加权因子的方法对数据进行抽取加权计算，根据自定义颜色-幅度转换比例将频谱数据转换成颜色数据，采用数据缓存方式实时读取数据在显示屏上显示出来。具有实时性高、不遗漏频谱的重要信息特征等特点。
本发明的目的是通过以下步骤来实现的：
一种频谱分析实时瀑布图的实现方法，其特征在于，包括以下步骤：
步骤1：模数转换器将接收的模拟信号转换成数字信号输出至FPGA，FPGA将数字信号进行下变频转换和快速傅立叶变化，形成频谱数据；
步骤2：频谱数据生成单元通过DMA模式读取FPGA产生的频谱数据，采用乒乓操作方式将频谱数据循环存储在频谱数据生成单元内两个数据区；
步骤3：频谱数据存储单元将频谱数据生成单元两个数据区中的频谱数据采用乒乓操作方式存储到DDR存储器芯片的第一数据区域；
步骤4：频谱数据处理单元设定频谱数据的抽取帧数，采用加权因子的方法将设定抽取帧数中每帧频谱数据相同位置的数据进行加权计算，生成加权因子权衡的频谱数据并存储到DDR存储器芯片的第二数据区域；
步骤5：颜色数据生成单元按照存储顺序依次读取第二数据区域中的加权因子权衡的频谱数据，将加权因子权衡的频谱数据的幅度根据自定义的颜色-幅度比例转换成颜色数据并显示到显示屏上。
其中，步骤4中采用加权因子的方法将每帧频谱数据中相同位置的数据进行加权计算，具体包括以下步骤：
(1)依次计算设定抽取帧数中每一帧频谱数据中每个位置的频谱数据和Ti(k)和加权因子和Si(k)：
Ti(k)＝Ti-1(k)+WiPi(k)，
Si(k)＝Si-1(k)+Wi，Wi＝m+1-i；
式中，Pi(k)为第i帧第k个位置的频谱数据，k为频谱数据的位置，i为当前接收到频谱数据的帧数，Wi为加权因子，m为设定的频谱数据的抽取帧数，T0(k)＝0，S0(k)＝0；
(2)依次将对应位置的频谱数据和Tm(k)和加权因子和Sm(k)相除，得到每个位置的加权因子权衡的频谱数据Q(k)，计算公式为：
本发明与背景技术相比具有如下优点：
1.本发明采用DMA模式乒乓操作方式读取频谱数据，能够实时采集、实时分析、不丢失频谱数据的特点。
2.本发明采用加权因子方式对频谱数据进行实时分析加权运算，能够不丢失频谱的重要信息特征的特点。
3.采用数据缓存方式对频谱数据进行处理，有效的避免了硬件之间处理速度的差异造成的丢帧问题。
4.本发明制造的设备具有集成化程度高、结构简单、可靠性高、人机交互界面友好、使用方便等特点。
附图说明
图1是本发明的瀑布图方法程序流程框图。
图2是本发明的频谱数据处理程序流程框图。
具体实施方式
本发明的实施步骤：如图1所示
步骤1：FPGA控制处理
模数转换器将接收的模拟信号转换成数字信号输出至FPGA，FPGA将数字信号进行下变频转换和快速傅立叶变化，形成频谱数据；
步骤2：频谱数据生成
频谱数据生成单元通过DMA模式读取FPGA产生的频谱数据，采用乒乓操作方式将频谱数据循环存储在频谱数据生成单元内两个数据区；
频谱数据生成单元在DSP芯片的内部缓存区开辟两个数据区A和B，每个区域大小为615点数据，每个数据为16比特。采用乒乓操作方式将频谱数据接收存储。设定当前接收频谱数据到数据区A，同时对数据区B的频谱数据进行存储和抽取计算处理，数据区B处理完成后，等待数据区A接收完成，一旦接收完成，数据区A转为存储、抽取计算处理流程，同时数据区B转为接收新的数据流程，依次类推，如此数据区A、B轮流循环接收数据。
步骤3：频谱数据存储
频谱数据存储单元将频谱数据生成单元两个数据区中的频谱数据采用乒乓操作方式存储到DDR存储器芯片的第一数据区域；
频谱数据存储单元在DDR2芯片存储区中开辟两个数据区域，第一数据区域C作为实时数据存储区，将频谱数据生成单元两个数据区中的频谱数据都实时存储在第一数据区域C中。当数据区A正在接收频谱数据时，将数据区B的频谱数据存储到数据区域C的地址i，当数据区域A的频谱数据接收完成后，将其存储在数据区域C的下一地址i+1，同时数据区B开始接收数据，依次类推，如此数据区A、B轮流循环存储到第一数据区域C中。
步骤4：频谱数据处理
频谱数据处理单元设定频谱数据的抽取帧数，采用加权因子的方法将设定抽取帧数中每帧频谱数据相同位置的数据进行加权计算，生成加权因子权衡的频谱数据并存储到DDR存储器芯片的第二数据区域；
DDR存储器芯片的第二数据区D为瀑布图模式数据的缓冲区，由于显示屏的显示速率小于频谱数据的生成速率，需要将第一数据区域C中频谱数据抽取采样处理并缓存在第二数据区D中，这样避免了硬件之间处理速度的差异造成的丢帧问题，而且方便显示屏显示。
为了不丢失频谱数据的信息，采用加权因子的方法将每帧频谱的相同位置数据进行加权计算。由于不同带宽下频谱生成的速率不同，所以不同带宽下频谱抽取帧数不同。越靠近抽取帧数的频谱数据，其加权因子的权值越大。
如图2所示。采用加权因子的方法将设定抽取帧数中每帧频谱数据相同位置的数据进行加权计算具体包括以下步骤：
(1)依次按下述公式计算设定抽取帧数中每一帧频谱数据中每个位置的频谱数据和Ti(k)和加权因子和Si(k)：
每个位置的频谱数据和Ti(k)为：
Ti(k)＝Ti-1(k)+WiPi(k)
加权因子和Si(k)为：
Si(k)＝Si-1(k)+Wi，其中Wi＝m+1-i；
其中，Pi(k)为第i帧第k个位置的频谱数据，k为频谱数据的位置，n为频谱数据的位置的位置个数，k∈[0,614]，i为当前接收到频谱数据的帧数i∈[1,m]；Wi为加权因子，m为设定的频谱数据的抽取帧数，T0(k)＝0，S0(k)＝0；
每一帧频谱数据的位置从0开始到614，初始i＝1，k＝0，计算第一帧频谱数据中每个位置的频谱数据和Ti(k)和加权因子和Si(k)，然后开始计算下一帧的频谱数据，直至设定抽取帧数m计算完成，得到每个位置的频谱数据和Tm(k)和加权因子和Sm(k)；
(2)k从0开始到614，依次将615个对应的k值相同的频谱数据和Tm(k)和加权因子和Sm(k)相除，得到加权因子权衡的频谱数据Q(k)，计算方法为：
计算完成后按顺序依次存储到DDR存储器芯片的第二数据区域。
步骤5：颜色数据生成并显示
颜色数据生成单元按照存储顺序依次读取第二数据区域中的加权因子权衡的频谱数据，将加权因子权衡的频谱数据的幅度根据自定义的颜色-幅度比例转换成颜色数据并显示到显示屏上。瀑布图的流动方向代表时间，宽度代表频率，各点的颜色代表信号的大小。瀑布图从下往上流，越往下越接近起始时间，越往上越接近当前时间。

Claims (2)

1.一种频谱分析实时瀑布图的实现方法，其特征在于，包括以下步骤：
步骤1：模数转换器将接收的模拟信号转换成数字信号输出至FPGA，FPGA将数字信号进行下变频转换和快速傅立叶变化，形成频谱数据；
步骤2：频谱数据生成单元通过DMA模式读取FPGA产生的频谱数据，采用乒乓操作方式将频谱数据循环存储在频谱数据生成单元内两个数据区；
步骤3：频谱数据存储单元将频谱数据生成单元两个数据区中的频谱数据采用乒乓操作方式存储到DDR存储器芯片的第一数据区域；
步骤4：频谱数据处理单元设定频谱数据的抽取帧数，采用加权因子的方法将设定抽取帧数中每帧频谱数据相同位置的数据进行加权计算，生成加权因子权衡的频谱数据并存储到DDR存储器芯片的第二数据区域；
步骤5：颜色数据生成单元按照存储顺序依次读取第二数据区域中的加权因子权衡的频谱数据，将加权因子权衡的频谱数据的幅度根据自定义的颜色-幅度比例转换成颜色数据并显示到显示屏上。
2.根据权利要求1所述的一种频谱分析实时瀑布图的实现方法，其特征在于，步骤4中采用加权因子的方法将每帧频谱数据中相同位置的数据进行加权计算，具体包括以下步骤：
(1)依次计算设定抽取帧数中每一帧频谱数据中每个位置的频谱数据和Ti(k)与加权因子和Si(k)：
Ti(k)＝Ti-1(k)+WiPi(k)，
Si(k)＝Si-1(k)+Wi，Wi＝m+1-i；
式中，Pi(k)为第i帧第k个位置的频谱数据，k为频谱数据的位置，i为当前接收到频谱数据的帧数，Wi为加权因子，m为设定的频谱数据的抽取帧数，T0(k)＝0，S0(k)＝0；
(2)依次将对应位置的频谱数据和Tm(k)与加权因子和Sm(k)相除，得到每个位置的加权因子权衡的频谱数据Q(k)，计算公式为：
CN201611243290.8A 2016-12-29 2016-12-29 一种频谱分析实时瀑布图的实现方法 Active CN106771589B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611243290.8A CN106771589B (zh) 2016-12-29 2016-12-29 一种频谱分析实时瀑布图的实现方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611243290.8A CN106771589B (zh) 2016-12-29 2016-12-29 一种频谱分析实时瀑布图的实现方法

Publications (2)

Publication Number Publication Date
CN106771589A CN106771589A (zh) 2017-05-31
CN106771589B true CN106771589B (zh) 2019-02-26

Family

ID=58925503

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611243290.8A Active CN106771589B (zh) 2016-12-29 2016-12-29 一种频谱分析实时瀑布图的实现方法

Country Status (1)

Country Link
CN (1) CN106771589B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107576848B (zh) * 2017-09-27 2019-07-26 中国电子科技集团公司第五十四研究所 一种基于频谱分析的模板设定与模板检测方法
CN108387779B (zh) * 2018-01-17 2020-10-23 成都中安频谱科技有限公司 频谱显示方法及装置
CN111026345A (zh) * 2019-12-09 2020-04-17 杭州电子科技大学 一种短波数据显示控制的装置及方法
CN112198365A (zh) * 2020-09-29 2021-01-08 中国船舶重工集团公司第七二四研究所 一种实时频谱监测设备的实现方法
CN112836169A (zh) * 2021-03-04 2021-05-25 天津光电通信技术有限公司 一种基于高速采集的fft频谱数据处理方法及装置
CN114629575A (zh) * 2022-02-10 2022-06-14 上海卫星工程研究所 卫星射频通道频谱监视与记录系统及装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS59116057A (ja) * 1982-12-22 1984-07-04 Victor Co Of Japan Ltd オ−デイオ信号のスペクトル表示装置
CN102455859A (zh) * 2010-11-03 2012-05-16 北京普源精电科技有限公司 一种频谱数据显示方法和装置
CN102928638A (zh) * 2012-11-01 2013-02-13 南京国睿安泰信科技股份有限公司 一种频谱仪瀑布图显示系统及方法
CN103389413A (zh) * 2013-07-09 2013-11-13 中国电子科技集团公司第四十一研究所 一种频谱直方图的实时统计方法
CN103487653A (zh) * 2013-09-18 2014-01-01 中国电子科技集团公司第四十一研究所 一种用于实时频谱轨迹处理的装置及方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS59116057A (ja) * 1982-12-22 1984-07-04 Victor Co Of Japan Ltd オ−デイオ信号のスペクトル表示装置
CN102455859A (zh) * 2010-11-03 2012-05-16 北京普源精电科技有限公司 一种频谱数据显示方法和装置
CN102928638A (zh) * 2012-11-01 2013-02-13 南京国睿安泰信科技股份有限公司 一种频谱仪瀑布图显示系统及方法
CN103389413A (zh) * 2013-07-09 2013-11-13 中国电子科技集团公司第四十一研究所 一种频谱直方图的实时统计方法
CN103487653A (zh) * 2013-09-18 2014-01-01 中国电子科技集团公司第四十一研究所 一种用于实时频谱轨迹处理的装置及方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
设备实时监测和故障诊断系统中三维瀑布图控件的实现;白登峰 等;《软件导刊》;20060915(第17期);83-85

Also Published As

Publication number Publication date
CN106771589A (zh) 2017-05-31

Similar Documents

Publication Publication Date Title
CN106771589B (zh) 一种频谱分析实时瀑布图的实现方法
CN100487720C (zh) 人脸对照装置
CN105841961A (zh) 一种基于Morlet小波变换和卷积神经网络的轴承故障诊断方法
CN109459130B (zh) 一种基于智能手机的建筑结构舒适度测试系统
CN106251249A (zh) 一种建筑施工信息化动态监测系统
CN110148425A (zh) 一种基于完整局部二进制模式的伪装语音检测方法
CN110222769A (zh) 一种基于YOLOV3-tiny的改进目标检测方法
CN110718235A (zh) 异常声音检测的方法、电子设备及存储介质
CN109725198A (zh) 一种频谱仪的余辉生成方法及频谱仪
CN112183643A (zh) 基于声发射的硬岩拉剪破裂识别方法及装置
CN111916089B (zh) 基于声信号特征分析的冰雹检测方法和装置
CN106205637A (zh) 音频信号的噪声检测方法与装置
CN108919241A (zh) 一种基于恒虚警检测的水下信号时频端点参数估计方法
CN108387757A (zh) 用于检测可移动设备的移动状态的方法和装置
CN110288026A (zh) 一种基于度量关系图学习的图像分割方法及装置
CN103884891A (zh) 一种具有高波形刷新率的数字示波器
CN102300116B (zh) 视频分辨率的快速检测方法及装置
CN112881812A (zh) 基于机器学习编码的全闪实时定位方法及装置
CN102750547B (zh) 一种基于压缩传感的水果大小分级方法
CN112241001A (zh) 雷达人体动作识别方法、装置、电子设备和存储介质
CN107991633A (zh) 地磁信号处理方法及装置
CN107894837A (zh) 动态情感分析模型样本处理方法及装置
CN111025100A (zh) 变压器特高频局部放电信号模式识别方法与装置
CN106618552B (zh) 一种检测模型获取方法及装置
CN102890822A (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