CN108776264B - 数字示波器的fft分析装置 - Google Patents

数字示波器的fft分析装置 Download PDF

Info

Publication number
CN108776264B
CN108776264B CN201810835993.2A CN201810835993A CN108776264B CN 108776264 B CN108776264 B CN 108776264B CN 201810835993 A CN201810835993 A CN 201810835993A CN 108776264 B CN108776264 B CN 108776264B
Authority
CN
China
Prior art keywords
frequency
module
data
signal
channel
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
CN201810835993.2A
Other languages
English (en)
Other versions
CN108776264A (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201810835993.2A priority Critical patent/CN108776264B/zh
Publication of CN108776264A publication Critical patent/CN108776264A/zh
Application granted granted Critical
Publication of CN108776264B publication Critical patent/CN108776264B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/16Spectrum analysis; Fourier analysis

Landscapes

  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Measuring Frequencies, Analyzing Spectra (AREA)

Abstract

本发明公开了一种数字示波器的FFT分析装置,对多通道ADC采集的并行数据进行缓存,然后测量得到每个通道的信号频率,根据抽取系数计算得到倍频参数,对每个通道数据的频率信号进行倍频得到倍频信号,同步抽取模块从缓存的采集数据中读取各个通道的并行数据,根据倍频信号进行数据抽取,对单个通道的并行抽点数据进行平均得到每个通道的平均数据,对该平均数据进行FFT运算和缓存,由上位机读取进行显示。本发明通过对多通道多路并行ADC采集信号进行同步抽取后取平均值,降低通道噪声,使输入信号的频谱更利于观测。

Description

数字示波器的FFT分析装置
技术领域
本发明属于数字示波器技术领域,更为具体地讲,涉及一种数字示波器的FFT分析装置。
背景技术
在传统数字示波器中,FFT分析方法是:将输入的模拟量电流电压通过AD采样环节变换成离散的数字量,然后进行快速傅立叶变换,计算获得基波和各次谐波的复数值。现有数字示波器采样率高,传统的FFT用软件实现,串行处理速度慢,不能满足高速数字示波器的实时FFT需求;传统示波器FFT分析数据点抽取信号不能与原始采集信号同步,并且由于不能满足整周期固定点采样导致存在截断误差,因此实际精度不高;传统FFT分析方法噪声频谱分量过大,影响观察被测信号的频谱,具体分析如下:
(1)传统FFT分析方法耗时较长,导致波形捕获率降低:
传统FFT分析是由软件实现的,高速数字示波器传输速率快且为多通道同步采样,由于软件是串行处理导致处理速度慢,信号处理需要消耗更多的时间,由数据处理时间与波形捕获率成反比,数字示波器的波形捕获率将降低。
(2)传统FFT分析方法精度较低:
传统数字示波器中基于软件的FFT方法,普遍采用根据适当的抽点系数抽点后将数据送入软件中,软件根据相关的FFT算法计算被测信号的频谱。该抽点算法无法实现被测信号基波频率f0的2q,q∈[1,2,3…]倍抽点,并且抽点使能信号与被测信号基波频率f0相位随机,从而导致数据截断误差、相位误差降低FFT的分析精度。
(3)传统FFT分析方法不能有效降低通道的底噪:
传统FFT分析无法抑制通道噪声,导致低信噪比的被测信号的频谱容易被噪声的频谱淹没,导致分析结果不精确,不利于复现原始信号的频谱,给原始采集信号频谱的观察带来很大的干扰。
发明内容
本发明的目的在于克服现有技术的不足,提供一种数字示波器的FFT分析装置,通过对多通道多路并行ADC采集信号进行同步固定点整数周期抽取后取平均值,提高了FFT分析的精度,降低通道噪声,使输入信号的频谱更利于观测。
为了实现上述发明目的,本发明数字示波器的FFT分析装置包括时钟模块、分频器模块、比较器模块、信号采集模块、相位校正模块、数据缓存模块、频率测量模块、倍频模块、同步抽取模块、FFT分析模块和上位机,其中相位校正模块、数据缓存模块、频率测量模块、倍频模块、同步抽取模块和FFT分析模块在FPGA中实现。
时钟模块用于生成N路采集时钟ad_clkn、1路频率测量参考时钟clks和1路倍频参考时钟clk0,其中n=1,2,…,N,N表示采集模块的通道数量,将N路采集时钟ad_clkn分别输入信号采集模块的每片ADC进行时钟控制,以保证N通道信号同步采集,频率测量参考时钟clks和倍频参考时钟clk0分别发送给频率测量模块和倍频模块;
比较器模块用于接收输入的N个通道模拟信号进行阈值比较,得到与N路模拟信号频率相同、占空比为50%的方波信号Wn,将N路方波信号Wn发送至分频器模块;
分频器模块用于将接收到的N路方波信号Wn进行固定K分频,K的值需要令分频后信号的频率符合FPGA信号输入的频率范围,将分频后得到的方波信号
Figure GDA0002255825570000021
发送至FPGA中的频率测量模块用于频率测量;
信号采集模块用于对输入的N个通道模拟信号进行采集,信号采集模块包括N个ADC模块,每个ADC模块包含M片交替采集的ADC子模块,每个通道的模拟信号由对应ADC按照接收到的采集时钟ad_clkn进行信号采集,每个ADC模块并行输出M路数据和1路采集时钟ad_clkn至数据缓存模块,记第n个通道所采集的数据datan为Un1,Un2…,UnM,其中Unm表示第n个通道中第m路数据;
相位校正模块用于对信号采集模块的每片ADC进行相位校正,使每个ADC模块输出相位相同;
数据缓存模块包括N个FIFO,第n个FIFO对第n个通道的数据datan进行降速缓存,根据ad_clkn实现N个通道的数据datan同步输出;
频率测量模块用于从时钟控制模块接收频率测量参考时钟clks,从分频器模块接收方波信号
Figure GDA0002255825570000031
进行频率测量,得到N个通道数据对应的K分频后的频率
Figure GDA0002255825570000032
然后计算得到N个通道数据对应的频率
Figure GDA0002255825570000033
将N个频率
Figure GDA0002255825570000034
发送给同步抽取模块,将方波信号
Figure GDA0002255825570000035
转发至倍频模块;
倍频模块用于接收倍频参考时钟clk0和N个方波信号
Figure GDA0002255825570000036
分别对N个方波信号
Figure GDA0002255825570000037
进行K×Dn倍倍频得到倍频信号
Figure GDA0002255825570000038
将N个倍频信号
Figure GDA0002255825570000039
和对应的分频系数An输出至同步抽取模块,Dn表示从上位机接收的对第n个通道的第一倍频系数;
同步抽取模块用于从上位机接收FFT分析的信号周期数P,从倍频模块接收N个倍频信号
Figure GDA00022558255700000310
然后从数据缓存模块并行读取N个通道的数据datan,分别对每个通道的M路数据Unm根据对应的倍频信号
Figure GDA00022558255700000311
进行抽取,抽取方法如下:
设置每个通道的M路数据Unm的P个周期对应的闸门时间tn
Figure GDA00022558255700000312
倍频信号
Figure GDA00022558255700000313
与闸门时间tn的闸门信号边沿对齐,在闸门时间tn内对数据Unm根据倍频信号
Figure GDA00022558255700000314
进行抽取,将倍频信号
Figure GDA00022558255700000315
每两个上升沿之间的数据作为一个数据序列提取出来,在闸门时间tn共计提取出Dsum个数据序列,Dsum=P×K×Dn,记数据Unm中提取出的第d个数据序列为datanm[d],d=1,2,…,Dsum,记数据序列datanm[d]对应的分频系数为An[d];
将每个通道提取出的数据按照以下计算公式进行平均,得到第n个通道的平均数据data′n
Figure GDA00022558255700000316
其中data′n[d]表示平均数据data′n的第d个数据点,datanm[d][i]表示数据序列datanm[d]中的第i个数据,i=1,2,…,An
将得到的每个通道的平均数据data′n发送给FFT分析模块;
FFT分析模块包括N组FFT IP核和平均数据FIFO缓存模块,第n个FFT IP核用于接收第n个通道的平均数据data′n,对其进行FFT运算,得到FFT运算结果Fn;第n个平均数据FIFO缓存模块用于对第n个通道的平均数据data′n进行缓存,缓存之后可以发送至上位机进行波形显示与FFT运算结果进行对比;
上位机接收用户设置的FFT分析的信号周期数P和第一倍频系数Dn,将FFT分析的信号周期数P发送给同步抽取模块,将第一倍频系数Dn发送给倍频模块;从FFT分析模块中读取每个通道的FFT运算结果Fn和缓存的平均数据data′n进行频谱显示。
本发明示波器的FFT分析装置,对多通道ADC采集的并行数据进行缓存,然后测量得到每个通道分频后的信号频率,对每个通道分频后的频率信号进行倍频得到倍频信号,对通道数据的频率进行还原,同步抽取模块从缓存的采集数据中读取各个通道的并行数据,根据倍频信号进行数据抽取,对单个通道的并行抽点数据进行平均得到每个通道的平均数据,对该平均数据进行FFT运算和缓存,由上位机读取进行对应通道的频谱显示。
本发明具有以下有益效果:
1)本发明是由FPGA实现,FPGA对数据处理速度比软件数据处理速度快,因此能够满足高速数字示波器实时的FFT功能的要求,同时极大提高了波形捕获率;
2)本发明可以实现整数周期固定点抽取,倍频模块产生的倍频信号为非等间隔信号,使用此倍频信号进行同步抽取可以实现对采集信号整数周期的固定点抽取,并且抽点信号能够与被测信号同步,减小了抽点时产生的截断误差,极大提高了FFT分析的精度,使得FFT分析后的频谱更接近被测信号的真实频谱;
3)本发明通过对多路并行ADC采集信号抽取后的数据取平均值能够有效抑制通道噪声和毛刺信号,提高了FFT分析的精度,更有利于FFT分析后得到信号源的真实频谱信息。
附图说明
图1是本发明数字示波器的FFT分析装置的具体实施方式结构图;
图2是本实施例中倍频模块结构图;
图3是本实施例中现有示波器FFT分析装置的输出频谱图;
图4是本实施例中本发明的输出频谱图。
具体实施方式
下面结合附图对本发明的具体实施方式进行描述,以便本领域的技术人员更好地理解本发明。需要特别提醒注意的是,在以下的描述中,当已知功能和设计的详细描述也许会淡化本发明的主要内容时,这些描述在这里将被忽略。
实施例
图1是本发明数字示波器的FFT分析装置的具体实施方式结构图。如图1所示,本发明数字示波器的FFT分析装置包括时钟模块1、分频器模块2、比较器模块3、信号采集模块4、相位校正模块5、数据缓存模块6、频率测量模块7、倍频模块8、同步抽取模块9、FFT分析模块10和上位机11,其中相位校正模块5、数据缓存模块6、频率测量模块7、倍频模块8、同步抽取模块9和FFT分析模块10在FPGA(Field-Programmable Gate Array,现场可编程门阵列)中实现,接下来分别对每个模块进行详细说明。
时钟模块1用于生成N路采集时钟ad_clkn、1路频率测量参考时钟clks和1路倍频参考时钟clk0,其中n=1,2,…,N,N表示信号采集模块2的通道数量,将N路采集时钟ad_clkn分别输入信号采集模块2的每片ADC进行时钟控制,以保证N通道信号同步采集,频率测量参考时钟clks和倍频参考时钟clk0分别发送给频率测量模块5和倍频模块6。本实施例中采集通道数量N=2,设置采集时钟ad_clkn频率为2.5GHz,频率测量参考时钟clks频率为100MHz,倍频参考时钟clk0频率为10MHz。
比较器模块3用于输入的N个通道模拟信号进行阈值比较,输出与比较器模块输入信号频率相同占空比为50%的方波信号Wn,将方波信号Wn发送至分频器模块。
分频器模块2用于将接收到的N路方波信号Wn进行固定K分频,K的值需要令分频后信号的频率符合FPGA信号输入的频率范围,将分频器模块输出的方波信号
Figure GDA0002255825570000051
发送至FPGA中的频率测量模块用于频率测量。一般来说K设置为奇数,本实施例中K=9。
信号采集模块4用于对输入的N个通道模拟信号进行采集,信号采集模块2包括N个ADC模块,每个ADC模块包含M片交替采集的ADC子模块,每个通道的模拟信号由对应ADC按照接收到的采集时钟ad_clkn进行信号采集,每个ADC模块并行输出M路数据和1路采集时钟ad_clkn至数据缓存模块4,记第n个通道所采集的M路数据datan为Un1,Un2…,UnM,其中Unm表示第n个通道中第m路数据。
本实施例中信号采集模块2采用2个数据采集通道包含两个ADC模块,本实施例中采用的ADC芯片有3种采集模式,单通道1.25Gsps采样率采集(使用1个ADC子模块)、双通道2.5Gsps采样率交替采集(使用2个ADC子模块)、4通道5Gsps采样率交替采集(使用4个ADC子模块),此处将ADC配置为4通道5Gsps采样率交替采集模式,即每个通道有4路并行数据,在本实施例中采用FFT分析装置测量开关电源的频谱。
相位校正模块5用于对信号采集模块的每片ADC进行相位校正,使每个ADC模块输出相位相同。本实施例中相位校正模块3的相位校正方法如下:
1)N个ADC模块同时对一个正弦信号源进行采集。
2)判断N个ADC模块的输出数据是否相同,如果相同则记载当前对每个ADC模块的相位调整值,相位校正完全,否则进入步骤3)。
3)以某一个ADC模块(一般选择第1个)的输出数据为基准,控制另外N-1个ADC模块的相位调整值,调节输出信号相位,返回步骤2)。
本实施例中N=2。
数据缓存模块6包括N个FIFO,第n个FIFO对第n个通道的数据datan进行降速缓存,根据ad_clkn实现N个通道的数据datan同步输出。本实施例中包括2个FIFO,由于每个数据datan为4路并行数据,那么经过FIFO缓存后仍为4路并行数据。
频率测量模块7用于从时钟控制模块1接收频率测量参考时钟clks,从分频器模块接收方波信号
Figure GDA0002255825570000061
进行频率测量,得到N个通道数据对应的K分频后的频率
Figure GDA0002255825570000062
然后计算得到N个通道数据对应的频率
Figure GDA0002255825570000063
将N个频率
Figure GDA0002255825570000064
发送给同步抽取模块9,将方波信号
Figure GDA0002255825570000065
转发至倍频模块8。本实施例中采用的频率测量方法如下:
由频率测量参考时钟clks计数产生一个预设闸门信号,然后用方波信号
Figure GDA0002255825570000066
同步预设闸门信号产生实际闸门信号,频率测量参考时钟clks和方波信号
Figure GDA0002255825570000071
在实际闸门内计数,记频率测量参考时钟clks的频率为fs,实际闸门信号的时间内对标准时钟计数为Nums,对方波信号
Figure GDA0002255825570000072
的计数值为
Figure GDA0002255825570000073
则方波信号频率
Figure GDA0002255825570000074
的计算公式如下:
Figure GDA0002255825570000075
倍频模块8用于接收倍频参考时钟clk0和N个方波信号
Figure GDA0002255825570000076
分别对N个方波信号
Figure GDA0002255825570000077
进行K×Dn倍倍频得到倍频信号
Figure GDA0002255825570000078
将N个倍频信号
Figure GDA0002255825570000079
和对应的分频系数An输出至同步抽取模块9,其中Dn表示从上位机11接收的对第n个通道的第一倍频系数,即每个基波周期抽取的数据点的个数。
图2是本实施例中倍频模块结构图。如图2所示,本实施例中倍频模块8包括N组周期计数器81和小数分频器82,其中第n个周期计数器81用于采用倍频参考时钟clk0对方波信号
Figure GDA00022558255700000710
进行周期计数,记方波信号
Figure GDA00022558255700000711
一个周期内对倍频参考时钟clk0的计数值为Cn,即
Figure GDA00022558255700000712
f0表示倍频参考时钟clk0的频率。周期计数器81将得到的计数值Cn发送给第n个小数分频器。
每个小数分频器82包括除法器821、加法器822、比较器823和分频器824,其中第n个小数分频器82中的除法器821用于接收计数值Cn,计算Cn/(KDn),记商为Sn、余数为Yn,将商Sn发送给分频器824,将余数Yn发送给加法器822;
加法器822用于计算和值Sumn=Sumn+Yn,Sumn的初始值为0,将得到的和值Sumn发送给比较器823;
比较器823用于将和值Sumn与第二倍频系数KDn进行比较,当Sumn>KDn时,令Sumn=Sumn-KDn,返回给加法器822,且令标识flagn=1,当Sumn≤KDn时,令标识flagn=0;比较器823将标识flagn发送给分频器824。
分频器824用于接收商Sn和标识flagn,根据标识flagn确定分频系数,当flagn=1时分频系数为Sn-1,当flagn=0时分频系数为Sn,根据分频系数对倍频参考时钟clk0进行分频,得到倍频信号
Figure GDA00022558255700000713
和分频系数An,可见,分频系数An可以采用如下表达式表示:
Figure GDA00022558255700000714
由于
Figure GDA00022558255700000715
倍频信号的频率
Figure GDA00022558255700000716
由此可以实现将方波信号
Figure GDA0002255825570000081
进行KDn倍倍频。根据分频器824的操作可知,本实施例中倍频模块8输出的倍频信号为非等间隔周期信号,同步抽取模块9基于此倍频信号进行抽取。
同步抽取模块9用于从上位机11接收FFT分析的信号周期数P,从倍频模块8接收N个倍频信号
Figure GDA0002255825570000082
然后从数据缓存模块6并行读取N个通道的数据datan,分别对每个通道的M路数据Unm根据对应的倍频信号
Figure GDA0002255825570000083
进行抽取,抽取方法如下:
设置每个通道的M路数据Unm的P个周期对应的闸门时间tn
Figure GDA0002255825570000084
倍频信号
Figure GDA0002255825570000085
与闸门时间tn的闸门信号边沿对齐,在闸门时间tn内对数据Unm根据倍频信号
Figure GDA0002255825570000086
进行抽取,将倍频信号
Figure GDA0002255825570000087
每两个上升沿之间的数据作为一个数据序列提取出来。根据倍频信号
Figure GDA0002255825570000088
的生成过程可知,对于同一通道的每路数据,在闸门时间tn共计提取出Dsum个数据序列,Dsum=P×Dn,记数据Unm中提取出的第d个数据序列为datanm[d],d=1,2,…,Dsum。根据倍频信号
Figure GDA0002255825570000089
的生成过程和抽点过程可知,在两个上升沿之间的数据点数即为其分频系数大小,因此记数据序列datanm[d]对应的分频系数为An[d],An[d]即为数据序列datanm[d]中的数据点数。
将每个通道提取出的数据按照以下计算公式进行平均,得到第n个通道的平均数据data′n
Figure GDA00022558255700000810
其中data′n[d]表示平均数据data′n的第d个数据点,datanm[d][i]表示数据序列datanm[d]中的第i个数据,i=1,2,…,An
将得到的每个通道的平均数据data′n发送给FFT分析模块。将得到的每个通道的平均数据data′n发送给FFT分析模块10。
FFT分析模块10包括N组FFT IP核101和平均数据FIFO缓存模块102,第n个FFT IP核101用于接收第n个通道的平均数据data′n,对其进行FFT运算,得到FFT运算结果Fn;第n个平均数据FIFO缓存模块102用于对第n个通道的平均数据data′n进行缓存,缓存之后可以发送至上位机进行波形显示与FFT运算结果进行对比。
本实施例中信号周期数P=8,第一倍频系数Dn=1024,因此最终抽取得到的数据点数Dsum=8192。
上位机11将FFT分析的信号周期数P和数据点数Q发送给同步抽取模块9,将抽取系数Dn发送给倍频模块8;从FFT分析模块10中读取每个通道的FFT运算结果Fn和缓存的平均数据data′n并进行频谱显示。
为了更好地说明本发明的技术效果,采用matlab软件对现有示波器FFT分析装置和本发明进行仿真对比。本次仿真对比的输入信号为:
U=1+sin(2π×104×t)+0.7×sin(2π×2×104×t)+0.6×sin(2π×3×104×t)
+0.5×sin(2π×4×104×t)+0.4×sin(2π×5×104×t)+0.3×sin(2π×6×104×t)
在信号U中叠加信噪比为1的高斯噪声。
图3是本实施例中现有示波器FFT分析装置的输出频谱图。图4是本实施例中本发明的输出频谱图。对比图3和图4可知,本发明所得到的输出频谱中信号频率更加精确,并且对噪声有很好的抑制作用,使输入信号的频谱更利于观测。
尽管上面对本发明说明性的具体实施方式进行了描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。

Claims (4)

1.一种数字示波器的FFT分析装置,其特征在于,包括时钟模块、分频器模块、比较器模块、信号采集模块、相位校正模块、数据缓存模块、频率测量模块、倍频模块、同步抽取模块、FFT分析模块和上位机,其中相位校正模块、数据缓存模块、频率测量模块、倍频模块、同步抽取模块和FFT分析模块在FPGA中实现;
时钟模块用于生成N路采集时钟ad_clkn、1路频率测量参考时钟clks和1路倍频参考时钟clk0,其中n=1,2,…,N,N表示采集模块的通道数量,将N路采集时钟ad_clkn分别输入信号采集模块的每片ADC进行时钟控制,以保证N通道信号同步采集,频率测量参考时钟clks和倍频参考时钟clk0分别发送给频率测量模块和倍频模块;
比较器模块用于接收输入的N个通道模拟信号进行阈值比较,得到与N路模拟信号频率相同、占空比为50%的方波信号Wn,将N路方波信号Wn发送至分频器模块;
分频器模块用于将接收到的N路方波信号Wn进行固定K分频,K的值需要令分频后信号的频率符合FPGA信号输入的频率范围,将分频后得到的方波信号
Figure FDA0002255825560000015
发送至FPGA中的频率测量模块用于频率测量;
信号采集模块用于对输入的N个通道模拟信号进行采集,信号采集模块包括N个ADC模块,每个ADC模块包含M片交替采集的ADC子模块,每个通道的模拟信号由对应ADC按照接收到的采集时钟ad_clkn进行信号采集,每个ADC模块并行输出M路数据和1路采集时钟ad_clkn至数据缓存模块,记第n个通道所采集的数据datan为Un1,Un2…,UnM,其中Unm表示第n个通道中第m路数据;
相位校正模块用于对信号采集模块的每片ADC进行相位校正,使每个ADC模块输出相位相同;
数据缓存模块包括N个FIFO,第n个FIFO对第n个通道的数据datan进行降速缓存,根据ad_clkn实现N个通道的数据datan同步输出;
频率测量模块用于从时钟控制模块接收频率测量参考时钟clks,从分频器模块接收方波信号
Figure FDA0002255825560000011
进行频率测量,得到N个通道数据对应的K分频后的频率
Figure FDA0002255825560000012
然后计算得到N个通道数据对应的频率
Figure FDA0002255825560000013
将N个频率
Figure FDA0002255825560000014
发送给同步抽取模块,将方波信号
Figure FDA0002255825560000021
转发至倍频模块;
倍频模块用于接收倍频参考时钟clk0和N个方波信号
Figure FDA0002255825560000022
分别对N个方波信号
Figure FDA0002255825560000023
进行K×Dn倍倍频得到倍频信号
Figure FDA0002255825560000024
将N个倍频信号
Figure FDA0002255825560000025
和对应的分频系数An输出至同步抽取模块,Dn表示从上位机接收的对第n个通道的第一倍频系数;
同步抽取模块用于从上位机接收FFT分析的信号周期数P,从倍频模块接收N个倍频信号
Figure FDA0002255825560000026
然后从数据缓存模块并行读取N个通道的数据datan,分别对每个通道的M路数据Unm根据对应的倍频信号
Figure FDA0002255825560000027
进行抽取,抽取方法如下:
设置每个通道的M路数据Unm的P个周期对应的闸门时间tn
Figure FDA0002255825560000028
倍频信号
Figure FDA0002255825560000029
与闸门时间tn的闸门信号边沿对齐,在闸门时间tn内对数据Unm根据倍频信号
Figure FDA00022558255600000210
进行抽取,将倍频信号
Figure FDA00022558255600000211
每两个上升沿之间的数据作为一个数据序列提取出来,在闸门时间tn共计提取出Dsum个数据序列,Dsum=P×K×Dn,记数据Unm中提取出的第d个数据序列为datanm[d],d=1,2,…,Dsum,记数据序列datanm[d]对应的分频系数为An[d];
将每个通道提取出的数据按照以下计算公式进行平均,得到第n个通道的平均数据data′n
Figure FDA00022558255600000212
其中data′n[d]表示平均数据data′n的第d个数据点,datanm[d][i]表示数据序列datanm[d]中的第i个数据,i=1,2,…,An
将得到的每个通道的平均数据data′n发送给FFT分析模块;
FFT分析模块包括N组FFTIP核和平均数据FIFO缓存模块,第n个FFTIP核用于接收第n个通道的平均数据data′n,对其进行FFT运算,得到FFT运算结果Fn;第n个平均数据FIFO缓存模块用于对第n个通道的平均数据data′n进行缓存,缓存之后可以发送至上位机进行波形显示与FFT运算结果进行对比;
上位机接收用户设置的FFT分析的信号周期数P和第一倍频系数Dn,将FFT分析的信号周期数P发送给同步抽取模块,将第一倍频系数Dn发送给倍频模块;从FFT分析模块中读取每个通道的FFT运算结果Fn和缓存的平均数据data′n进行频谱显示。
2.根据权利要求1所述的FFT分析装置,其特征在于,所述相位校正模块的相位校正方法如下:
1)N个ADC模块同时对一个正弦信号源进行采集;
2)判断N个ADC模块的输出数据是否相同,如果相同则记载当前对每个ADC模块的相位调整值,相位校正完全,否则进入步骤3);
3)以某一个ADC模块的输出数据为基准,控制另外N-1个ADC模块的相位调整值,调节输出信号相位,返回步骤2)。
3.根据权利要求1所述的FFT分析装置,其特征在于,所述频率测量模块的频率测量方法如下:
由频率测量参考时钟clks计数产生一个预设闸门信号,然后用方波信号
Figure FDA0002255825560000031
同步预设闸门信号产生实际闸门信号,频率测量参考时钟clks和方波信号
Figure FDA0002255825560000032
在实际闸门内计数,记频率测量参考时钟clks的频率为fs,实际闸门信号的时间内对标准时钟计数为Nums,对方波信号
Figure FDA0002255825560000033
的计数值为
Figure FDA0002255825560000034
则方波信号频率
Figure FDA0002255825560000035
的计算公式如下:
Figure FDA0002255825560000036
4.根据权利要求1所述的FFT分析装置,其特征在于,所述倍频模块包括N组周期计数器和小数分频器,其中第n个周期计数器用于采用倍频参考时钟clk0对方波信号
Figure FDA0002255825560000037
进行周期计数,记方波信号
Figure FDA0002255825560000038
一个周期内对倍频参考时钟clk0的计数值为Cn,将得到的计数值Cn发送给第n个小数分频器;
每个小数分频器包括除法器、加法器、比较器和分频器,其中第n个小数分频器中的除法器用于接收计数值Cn,计算Cn/(KDn),记商为Sn、余数为Yn,将商Sn发送给分频器,将余数Yn发送给加法器;
加法器用于计算和值Sumn=Sumn+Yn,Sumn的初始值为0,将得到的和值Sumn发送给比较器;
比较器用于将和值Sumn与第二倍频系数KDn进行比较,当Sumn>KDn时,令Sumn=Sumn-KDn,返回给加法器,且令标识flagn=1,当Sumn≤KDn时,令标识flagn=0;比较器将标识flagn发送给分频器;
分频器用于接收商Sn和标识flagn,根据标识flagn确定分频系数,当flagn=1时分频系数为Sn-1,当flagn=0时分频系数为Sn,根据分频系数对倍频参考时钟clk0进行分频,得到倍频信号
Figure FDA0002255825560000041
和分频系数An,分频系数An的表达式:
Figure FDA0002255825560000042
CN201810835993.2A 2018-07-26 2018-07-26 数字示波器的fft分析装置 Active CN108776264B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810835993.2A CN108776264B (zh) 2018-07-26 2018-07-26 数字示波器的fft分析装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810835993.2A CN108776264B (zh) 2018-07-26 2018-07-26 数字示波器的fft分析装置

Publications (2)

Publication Number Publication Date
CN108776264A CN108776264A (zh) 2018-11-09
CN108776264B true CN108776264B (zh) 2020-04-14

Family

ID=64030190

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810835993.2A Active CN108776264B (zh) 2018-07-26 2018-07-26 数字示波器的fft分析装置

Country Status (1)

Country Link
CN (1) CN108776264B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111796132B (zh) * 2020-09-08 2020-12-22 深圳市鼎阳科技股份有限公司 一种集成实时频谱分析功能的数字示波器
CN112067870B (zh) * 2020-09-15 2023-04-28 中电科思仪科技股份有限公司 一种基于fpga的示波器参数自动测量装置及方法
CN113486284A (zh) * 2021-05-24 2021-10-08 广州致远电子有限公司 支持外部数据源的示波器、示波方法及装置

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101126772A (zh) * 2007-09-11 2008-02-20 电子科技大学 一种示波器的高速信号重构方法
CN201053983Y (zh) * 2006-12-29 2008-04-30 徐红启 一种便携式通用数字存储示波器
JP2008301528A (ja) * 2008-09-01 2008-12-11 Kddi R & D Laboratories Inc Ofdm信号復調装置
CN203069643U (zh) * 2012-12-21 2013-07-17 优利德科技(成都)有限公司 多可变带宽数字存储示波器
CN103869157A (zh) * 2012-12-13 2014-06-18 特克特朗尼克公司 在测试和测量仪器中的自动中心频率和量程设定
CN104459315A (zh) * 2014-11-06 2015-03-25 电子科技大学 基于非基2fft变换的间谐波检测方法
CN106918731A (zh) * 2017-02-09 2017-07-04 深圳市鼎阳科技有限公司 一种数字示波器及其信号频率测量方法
CN106932640A (zh) * 2017-04-28 2017-07-07 金陵科技学院 基于fpga的脉冲信号参数测量方法及系统

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050165568A1 (en) * 2004-01-28 2005-07-28 Sullivan Steven K. Alias detection when displaying FFTS
CN105406859A (zh) * 2015-12-10 2016-03-16 武汉理工大学 单片全数字锁相环

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201053983Y (zh) * 2006-12-29 2008-04-30 徐红启 一种便携式通用数字存储示波器
CN101126772A (zh) * 2007-09-11 2008-02-20 电子科技大学 一种示波器的高速信号重构方法
JP2008301528A (ja) * 2008-09-01 2008-12-11 Kddi R & D Laboratories Inc Ofdm信号復調装置
CN103869157A (zh) * 2012-12-13 2014-06-18 特克特朗尼克公司 在测试和测量仪器中的自动中心频率和量程设定
CN203069643U (zh) * 2012-12-21 2013-07-17 优利德科技(成都)有限公司 多可变带宽数字存储示波器
CN104459315A (zh) * 2014-11-06 2015-03-25 电子科技大学 基于非基2fft变换的间谐波检测方法
CN106918731A (zh) * 2017-02-09 2017-07-04 深圳市鼎阳科技有限公司 一种数字示波器及其信号频率测量方法
CN106932640A (zh) * 2017-04-28 2017-07-07 金陵科技学院 基于fpga的脉冲信号参数测量方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《Improving and Qualifying Spectrum Analysis Made by Digital Scopes》;Consolatina Liguori etc.;《IEEE Transactions on Instrumentation and Measurement》;20071206;第56卷(第6期);第2411-2419页 *
《基于 SOPC 的便携式数字示波器设计与实现》;申礼斌等;《电子设计工程》;20140930;第22卷(第18期);第65-73页 *

Also Published As

Publication number Publication date
CN108776264A (zh) 2018-11-09

Similar Documents

Publication Publication Date Title
CN108776264B (zh) 数字示波器的fft分析装置
Chang et al. Measuring power system harmonics and interharmonics by an improved fast Fourier transform-based algorithm
Tiwari et al. Hardware implementation of polyphase-decomposition-based wavelet filters for power system harmonics estimation
CN102288807A (zh) 一种测量电网电压闪变的方法
CN103575984A (zh) 基于凯塞窗双谱线插值fft的谐波分析方法
CN110836992B (zh) 基于fpga的示波功率仪采集系统
CN108519511B (zh) 一种线性调频信号频率特征参数的时域测量方法
CN103983849A (zh) 一种实时高精度的电力谐波分析方法
CN110837000B (zh) 基于fpga的频率测量系统
CN104181391A (zh) 数字功率计谐波检测的方法
CN103969508A (zh) 一种实时高精密的电力谐波分析方法及装置
CN114252837B (zh) 用于数字功率计谐波测量的高精度均等取样系统及方法
CN101937017A (zh) 智能电表的交流采样动态去直流方法
CN101556325A (zh) 快速电能误差检定方法
CN113533848A (zh) 一种数字信号的测量方法和装置
CN110780251B (zh) 一种基于逐点计算delta的电力波形比对系统及方法
CN102072987B (zh) 短区间正弦信号的相位估计法及其实验装置
JP2000180484A (ja) 高調波測定装置
JP5553334B2 (ja) 正弦波パラメータ推定方法
Kwiatkowski et al. Hardware Computational Module for Frequency Ω-Counter
JP2010160093A (ja) 波形測定装置
CN110794208A (zh) 一种基于数字相位推算的宽带快速调制域测频方法
CN111092600B (zh) 基于相位叠加方法的fpga倍频方法
CN117169590B (zh) 一种基于软件变采样率的电力谐波分析的方法和装置
Kwiatkowski et al. Hardware computing module for frequency Ω-counter

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