CN106771589A - 一种频谱分析实时瀑布图的实现方法 - Google Patents
一种频谱分析实时瀑布图的实现方法 Download PDFInfo
- Publication number
- CN106771589A CN106771589A CN201611243290.8A CN201611243290A CN106771589A CN 106771589 A CN106771589 A CN 106771589A CN 201611243290 A CN201611243290 A CN 201611243290A CN 106771589 A CN106771589 A CN 106771589A
- Authority
- CN
- China
- Prior art keywords
- data
- frequency spectrum
- spectrum data
- real
- frequency
- 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
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),计算公式为:
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 true CN106771589A (zh) | 2017-05-31 |
CN106771589B 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) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107576848A (zh) * | 2017-09-27 | 2018-01-12 | 中国电子科技集团公司第五十四研究所 | 一种基于频谱分析的模板设定与模板检测方法 |
CN108387779A (zh) * | 2018-01-17 | 2018-08-10 | 成都中安频谱科技有限公司 | 频谱显示方法及装置 |
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)
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 | 中国电子科技集团公司第四十一研究所 | 一种用于实时频谱轨迹处理的装置及方法 |
-
2016
- 2016-12-29 CN CN201611243290.8A patent/CN106771589B/zh active Active
Patent Citations (5)
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)
Title |
---|
白登峰 等: "设备实时监测和故障诊断***中三维瀑布图控件的实现", 《软件导刊》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107576848A (zh) * | 2017-09-27 | 2018-01-12 | 中国电子科技集团公司第五十四研究所 | 一种基于频谱分析的模板设定与模板检测方法 |
CN107576848B (zh) * | 2017-09-27 | 2019-07-26 | 中国电子科技集团公司第五十四研究所 | 一种基于频谱分析的模板设定与模板检测方法 |
CN108387779A (zh) * | 2018-01-17 | 2018-08-10 | 成都中安频谱科技有限公司 | 频谱显示方法及装置 |
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 | 上海卫星工程研究所 | 卫星射频通道频谱监视与记录***及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN106771589B (zh) | 2019-02-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106771589A (zh) | 一种频谱分析实时瀑布图的实现方法 | |
CN105841961A (zh) | 一种基于Morlet小波变换和卷积神经网络的轴承故障诊断方法 | |
CN109445882B (zh) | 一种频谱三维显示装置、方法及计算机可读存储介质 | |
CN109459130B (zh) | 一种基于智能手机的建筑结构舒适度测试*** | |
CN110148425A (zh) | 一种基于完整局部二进制模式的伪装语音检测方法 | |
CN108042132A (zh) | 基于dwt和emd融合csp的脑电特征提取方法 | |
CN109725198B (zh) | 一种频谱仪的余辉生成方法及频谱仪 | |
US20170296081A1 (en) | Frame based spike detection module | |
WO2023138549A1 (zh) | 图像处理方法、装置、电子设备及存储介质 | |
CN106571146A (zh) | 噪音信号确定方法、语音去噪方法及装置 | |
CN107478418A (zh) | 一种旋转机械故障特征自动提取方法 | |
CN108919241A (zh) | 一种基于恒虚警检测的水下信号时频端点参数估计方法 | |
CN109058089A (zh) | 一种基于声发射信号的真空泵过载故障检测的方法 | |
CN110288026A (zh) | 一种基于度量关系图学习的图像分割方法及装置 | |
CN109343043A (zh) | 一种基于选择主成分分析的雷达hrrp目标识别方法 | |
CN103487653B (zh) | 一种用于实时频谱轨迹处理的装置及方法 | |
CN105105757A (zh) | 一种可穿戴人体运动姿态轨迹记录与测评装置 | |
CN104933178B (zh) | 官方网站确定方法和***及官方网站的排序方法 | |
CN103106412B (zh) | 薄片类介质识别方法和识别装置 | |
KR20120094415A (ko) | 파형의 정밀 측정 | |
CN109300484A (zh) | 音频对齐方法、装置、计算机设备以及可读存储介质 | |
CN106488554A (zh) | 一种指纹数据库建立方法和*** | |
CN112241001A (zh) | 雷达人体动作识别方法、装置、电子设备和存储介质 | |
Ma et al. | Complex PQD Classification Using Time–Frequency Analysis and Multiscale Parallel Attention Residual Network | |
CN205163106U (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 |