CN102739238A - 计数装置以及计数方法 - Google Patents
计数装置以及计数方法 Download PDFInfo
- Publication number
- CN102739238A CN102739238A CN2012100922645A CN201210092264A CN102739238A CN 102739238 A CN102739238 A CN 102739238A CN 2012100922645 A CN2012100922645 A CN 2012100922645A CN 201210092264 A CN201210092264 A CN 201210092264A CN 102739238 A CN102739238 A CN 102739238A
- Authority
- CN
- China
- Prior art keywords
- service cycle
- counting
- value
- signal
- less
- 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
- 238000000034 method Methods 0.000 title claims description 12
- 238000009826 distribution Methods 0.000 claims abstract description 45
- 230000000630 rising effect Effects 0.000 claims abstract description 23
- 238000004364 calculation method Methods 0.000 claims description 21
- 230000007423 decrease Effects 0.000 claims description 11
- 229910052739 hydrogen Inorganic materials 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims 2
- 230000010355 oscillation Effects 0.000 description 21
- 239000004065 semiconductor Substances 0.000 description 15
- 238000003860 storage Methods 0.000 description 8
- 230000015572 biosynthetic process Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 4
- 230000008034 disappearance Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000033228 biological regulation Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000036962 time dependent Effects 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000001915 proofreading effect Effects 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03K—PULSE TECHNIQUE
- H03K21/00—Details of pulse counters or frequency dividers
- H03K21/02—Input circuits
- H03K21/023—Input circuits comprising pulse shaping or differentiating circuits
Landscapes
- Optical Radar Systems And Details Thereof (AREA)
- Instruments For Measurement Of Length By Optical Means (AREA)
Abstract
本发明提供一种即使输入信号在时间上非对称也可以对计数结果进行高精度地校正的计数装置。计数器(13)能对二值化信号的运转周期进行计数。计数结果校正部(14)分别对从信号的上升到下降为止的第1运转周期及从信号的下降到上升为止的第2运转周期,制作运转周期的频数分布,求出第1运转周期的代表值TH的0倍以上且不到1倍的长度的第1运转周期的数量的总和、第2运转周期的代表值TL的0倍以上且不到1倍的长度的第2运转周期的数量的总和、{TH+(n-0.5)×(TH+TL)}以上且不到{TH+(n+0.5)×(TH+TL)}的长度的第1运转周期的数量的总和、以及{TL+(n-0.5)×(TH+TL)}以上且不到{TL+(n+0.5)×(TH+TL)}的长度的第2运转周期的数量的总和,从而对计数结果进行校正。
Description
技术领域
本发明涉及一种对信号的数量进行计数的计数装置以及计数方法。
背景技术
以往,提出了一种使用半导体激光器的自混合效应的波长调制型激光测量器(参照专利文献1)。该激光测量器的构成如图9所示。图9的激光测量器包括:半导体激光器201,其向物体210发射激光;光电二极管202,其将半导体激光器201的光输出变换为电信号;透镜203,其将来自半导体激光器201的光聚光并照射到物体210,并将从物体210返回的光聚光使其入射到半导体激光器201中;激光驱动器204,其使半导体激光器201的振荡波长连续地增加的第1振荡期间与振荡波长连续地减少的第2振荡期间交替反复;电流-电压变换放大部205,其将光电二极管202的输出电流变换为电压并放大;信号提取电路206,其将电流-电压变换放大部205的输出电压进行2次微分;计数装置207,其对信号提取电路206的输出电压中所含有的跳模脉冲(以下,记为MHP)的数量计数;运算装置208,其计算出与物体210的距离以及物体210的速度;以及,显示装置209,显示运算装置208的算出结果。
激光驱动器204将相对于时间以一定的变化率反复增减的三角波驱动电流作为注入电流提供给半导体激光器201。由此,半导体激光器201被驱动为振荡波长以一定的变化率连续地增加的第1振荡期间与振荡波长以一定的变化率连续地减少的第2振荡期间交替反复。图10是示出半导体激光器201的振荡波长随时间变化的图。在图10中,P1是第1振荡期间,P2是第2振荡期间,λa是各期间的振荡波长的最小值,λb是各期间的振荡波长的最大值,Tt是三角波的周期。
从半导体激光器201出射的激光通过透镜203被聚光之后入射到物体210中。由物体210反射的光通过透镜203被聚光之后入射到半导体激光器201中。光电二极管202将半导体激光器201的光输出变换为电流。电流-电压变换放大部205将光电二极管202的输出电流变换为电压并放大,信号提取电路206对电流-电压变换放大部205的输出电压进行2次微分。计数装置207对信号提取电路206的输出电压中所含有的MHP的数量分别就第1振荡期间P1和第2振荡期间P2进行计数。基于半导体激光器1的最小振荡波长λa、最大振荡波长λb、第1振荡期间P1的MHP的数量以及第2振荡期间P2的MHP的数量,运算装置208计算出与物体210的距离以及物体210的速度。利用这样地自混合型的激光测量器的技术对MHP的数量进行测定的话,能够从该MHP数量计算出物体的振动频率。
采用上文所述的激光测量器的话,会将例如干扰光等的噪音作为MHP进行计数、或者因为信号的缺失而导致MHP漏记,存在着计数装置所计数的MHP的数量产生误差,所计算出的距离、振动频率等的物理量也产生误差这样的问题。
鉴于此,发明者提出了一种计数装置,其测定计数期间中的MHP的周期,并根据测定结果做成计数期间中的周期的频数分布,根据频数分布计算出MHP的周期的代表值,根据频数分布求出代表值的第1规定数倍以下的等级的频数的总和Ns和代表值的第2规定数倍以上的等级的频数的总和Nw,基于这些频数Ns和Nw对MHP的计数结果进行校正,由此可以去除计数时的缺漏和过剩的计数的影响(参照专利文献2)。
根据专利文献2所公开的计数装置,只是在SN(Signal to Noise ratio)没有极端下降的情况下可以进行大致良好的校正。
但是,采用专利文献2中所公开的计数装置的话,在短距离测定中,信号强度和滞后宽度(ヒステリシス幅)相比非常强的情况下,由于频率高于MHP的噪音的影响,输入到计数装置的信号会在二值化的阈值附近产生颤振(チヤタリング),短周期的信号、MHP的原来周期的一半左右的周期的信号会多有发生。此时,由于比MHP的原来的周期短的周期成为周期的分布的代表值,因此,存在着不能对MHP的计数结果进行正确地校正,MHP的计数结果是原来的值的例如数倍大这样的问题。
这里,进一步提出了即使在输入信号中连续地产生有高频率的噪音的情况下,也可以对计数误差进行校正的计数装置的方案(参照专利文献3)。专利文献3所揭示的计数装置是如下所述的计数装置:能对计数期间中的输入信号的运转周期(ランレングス)的数量进行计数,测定计数期间中的输入信号的运转周期,根据该测定结构制作计数期间中的输入信号的运转周期的频数分布,根据该频数分布计算出输入信号的运转周期的分布的代表值,求出不到代表值的0.5倍的运转周期的数量的总和Ns和代表值的2n倍以上且不到(2n+2)倍(n为1以上的自然数)的运转周期的数量的总和Nwn,基于这些频数Ns和Nwn,对MHP的计数结果进行校正。
这里,MHP等的干涉波的波形根据传输波去除电路的特性、对象物的状态在时间上呈非对称的形状(参照专利文献4)。图11(A)是表示这样的非对称的干涉波形的图,图11(B)是表示将图11(A)的波形二值化后的结果的图。图11(A)中的TH1、TH2是用于二值化的阈值。干涉波形像这样在时间上非对称的话,则二值化信号的占空比不为0.5。因此,采用专利文献3中所揭示的计数装置,存在计数校正的精度降低这样的问题。
现有技术文献
专利文献
专利文献1日本特开2006-313080号公报
专利文献2日本特开2009-47676号公报
专利文献3日本特开2011-33525号公报
专利文献4日本专利第3282746号公报
发明内容
发明要解决的问题
根据专利文献3所揭示的计数装置,即使在输入信号连续地产生高频率的噪声的情况下,也可以对计数误差进行校正。
但是,采用专利文献3中所揭示的计数装置的话,存在干扰波形在时间上非对称的话,将干扰波形二值化了的的信号的占空比不为0.5,计数校正的精度降低这样的问题。
本发明为了解决上述课题而作出,其目的是提供即便在被输入至计数装置的信号在时间上非对称的情况下,也可以高精度地对计数结果进行校正的计数装置以及计数方法。
解决课题的手段
本发明提供一种计数装置,其在特定的物理量和信号的数量具有线性关系、且所述特定的物理量为一定的情况下,对为大致单一频率的所述信号进行计数,其特征在于,该计数装置包括:二值化单元,其将输入信号进行二值化;信号计数单元,其在一定的计数期间中,对从所述二值化单元输出的二值化信号的运转周期(行程长度)的数量进行计数;运转周期测定单元,其在每次输入信号的运转周期部分时,对所述计数期间中的所述二值化信号的运转周期进行测定;频数分布制作单元,其根据所述运转周期测定单元的测定结果分别对从二值化信号的上升到接下来的下降为止的第1运转周期及从二值化信号的下降到接下来的上升为止的第2运转周期,制作所述计数期间中的二值化信号的运转周期的频数分布;代表值算出单元,其根据所述第1运转周期的频数分布计算出所述第1运转周期的分布的代表值TH,并根据所述第2运转周期的频数分布计算出所述第2运转周期的分布的代表值TL;校正值算出单元,其根据所述运转周期测定单元的测定结果和所述代表值算出单元的计算结果,求出所述代表值TH的0倍以上且不到1倍的长度的第1运转周期的数量的总和NsH、所述代表值TL的0倍以上且不到1倍的长度的第2运转周期的数量的总和NsL、{TH+(n-0.5)×(TH+TL)}以上且不到{TH+(n+0.5)×(TH+TL)}的长度的第1运转周期的数量的总和NwnH和{TL+(n-0.5)×(TH+TL)}以上且不到{TL+(n+0.5)×(TH+TL)}的长度的第2运转周期的数量的总和NwnL,基于该些频数NsH、NsL、NwnH、NwnL对所述计数单元的计数结果进行校正,以计算出所述输入信号的数量,其中,n为1以上的自然数。
又,在本发明的计数装置的1个构成例中,将所述信号计数单元的计数结果设为N,将所述第1运转周期与所述第2运转周期的和能取得的最大值设为(TH+TL)max时,所述校正值算出单元按照下式
求出校正后的计数结果N’。
又,在本发明的计数装置的1个构成例中,所述代表值TH、TL是中央值、最频值、平均值、等级值和频数的积为最大时的等级值、等级值的a次幂与频数的积为最大时的等级值中的任意一个,其中,0<a<1。
又,在本发明的计数装置的1个构成例中,求出所述第1运转周期的数量的总和NsH的阈值是所述代表值TH的0倍以上且不到0.5倍,或者是所述代表值TH的0倍以上且不到(TH+TL)/4;求出所述第2运转周期的数量的总和NsL的阈值是所述代表值TL的0倍以上且不到0.5倍,或者是所述代表值TL的0倍以上且不到(TH+TL)/4。
本发明提供一种计数方法,其在特定的物理量和信号的数量具有线性关系、且所述特定的物理量为一定的情况下,对为大致单一频率的所述信号进行计数,其特征在于,该计数方法包括:二值化步骤,其将输入信号进行二值化;信号计数步骤,其在一定的计数期间中,对在所述二值化步骤所得到的二值化信号的运转周期的数量进行计数;运转周期测定步骤,其在每次输入信号的运转周期部分时,对所述计数期间中的所述二值化信号的运转周期进行测定;频数分布制作步骤,其根据所述运转周期测定步骤的测定结果分别对从二值化信号的上升到接下来的下降为止的第1运转周期及从二值化信号的下降到接下来的上升为止的第2运转周期,制作所述计数期间中的二值化信号的运转周期的频数分布;代表值算出步骤,其根据所述第1运转周期的频数分布计算出所述第1运转周期的分布的代表值TH,并根据所述第2运转周期的频数分布计算出所述第2运转周期的分布的代表值TL;校正值算出步骤,其根据所述运转周期测定步骤的测定结果和所述代表值算出步骤的计算结果,求出所述代表值TH的0倍以上且不到1倍的长度的第1运转周期的数量的总和NsH、所述代表值TL的0倍以上且不到1倍的长度的第2运转周期的数量的总和NsL、{TH+(n-0.5)×(TH+TL)}以上且不到{TH+(n+0.5)×(TH+TL)}的长度的第1运转周期的数量的总和NwnH和{TL+(n-0.5)×(TH+TL)}以上且不到{TL+(n+0.5)×(TH+TL)}的长度的第2运转周期的数量的总和NwnL,基于该些频数NsH、NsL、NwnH、NwnL对所述计数步骤的计数结果进行校正,以计算出所述输入信号的数量,其中,n为1以上的自然数。
又,在本发明的计数方法的1个构成例中,将所述信号计数步骤的计数结果设为N,将所述第1运转周期与所述第2运转周期的和能取得的最大值设为(TH+TL)max时,所述校正值算出步骤按照下式
求出校正后的计数结果N’。
发明效果
根据本发明,将输入信号二值化,能在一定的计数期间中对从二值化单元所输出的二值化信号的运转周期的数量进行计数,测定计数期间中的二值化信号的运转周期,根据该测定结果分别对从二值化信号的上升到接下来的下降为止的第1运转周期及从二值化信号的下降到接下来的上升为止的第2运转周期,制作计数期间中的二值化信号的运转周期的频数分布,从第1运转周期的频数分布计算出第1运转周期的代表值TH,并从第2运转周期的频数分布计算出第2运转周期的分布的代表值TL,求出代表值TH的0倍以上且不到1倍的长度的第1运转周期的数量的总和NsH、代表值TL的0倍以上且不到1倍的长度的第2运转周期的数量的总和NsL、{TH+(n-0.5)×(TH+TL)}以上且不到{TH+(n+0.5)×(TH+TL)}的长度的第1运转周期的数量的总和NwnH、以及{TL+(n-0.5)×(TH+TL)}以上且不到{TL+(n+0.5)×(TH+TL)}的长度的第2运转周期的数量的总和NwnL,通过基于该些频数NsH、NsL、NwnH、NwnL对信号计数单元的计数结果进行校正,即使在被输入至计数装置的信号连续地产生高频率的噪音,且输入信号的波形在时间上不对称的情况下,也可以对计数误差进行高精度地校正。
附图说明
图1是表示本发明的实施形态所涉及的计数装置的构成的框图。
图2是表示本发明的实施形态所涉及的计数装置的动作的流程图。
图3是表示本发明的实施形态所涉及的计数装置的计数结果校正部的构成的一个例子的框图。
图4是对本发明的实施形态所涉及的计数装置的计数的动作进行说明的图。
图5是对本发明的实施形态所涉及的计数装置的运转周期测定部的动作进行说明的图。
图6是表示在被输入至计数装置的信号中混入了高频率噪音时的MHP的周期的频数分布的一个例子的图。
图7是表示MHP的运转周期的频数分布的一个例子的图。
图8是对本发明的实施形态所涉及的计数装置的计数结果的校正原理进行说明的图。
图9是表示以往的激光计测器的构成的框图。
图10是表示在图9的激光计测器中的半导体激光器的振荡波长随时间变化的一个例子的图。
图11是表示时间上非对称的干扰波形的图。
符号的说明
1…计数装置、11…二值化部、12…逻辑乘运算部、13…计数器、14…计数结果校正部、15…存储部、140…运转周期测定部、141…频数分布制作部、142…代表制算出部、143…校正值算出部。
具体实施方式
第1实施形态
以下,参照附图对本发明的实施形态进行说明。图1是表示本发明的实施形态所涉及的计数装置的构成的框图,图2是表示计数装置的动作的流程图。计数装置1包括:二值化部11、逻辑乘运算部(AND)12、计数器13、计数结果校正部14以及存储部15。计数器13构成信号计数单元。
图3是示出计数结果校正部14的构成的一个实例的框图。计数结果校正部14包括:运转周期测定部140、频数分布生成部141、代表值算出部142以及校正值算出部143。
以下,在本实施形态中,将计数装置1适用于图9中示出那样的自混合型激光计测器,以对作为自混合信号的跳模脉冲(MHP)的数量进行计数的情况为例子进行说明。
图4(A)~图4(D)是对计数器13的动作进行说明的图,图4(A)是示意性地示出向计数装置1的输入信号的波形、即MHP的波形的图,图4(B)是示出与图4(A)对应的二值化部11的输出的图,图4(C)是示出输入到计数装置1中的门信号GS的图,图4(D)是示出与图4(B)对应的计数器13的计数结果的图。
首先,计数装置1的二值化部11对图4(A)所示的输入信号是高电平(H)还是低电平(L)进行判断,并输出如图4(B)所示的判断结果。此时,在输入信号的电压上升达到阈值TH1以上时,二值化部11判断为高电平,输入信号的电压下降达到阈值TH2(TH2<TH1)以下时,二值化部11判断为低电平,由此将输入信号二值化。
AND12输出二值化部11的输出与如图4(C)所示的门信号GS的逻辑乘运算的结果,计数器13对AND12的输出(二值化信号)的上升和下降进行计数(图4(D))。这里,门信号GS是在计数期间(例如,计数装置1适用于自混合型激光计测器的情况下,第1振荡期间P1或者第2振荡期间P2)的起点上升、在计数期间的终点下降的信号。因此,计数器13计数的是计数期间中的AND12的输出的上升沿的数量和下降沿的数量(即MHP的运转周期的数量)(图2的步骤S100)。
图5是对计数结果校正部14的运转周期测定部140的动作进行说明的图。运转周期测定部140对计数期间中的MHP的运转周期进行测定(图2的步骤S101)。即,运转周期测定部140通过将计数期间中的AND12的输出与阈值TH3进行比较,检测AND12的输出的上升,并通过将AND12的输出与阈值TH4进行比较,检测AND12的输出的下降。而且,运转周期测定部140通过对从AND12的输出的上升到接下来的下降为止的时间tud及从AND12的输出的下降到接下来的上升为止的时间tdu进行测定,来对计数期间中的AND12的输出的运转周期(即,MHP的运转周期)进行测定。这样地,MHP的运转周期是时间tud或者tdu。运转周期测定部140在每当检测出AND12的输出的上升或者下降时都进行上述的测定。
存储部15存储计数器13的计数结果和运转周期测定部140的测定结果。
门信号GS下降,计数期间结束之后,计数结果校正部14的频数分布制作部141根据存储于存储部15中的运转周期测定部140的测定结果制作计数期间中的MHP的运转周期tud、tdu的频数(度数)分布(图2步骤S102)。这时,频数分布制作部141分别对从MHP的上升到接下来的下降为止的第1运转周期tud及从MHP的下降到接下来的上升为止的第2运转周期tdu,制作计数期间中的MHP的运转周期的频数分布。
然后,计数结果校正部14的代表值算出部142根据频数分布制作部141制作了的第1运转周期tud的频数分布计算出第1运转周期tud的代表值TH,并根据频数分布制作部141制作了的第2运转周期tdu的频数分布计算出第2运转周期tdu的代表值TL(图2的步骤S103)。这里,代表值算出部142可以将第1运转周期tud的最频值、中央值或者平均值作为代表值TH,也可以将等级值和频数的积为最大时的等级值作为代表值TH,也可以将等级值的a次方(0<a<1)与频数的积为最大时的等级值作为代表值TH。关于第2运转周期tdu的代表制TL可以与TH相同地求得。表1示出频数分布的数值例以及该数值例中的等级值和频数的积。
表1
频数分布的数值例
等级值 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 |
频数 | 11 | 2 | 0 | 3 | 7 | 10 | 6 | 2 | 3 | 1 |
积 | 11 | 4 | 0 | 12 | 35 | 60 | 42 | 16 | 27 | 10 |
在表1的例子中,频数最大的最频值(等级值)是1。相对于此,等级值和频数的积为最大时的等级值是6,是与最频值不同的值。由于将等级值和频数的积为最大时的等级值作为代表值TH、TL的理由在专利文献3中被揭示,故省略其说明。代表值算出部142所计算出的代表值TH、TL存储于存储部15中。代表值算出部142在每次通过频数分布制作部141制作频数分布时进行这样的代表值TH、TL的计算。
计数结果校正部14的校正值算出部143根据运转周期测定部140的测定结果和代表值算出部142的计算结果,求出代表值TH的0倍以上且不到1倍的长度的第1运转周期tud的数量的总和NsH、在代表值TL的0倍以上且不到1倍的长度的第2运转周期tdu的数量的总和NsL、在{TH+(n-0.5)×(TH+TL)}以上且不到{TH+(n+0.5)×(TH+TL)}(n是1以上nmax以下的自然数)的长度的第1运转周期tud的数量的总和NwnH、和在{TL+(n-0.5)×(TH+TL)}以上且不到{TL+(n+0.5)×(TH+TL)}的长度的第2运转周期tdu的数量的总和NwnL,并按照下式对计数器13的计数结果进行校正(图2的步骤S104)。
【数式3】
在式(1)中,N是作为计数器13的计数结果的MHP的运转周期的数量,N’是校正后得到的MHP的数量,(TH+TL)max是第1运转周期tud和第2运转周期tdu的和能取得的最大值。求得第1运转周期tud的数量的总和NsH的阈值可以是代表值TH的0倍以上且不到0.5倍,也可以是代表值TH的0倍以上且不到(TH+TL)/4。同样地,求得第2运转周期tdu的数量的总和NsL的阈值可以是代表值TL的0倍以上且不到0.5倍,也可以是代表值TL的0倍以上且不到(TH+TL)/4。计数装置1在每个计数期间进行上述那样的处理。
接着,对计数装置1的计数器13的计数结果的校正原理进行说明。式(1)所示的计数结果的校正的基本原理与专利文献2所公开的计数结果的校正原理相同。但是,根据专利文献2所公开的校正原理,在输入到计数装置的信号中混入了频率高于MHP的突发噪音的情况下,无法正确地对计数器13的计数结果进行校正。
图6是表示在被输入至计数装置的信号中混入了高频率噪音的情况下,MHP的周期的频数分布的一个例子的图。在输入信号中混入了高频率的噪音的情况下,如图6所示,MHP的周期的频数分布中,除了在MHP的原来的周期Ta处出现具有频数的极大值的分布170之外,在周期Ta的大约一半的周期处还出现了具有频数的极大值的分布171、噪音的短周期172。而且,由于混入了高频的噪音,这些频数的极大值会有略微向时间短的一方移动的倾向。进一步地,有时高频的噪音会连续地混入。采用专利文献2所公开的现有的计数装置的话,混入这样的高频的连续噪音时,就不能够充分地对MHP的计数结果进行校正。关于这样的问题,在专利文献3中有详细地说明。
因此,在本实施形态中,不是使用MHP的周期的代表值,而是与专利文献3同样地使用MHP的运转周期的代表值T0对计数结果进行校正。图7示出MHP的运转周期的频数分布的一个实例。从图7可以清楚地看出,求出MHP的运转周期的频数分布的话,即使在被输入到计数装置1的信号中混入高频的噪音,在0.5T0附近也不会出现频数的极大值。即,由于求出运转周期的数量的总和NsH、NsL的阈值附近的频数的极大值消失了,因此可以正确地求出所述NsH、NsL,从而可以抑制校正的误差。
但是,采用专利文献3所揭示的计数装置的话,MHP的波形在时间上非对称的话,MHP的运转周期的频数分布不是图7、图8(A)中所示那样的形状,而是变成如图8(B)所示的在TH具有最频值的形状和如图8(C)所示的在TL具有最频值的形状。MHP的波形是时间上对称的情况下,由于噪音产生MHP的波形缺失的话,产生T0的奇数倍的运转周期。因此,在图8(A)的例子中,在3T0、5T0处产生频数的极大值。
一方面,在如图11(A)所示MHP的波形是时间上非对称的情况下,在MHP的波形产生缺失的话,是就产生在TH加上(TH+TL)的整数倍的值的长度的运转周期和在TL加上(TH+TL)的整数倍的值的长度的运转周期,而不是T0的奇数倍这样的运转周期。在MHP的波形产生缺失,2n+1个运转周期变为1个时的运转周期,在第1运转周期tud的情况下为TH+n×(TH+TL),在第2运转周期tdu的情况下为TL+n×(TH+TL)。由于在MHP中重叠有各种频率的噪音,因此运转周期成为以TH+n×(TH+TL)为中心的高斯分布和以TL+n×(TH+TL)为中心的高斯分布。因此,在图8(B)的例子中,在TH+(TH+TL)和TH+2(TH+TL)处产生频数的极大值,在图8(C)的例子中,在TL+(TH+TL)和TL+2(TH+TL)处产生频数的极大值。
这里,在本实施形态中,求出第1运转周期tud的数量的总和NwnH的阈值为{TH+(n-0.5)×(TH+TL)}以上且不到{TH+(n+0.5)×(TH+TL)},求出第2运转周期tdu的数量的总和NwnL的阈值为TL+(n-0.5)×(TH+TL)}以上且不到{TL+(n+0.5)×(TH+TL)}。因此,可以正确地求出(NwnH+NwnL),从而可以抑制校正的误差。
在图8(D)的例子中示出代表值TH的0倍以上且不到1倍的第1运转周期tud的数量的总和NsH和{TH+0.5×(TH+TL)}以上且不到{TH+1.5×(TH+TL)的第2运转周期tdu的数量的总和Nw1H;在图8(E)的例子中示出代表值TL的0倍以上且不到1倍的第2运转周期tdu的数量的总和NsL和{TL+0.5×(TH+TL)}以上且不到{TL+1.5×(TH+TL)}的第2运转周期tdu的数量的总和Nw1L。以上就是数式(1)所示的计数结果的校正原理。又,对数式(1)的右边乘以1/2倍的理由是为了将MHP的运转周期的数量变换为MHP的数量。
如上所述,在本实施形态中,通过计数器13中对计数期间中的MHP的运转周期的数量进行计数,测定计数期间中的MHP的运转周期,根据该测定结果分别对MHP的上升到接下来的下降为止的第1运转周期及从MHP的下降到接下来的上升为止的第2运转周期制作计数期间中的MHP的运转周期的频数分布,根据第1运转周期的频数分布计算出第1运转周期的分布的代表值TH,并根据第2运转周期的频数分布计算出第2运转周期的分布的代表值TL,求出代表值TH的0倍以上且不到1倍的长度的第1运转周期的数量的总和NsH、代表值TL的0倍以上且不到1倍的长度的第2运转周期的数量的总和NsL、{TH+(n-0.5)×(TH+TL)}以上且不到{TH+(n+0.5)×(TH+TL)}的长度的第1运转周期的数量的总和NwnH和{TL+(n-0.5)×(TH+TL)}以上且不到{TL+(n+0.5)×(TH+TL)}的长度的第2运转周期的数量的总和NwnL,通过基于该些频数NsH、NsL、NwnH、NwnL对计数器13的计数结果进行校正,即使在被输入至计数装置的信号中连续地产生比MHP高频率的噪音,且MHP的波形在时间上是不对称的情况下,也可以对MHP的计数误差进行高精度地校正。
另外,在本实施形态中,计数装置1可以通过包括例如CPU、存储装置以及接口的计算机和控制这些硬件资源的程序来实现。使这样的计算机动作的程序以存储在软盘、CD-ROM、DVD-ROM以及存储器卡等的存储介质中的状态被提供。CPU将读取的程序写入存储装置中,并根据该程序进行本实施形态所说明的处理。
另外,在本实施形态中,对将本发明的计数装置应用到自混合型的激光计测器的情况进行了说明,但并不限定于此,本发明的计数装置也可以应用于光电传感器等其他的领域。本发明的计数装置有效的情况下,作为计数对象的信号的数量与特定的物理量(在本实施形态的情况中,半导体激光器与物体的距离、及物体的位移)具有线性的关系,特定的物理量为一定的情形是信号为大约单一频率的情形。
而且,即使信号不是单一的频率,在特定的物理量和计数期间相比为足够低的频率、且如例如以1/10以下的频率振动的对象物的速度那样,周期分布的范围较小的情况下,则也可以作为大约单一的频率,这对本发明的计数装置来说是有效的。
在本实施形态中,作为适用计测装置的物理量传感器的实例,举出了如专利文献1所公开那样根据计数装置的计数结果计算出距物体的距离及物体的速度的激光计测器的实例进行了说明,但并不限定于此,本发明也可以应用于其他的物理量传感器。即,既可以根据计数装置的计数结果计算出物体的张力,也可以根据计数装置的计数结果计算出和物体的振动频率。物理量传感器所计算出的物理量是各种各样的,可见所述特定的物理量和物理量传感器计算出的物理量有相同的情况,也有不同的情况。
又,本发明中的输入信号是指连续的变化量(自混合的情况下,自混合信号)上的事件(イベント)和波动(自混合情况下的干扰条纹)。
产业上的可利用性
本发明可以应用于对信号的数量进行计数的计数装置。
Claims (8)
1.一种计数装置,其在特定的物理量和信号的数量具有线性关系、且所述特定的物理量为一定的情况下,对为大致单一频率的所述信号进行计数,其特征在于,该计数装置包括:
二值化单元,其将输入信号进行二值化;
信号计数单元,其在一定的计数期间中,对从所述二值化单元输出的二值化信号的运转周期(行程长度)的数量进行计数;
运转周期测定单元,其在每次输入信号的运转周期部分时,对所述计数期间中的所述二值化信号的运转周期进行测定;
频数分布制作单元,其根据所述运转周期测定单元的测定结果分别对从二值化信号的上升到接下来的下降为止的第1运转周期及从二值化信号的下降到接下来的上升为止的第2运转周期,制作所述计数期间中的二值化信号的运转周期的频数分布;
代表值算出单元,其根据所述第1运转周期的频数分布计算出所述第1运转周期的分布的代表值TH,并根据所述第2运转周期的频数分布计算出所述第2运转周期的分布的代表值TL;
校正值算出单元,其根据所述运转周期测定单元的测定结果和所述代表值算出单元的计算结果,求出所述代表值TH的0倍以上且不到1倍的长度的第1运转周期的数量的总和NsH、所述代表值TL的0倍以上且不到1倍的长度的第2运转周期的数量的总和NsL、{TH+(n-0.5)×(TH+TL)}以上且不到{TH+(n+0.5)×(TH+TL)}的长度的第1运转周期的数量的总和NwnH和{TL+(n-0.5)×(TH+TL)}以上且不到{TL+(n+0.5)×(TH+TL)}的长度的第2运转周期的数量的总和NwnL,基于该些频数NsH、NsL、NwnH、NwnL对所述计数单元的计数结果进行校正,以计算出所述输入信号的数量,其中,n为1以上的自然数。
2.如权利要求1所述的计数装置,其特征在于,
将所述信号计数单元的计数结果设为N,将所述第1运转周期与所述第2运转周期的和能取得的最大值设为(TH+TL)max时,所述校正值算出单元按照下式
求出校正后的计数结果N’。
3.如权利要求1或2所述的计数装置,其特征在于,
所述代表值TH、TL是中央值、最频值、平均值、等级值和频数的积为最大时的等级值、等级值的a次幂与频数的积为最大时的等级值中的任意一个,其中,0<a<1。
4.如权利要求1-3中的任一项所述的计数装置,其特征在于,
求出所述第1运转周期的数量的总和NsH的阈值是所述代表值TH的0倍以上且不到0.5倍,或者是所述代表值TH的0倍以上且不到(TH+TL)/4;
求出所述第2运转周期的数量的总和NsL的阈值是所述代表值TL的0倍以上且不到0.5倍,或者是所述代表值TL的0倍以上且不到(TH+TL)/4。
5.一种计数方法,其在特定的物理量和信号的数量具有线性关系、且所述特定的物理量为一定的情况下,对为大致单一频率的所述信号进行计数,其特征在于,该计数方法包括:
二值化步骤,其将输入信号进行二值化;
信号计数步骤,其在一定的计数期间中,对在所述二值化步骤所得到的二值化信号的运转周期的数量进行计数;
运转周期测定步骤,其在每次输入信号的运转周期部分时,对所述计数期间中的所述二值化信号的运转周期进行测定;
频数分布制作步骤,其根据所述运转周期测定步骤的测定结果分别对从二值化信号的上升到接下来的下降为止的第1运转周期及从二值化信号的下降到接下来的上升为止的第2运转周期,制作所述计数期间中的二值化信号的运转周期的频数分布;
代表值算出步骤,其根据所述第1运转周期的频数分布计算出所述第1运转周期的分布的代表值TH,并根据所述第2运转周期的频数分布计算出所述第2运转周期的分布的代表值TL;
校正值算出步骤,其根据所述运转周期测定步骤的测定结果和所述代表值算出步骤的计算结果,求出所述代表值TH的0倍以上且不到1倍的长度的第1运转周期的数量的总和NsH、所述代表值TL的0倍以上且不到1倍的长度的第2运转周期的数量的总和NsL、{TH+(n-0.5)×(TH+TL)}以上且不到{TH+(n+0.5)×(TH+TL)}的长度的第1运转周期的数量的总和NwnH和{TL+(n-0.5)×(TH+TL)}以上且不到{TL+(n+0.5)×(TH+TL)}的长度的第2运转周期的数量的总和NwnL,基于该些频数NsH、NsL、NwnH、NwnL对所述计数步骤的计数结果进行校正,以计算出所述输入信号的数量,其中,n为1以上的自然数。
6.如权利要求5所述的计数方法,其特征在于,
将所述信号计数步骤的计数结果设为N,将所述第1运转周期与所述第2运转周期的和能取得的最大值设为(TH+TL)max时,所述校正值算出步骤按照下式
求出校正后的计数结果N’。
7.如权利要求5或6所述的计数方法,其特征在于,
所述代表值TH、TL是中央值、最频值、平均值、等级值和频数的积为最大时的等级值、等级值的a次幂和频数的积为最大时的等级值中的任意一个,其中,0<a<1。
8.如权利要求5-7中的任一项所述的计数方法,其特征在于,
求出所述第1运转周期的数量的总和NsH的阈值是所述代表值TH的0倍以上且不到0.5倍,或者是所述代表值TH的0倍以上且不到(TH+TL)/4;
求出所述第2运转周期的数量的总和NsL的阈值是所述代表值TL的0倍以上且不到0.5倍,或者是所述代表值TL的0倍以上且不到(TH+TL)/4。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011083425A JP5818485B2 (ja) | 2011-04-05 | 2011-04-05 | 計数装置および計数方法 |
JP2011-083425 | 2011-04-05 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102739238A true CN102739238A (zh) | 2012-10-17 |
CN102739238B CN102739238B (zh) | 2015-06-24 |
Family
ID=46966138
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210092264.5A Expired - Fee Related CN102739238B (zh) | 2011-04-05 | 2012-03-31 | 计数装置以及计数方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US8599996B2 (zh) |
JP (1) | JP5818485B2 (zh) |
CN (1) | CN102739238B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105738074B (zh) * | 2016-03-07 | 2018-07-03 | 山信软件股份有限公司 | 一种光电管异常状态的检测方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1175129A (zh) * | 1996-08-19 | 1998-03-04 | 日本电气株式会社 | 二进制计数器读出电路 |
US20030034810A1 (en) * | 2001-08-14 | 2003-02-20 | Trivedi Pradeep R. | Non-integer division of frequency |
US20090084941A1 (en) * | 2004-04-26 | 2009-04-02 | Sony Corporation | Counter circuit, ad conversion method, ad converter, semiconductor device for detecting distribution of physical quantities, and electronic apparatus |
CN101988851A (zh) * | 2009-08-04 | 2011-03-23 | 株式会社山武 | 计数装置、物理量传感器、计数方法以及物理量测量方法 |
JP2011509568A (ja) * | 2007-12-20 | 2011-03-24 | クゥアルコム・インコーポレイテッド | 1つ又は複数のサイクル吸収クロック信号を生成又は利用する方法及び装置 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3282746B2 (ja) | 1993-04-09 | 2002-05-20 | 科学技術振興事業団 | 半導体レーザデジタル振動変位計測装置 |
JP5172077B2 (ja) | 2005-05-06 | 2013-03-27 | アズビル株式会社 | 距離・速度計および距離・速度計測方法 |
EP1953567A3 (en) * | 2007-01-25 | 2010-09-22 | Yamatake Corporation | Counting device, distance meter, counting method, and distance measuring method |
JP5545915B2 (ja) | 2007-01-25 | 2014-07-09 | アズビル株式会社 | 計数装置、距離計、計数方法および距離計測方法 |
JP5545916B2 (ja) * | 2008-01-31 | 2014-07-09 | アズビル株式会社 | 物理量センサおよび物理量計測方法 |
JP5663148B2 (ja) * | 2009-06-29 | 2015-02-04 | アズビル株式会社 | 計数装置、物理量センサ、計数方法および物理量計測方法 |
JP5548048B2 (ja) * | 2010-06-30 | 2014-07-16 | アズビル株式会社 | 信号判定装置および信号判定方法 |
JP5548055B2 (ja) * | 2010-07-07 | 2014-07-16 | アズビル株式会社 | 信号判定装置および信号判定方法 |
-
2011
- 2011-04-05 JP JP2011083425A patent/JP5818485B2/ja not_active Expired - Fee Related
-
2012
- 2012-03-15 US US13/420,683 patent/US8599996B2/en not_active Expired - Fee Related
- 2012-03-31 CN CN201210092264.5A patent/CN102739238B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1175129A (zh) * | 1996-08-19 | 1998-03-04 | 日本电气株式会社 | 二进制计数器读出电路 |
US20030034810A1 (en) * | 2001-08-14 | 2003-02-20 | Trivedi Pradeep R. | Non-integer division of frequency |
US20090084941A1 (en) * | 2004-04-26 | 2009-04-02 | Sony Corporation | Counter circuit, ad conversion method, ad converter, semiconductor device for detecting distribution of physical quantities, and electronic apparatus |
JP2011509568A (ja) * | 2007-12-20 | 2011-03-24 | クゥアルコム・インコーポレイテッド | 1つ又は複数のサイクル吸収クロック信号を生成又は利用する方法及び装置 |
CN101988851A (zh) * | 2009-08-04 | 2011-03-23 | 株式会社山武 | 计数装置、物理量传感器、计数方法以及物理量测量方法 |
Also Published As
Publication number | Publication date |
---|---|
JP5818485B2 (ja) | 2015-11-18 |
US8599996B2 (en) | 2013-12-03 |
CN102739238B (zh) | 2015-06-24 |
JP2012220234A (ja) | 2012-11-12 |
US20120257708A1 (en) | 2012-10-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101988851B (zh) | 计数装置、物理量传感器、计数方法以及物理量测量方法 | |
CN101936766B (zh) | 计数装置、物理量传感器、计数方法以及物理量测量方法 | |
CN102116775A (zh) | 速度测量装置以及方法 | |
EP3112799A1 (en) | Length measuring apparatus | |
CN102739238A (zh) | 计数装置以及计数方法 | |
CN102192707B (zh) | 物理量传感器和物理量计测方法 | |
US20140303936A1 (en) | Signal evaluating device and signal evaluating method | |
CN102944696A (zh) | 啁啾光纤光栅传感器及处理方法 | |
CN102062639B (zh) | 基于频数直方图测量脉冲激光全脉宽的方法 | |
EP2402782A2 (en) | Signal evaluating device and signal evaluating method | |
CN101398320B (zh) | 超声波流量计 | |
RU2359226C1 (ru) | Способ некогерентного накопления светолокационных сигналов | |
CN110045592A (zh) | 时间校正方法、装置、系统及计算机存储介质 | |
CN103023486A (zh) | 一种高速脉冲的计数方法 | |
JP4821240B2 (ja) | 流体の流れ計測装置 | |
JP5718410B2 (ja) | 計数装置、物理量センサ、計数方法および物理量計測方法 | |
JP5551378B2 (ja) | 張力・速度計測装置および方法 | |
RU2659466C1 (ru) | Преобразователь частоты следования импульсов в код | |
RU1840978C (ru) | Измеритель девиации | |
WO2018078727A1 (ja) | 流体測定装置 | |
Kandiros | Moment and Chernoff bounds | |
JP2020071079A (ja) | 超音波流量計および超音波流量計におけるピーク高さ情報取得方法 | |
JP5969198B2 (ja) | 超音波式流量計 | |
JP5097293B2 (ja) | 超音波式メータ装置 | |
CN102478610A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150624 Termination date: 20170331 |