CN103606004A - 一种倒数计数器频率分辨率增强方法 - Google Patents

一种倒数计数器频率分辨率增强方法 Download PDF

Info

Publication number
CN103606004A
CN103606004A CN201310576320.7A CN201310576320A CN103606004A CN 103606004 A CN103606004 A CN 103606004A CN 201310576320 A CN201310576320 A CN 201310576320A CN 103606004 A CN103606004 A CN 103606004A
Authority
CN
China
Prior art keywords
gate
arithmetic mean
frequency
time
frequency resolution
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
Application number
CN201310576320.7A
Other languages
English (en)
Other versions
CN103606004B (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 CN201310576320.7A priority Critical patent/CN103606004B/zh
Publication of CN103606004A publication Critical patent/CN103606004A/zh
Application granted granted Critical
Publication of CN103606004B publication Critical patent/CN103606004B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measuring Frequencies, Analyzing Spectra (AREA)
  • Measurement Of Unknown Time Intervals (AREA)

Abstract

本发明公开了一种倒数计数器频率分辨率增强方法,特征是利用高频率分辨率倒数计数器,通过FPGA发出周期闸门,在检测输入信号频率不发生改变的前提下,进行连续测量,对N次测量频率值求取算术平均值,显示结果为N次测量值的算术平均值,以达到频率分辨率增强的目的。本发明具有以下优点:1)可以在以往的倒数计数器高频率分辨率的基础上,进一步提高频率分辨率,频率分辨率一般能够增加2~3个有效位;2)通过连续测量,求取算术平均值,大大减少了随机误差,进一步提高了测量精度;3)通过FPGA控制实现与进行数据处理,易于实现而且测量速度快。

Description

一种倒数计数器频率分辨率增强方法
技术领域
本发明涉及一种倒数计数器频率分辨率增强方法。
背景技术
计数器的频率分辨率是指计数器能够在相近频率中检测的最小变化量。倒数计数器的实质是测量时间间隔,通过测量同步闸门的时间间隔与闸门内信号的脉冲数,计算出频率。单次时间间隔分辨率为100ps的倒数计数器的频率分辨率为10位/秒,即闸门时间为1s时,倒数计数器显示的有效位为10位。例如,测量1kHz的信号,频率分辨率为1uHz(1uHz=10-6Hz)。随机误差是影响分辨率最重要的因素,随机误差包括计数误差与触发误差。倒数计数器不可避免地出现±1个计数误差与触发误差。
倒数计数器的频率分辨率由给定的闸门时间决定,一般是采用显示的有效位数表示。在相同的闸门时间内,倒数计数器显示同样的分辨率位数,且与输入频率无关。现有技术中倒数计数器通过内插技术,一般能够实现单次时间间隔分辨率为100ps,频率分辨率由以下公式计算得到:Log10(闸门时间/100ps),若闸门时间为1s,频率分辨率为10位/秒,即显示的有效位为10位。
现有倒数计数器的倒数测量波形原理如图1所示。在现有倒数计数器方案中,信号为经过比较器整形后的方波,先与FPGA发出的闸门同步,得到与信号同步的闸门。同步闸门控制一个计数器对信号的脉冲计数,假定计的脉冲数为N1。同步闸门的时间间隔由另一计数器与内插电路测量得到。该计数器对时钟的脉冲计数,假定计的脉冲数为n1(计数时钟由系统内部提供,周期已知);因为计数时钟沿与同步闸门的前后沿不一定完全同步,通过内插电路分别测量出时钟沿与同步闸门前沿与后沿的时间间隔为△t1、△t2。△t1、△t2可通过基于延迟线原理的TDC芯片测量得出,如今采用TDC芯片能够实现的时间分辨率可达100ps。
测量出同步闸门的时间间隔t为:
t=n1×T0+Δt1-Δt2,T0为系统时钟的周期。
闸门时间内的频率的F可由以下公式计算:
1 s F = t N 1 = n 1 × T 0 + Δt 1 - Δt 2 N 1 ⇒ F = N 1 t = N 1 n 1 × T 0 + Δt 1 - Δt 2
在倒数计数器中,通过精确测量时间间隔,求取信号频率。由上式可以看出,通过内插技术,±1计数时钟误差得到了进一步的降低,测量值的计数误差转移到△t1、△t2中。此时系统的计数误差变为±tss(tss为TDC芯片的时间分辨率),若TDC芯片能够实现的时间分辨率:tss=100ps,则系统的±1个计数误差:
Figure BDA0000416117290000021
只进行单次倒数测量,闸门时间为1s时,频率分辨率即为10位/秒,显示的有效位数为10位。
我们知道,倒数计数器的测量精度由系统误差、时基误差与随机误差确定。系统误差是系统特有的,为固定值,一般是通过校准降低;时基误差由晶振特性决定,只能通过选择高精度、高性能的晶体来减少。系统的随机误差,包括量计数误差与触发误差两部分。通常所说的±1个计数误差是内部时钟频率和输入信号间的非同步性造成的误差。触发误差输入信号的噪声或来自计数器输入通道的噪声都可能引起噪声毛刺触发。倒数计数器单次倒数测量求值时,随机误差是频率分辨率的最大影响因素。
综上所述,现有技术中的倒数计数器频率分辨率还有待于进一步增强。
发明内容
本发明的任务在于提供一种倒数计数器频率分辨率增强方法。
其技术解决方案是:
一种倒数计数器频率分辨率增强方法,利用高频率分辨率倒数计数器,通过FPGA发出周期闸门,在检测输入信号频率不发生改变的前提下,进行连续测量,对N次测量频率值求取算术平均值,显示结果为N次测量值的算术平均值,以达到频率分辨率增强的目的。
上述倒数计数器频率分辨率增强方法,具体包括如下步骤:
a待测频率为输入信号经过比较器整形后的方波,先与FPGA发出的周期闸门同步,得到与信号同步的周期闸门;
b在第i次测量中,由同步闸门控制第一计数器对信号的脉冲计数,假设获得的脉冲数为Ni;同步闸门的时间间隔由第二计数器与内插电路测量得到,由第二计数器对时钟的脉冲计数,假定计的脉冲数为ni
c假设通过内插电路分别测量出时钟沿与同步闸门前沿与后沿的时间间隔为△t1、△t2;△t1、△t2通过基于延迟线原理的TDC芯片测量得出,该TDC芯片能够实现的时间分辨率达到100ps或100ps以上;
计算出第i次同步闸门的时间间隔ti为:
ti=ni×T0+Δt1-Δt2,T0为系统时钟的周期;
并且,第i次闸门时间内的频率的Fi由以下公式计算:
1 s F i = t i N i = n i × T 0 + Δt 1 - Δt 2 N i ⇒ F i = N i t i = N i n i × T 0 + Δt 1 - Δt 2
连续测量时,假设计算算术平均值的最大次数为N,判断输入信号频率没有改变的前提下,在第i次测量时,若i≤N,显示测量结果为前i次测量的算术平均值;若i>N时,显示的测量结果则为当前N次测量的算术平均值;在第N个闸门周期后,显示的测量结果都是N次测量的算术平均值;此时,随机误差的标准偏差公式变为:
σ ( x ‾ ) = σ ( X ) N ;
其中,X为有限次等精度测量的频率值,
Figure BDA0000416117290000033
为N次测量的算术平均值。
依据该公式可知,进行N次连续测量后,算术平均值的随机误差为单次测量值的随机误差的±
Figure BDA0000416117290000034
即算术平均值的随机误差比单次测量值的随机误差减小
Figure BDA0000416117290000035
倍,时间分辨率变为±1计数/
Figure BDA0000416117290000036
利用连续测量求取算术平均值与显示的有效位数的关系,即10次算术平均值能增加一个有效位,假设算术平均值的最大次数为100次,在判断输入信号频率不发生改变的前提下,在1s闸门时间时、100个闸门周期后,实现12位的频率分辨率。
本发明具有以下有益技术效果:
1)可以在以往的倒数计数器高频率分辨率的基础上,进一步提高频率分辨率,频率分辨率一般能够增加2~3个有效位。
2)通过连续测量,求取算术平均值,大大减少了随机误差,进一步提高了测量精度。
3)通过FPGA控制实现与进行数据处理,易于实现而且测量速度快。
附图说明
下面结合附图与具体实施方式对本发明作更进一步的说明:
图1为现有技术中倒数计数器倒数测量的波形原理示意图。
图2为本发明中倒数计数器的连续测量原理示意框图。
图3为本发明中倒数计数器连续测量的波形原理示意图。
图4为本发明中连续测量FPGA控制流程示意图。
具体实施方式
结合图2与图3,一种倒数计数器频率分辨率增强方法,利用高频率分辨率倒数计数器,即建立在倒数计数器高频率分辨率的基础上,通过FPGA发出周期闸门,在检测输入信号频率不发生改变的前提下,进行连续测量,对N次测量频率值求取算术平均值,显示结果为N次测量值的算术平均值,以达到频率分辨率增强的目的。
上述倒数计数器频率分辨率增强方法,具体包括如下步骤:
a待测频率为输入信号经过比较器整形后的方波,先与FPGA发出的周期闸门同步,得到与信号同步的周期闸门。
b在第i次测量中,由同步闸门控制第一计数器对信号的脉冲计数,假设获得的脉冲数为Ni;同步闸门的时间间隔由第二计数器与内插电路测量得到,由第二计数器对时钟的脉冲计数,假定计的脉冲数为ni,(计数时钟由系统内部提供,周期已知)。
c因为计数时钟沿与同步闸门的前后沿不一定完全同步,假设通过内插电路分别测量出时钟沿与同步闸门前沿与后沿的时间间隔为△t1、△t2;△t1、△t2通过基于延迟线原理的TDC芯片测量得出,该TDC芯片能够实现的时间分辨率达到100ps或100ps以上;
计算出第i次同步闸门的时间间隔ti为:
ti=ni×T0+Δt1-Δt2,T0为系统时钟的周期;并且,
第i次闸门时间内的频率的Fi由以下公式计算:
1 s F i = t i N i = n i × T 0 + Δt 1 - Δt 2 N i ⇒ F i = N i t i = N i n i × T 0 + Δt 1 - Δt 2
连续测量时,假设计算算术平均值的最大次数为N,判断输入信号频率没有改变的前提下,在第i次测量时,若i≤N,显示测量结果为前i次测量的算术平均值;若i>N时,显示的测量结果则为当前N次测量的算术平均值;在第N个闸门周期后,显示的测量结果都是N次测量的算术平均值;此时,随机误差的标准偏差公式变为:
σ ( x ‾ ) = σ ( X ) N ;
其中,X为有限次等精度测量的频率值,
Figure BDA0000416117290000043
为N次测量的算术平均值。
由该公式可知,进行N次连续测量后,算术平均值的随机误差为单次测量值的随机误差±
Figure BDA0000416117290000044
即算术平均值的随机误差比单次测量值的随机误差减小
Figure BDA0000416117290000045
倍。时间分辨率变为:±1计数/
Figure BDA0000416117290000046
而连续测量求取算术平均值与显示的有效位数的关系则是:10次算术平均值能增加一个有效位,10000次算术平均值能够增加4个有效位。当然,时间分辨率与频率分辨率不可能经过无限次平均后的无限提高,在进行一定数目的平均后,算术平均值会趋近于理论真值。
本发明中,假设算术平均的最大次数为100次。在判断输入信号频率不发生改变的前提下,在1s闸门时间时,100个闸门周期后,可以实现12位的频率分辨率。
参看图4,在本发明中,FPGA发出周期闸门,进行连续测量,算术平均的数目为100次。
F为测量的算术平均值,Fss为100ps时间分辨率所对应的频率分辨率,a为检测输入频率变化因子(理论上a=1,实际上,因为其它误差因素影响,适当增大)。判断关系式为输入信号频率发生变化的条件式。若判断关系式为真,说明输入信号频率发生改变,则测量重新开始;若判断关系式为假,说明输入信号频率没发生变化,连续测量继续。第i次测量时,i≤100时,则显示测量结果为当前i次测量算术平均值;i>100时,则显示测量结果为当前100次测量算术平均值,此时算术平均的计算公式如下:
Figure BDA0000416117290000051
从上式中,
Figure BDA0000416117290000052
为第i次闸门前N次频率测量的算术平均值,
Figure BDA0000416117290000053
为第i-1次闸门前N次频率测量的算术平均值,F(i)为第i次闸门频率测量值,F(i-N)为第i-N次闸门频率测量值,F(i)与F(i-N)的差实质为频率分辨率Fss。即为:
F ( i ) - F ( i - N ) = ± Fss = ± 10 - 10 ⇒ F ( i ) - F ( i - N ) 100 = ± 10 - 12
即在100个闸门周期后,频率显示的有效位数能增强2位。
对于1s的闸门时间,计数器的频率分辨率最大能达到12位,即12位/秒。例如测量1KHz的信号,设定1s的闸门时间,采用连续测量,在100个闸门周期(200s,闸门高有效,闸门周期为2s的方波)时,FPGA已经完成100个数据的测量,此后,显示的测量结果为前99个测量值与当前测量值的算术平均值,显示结果的有效位为12位,此时频率分辨率为10nHz(1nHz=10-9Hz),大大提高了频率分辨率。
上述方式中未述及的有关技术内容采取或借鉴已有技术即可实现。
需要说明的是,在本说明书的教导下本领域技术人员还可以作出这样或那样的容易变化方式,诸如等同方式,或明显变形方式。上述的变化方式均应在本发明的保护范围之内。

Claims (2)

1.一种倒数计数器频率分辨率增强方法,其特征在于利用高频率分辨率倒数计数器,通过FPGA发出周期闸门,在检测输入信号频率不发生改变的前提下,进行连续测量,对N次测量频率值求取算术平均值,显示结果为N次测量值的算术平均值,以达到频率分辨率增强的目的。
2.根据权利要求1所述的倒数计数器频率分辨率增强方法,其特征在于具体包括如下步骤:
a待测频率为输入信号经过比较器整形后的方波,先与FPGA发出的周期闸门同步,得到与信号同步的周期闸门;
b在第i次测量中,由同步闸门控制第一计数器对信号的脉冲计数,假设获得的脉冲数为Ni;同步闸门的时间间隔由第二计数器与内插电路测量得到,由第二计数器对时钟的脉冲计数,假定计的脉冲数为ni
c假设通过内插电路分别测量出时钟沿与同步闸门前沿与后沿的时间间隔为△t1、△t2;△t1、△t2通过基于延迟线原理的TDC芯片测量得出,该TDC芯片能够实现的时间分辨率达到100ps或100ps以上;
计算出第i次同步闸门的时间间隔ti为:
ti=ni×T0+Δt1-Δt2,T0为系统时钟的周期;
并且,第i次闸门时间内的频率的Fi由以下公式计算:
1 s F i = t i N i = n i × T 0 + Δt 1 - Δt 2 N i ⇒ F i = N i t i = N i n i × T 0 + Δt 1 - Δt 2
连续测量时,假设计算算术平均值的最大次数为N,判断输入信号频率没有改变的前提下,在第i次测量时,若i≤N,显示测量结果为前i次测量的算术平均值;若i>N时,显示的测量结果则为当前N次测量的算术平均值;在第N个闸门周期后,显示的测量结果都是N次测量的算术平均值;此时,随机误差的标准偏差公式变为:
σ ( x ‾ ) = σ ( X ) N ;
其中,X为有限次等精度测量的频率值,
Figure FDA0000416117280000016
为N次测量的算术平均值。
依据该公式可知,进行N次连续测量后,算术平均值的随机误差为单次测量值的随机误差的±
Figure FDA0000416117280000013
即算术平均值的随机误差比单次测量值的随机误差减小
Figure FDA0000416117280000014
倍,时间分辨率变为±1计数/
Figure FDA0000416117280000015
利用连续测量求取算术平均值与显示的有效位数的关系,即10次算术平均值能增加一个有效位,假设算术平均值的最大次数为100次,在判断输入信号频率不发生改变的前提下,在1s闸门时间时、100个闸门周期后,实现12位的频率分辨率。
CN201310576320.7A 2013-11-18 2013-11-18 一种倒数计数器频率分辨率增强方法 Expired - Fee Related CN103606004B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310576320.7A CN103606004B (zh) 2013-11-18 2013-11-18 一种倒数计数器频率分辨率增强方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310576320.7A CN103606004B (zh) 2013-11-18 2013-11-18 一种倒数计数器频率分辨率增强方法

Publications (2)

Publication Number Publication Date
CN103606004A true CN103606004A (zh) 2014-02-26
CN103606004B CN103606004B (zh) 2016-08-17

Family

ID=50124224

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310576320.7A Expired - Fee Related CN103606004B (zh) 2013-11-18 2013-11-18 一种倒数计数器频率分辨率增强方法

Country Status (1)

Country Link
CN (1) CN103606004B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1997013073A1 (en) * 1995-10-02 1997-04-10 Pos-Line Ab Method and valve apparatus for controlling a reciprocatable fuid actated power machine
CN101819231A (zh) * 2010-03-18 2010-09-01 天津农学院 连续计数间隔标记的高分辨率频率测量方法
CN103187967A (zh) * 2011-12-29 2013-07-03 深圳市汇川控制技术有限公司 基于fpga的plc高速脉冲计数实现系统及方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1997013073A1 (en) * 1995-10-02 1997-04-10 Pos-Line Ab Method and valve apparatus for controlling a reciprocatable fuid actated power machine
CN101819231A (zh) * 2010-03-18 2010-09-01 天津农学院 连续计数间隔标记的高分辨率频率测量方法
CN103187967A (zh) * 2011-12-29 2013-07-03 深圳市汇川控制技术有限公司 基于fpga的plc高速脉冲计数实现系统及方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王厚军 等: "《现代测试技术》", 31 August 2002, 电子科技大学出版社 *

Also Published As

Publication number Publication date
CN103606004B (zh) 2016-08-17

Similar Documents

Publication Publication Date Title
CN103837741B (zh) 基于fpga的等精度频率测试系统及其设计方法
CN100575965C (zh) 一种单片机测量频率的方法
CN104460304B (zh) 一种具有自动校正功能的高分辨率时间间隔测量装置
CN103487649B (zh) 一种兼容连续波和脉冲调制载波频率测量的方法及装置
CN101976037B (zh) 一种多次同步模拟内插的时间间隔测量方法和装置
CN102998970B (zh) 一种基于crio平台的高精度对时同步方法
TW201303315A (zh) 頻率量測方法及系統
CN106597097A (zh) 一种高精度频率测量方法
KR101240798B1 (ko) 리얼타임클럭 주파수 오프셋 검출장치 및 그 방법
CN106645952A (zh) 一种信号相位差的检测方法及系统
CN103529293A (zh) 基于边沿效应的并行的频率和周期性信号参数测量方法
KR100735407B1 (ko) Toa를 이용한 거리 측정 장치 및 방법
CN101556325A (zh) 快速电能误差检定方法
CN203929885U (zh) 基于fpga的等精度频率测试系统
CN103412474B (zh) 基于fpga的tdc-gp2测时范围高精度扩展电路
CN102495283B (zh) 自适应等精度测频方法
CN104635046A (zh) 一种调制域频率计数器及其连续测频方法
CN107908097B (zh) 采用混合内插级联结构的时间间隔测量系统及测量方法
CN103606004A (zh) 一种倒数计数器频率分辨率增强方法
JP2013024854A (ja) 距離測定方法及びそのシステム
JP5787096B2 (ja) 物理量測定装置、物理量測定方法
US7649969B2 (en) Timing device with coarse-duration and fine-phase measurement
KR20150015200A (ko) 제로-크로싱 방식의 전원 주파수 검출장치 및 제로-크로싱 방식의 전원 주파수 검출방법
CN106501621A (zh) 脉冲宽度采集方法
CN106443180B (zh) 测量闸门自适应的高精度频率测量方法及测量装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for 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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160817

Termination date: 20201118