CN100520419C - 改进的频率确定 - Google Patents

改进的频率确定 Download PDF

Info

Publication number
CN100520419C
CN100520419C CNB2004800197387A CN200480019738A CN100520419C CN 100520419 C CN100520419 C CN 100520419C CN B2004800197387 A CNB2004800197387 A CN B2004800197387A CN 200480019738 A CN200480019738 A CN 200480019738A CN 100520419 C CN100520419 C CN 100520419C
Authority
CN
China
Prior art keywords
count value
counter
intermediate count
ratio
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.)
Expired - Fee Related
Application number
CNB2004800197387A
Other languages
English (en)
Other versions
CN1820203A (zh
Inventor
P·博德
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.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips Electronics NV
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 Koninklijke Philips Electronics NV filed Critical Koninklijke Philips Electronics NV
Publication of CN1820203A publication Critical patent/CN1820203A/zh
Application granted granted Critical
Publication of CN100520419C publication Critical patent/CN100520419C/zh
Expired - Fee Related 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/02Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage
    • G01R23/10Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage by converting frequency into a train of pulses, which are then counted, i.e. converting the signal into a square wave

Abstract

简单的传统频率测量电路包含作为必要单元的两个计数器。在这些计数器中可能出现这样的问题:与测量时间间隔相应,计数值有±1/N的相对误差。按照本发明,使用这些计数器的中间计数值对来估算频率比例。按照本发明的一个方面,这可以基于根据计数器的中间计数值对确定的回归线的陡度来进行。而且,按照本发明的一个方面,在取样中间取样值的同时,两个计数器继续进行计数。

Description

改进的频率确定
技术领域
本发明涉及测量或者确定频率的领域。具体来说,本发明涉及一种用于确定第一信号的第一频率和第二信号的第二频率的比值的电路、一种确定第一信号的第一频率和第二信号的第二频率的比值的方法以及一种包含计算程序装置的计算机程序产品。
背景技术
用于频率测量的电子电路在电机和电子系统的广阔领域中得到具体表现,比如CD/硬盘驱动器中的旋转控制器或者用于通过独立的时钟振荡器进行操作的电子子系统的同步的接口电路。
针对后者的典型实例是具有GPS能力的GSM移动电话,包括多达三个独立的晶体振荡器(XO):分别在工作时间和空闲时间供应GSM子系统的26MHz XO和32kHz XO以及供应GPS子系统的20MHz X0。在GSM工作模式下,需要26MHz XO作为GSM RF合成器的归零基准,低功率32kHz XO在GSM空闲模式下有助于减少功率消耗。在GSM空闲模式下,GSM系统z中的大部分,包括26MHz XO,将会断电。只有32kHzX0和一个计数器(称为32kHz“睡眠计数器”)工作。该睡眠计数器用作在由基站发送的新传呼消息到达之前及时对GSM系统加电的定时器。由于32kHz X0的典型公差,GSM移动电话可能丢失一部分传呼消息。这就是为什么GSM移动电话要包括一个能够实现对32kHz XO的频率和26MHz XO的频率之间的比例进行评估的频率测量单元的原因。近来随着具有GPS性能的GSM移动电话的出现,激起了对频率测量的同样的需求。这是因为GPS接收机IC传统上使用另一种参考频率,而不是GSM的26MHz,因此需要一个专用的GPS XO。
本申请的要点是由精度要求决定的测量反应时间。可以将相对误差随着测量间隔的增加而减少看作基本原则。传统的解决办法在两种情况下都产生几秒的测量时间间隔。主流观点是基于计数器的传统电路是最适宜的,在简单数字硬件实现的约束下,没有给改进留出空间。另一方面,众所周知的是,如果更复杂的硬件可被接受,不兼顾精度的话,测量间隔可被减小。我们例如可以想象到下面的方法:a)使数字时钟通过一个消除谐波的滤波器,b)用ADC对结果得到的正弦波形进行数字化,c)应用某些DSP算法来进行频率评估,类似于GSM接收机中使用的算法。
两个数字时钟信号的频率的比例的测量通常是借助相对简单的数字电路来进行的。这些电路通常包含两个由时钟信号触发的计数器。为了确定两个频率的比例,在某一测量时间之后,两个计数值的比例就是两个时钟信号的频率的频率比例的示值。
典型地将测量时间选择成覆盖多个时钟周期。测量时间的选择通常是借助计数器进行的,从而定义出一个测量窗口。因此,这个定义测量窗口的计数器的变化是预设的。因此,仅对各个其它计数器的计数器值的变化进行取样。取决于两个时钟信号的相位,从计数器读取的计数值可依据一个脉冲或者周期而变化。这样,测量结果可能有±1/N的最大相对误差,其中N是从计数器读取的计数值的变化。该相对误差可以通过选择更长的测量时间(因此有更大的N)来减少。不过,对于更大的应用范围,比如移动电信(如上所述)、例如用于CD播放器、用于计算机的硬盘或汽车电子中的应用的回转或旋转速度控制或者用于实验室或用于制造工艺的电子测量设备,越来越希望得到精确且快速的频率测量,并且在某些应用中必须要的得到这种精确且快速的频率测量。
US 4,350,950公开了一种频率测量方法,其中借助于振荡器中生成的参考信号来测量脉冲信号的频率。通过分离的计数器对被测脉冲频率和参考脉冲频率二者进行计数,并且计数值或计数结果被存储在分离的存储器中。触发器被置位和复位,以确定每个测量时段。被测频率计数器的两个连续计数之间的差对应于在测量时段期间被计数的脉冲或周期。参考频率计数器的两连续计数之间的差对应于测量时段的长度,从而通过参频率确定分辨率。微处理器划分测量时段计数得到的周期数,以确定被测频率。
US2003/0117181A1公开了用于确定脉冲输入信号的频率的方法和系统。在预定测量间隔的开始和末尾捕捉脉冲计数和定时器计数,以获得开始脉冲计数和末尾脉冲计数以及开始脉冲时刻和末尾脉冲时刻。通过末尾脉冲计数减去开始脉冲计数除以末尾脉冲时刻减去开始脉冲时刻来确定脉冲频率。
发明内容
本发明的一个目的是提供一种快速且精确的频率测量。
按照权利要求1所述的本发明的示范性实施方式,上述发明目的可以通过权利要求1中提出的用于确定第一信号的第一频率和第二信号的第二频率的比例的电路来解决。按照本发明地示范性实施方式的一个方面,成对的取样中间计数值的结果序列对可在笛卡尔坐标系统中用圆点表示。因此,取决于测量频率的稳定性,中间计数值近似于一条直线。按照本发明的一个方面,相应回归线的陡度可被确定,并且可被看作想要得到的频率比例的估算。
具有优势地,当保持测量间隔与上述的传统方法的测量间隔相同时,测量误差明显减少。根据本发明的一个方面,发现与上述传统方法的测量窗口相同时,测量误差可以减少,如果25个取样,测量误差减少到二分之一,也就是,第一和第二计数器的成对中间计数值被考虑。换句话说,成对中间计数值的数目的增加使测量误差减少。
根据本发明的一个方面,取样数目另外增加到n倍,将导致测量误差另外减少
Figure C200480019738D0007100313QIETU
另一方面,当传统方法的精度足够时,在保持相同测量精度的情况下,上述电路可以减少测量时间。
根据本发明另一示范性实施方式,第一计数器被第一信号的上升边缘和下降边缘之一触发,第二计数器被第二信号的上升边缘和下降边缘之一触发,第二计数器的时钟信号是第一和第二信号之一。第一计数器被取样时,第二计数器的第二中间计数值被预置在寄存器里。
根据本发明另一示范性实施方式,提供一个存储器,包括第一和第二存储器。第一存储器用于存储第一计数器的第一中间计数值,以致于第一计数器的第一中间计数值的序列被提供,第二存储器用于存储第二计数器的第二中间计数值,以致于第二计数器的第二中间计数值的序列被提供。
根据本发明另一示范性实施方式,提供一种方法,该方法用于确定第一信号的第一频率和第二信号的第二频率的比例。
根据本发明另一示范性实施方式,提供一个计算机程序产品,该计算机程序产品包含计算机程序代码装置。根据本发明示范性实施方式的一个方面,计算机程序产品可以是一个计算机可读介质,比如CD-ROM。计算机程序代码装置涉及一个计算机程序,当计算机程序代码装置在一个处理器执行时,使得处理器执行相应于本发明的方法的操作,计算机程序代码装置可用任何适合的程序语言书写,比如C++。取代被存储在计算机程序产品中,计算机程序代码装置(比如计算机程序)也可从网络(比如万维网)获得,可从网络下载到计算机的内部存储器,处理器或者其它适当的设备。
下述可以看作本发明的示范性实施方式的要点,在计数器保持运行的同时,成对中间计数值被取样和存储。随后,基于这些成对中间计数值,希望频率比例的估算被确定。根据本发明,发现当成对中间计数值的结果序列如笛卡尔坐标中的圆点所示时,其近似于一条直线。根据本发明的一个方面,相应回归线的陡度可被计算,并且可被看作希望频率比例的估算。根据本发明,一部分电路可以在时钟域操作,另一部分电路可以在选通域操作。
附图说明
参照下面的示范性实施方式,本发明的这些和其它方面将会变得显而易见。
本发明的示范性实施方式将随后描述,参考下面的附图:
附图1示出了根据本发明的电路的一个示范性实施方式的简化电路图;
附图2是操作附图1所示的电路的方法的流程图;
附图3是进一步解释本发明的采集时钟的标准相位和选通时钟的标准相位的比较图;
附图4是根据本发明示范性实施方式的采集时钟的标准相位和选通时钟的标准相位的另一比较图;
附图5示出了根据本发明的电路的第二示范性实施方式的简化电路图;
附图6示出了附图5的电路中出现的信号的时间图;
附图7示出了根据本发明电路的第三示范性实施方式的简化电路图;
附图8示出了附图7的电路中出现的信号的时间图。
具体实施方式
附图1表示用于确定第一数字信号S1的第一频率f1和第二数字信号S2的第一频率f2的比例的电路的一种示范性实施方式的简化电路图。从附图1可以看出,第一信号S1输入到第一计数器2,第二信号S2输入到第二计数器4。第一计数器2在第一信号S1的各个上升边缘或者下降边缘处受到触发或者进行递增。
从附图1可以看出,第一信号S1也作为时钟信号输入到第二计数器4。第二计数器4按照第一信号S1的上升边缘或者下降边缘计时。随后,计数器4的计数值根据第二信号S2的上升边缘或者下降边缘变化。
附图标记6表示连接到第一计数器2的第一寄存器。第一寄存器设置为用于在第一计数器2的计数期间存储第一计数器的中间计数值。换句话说,在第一计数器2继续计数的同时,第一寄存器6可以取样第一计数器2的中间计数值。
而且,提供了连接到第二计数器4的第二寄存器8。从附图1可以看出,第二寄存器8也连接到第一寄存器6。可将第二寄存器8设置成,在第二计数器4的预设计数值下,第二寄存器8向第一寄存器6输出一个触发信号。随后,当第一寄存器6接收到来自第二寄存器8的触发信号时,第一寄存器6取样第一计数器2的中间计数值。据此,第二寄存器8可以定义第二计数器4的中间计数值,以致当第二计数器4达到这些预设中间计数值时,第二寄存器8触发第一寄存器6,使得第一寄存器6取样第一寄存器6的中间计数值。
第一寄存器6和第二寄存器8均连接到同步单元10。同步单元10适用于控制第二计数器4到达预设计数值与第一寄存器6对第一计数器2进行取样之间的时间段。换句话说,同步单元10控制第二计数器4到达预置计数值和取样第一计数器2的中间计数值之间的时间关系。
第一寄存器6连接到第一存储器12。第一存储器12起到第一寄存器的扩展的作用。因此,在附图1示出的电路的操作期间,第一计数器2的中间计数值的序列存储在第一存储器12中。
计数器2和4可以是有限状态机(FSMs)。而且,包含在虚线框中的所有单元(包括计数器2和4)可借助FPGAs、PLDs、EPLDs、ASICs或者适合的Ics来实现。
还有第二存储器14,它与第二寄存器8相连。在第二存储器14中,可以预设要对第一计数器2进行取样时的计数值或者时间点。
第一和第二存储器12和14连接到计算单元16。计算单元适用于,基于存储在第一和第二存储器12和14中的第一和第二计数器2和4的计数值序列,计算两个频率S1和S2的比例的估算值。在频率比例确定之后,计算单元16将确定或者测量结果输出给输出单元18。
同步单元10和计算单元16可通过适当的硬件实现。它们也可借助有限状态机(FSMs)或者FPGAs、PLDs、EPLDs、ASICs或者适当的Ics来实现。但是,如围绕着第一和第二寄存器6和8、同步单元10、第一和第二存储器12和14、计算单元16的虚线所示,这些组成部分也可用适当的处理器系统实现。然后,上述系统的操作由适当的程序控制。这一程序可被存储在适当的机器可读介质中,比如CD-ROM。上述程序可用任何适当的语言编写,比如汇编语言或者C++。寄存器6和8、存储器12和14可由处理器的内部存储器实现。对这种存储区域的读访问和写访问可借助中断服务程序实现。
两个频率的比例f1/f2的估算值的求算结果在计算单元16中进行求算,可描述如下:
相应的中间计数值对,也就是在相应时间点确定的第一和第二计数器2和4的中间计数值,在笛卡尔坐标系统中用点表示。然后,按照本发明的一个方面,这些点成一条直线。按照本发明的一个方面,计算相应回归线的陡度,该陡度相当于想得到的频率比例的估算值。按照本发明的一个方面,借助线性回归确定该直线。
按照本发明的这种示范性实施方式的变形方式,计算单元16适用于计算随时间变化的调制率,也就是两个频率f1/f2的比例。
频率测量单元的任务是测量两个无关数字时钟信号的频率比例,上述两个独立数字时钟信号分别用选通时钟(gating clock)和采集时钟(collecting clock)表示。相应的时钟频率和时钟周期可表示为fg,fc,Tg=1/fg和Tc=1/fc。如上所述,简单的传统频率测量电路包含两个计数器作为基本单元:选通计数器通过选通时钟定时,采集计数器通过采集时钟定时。与某些比较器逻辑组合的选通计数器,通过使得采集计数器对于一定数目的选通计数器周期Ng有效,来确定测量时间间隔。采集计数器在测量时间间隔期间采集的周期数用Nc表示。想要得到的频率比例fc/fg可如下确定:
用Tg和Tc的多项式表示测量时间间隔
NgTg=NcTc+dTc                      (1)
其中,NgTg是几秒的测量时间间隔,项dTc(-1<d<1)针对的是测量时间间隔可能不是Tc的整数倍而采集计数器只能处理Tc的整数倍的情况。d的性质在后文中给出。(1)也可写为:
f c f g = N c N g &CenterDot; ( 1 + d N c )
其中项ε=d/Nc表示相对测量误差。举一个例子,在最坏的情况下,d=1,Nc=26·106,相对误差ε=0.04·10-6。这反映了具有测量间隔NcTc=1秒和Tc=1/fc=1/fGSM=1/26MHz的典型的GSM/GPS测量的情况。
为了促成新的频率测量方法,需要进一步深入了解问题。为此,引入了时钟信号的归一化相位,并且表明计数器的当前值表示它的量化形式。
矩形波时钟信号可与
Figure C200480019738D00112
(t)相关联,
Figure C200480019738D00113
(t)是该信号的傅立叶基音(tone)的相位。
Figure C200480019738D00114
(t)是线性倾斜的,它可由初始相位φ和频率f表示:
Figure C200480019738D00115
注意,在本文中通篇将(t)简单地理解为是经过展开的,意思就是说不局限于区间[0,2π]。
因此,在计数器的当前值和时钟信号的相位斜率的当前值之间存在关系。为了使该关系更加明显,量化和归一化后的相位
Figure C200480019738D0011100527QIETU
可如下定义:
Figure C200480019738D00117
其中,[x]是等于或者小于x的最接近的整数,以致0≤x-[x]<1。量化导致相位倾斜变成分段式的,同时归一化使2π区间成1的周期。因此,可以将计数器看成提供量化和归一化的时钟相位
Figure C200480019738D0011100527QIETU
(t),假设时钟信号的触发边缘对应于时钟相位
Figure C200480019738D00121
(以2π为模)。
假设ηg=fgt,ηc=φc/2π+fct分别对应选通时钟和采集时钟的非量化归一化相位。将选通时钟的初始相位假设为φg=0并不意味着丧失一般性,因为所关注的ηg和ηc之间的相位关系已由φc覆盖了。注意的是,由于fg=1/Tg,同样ηg可被看作归一化时间ηg=t/Tg。代入得到
ηg=φc/2π+fcηg/fg
附图3示出了ηcg)的实例,以及量化后的
Figure C200480019738D00122
附图3示出了采集时钟的归一化相位和选通时钟的归一化相位的函数关系。采集计数器的当前值可用
Figure C200480019738D00123
表示。传统频率测量电路依赖于来自采集计数器的两个取样
Figure C200480019738D00124
Figure C200480019738D00125
虚曲线表示φc的影响,恒定相位偏移ηc:在黑曲线的情况下,将φc选择成使得φc的任何进一步减少将导致从5变化为4。在虚曲线的情况下,将φc被选择成使得φc的任何进一步减少将导致
Figure C200480019738D00127
从16变化为15。注意的是,在上述两种情况下, &eta; ^ c [ k 0 ] = 5 适用于这两种情况。
因此,频率测量转变为估算相位斜度ηcg)的陡度的问题。如果使用非量化归一化相位,则正确的解为:
fc/fg≈(η[k1]-ηc[k0])/(k1-k0)
其中,k0=ηg(t0)和k1=ηg(t1)可任意选择。按照本发明,进行的是近似计算。
直接的方法是传统方法,也就是从采集计数器中取得取样
Figure C200480019738D00129
Figure C200480019738D001210
作为ηc[k0]和ηc[k1]的近似值。这里,k0和k1是跨越测量间隔的选通计数器的第一个值和最后一个值,是第k个选通时钟周期起始点处的采集计数器的值。假设
Figure C200480019738D001212
可以通过在选通时钟的上升边缘取样采集计数器的当前值获得就足够了。
Figure C200480019738D001213
Figure C200480019738D001214
定义连接线ηc,该连接线可被认为ηc的近似值。因此,频率比例可估算为:
f c / f g &ap; ( &eta; ^ c [ k 1 ] - &eta; ^ c [ k 0 ] ) / ( k 1 - k 0 ) = N c / N g
附图3描述对应于选择k0=3和k1=7的
Figure C200480019738D001216
Figure C200480019738D001217
黑曲线和虚曲线表示ηc(t)=φc/2π+fct的初始相位φc如何影响估算值ηc(t)的陡度。取决于φc,Nc按照频率比例fc/fg的估算值太小或者太大而取15-5=10或者16-5=11的值。适用于任何频率比例fc/fg的总的规则是,Nc不能超过两个值,并且取值依赖于采集时钟和选通时钟之间的相位关系。
先前的讨论提出频率估算问题作为基于量化形式的ηc[k]估算采集时钟的非量化相位斜度ηc[k]的陡度的任务。在这个前提下,传统方法表现得不是最优,因为它只利用了来自
Figure C200480019738D00131
的两个取样,但是忽略了来自测量时间间隔内的信息:假设测量时间间隔持续Ng=k1-k0个选通时钟周期,我们可以利用选通时钟的每个上升边缘从采集计数器中得到总数为Ng+1个取样。
这些取样
Figure C200480019738D00132
在附图5中用黑点标记。它们分布在黑色回归线上。来自回归线的取样
Figure C200480019738D00134
用空点标记。按照本发明,就陡度而言,与相比,
Figure C200480019738D00136
看起来是更好的估算值。
由于实际原因,可能希望考虑来自的Ng+1个取样的子集。例如,如果Ng具有106的数量级,可能不希望基于如此大量的取样值计算回归线。如上所述,k表示选通时钟ηg的连续归一化相位的整数值。来自这些整数值的子集可表示为
ηg[m]
其中m是大小M的子集内的系数,1≤m<M。
子集的实例是等距栅格
ηg[m]=m[Ng/M]
其中整数[Ng/M]是栅格间距。
子集的另一个实例是均一抖动栅格
ηg[m]=m[Ng/M]+rnd[m]
此时,rnd[m]是R0≤rnd[m]≤R1范围内的均一分布的整数随机变量,其中R0,R1是整数。该仿真结果表明抖动栅格可能有益于某些格局。
假设取样
Figure C200480019738D00138
对应于某一子集ηg[m],回归线可如下确定:下面的前提是假设
&eta; &RightArrow; c [ &eta; g [ m ] ] = a 0 + a 1 &eta; g [ m ]
将恒定归一化相位偏移a0和陡度a1选择成使得误差信号
Figure C200480019738D001310
Figure C200480019738D001311
的能量E最小化。未知的a0和a1通过求解下面的等式而得出:
&PartialD; &PartialD; a 0 E ( a 0 , a 1 ) = &PartialD; &PartialD; a 0 &Sigma; m ( &eta; ^ c [ &eta; g [ m ] ] - a 0 - a 1 &eta; g [ m ] ) 2 = 0 - - - ( 3 )
&PartialD; &PartialD; a 1 E ( a 0 , a 1 ) &ap; &PartialD; &PartialD; a 1 &Sigma; m ( &eta; ^ c [ &eta; g [ m ] ] - a 0 - a 1 &eta; g [ m ] ) 2 = 0 - - - ( 4 )
因此产生线性等式系统:
B 1 B 2 = A 11 A 12 A 21 A 22 a 0 a 1 - - - ( 5 )
其中数据相关系数为:
B 1 = &Sigma; m &eta; ^ c [ &eta; g [ m ] ] - - - ( 6 )
B 2 = &Sigma; m &eta; ^ c [ &eta; g [ m ] ] &eta; g [ m ] - - - ( 7 )
和数据无关系数为:
A11=M(8)
A 12 = A 21 = &Sigma; m &eta; g [ m ] - - - ( 9 )
A 22 = &Sigma; m ( &eta; g [ m ] ) 2
首先计算等式(6)到(10),随后求解等式(5)得到a0和a1。由于系数A11、A12、A21、A22无关于(数据)取样
Figure C200480019738D00147
只要选择了某些适当的值M和ηg[m]之后,就足以将它们计算出来。仅有B1、B2是数据相关的,并且对于一组新的取样ηg[m],需要重新计算。
计算量与取样数M是成比例的,如果M很大,计算B1、B2的工作量支配求解等式系统(5)的工作量。
由于未知的恒定相位a0不是我们关注的,不要求对其进行明确计算。
早期精度较低的频率估算值可通过首先对M0<M个取样计算a1获得。当更多的取样到达时,由于(6)到(10)的累积特性,用较低的工作量就可计算出a1的更新值。
由于等式系统(5)的简单特性,可能简单地引用解a1。结果,可以将这个简单的解理解为与有限脉冲响应(FIR)滤波器的单一输出取样密切相关。(该滤波器有斜坡形状的脉冲响应,并且可看作通信理论中公知的匹配滤波器)。因此,可以选择某些适当的FIR滤波器实现方式。
如果要测量的频率比例随时间变化,人们可能希望对其进行重复测量。取决于更新率,一组新的M个取样可能或者可能不与前一组取样重叠。如果重叠,则在一定程度上降低计算量就是可能的。
可能已知频率比例在取得M个取样的测量时间间隔内变化。这可能是例如由于在启动各个晶体振荡器之后,两个频率之一的渐近指数稳定造成的。在这种情况下,我们可以估算更复杂回归线的参数。
附图5示出了按照本发明的电子电路的第二种示范性实施方式的简化电路图,该电路按照上面介绍的按照本发明的一个方面的原理进行操作。
从附图5可以看出,提供了一个计数器30和一个锁存器40。从附图5可以看出,将信号uc输入到计数器30的时钟输入端。然后,将输出信号xc1输入到锁存器40,锁存器40的时钟输入端接收信号ug。锁存器40的输出信号是信号xg1
而且,提供了计数器32和比较器34。计数器32的时钟输入端接收选通时钟信号ug。计数器32的输出信号输入到比较器34,每次计数器32输出的计数值达到n,则该比较器向锁存器40输出一个使能信号。
灰线36指示时钟域的过渡,即由采集时钟uc和选通时钟ug操作的电路单元之间的边界。
附图6示出了附图5的电路中出现的各个信号的时序图。从附图6可以看出,在理想条件的情况下,在附图5的电路中可能会出现问题。不过,由于计数器输出信号xc1的位不正确地同步变化,因此会出现问题。因此,如果选通时钟uc的取样边缘在计数器位变化的同时时出现,锁存器40的输出信号,即信号xg1,在某些情况下可能是错误的。
附图7示出了按照本发明的电路的第三种示范性实施方式的简化电路图。通过本发明的第三种示范性实施方式的电路,参考附图5和6描述的同步问题是可以避免的。
从附图7可以看出,将信号uc输入到计数器50的时钟端口,输出信号xc1输入到锁存器2 52。锁存器2 52的时钟输入端也接收信号uc。锁存器2 52的输出信号xc2输入到锁存器3 54,锁存器3 54的时钟输入端也与信号uc相接。
输入信号ug输入到AND门56的一个输入端,AND门56的输出信号输入到锁存器1 58。锁存器1 58的时钟端口也连接到输入信号uc。锁存器1 58的输出信号xc4通过一个反相器60返回到AND门56的另一个输入端。进一步的,信号xc4输入到另一个计数器62的使能端口,该计数器62由采集时钟uc进行时钟控制。计数器62的计数信号输入到比较器64,其中,计数信号与比较值n进行比较,以致于每次计数信号62的计数值达到n时,比较器64输出一个输出信号,该信号作为使能信号输入到锁存器354。时钟域过渡处于锁存器1 58处。
附图8示出了附图7的电路中出现的信号的时序图。附图8的第一个时序图描绘出随时间变化的信号uc。第二个时序图描绘出随时间变化的信号xc1。第三个时序图描绘出随时间变化的信号xc2。第四个时序图描绘出随时间变化的信号ug。第五个时序图描绘出随时间变化的信号xc4。第六个时序图描绘出随时间变化的信号xc3
从附图8可以看出,在附图5所示的电路中出现的同步问题,可通过附图7所示的电路得到避免,在附图7中,选通时钟ug被采集时钟uc取样,并且所有进一步的处理发生在采集时钟域中。“选通时钟检测”信号xc4在该电路中起了重要的角色。如果上升边缘出现在第n个采集时钟周期期间的选通时钟信号ug中,将导致xc4在采集时钟第n+1个周期期间变高,并且在第n+2个周期期间再次变低。因此,xc4可用于使能锁存器3 54,以接替xc2,其中xc2是xc1的延迟形式。
不过,应该注意,从信号xg1和xc3表示来自采集计数器的从头至尾轮换相同的计数器值序列的取样这个意义上说,附图5和7的电路是相同的。这是因为a)按照选通时钟ug上升边缘在第n个采集时钟周期期间的某一时刻在xc1的控制下取样采集时钟与b)在第n个采集时钟周期期间的某一时刻检测选通时钟ug的上升边缘并且之后检查时钟周期数目是相同的。
按照本发明的这些示范性实施方式的变形方式,为了定义所有可能的取样的子集(取样栅格),电路可通过比较器和保持km=ηg[m]、将会触发采样的选通时钟周期的数目的寄存器得到扩展。一旦触发条件出现并且取样可以获得,取样必须被存储起来,并且在下一个触发的时刻将其加载。在中断服务程序的控制下,可以存储和加载到FIFOs中和从FIFOs存储和加载或者可以存储和加载到处理器的存储器中和从处理器的存储器存储和加载。
很有优势地,按照本发明,测量误差明显减少,但不增加测量时间间隔。同样,在保持与现有技术已知的传统方法相同测量误差的情况下,测量间隔可减少。用于典型应用的仿真示出,如果可获得25个取样,测量误差可被减少到二分之一。取样数目增加到n倍,测量误差另外减少
Figure C200480019738D00171

Claims (7)

1.用于确定第一信号(S1)的第一频率和第二信号(S2)的第二频率的比例的电路,该电路包括:
第一计数器(2)和第二计数器(4);
取样装置(6),用于当第二计数器(4)达到预设的第二中间计数值时,对第一计数器(2)的第一中间计数值进行取样,从而第一计数器(2)是在第二计数器(4)的控制下被取样的;
其中第一和第二中间计数值形成多个第一和第二计数器(2,4)的中间计数值对;
其中,在第一中间计数值的取样期间,第一和第二计数器(2,4)继续计数;
计算单元(16),基于多个中间计数值对,确定第一和第二频率的比例,
其中计算单元(16)使用多于两个中间计数值对来确定第一和第二频率的比例;
其特征在于通过从与不同数量的取样相对应的中间计数值对进行回归运算来计算第一和第二频率的比例。
2.根据权利要求1的电路,其中第一计数器(2)是由第一信号(S1)的上升边缘和下降边缘之一触发的;其中第二计数器(4)是由第二信号(S2)的上升边缘和下降边缘之一触发的;其中第二计数器(4)的时钟信号是第一和第二信号(S1,S2)之一;其中在其值对第一计数器(2)进行取样的第二计数器(4)的第二中间计数值是预设在寄存器里的。
3.根据权利要求1的电路,进一步包括:一个存储器(12,14);其中该存储器包括第一和第二存储装置;其中第一存储装置(12)用于存储第一计数器(2)的第一中间计数值,以致提供第一计数器(2)的第一中间计数值的序列,而第二存储装置(14)用于存储第二计数器(4)的第二中间计数值,以致提供第二计数器(4)的第二中间计数值的序列。
4.根据权利要求1的电路,
其中计算单元(16)由处理器实现;
其中多个中间计数值对存储在处理器的工作存储器中;并且
其中针对读取和写入多个中间计数值对之一通过中断程序对工作存储器进行访问。
5.根据权利要求1的电路,其中计算单元(16)确定频率比例随时间的变化。
6.确定第一信号(S1)的第一频率和第二信号(S2)的第二频率的比例的方法,该方法包括以下步骤:
当第二计数器(4)达到预设的第二中间计数值时,对第一计数器(2)的第一中间计数值进行取样,从而第一计数器(2)是在第二计数器(4)的控制下得以取样的;
其中第一和第二中间计数值形成多个第一和第二计数器(2,4)的中间计数值对;
其中,在第一中间计数值的取样期间,第一和第二计数器(2,4)继续计数;并且
基于多个中间计数值对,确定第一和第二频率的比例;
其中使用多于两个中间计数值对来确定第一和第二频率的比例;
其特征在于通过从与不同数量的取样相对应的中间计数值对进行回归运算来计算第一和第二频率的比例。
7.包括计算机程序代码装置的计算机程序产品,其中当在处理器上运行该计算机程序代码装置时,该计算机程序代码装置促使处理器执行如下操作:
当第二计数器(4)达到预设的第二中间计数值时,对第一计数器(2)的第一中间计数值进行取样,从而第一计数器(2)是在第二计数器(4)的控制下得以取样的;
其中第一和第二中间计数值形成多个第一和第二计数器(2,4)的中间计数值对;
其中,在第一中间计数值的取样期间,第一和第二计数器(2,4)继续计数;并且
基于多个中间计数值对,确定第一和第二频率的比例;
其中使用多于两个中间计数值对来确定第一和第二频率的比例;
其特征在于通过从与不同数量的取样相对应的中间计数值对进行回归运算来计算第一和第二频率的比例。
CNB2004800197387A 2003-07-11 2004-07-05 改进的频率确定 Expired - Fee Related CN100520419C (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
EP03102128.0 2003-07-11
EP03102128 2003-07-11

Publications (2)

Publication Number Publication Date
CN1820203A CN1820203A (zh) 2006-08-16
CN100520419C true CN100520419C (zh) 2009-07-29

Family

ID=34042950

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2004800197387A Expired - Fee Related CN100520419C (zh) 2003-07-11 2004-07-05 改进的频率确定

Country Status (5)

Country Link
US (1) US20060155491A1 (zh)
EP (1) EP1646881A1 (zh)
JP (1) JP2007530917A (zh)
CN (1) CN100520419C (zh)
WO (1) WO2005005999A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104614586A (zh) * 2015-01-13 2015-05-13 上海惠安系统控制有限公司 一种基于gps的电气量信号自适应采样方法

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005074918A1 (en) * 2004-02-06 2005-08-18 Benzstrasse 1 D-61352 Bad Homburd The combination of anticholinergics and glucocorticoids for the long-term treatment of asthma and copd
AU2005210084B2 (en) 2004-02-06 2010-06-24 Meda Pharma Gmbh & Co. Kg Combination of anticholinergics and inhibitors of phosphodiesterase type 4 for the treatment of respiratory diseases
JP2008516253A (ja) * 2004-10-11 2008-05-15 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 非線形周波数及び位相測定方式
PL1971369T3 (pl) 2005-12-21 2010-01-29 Meda Pharma Gmbh & Co Kg Połączenie R,R-glikopirolanu, rolipramu i budesonidu dla potrzeb leczenia chorób zapalnych
US7936162B2 (en) * 2008-01-23 2011-05-03 Anritsu Corporation Measured-signal repetition frequency detection method, and sampling apparatus and waveform observation system using the method
US7881895B2 (en) * 2008-05-27 2011-02-01 Sony Ericsson Mobile Communications Ab Methods of calibrating a clock using multiple clock periods with a single counter and related devices and methods
US8560875B2 (en) * 2009-09-17 2013-10-15 Avago Technologies General Ip (Singapore) Pte. Ltd. Apparatus for clock calibrating a less precise second clock signal with a more precise first clock signal wherein the first clock signal is inactive during a sniff mode and the second clock signal is active during a sniff mode
CN101968510B (zh) * 2010-09-21 2013-04-17 天津农学院 谐振式传感器输出信号频率的高分辨率测量电路及其方法
CN102901870B (zh) * 2012-08-30 2014-12-17 河南科技大学 基于单片机的频率测试系统及使用上述系统的测试方法
DE102013109818B4 (de) * 2013-09-09 2015-05-13 Intel IP Corporation Datenverarbeitungsvorrichtung und Verfahren zum Bestimmen eines Taktverhältnisses
CN112732445A (zh) * 2021-01-14 2021-04-30 联想(北京)有限公司 一种数据获取方法、装置及电子设备

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE2921899A1 (de) * 1979-05-30 1980-12-04 Nord Micro Elektronik Feinmech Verfahren zur frequenzmessung
US5107439A (en) * 1990-11-09 1992-04-21 Hewlett-Packard Company Continuous overlapping frequency measurement
JPH04344476A (ja) * 1991-05-21 1992-12-01 Tokyo Electric Power Co Inc:The 周波数比測定回路
JPH0677837A (ja) * 1992-08-24 1994-03-18 Hitachi Cable Ltd サンプリング周波数変換に基づくデータ出力装置における変換周波数比演算装置
JP3252581B2 (ja) * 1994-01-26 2002-02-04 ソニー株式会社 標本化周波数変換装置
JP3464913B2 (ja) * 1998-05-25 2003-11-10 松下電器産業株式会社 サンプルレートコンバータ
JP2002162427A (ja) * 2000-11-24 2002-06-07 Yamaha Corp 周波数測定回路及びそれを用いたデジタル信号処理装置
JP2002228737A (ja) * 2001-02-01 2002-08-14 Mitsubishi Electric Corp 測位情報受信装置および測位情報受信方法
US6597205B2 (en) * 2001-12-21 2003-07-22 Honeywell International Inc. High accuracy method for determining the frequency of a pulse input signal over a wide frequency range
US7130368B1 (en) * 2002-09-19 2006-10-31 Nortel Network Limited Clock recovery using a direct smoothing process

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104614586A (zh) * 2015-01-13 2015-05-13 上海惠安系统控制有限公司 一种基于gps的电气量信号自适应采样方法
CN104614586B (zh) * 2015-01-13 2019-01-18 上海惠安系统控制有限公司 一种基于gps的电气量信号自适应采样方法

Also Published As

Publication number Publication date
EP1646881A1 (en) 2006-04-19
JP2007530917A (ja) 2007-11-01
WO2005005999A1 (en) 2005-01-20
CN1820203A (zh) 2006-08-16
US20060155491A1 (en) 2006-07-13

Similar Documents

Publication Publication Date Title
CN100451666C (zh) 改进的时域信号频率稳定度测量方法和装置
CN100520419C (zh) 改进的频率确定
US4875201A (en) Electronic pulse time measurement apparatus
CN101976037B (zh) 一种多次同步模拟内插的时间间隔测量方法和装置
CN102928677A (zh) 一种纳米级脉冲信号采集方法
EP0177557B1 (en) Counting apparatus and method for frequency sampling
US7792650B2 (en) Edge-aligned ratio counter
CN102073268A (zh) 一种高精度脉冲时间间隔测量电路
CN1123775C (zh) 一种频率、周期测量方法及装置
US9568889B1 (en) Time to digital converter with high resolution
CN109283354A (zh) 一种基于增量式光电编码器的变m/t测速方法
CN106597097A (zh) 一种高精度频率测量方法
EP2499741B1 (en) Time-to-digital converter with successive measurements
EP0484975A2 (en) Continuous overlapping frequency measurement
TW202101908A (zh) 用於產生具有可控脈衝寬度之脈衝輸出的電路
US6944099B1 (en) Precise time period measurement
CN104569582B (zh) 一种用于实现频率量测的方法及fpga电路
US4392749A (en) Instrument for determining coincidence and elapse time between independent sources of random sequential events
CN112152596B (zh) 用于产生脉冲输出的电路及方法
RU197391U1 (ru) Цифровой частотомер
RU2173857C1 (ru) Способ измерения частоты импульсов
SU1695328A1 (ru) Статистический анализатор
RU2055395C1 (ru) Знаковый коррелометр
SU1118933A1 (ru) Цифровой фазовый детектор
RU2517783C1 (ru) Цифровой измеритель частоты

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
ASS Succession or assignment of patent right

Owner name: NXP CO., LTD.

Free format text: FORMER OWNER: KONINKLIJKE PHILIPS ELECTRONICS N.V.

Effective date: 20071012

C41 Transfer of patent application or patent right or utility model
TA01 Transfer of patent application right

Effective date of registration: 20071012

Address after: Holland Ian Deho Finn

Applicant after: Koninkl Philips Electronics NV

Address before: Holland Ian Deho Finn

Applicant before: Koninklijke Philips Electronics N.V.

C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20090729

Termination date: 20120705