CN101639499B - 基波跟踪的小信号谐波失真的测量装置及方法 - Google Patents

基波跟踪的小信号谐波失真的测量装置及方法 Download PDF

Info

Publication number
CN101639499B
CN101639499B CN2009101700221A CN200910170022A CN101639499B CN 101639499 B CN101639499 B CN 101639499B CN 2009101700221 A CN2009101700221 A CN 2009101700221A CN 200910170022 A CN200910170022 A CN 200910170022A CN 101639499 B CN101639499 B CN 101639499B
Authority
CN
China
Prior art keywords
frequency
measured signal
sample
harmonic
signal
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
Application number
CN2009101700221A
Other languages
English (en)
Other versions
CN101639499A (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 41 Institute
Original Assignee
CETC 41 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 41 Institute filed Critical CETC 41 Institute
Priority to CN2009101700221A priority Critical patent/CN101639499B/zh
Publication of CN101639499A publication Critical patent/CN101639499A/zh
Application granted granted Critical
Publication of CN101639499B publication Critical patent/CN101639499B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明公开了一种基波跟踪的小信号谐波失真的测量装置及方法,其中装置包括:接口单元、自动增益单元、峰值检测单元、低通滤波单元、频率测量单元、AD转换单元、FIFO、DDS、FPGA和DSP信号处理单元;本发明首先对被测信号频率进行准确测量,然后调整AD转换器的采样频率,使AD转换器采样频率根据设定的采样样本点数和测量的基波频率成整数倍关系,最后对样本数据进行FFT变换,在频域提取被测信号基波幅度和相应的N次谐波幅度计算出被测信号的谐波失真。

Description

基波跟踪的小信号谐波失真的测量装置及方法
技术领域
本发明涉及信号测量技术领域,尤其涉及一种基波跟踪的小信号谐波失真的测量装置及方法。
背景技术
目前最常用的失真度测量方法是基波剔除法,附图1为基波剔除法失真度测量的简化框图。其中最主要功能块就是可调谐的陷波器,在工作时,该滤波器手动或自动调谐到正弦波的基波上,基波被很大衰减,对2次和高次谐波实现尽量低的衰减或无衰减,从而实现谐波基本上无衰减通过,宽带噪声和交流电源有关的哼声和任何其他处在陷波器频率上下的干扰信号也可无衰减的通过。从电路实现上看设计一个将基波衰减较大,而对2次及其以上各次谐波产生较小衰减的陷波器是可以实现的,但是实现难度较大,而且可调谐的陷波器特性随环境温度变化较大,校准比较困难。
发明内容
鉴于上述的分析,本发明旨在提供一种基波跟踪的小信号谐波失真的测量装置及方法,用以解决现有技术中存在的的问题。
本发明的目的主要是通过以下技术方案实现的:
本发明提供了一种基波跟踪的小信号谐波失真的测量装置,包括:接口单元、自动增益单元、峰值检测单元、低通滤波单元、频率测量单元、AD转换单元、先入先出存储器FIFO、直接数字合成器DDS、现场可编程逻辑阵列 FPGA和数字信号处理器DSP,其中,
所述接口单元,用于根据被测信号选择输入端口;
所述峰值检测单元,用于对被测信号进行峰值测量,当检测到被测信号峰值超过峰值上限或低于峰值下限时,把被测信号峰值超过峰值上限或低于峰值下限的超限信息传给DSP;由DSP根据获得的超限信息通过FPGA控制所述自动增益单元的增益幅度;
所述自动增益单元,用于根据DSP的控制调节增益幅度,使被测信号幅度回到限定幅度范围内;
所述低通滤波单元,用于对所述自动增益单元输出的被测信号进行过滤,过滤掉测量范围之外的高频信号;
所述频率测量单元,用于对经过所述低通滤波单元滤波处理后的被测信号的基波频率进行测量,根据测量得到的基波频率设定AD转换器的采样频率;
所述DDS,用于根据所述频率测量单元输出的被测信号的基波频率,设定DDS频率控制字值,输出AD转换器所需的采样频率;
所述AD转换器,用于根据DDS设定的采样频率对所述低通滤波器输出的被测量信息进行采样,并将采样获得被测量信号的样本信息输出给FIFO;
所述FIFO,用于保存所述AD转换器输出的样本信息,当FIFO中的数据存储快满时发送一个FIFO快满信号给DSP;
所述DSP,用于在收到FIFO快满信号后产生一个中断,从FIFO中把被测信号的样本信息读走,对被测信号的样本信息进行FFT运算,把被测信号由时域转换到频域,然后在频域提取相应的基波和N次谐波的幅度值,进而计算得到被测信号的谐波失真。
进一步地,输入端口包括平衡输入和非平衡输入两种选择电路,平衡输入主要用于差分信号的测量,非平衡输入用于非差分信号的测量。
进一步地,所述频率测量单元具体用于,对经过所述低通滤波单元滤波处理后的被测信号的基波频率进行测量,根据测量得到的基波频率设定AD转换器的采样频率,采用如下公式设定AD转换器的采样频率:
f sample = N n * 2 * FREQ ;
其中,fsample表示AD转换器的采样频率,FREQ表示基波频率,N表示FFT变换的样本点数,n表示n为基波在频域的点位置,2表示有效点数为原来的一半。
进一步地,根据AD转换器的采样频率值生成AD9850的频率控制字,AD9850生成的频率作为AD转换器的采样频率,
AD9850频率控制字=232*fsample/CLKIN,其中,fsample表示AD转换器的采样频率,CLKIN表示AD9850的参考时钟。
进一步地,所述DSP具体用于,在收到FIFO快满信号后产生一个中断,从FIFO中把被测信号的样本信息读走,对被测信号的样本信息进行FFT运算,把被测信号由时域转换到频域,然后在频域提取相应的基波和N次谐波的幅度值,并根据如下公式计算出被测信号的谐波失真:
V BASE = V BASER 2 + V BASEI 2 ;
V hum - N = V hum - NR 2 + V hum - NI 2 ;
K f = Σ V hum - N 2 V BASE ;
其中,Kf表示谐波失真,VBASE表示基波幅度,VBASER表示基波实部,VBASEI表示基波虚部,Vhum-N表示谐波幅度,Vhum-NR表示谐波实部,Vhum-NI表示谐波虚部。
本发明还提供了一种基波跟踪的小信号谐波失真的测量方法,包括:
步骤A:峰值检测单元对被测信号进行峰值测量,当检测到被测信号峰值超过峰值上限或低于峰值下限时,把被测信号峰值超过峰值上限或低于峰值下限的超限信息传给数字信号处理器DSP,由DSP根据获得的超限信息设置自动增益单元的增益幅度,使被测信号幅度回到限定幅度范围内;
步骤B:低通滤波单元对所述自动增益单元输出的被测信号进行过滤,过滤掉测量范围之外的高频信号;
步骤C:频率测量单元对经过所述低通滤波单元滤波处理后的被测信号的基波频率进行测量,根据测量得到的基波频率确定AD转换器的采样频率;
步骤D:直接数字合成器DDS根据确定的AD转换器的采样频率值确定输出频率,AD转换器根据DDS确定的采样频率对所述低通滤波器输出的被测量信息进行采样,并将采样获得被测量信号的样本信息输出给先入先出存储器FIFO;
步骤E:FIFO保存所述AD转换器输出的样本信息,当FIFO中的数据存储快满时发送一个FIFO快满信号给DSP;
步骤F:DSP在收到FIFO快满信号后产生一个中断,从FIFO中把被测信号的样本信息读走,对被测信号的样本信息进行FFT运算,把被测信号由时域转换到频域,然后在频域提取相应的基波和N次谐波的幅度值,进而计算得到被测信号的谐波失真。
进一步地,所述步骤C具体包括:
对经过所述低通滤波单元滤波处理后的被测信号的基波频率进行测量,根据测量得到的基波频率设定AD转换器的采样频率,采用如下公式设定AD转换器的采样频率:
f sample = N n * 2 * FREQ ;
其中,fsample表示AD转换器的采样频率,FREQ表示基波频率,N表示FFT变换的样本点数,n表示n为基波在频域的点位置,2表示有效点数为原来的一半。
进一步地,所述步骤F具体包括:
DSP在收到FIFO快满信号后产生一个中断,从FIFO中把被测信号的样本信息读走,对被测信号的样本信息进行FFT运算,把被测信号由时域转换到频域,然后在频域提取相应的基波和N次谐波的幅度值,并根据如下公式计算出被测信号的谐波失真:
V BASE = V BASER 2 + V BASEI 2 ;
V hum - N = V hum - NR 2 + V hum - NI 2 ;
K f = Σ V hum - N 2 V BASE ;
其中,Kf表示谐波失真,VBASE表示基波幅度,VBASER表示基波实部,VBASEI表示基波虚部,Vhum-N表示谐波幅度,Vhum-NR表示谐波实部,Vhum-NI表示谐波虚部。
本发明有益效果如下:
本发明使用基波跟踪和自动增益技术,通过FFT变换把信号由时域转到频域,在频域通过被测信号基波和谐波的频谱分析,实现对小信号谐波失真测量。
本发明的其他特征和优点将在随后的说明书中阐述,并且,部分的从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过在所写的说明书、权利要求书、以及附图中所特别指出的结构来实现和获得。
附图说明
图1为现有技术中,基波剔除法失真度测量的简化框图;
图2为本发明实施例所述装置的结构示意图;
图3为本发明实施例所述方法的流程示意图。
具体实施方式
本发明的目的在于,将被测量信号由时域转到频域,通过对信号在频域的频谱分析实现对被测信号谐波失真精确测量。本发明首先对被测信号频率进行准确测量,然后调整AD转换器的采样频率,使AD转换器采样频率根据设定的采样样本点数和测量的基波频率成整数倍关系,最后对样本数据进行FFT(快速傅里叶变换)变换,在频域提取被测信号基波幅度和相应的N次谐波幅度计算出被测信号的谐波失真。
下面结合附图来具体描述本发明的优先实施例,其中,附图构成本申请一部分,并与本发明的实施例一起用于阐释本发明的原理。为了清楚和简化目的,当其可能使本发明的主题模糊不清时,将省略本文所描述的器件中已知功能和结构的详细具体说明。
首先结合附图2对对发明实施例所述装置进行详细说明。
如图2所示,图2为本发明实施例所述装置的结构示意图,主要包括:接口单元、自动增益单元、峰值检测单元、低通滤波单元、频率测量单元、AD转换单元、FIFO(先入先出存储器)、DDS(直接数字合成器)和由FPGA(现场可编程逻辑阵列)和DSP(数字信号处理器)信号处理单元组成的控制单元组成;以下将对各个部分分别予以详细说明。
(一)接口单元,用于信号输入端口的选择,输入端口包括平衡输入和非 平衡输入两种选择电路,平衡输入主要用于差分信号的测量,非平衡输入用于非差分信号的测量,DSP根据被测信号类型发出控制信号选择输入端口类型。
(二)峰值检测单元和自动增益单元,峰值检测单元对被测信号进行峰值测量,峰值检测电路中设置了被测信号峰值上限和下限检测,当峰值检测电路检测到信号峰值超过峰值上限或低于峰值下限时,把被测信号峰值超过峰值上限或低于峰值下限信息传给DSP,DSP根据获得的超限信息设置自动增益电路的增益幅度,使被测信号幅度回到限定幅度范围内,即使被测信号幅度控制在AD转换器满量程的30%~90%之间,自动增益部分的增益幅度从-20dB~60dB,从而可实现1mV~10V幅度范围的信号分析测量。
(三)低通滤波单元,用来实现对被测信号测量范围之外的高频信号的过滤,防止测量范围之外的高频信号干扰测量结果。
(四)频率测量单元,用于测量被测信号的基波频率,只有对被测信号基波频率的准确测量才能实现后面AD转换器的采样频率对基波频率的准确跟踪。如果AD转换器采样频率不能准确跟踪被测信号的基波频率,结果将造成通过FFT变换转换到频域时的信号频谱泄漏。由频率检测电路获得被测信号的基波频率,从而确定了AD转换器的采样频率;具体地,频率检测电路对测量被测信号的基波频率进行测量,然后根据测量得到的频率并结合进行FFT变换的样本点数设定AD转换器的采样频率fsample,公式(1)如下:
f sample = N n * 2 * FREQ 公式(1)
其中,fsample---AD转换器的采样频率,FREQ---基波频率,N---表示FFT变换的样本点数,n---表示n为基波在频域的Bin位置,2----是因为FFT运算的对称性,有效点数为原来的一半。
(五)DDS,有了AD转换器的采样频率值,如何产生AD转换器所需的精准采样频率也是影响测量精度的一个重要因素,为有效地使用基于FFT的频谱分析方法,AD转换器采样频率有以下几个基本关系:
1)与任何采样处理一样,被测信号的频率上限是1/2的采样频率;
2)FFT处理将产生一个幅度和频率的数据组,其中点(Bin)的数目是幅度和时间记录中样本数的一半;
3)FFT是线性频率的处理,所有的Bin均有同样的宽度,其宽度从0Hz(DC)到1/2采样频率。
例如,假定采用48KHz采样频率捕获一个1024个样本构成的幅度和时间记录,对这一记录执行FFT变换将产生512个幅度值,每个幅度值代表了频率Bin中频谱成分的振幅。频谱范围从0Hz到1/2采样频率(奈奎斯特频率),对应于24KHz的奈奎斯特频率和512个线性的Bin,每个Bin的宽度就是24000/512=46.875Hz,第一个Bin的中心是0Hz(DC),范围-23.4375Hz~23.4375Hz。第二个Bin的幅度值代表的是23.4375Hz~70.3125Hz间的信号幅度,以此类推可知其它Bin的信号范围。如果被测信号的基波频率没有落在FFT变换后对应Bin频率区间的中心点处,那末被测信号的基波频率将在整个频谱范围扩散造成频谱泄漏,同样可知对应的N次谐波频率的也不会落在对应Bin频率区间的中心点处,也会造成谐波频率的频谱泄漏,从而致使测量结果不准确。对于频谱泄漏一般可采取的方法是使用窗函数法,对样本数据加上不同的窗,对样本数据加窗后,在频域被测信号的能量就不会在整个频谱上分布了,但会溢出到相邻的Bin上,使用不同的窗会得到不同的频谱幅度测量误差,并且溢出到邻近Bin的能量也不相同。所以采用加窗的办法只是减轻了信号在频域的频谱泄漏,也 不能很好的对低谐波失真信号进行准确测量。
如果采用DSP的定时器来实现对AD转换器的采样频率,由于定时器本身精度的限制可能无法实现对基波信号的准确跟踪,其结果将是基波频率没有落在FFT变换后对应Bin频率区间的中心点处,从而造成信号频谱的部分泄漏,从而影响对低谐波失真信号的准确测量;本发明实施例中采用DDS来产生AD转换器所需的采样频率,DDS可以产生频率精确度较高的输出信号,这里可以选用型号为AD9850的DDS来产生采样频率,AD9850的频率设置字采用公式(2)计算得出,AD9850的频率设置字用来设置AD9850输出时钟频率,该时钟频率就是AD转换器的采样频率。
频率控制字=232*fsample/CLKIN    公式(2)
fsample---AD转换器的采样频率
CLKIN---AD9850的参考时钟
(六)AD转换器根据DDS确定的采样频率采样获得被测量信号的样本信息,并将样本信息输出给FIFO;
(七)FIFO,保存样本信息,当FIFO中的数据存储快满时发送一个FIFO快满信号给DSP;
(八)DSP,收到FIFO快满信号后产生一个中断,从FIFO中把被测信号的样本信息读走,对被测信号的样本信息进行FFT运算,把被测信号由时域转换到频域,然后在频域提取相应的基波和N次谐波的幅度值,使用公式(3)~公式(5)计算出被测信号的谐波失真。
V BASE = V BASER 2 + V BASEI 2 公式(3)
V hum - N = V hum - NR 2 + V hum - NI 2 公式(4)
K f = Σ V hum - N 2 V BASE 公式(5)
其中,Kf---谐波失真,VBASE---基波幅度,VBASER---基波实部,VBASEI---基波虚部,Vhum-N---谐波幅度,Vhum-NR---谐波实部,Vhum-NI---谐波虚部。
(九)FPGA电路,使用FPGA可编程设计的好处是避免了大量使用分立的逻辑器件,通过FPGA编程设计实现以往使用分立器件实现的功能,提高了硬件电路设计的灵活性;通过FPGA电路的逻辑设计实现DSP对接口电路、自动增益电路、峰值检测电路、AD转换器、FIFO以及DDS的读写控制;通过在FPGA内部设计计数器实现对输入信号的频率测量。
下面结合附图3对利用本发明实施例所述装置,进行基波跟踪的小信号谐波失真的测量的方法进行详细说明。
如图3所示,图3为本发明实施例所述方法的逻辑流程示意图,包括但不限于如下步骤:
步骤301:峰值检测单元对被测信号进行峰值测量,当检测到被测信号峰值超过峰值上限或低于峰值下限时,把被测信号峰值超过峰值上限或低于峰值下限的超限信息传给DSP,由DSP根据获得的超限信息设置自动增益单元的增益幅度,使被测信号幅度回到限定幅度范围内;
步骤302:低通滤波单元对所述自动增益单元输出的被测信号进行过滤,过滤掉测量范围之外的高频信号,并将经过滤波处理的被测信号输出给AD转换器和频率测量单元;
步骤303:频率测量单元对经过所述低通滤波单元滤波处理后的被测信号的基波频率进行测量,根据测量得到的基波频率并结合进行FFT变换的样本点数设定AD转换器的采样频率;具体的说就是,频率测量单元对经过所述低通滤 波单元滤波处理后的被测信号的基波频率进行测量,根据测量得到的基波频率设定AD转换器的采样频率,采用如下公式设定AD转换器的采样频率:
f sample = N n * 2 * FREQ ;
其中,fsample表示AD转换器的采样频率,FREQ表示基波频率,N表示FFT变换的样本点数,n表示n为基波在频域的Bin位置,2表示有效点数为原来的一半。
步骤304:DDS根据频率测量单元确定的的AD转换器的采样频率值确定采样频率,AD转换器根据DDS确定的采样频率对所述低通滤波器输出的被测量信息进行采样,并将采样获得被测量信号的样本信息输出给FIFO;
步骤305:FIFO保存所述AD转换器输出的样本信息,当FIFO中的数据存储快满时发送一个FIFO快满信号给DSP;
步骤306:DSP在收到FIFO快满信号后产生一个中断,从FIFO中把被测信号的样本信息读走,对被测信号的样本信息进行FFT运算,把被测信号由时域转换到频域,然后在频域提取相应的基波和N次谐波的幅度值,进而计算得到被测信号的谐波失真。具体的说就是,DSP在收到FIFO快满信号后产生一个中断,从FIFO中把被测信号的样本信息读走,对被测信号的样本信息进行FFT运算,把被测信号由时域转换到频域,然后在频域提取相应的基波和N次谐波的幅度值,并根据如下公式计算出被测信号的谐波失真:
V BASE = V BASER 2 + V BASEI 2 ;
V hum - N = V hum - NR 2 + V hum - NI 2 ;
K f = Σ V hum - N 2 V BASE ;
其中,Kf表示谐波失真,VBASE表示基波幅度,VBASER表示基波实部,VBASEI表示基波虚部,Vhum-N表示谐波幅度,Vhum-NR表示谐波实部,Vhum-NI表示谐波虚部。
综上所述,本发明实施例提供了一种基波跟踪的小信号谐波失真的测量装置及方法,使用基波跟踪和自动增益技术,通过FFT变换把信号由时域转到频域,在频域通过被测信号基波和谐波的频谱分析,实现对小信号谐波失真测量。
显然,本领域的技术人员应该明白,上述的本发明的各模块或者各步骤可以用通用的计算装置来实现,他们可以集中在单个的计算装置上,或者分布在多个计算装置所组成的网络上,可选的,它们可以用计算装置可执行的程序代码来实现,从而将它们存储在存储装置中,由计算装置来执行。或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样本发明不限于任何特定的硬件和软件结合。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求书的保护范围为准。

Claims (8)

1.一种基波跟踪的小信号谐波失真的测量装置,其特征在于,包括:接口单元、自动增益单元、峰值检测单元、低通滤波单元、频率测量单元、AD转换单元、先入先出存储器FIFO、直接数字合成器DDS、现场可编程逻辑阵列FPGA和数字信号处理器DSP,其中,
所述接口单元,用于根据被测信号选择输入端口;
所述峰值检测单元,用于对被测信号进行峰值测量,当检测到被测信号峰值超过峰值上限或低于峰值下限时,把被测信号峰值超过峰值上限或低于峰值下限的超限信息传给DSP;由DSP根据获得的超限信息通过FPGA控制所述自动增益单元的增益幅度;
所述自动增益单元,用于根据DSP的控制调节增益幅度,使被测信号幅度回到限定幅度范围内;
所述低通滤波单元,用于对所述自动增益单元输出的被测信号进行过滤,过滤掉测量范围之外的高频信号;
所述频率测量单元,用于对经过所述低通滤波单元滤波处理后的被测信号的基波频率进行测量,根据测量得到的基波频率设定AD转换器的采样频率;
所述DDS,用于根据所述频率测量单元输出的被测信号的基波频率,设定DDS频率控制字值,输出AD转换器所需的采样频率;
所述AD转换器,用于根据DDS设定的采样频率对所述低通滤波器输出的被测量信息进行采样,并将采样获得被测量信号的样本信息输出给FIFO;
所述FIFO,用于保存所述AD转换器输出的样本信息,当FIFO中的数据存储快满时发送一个FIFO快满信号给DSP;
所述DSP,用于在收到FIFO快满信号后产生一个中断,从FIFO中把被测信号的样本信息读走,对被测信号的样本信息进行FFT运算,把被测信号由时域转换到频域,然后在频域提取相应的基波和N次谐波的幅度值,进而计算得到被测信号的谐波失真。
2.根据权利要求1所述的装置,其特征在于,输入端口包括平衡输入和非平衡输入两种选择电路,平衡输入主要用于差分信号的测量,非平衡输入用于非差分信号的测量。
3.根据权利要求1所述的装置,其特征在于,所述频率测量单元具体用于,对经过所述低通滤波单元滤波处理后的被测信号的基波频率进行测量,根据测量得到的基波频率设定AD转换器的采样频率,采用如下公式设定AD转换器的采样频率:
f sample = N n * 2 * FREQ ;
其中,fsample表示AD转换器的采样频率,FREQ表示基波频率,N表示FFT变换的样本点数,n表示n为基波在频域的点位置,2表示有效点数为原来的一半。
4.根据权利要求1所述的装置,其特征在于,根据AD转换器的采样频率值生成AD9850的频率控制字,AD9850生成的频率作为AD转换器的采样频率,
AD9850频率控制字=232*fsample/CLKIN,其中,fsample表示AD转换器的采样频率,CLKIN表示AD9850的参考时钟。
5.根据权利要求1到4中任意一项所述的装置,其特征在于,所述DSP具体用于,在收到FIFO快满信号后产生一个中断,从FIFO中把被测信号的样本信息读走,对被测信号的样本信息进行FFT运算,把被测信号由时域转换到频域,然后在频域提取相应的基波和N次谐波的幅度值,并根据如下公式计算出被测信号的谐波失真:
V BASE = V BASER 2 + V BASEI 2 ;
V hum - N = V hum - NR 2 + V hum - NI 2 ;
K f = Σ V hum - N 2 V BASE ;
其中,Kf表示谐波失真,VBASE表示基波幅度,VBASER表示基波实部,VBASEI表示基波虚部,Vhum-N表示谐波幅度,Vhum-NR表示谐波实部,Vhum-NI表示谐波虚部。
6.一种基波跟踪的小信号谐波失真的测量方法,其特征在于,包括:
步骤A:峰值检测单元对被测信号进行峰值测量,当检测到被测信号峰值超过峰值上限或低于峰值下限时,把被测信号峰值超过峰值上限或低于峰值下限的超限信息传给数字信号处理器DSP,由DSP根据获得的超限信息设置自动增益单元的增益幅度,使被测信号幅度回到限定幅度范围内;
步骤B:低通滤波单元对所述自动增益单元输出的被测信号进行过滤,过滤掉测量范围之外的高频信号;
步骤C:频率测量单元对经过所述低通滤波单元滤波处理后的被测信号的基波频率进行测量,根据测量得到的基波频率确定AD转换器的采样频率;
步骤D:直接数字合成器DDS根据确定的AD转换器的采样频率值确定输出频率,AD转换器根据DDS确定的采样频率对所述低通滤波器输出的被测量信息进行采样,并将采样获得被测量信号的样本信息输出给先入先出存储器FIFO;
步骤E:FIFO保存所述AD转换器输出的样本信息,当FIFO中的数据存储快满时发送一个FIFO快满信号给DSP;
步骤F:DSP在收到FIFO快满信号后产生一个中断,从FIFO中把被测信号的样本信息读走,对被测信号的样本信息进行FFT运算,把被测信号由时域转换到频域,然后在频域提取相应的基波和N次谐波的幅度值,进而计算得到被测信号的谐波失真。
7.根据权利要求6所述的方法,其特征在于,所述步骤C具体包括:
对经过所述低通滤波单元滤波处理后的被测信号的基波频率进行测量,根据测量得到的基波频率设定AD转换器的采样频率,采用如下公式设定AD转换器的采样频率:
f sample = N n * 2 * FREQ ;
其中,fsample表示AD转换器的采样频率,FREQ表示基波频率,N表示FFT变换的样本点数,n表示n为基波在频域的点位置,2表示有效点数为原来的一半。
8.根据权利要求6所述的方法,其特征在于,所述步骤F具体包括:
DSP在收到FIFO快满信号后产生一个中断,从FIFO中把被测信号的样本信息读走,对被测信号的样本信息进行FFT运算,把被测信号由时域转换到频域,然后在频域提取相应的基波和N次谐波的幅度值,并根据如下公式计算出被测信号的谐波失真:
V BASE = V BASER 2 + V BASEI 2 ;
V hum - N = V hum - NR 2 + V hum - NI 2 ;
K f = Σ V hum - N 2 V BASE ;
其中,Kf表示谐波失真,VBASE表示基波幅度,VBASER表示基波实部,VBADEI表示基波虚部,Vhum-N表示谐波幅度,Vhum-NR表示谐波实部,Vhum-NI表示谐波虚部。
CN2009101700221A 2009-09-01 2009-09-01 基波跟踪的小信号谐波失真的测量装置及方法 Expired - Fee Related CN101639499B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009101700221A CN101639499B (zh) 2009-09-01 2009-09-01 基波跟踪的小信号谐波失真的测量装置及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009101700221A CN101639499B (zh) 2009-09-01 2009-09-01 基波跟踪的小信号谐波失真的测量装置及方法

Publications (2)

Publication Number Publication Date
CN101639499A CN101639499A (zh) 2010-02-03
CN101639499B true CN101639499B (zh) 2011-09-14

Family

ID=41614585

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009101700221A Expired - Fee Related CN101639499B (zh) 2009-09-01 2009-09-01 基波跟踪的小信号谐波失真的测量装置及方法

Country Status (1)

Country Link
CN (1) CN101639499B (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102195725A (zh) * 2011-05-03 2011-09-21 牛刚 无线电台检测模块
CN102353525A (zh) * 2011-07-08 2012-02-15 电子科技大学 一种电光强度调制器谐波失真测量方法
CN105548707B (zh) * 2015-12-11 2018-04-06 中国电子科技集团公司第四十一研究所 一种电磁信号谐波失真自动测量方法及系统
CN105699770B (zh) * 2016-02-01 2019-01-08 中国科学院国家天文台 百万通道频谱仪的频点输出信号失真恢复方法
CN105891601A (zh) * 2016-06-23 2016-08-24 中国人民解放军63888部队 一种谐波分量辅助测试装置
CN106645942B (zh) * 2016-07-26 2020-02-18 北京东方计量测试研究所 一种低成本高精度嵌入式信号采集分析系统和方法
CN109283390A (zh) * 2018-10-10 2019-01-29 江苏理工学院 一种智能自动增益控制的频谱分析仪
CN109828145A (zh) * 2018-12-30 2019-05-31 南京传麒信息技术合伙企业(普通合伙) 一种宽范围自适应电压监测装置及方法
CN114546029B (zh) * 2019-12-30 2022-12-02 珠海极海半导体有限公司 控制芯片、mcu芯片、mpu芯片及dsp芯片
CN112462140B (zh) * 2021-02-02 2021-04-13 成都能通科技有限公司 一种提供电能参数分析的频率跟踪方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4902960A (en) * 1988-03-17 1990-02-20 Myron Zucker, Inc. Voltage input for harmonimeter, and methods of constructing and utilizing same
US4954785A (en) * 1989-04-12 1990-09-04 Sundstrand Corporation Auto tracking notch filter using switched capacitors to measure harmonic distortion and noise contained in a signal source

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4902960A (en) * 1988-03-17 1990-02-20 Myron Zucker, Inc. Voltage input for harmonimeter, and methods of constructing and utilizing same
US4954785A (en) * 1989-04-12 1990-09-04 Sundstrand Corporation Auto tracking notch filter using switched capacitors to measure harmonic distortion and noise contained in a signal source

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
付平等.直接数字合成正弦波谐波失真度测试方法.《测控技术》.2002,第21卷(第5期),46-48. *
余厚全等.数字正弦波发生器的谐波失真度与取样点数和样点长度的关系.《石油仪器》.1991,(第2期),76-80. *
倪陈强.基于虚拟仪器的正弦信号失真度分析仪.《微计算机信息》.2007,第23卷(第10-1期),118-119. *
李少东.DDS技术在VXI总线误码仪模块中的应用.《电子质量》.2006,(第6期),8-9. *
林琳等.基于LabVIEW 的谐波失真分析.《西华大学学报·自然科学版》.2006,第25卷(第5期),57-59,70. *
马月辉等.基于FFT 的低频谐波失真度测试仪.《自动化与仪器仪表》.2001,(第6期),30-32. *
魏燕定等.基于FFT 的超低频谐波失真度测试仪.《计量与测试技术》.1999,(第3期),2-3. *

Also Published As

Publication number Publication date
CN101639499A (zh) 2010-02-03

Similar Documents

Publication Publication Date Title
CN101639499B (zh) 基波跟踪的小信号谐波失真的测量装置及方法
US9182500B2 (en) Method and system for amplitude digitization of nuclear radiation pulses
CN106506101B (zh) 一种接收机幅频响应自动化校准方法及装置
CN101183050B (zh) 一种位移测量的电液伺服阀动态性能测试方法
CN102506994A (zh) 一种数字化声学检测系统
JP2013181766A (ja) 中性子測定装置およびその測定方法
CN102445600A (zh) 基于传统相位差测量的优化方法和电路
CN101320060A (zh) 快速相位计
CN103236825A (zh) 一种用于高精度数据采集系统的数据校正方法
CN1220189C (zh) 确定盘的品质的装置和方法
CN102495280B (zh) 一种抗噪音宽带频率测量方法及锁相频率计
CN106199389A (zh) 一种高效锁相放大器性能测试装置及方法
CN103543331B (zh) 一种计算电信号谐波和间谐波的方法
CN108332845B (zh) 噪声测量方法及噪声仪
CN112014810B (zh) 基于fpga的电子侦察信号参数高精度测量方法
Peterson Frequency detection and speech formants
Xu et al. Frequency-feature based antistrong-disturbance signal processing method and system for vortex flowmeter with single sensor
CN108181486A (zh) 加速度信号的处理方法及装置
CN103499386B (zh) 一种精密测量频率的方法和装置
CN104897250A (zh) 一种抗强谐波干扰的多台位流量计脉冲计数补偿方法
CN201917612U (zh) 一种相位差测量装置
CN101419090A (zh) 一种目标噪声测量中的阵列噪声信号的聚焦方法
CN105784103A (zh) 一种基于非线性调频激励的变信噪比的频率特性测量方法
CN104990616A (zh) 基于级联自适应陷波器的多路不同步脉冲计数补偿方法
CN116299643A (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: 20110914

Termination date: 20170901

CF01 Termination of patent right due to non-payment of annual fee