JP5545913B2 - Counting device, distance meter, counting method, and distance measuring method - Google Patents

Counting device, distance meter, counting method, and distance measuring method Download PDF

Info

Publication number
JP5545913B2
JP5545913B2 JP2007259549A JP2007259549A JP5545913B2 JP 5545913 B2 JP5545913 B2 JP 5545913B2 JP 2007259549 A JP2007259549 A JP 2007259549A JP 2007259549 A JP2007259549 A JP 2007259549A JP 5545913 B2 JP5545913 B2 JP 5545913B2
Authority
JP
Japan
Prior art keywords
period
counting
frequency
representative value
procedure
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
JP2007259549A
Other languages
Japanese (ja)
Other versions
JP2009085910A (en
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.)
Azbil Corp
Original Assignee
Azbil Corp
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 Azbil Corp filed Critical Azbil Corp
Priority to JP2007259549A priority Critical patent/JP5545913B2/en
Publication of JP2009085910A publication Critical patent/JP2009085910A/en
Application granted granted Critical
Publication of JP5545913B2 publication Critical patent/JP5545913B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Instruments For Measurement Of Length By Optical Means (AREA)
  • Length Measuring Devices By Optical Means (AREA)
  • Optical Radar Systems And Details Thereof (AREA)

Description

本発明は、信号の数を数える計数装置、および計数装置を用いて干渉波形の数を測定し測定対象との距離を求める干渉型の距離計に関するものである。   The present invention relates to a counting device that counts the number of signals, and an interference type distance meter that determines the distance from a measurement object by measuring the number of interference waveforms using the counting device.

レーザによる光の干渉を利用した距離計測は、非接触測定のため測定対象を乱すことなく、高精度の測定方法として古くから用いられている。最近では、半導体レーザは装置の小型化のため、光計測用光源として利用されようとしている。その代表的な例として、FMヘテロダイン干渉計を利用したものがある。これは、比較的長距離測定が可能で精度もよいが、半導体レーザの外部に干渉計を用いているため、光学系が複雑になるという欠点を有する。   Distance measurement using light interference by a laser has long been used as a highly accurate measurement method without disturbing the measurement object for non-contact measurement. Recently, a semiconductor laser is being used as a light source for optical measurement in order to reduce the size of the apparatus. A typical example is one using an FM heterodyne interferometer. This is capable of relatively long distance measurement and good accuracy, but has the disadvantage that the optical system becomes complicated because an interferometer is used outside the semiconductor laser.

これに対して、レーザの出力光と測定対象からの戻り光との半導体レーザ内部での干渉(自己結合効果)を利用した計測器が提案されている(例えば、非特許文献1、非特許文献2、非特許文献3参照)。このような自己結合型のレーザ計測器によれば、フォトダイオード内蔵の半導体レーザが発光、干渉、受光の各機能を兼ねているため、外部干渉光学系を大幅に簡略化することができる。したがって、センサ部が半導体レーザとレンズのみとなり、従来のものに比べて小型となる。また、三角測量法より距離測定範囲が広いという特徴を有する。   On the other hand, a measuring instrument using interference (self-coupling effect) in the semiconductor laser between the laser output light and the return light from the measurement object has been proposed (for example, Non-Patent Document 1, Non-Patent Document). 2, see Non-Patent Document 3). According to such a self-coupled laser measuring instrument, the semiconductor laser with a built-in photodiode serves as the functions of light emission, interference, and light reception, so that the external interference optical system can be greatly simplified. Therefore, the sensor unit is only a semiconductor laser and a lens, and is smaller than the conventional one. In addition, the distance measurement range is wider than the triangulation method.

FP型(ファブリペロー型)半導体レーザの複合共振器モデルを図25に示す。図25において、101は半導体レーザ、102は半導体結晶の壁開面、103はフォトダイオード、104は測定対象である。測定対象104からの反射光の一部が発振領域内に戻り易い。戻って来たわずかな光は、共振器101内のレーザ光と結合し、動作が不安定となり雑音(複合共振器ノイズまたは戻り光ノイズ)を生じる。戻り光による半導体レーザの特性の変化は、出力光に対する相対的な戻り光量が、極めてわずかであっても顕著に現れる。このような現象は、ファブリペロー型(以下、FP型)半導体レーザに限らず、Vertical Cavity Surface Emitting Laser型(以下、VCSEL型)、Distributed FeedBack laser型(以下、DFBレーザ型)など、他の種類の半導体レーザにおいても同様に現れる。   FIG. 25 shows a composite resonator model of an FP type (Fabry-Perot type) semiconductor laser. In FIG. 25, 101 is a semiconductor laser, 102 is a wall opening of a semiconductor crystal, 103 is a photodiode, and 104 is an object to be measured. Part of the reflected light from the measurement object 104 easily returns to the oscillation region. The small amount of light that has returned returns to the laser beam in the resonator 101, and the operation becomes unstable, causing noise (composite resonator noise or return light noise). The change in the characteristics of the semiconductor laser due to the return light appears remarkably even if the amount of return light relative to the output light is very small. Such a phenomenon is not limited to a Fabry-Perot type (hereinafter referred to as FP type) semiconductor laser, but also other types such as a vertical cavity surface emitting laser type (hereinafter referred to as a VCSEL type) and a distributed fed back laser type (hereinafter referred to as a DFB laser type). This also appears in the same semiconductor laser.

レーザの発振波長をλ、測定対象104に近い方の壁開面102から測定対象104までの距離をLとすると、以下の共振条件を満足するとき、戻り光と共振器101内のレーザ光は強め合い、レーザ出力がわずかに増加する。
L=nλ/2 ・・・(1)
式(1)において、nは整数である。この現象は、測定対象104からの散乱光が極めて微弱であっても、半導体レーザの共振器101内の見かけの反射率が増加することにより、増幅作用が生じ、十分観測できる。
If the oscillation wavelength of the laser is λ and the distance from the wall open surface 102 closer to the measurement target 104 to the measurement target 104 is L, the return light and the laser light in the resonator 101 are as follows when the following resonance condition is satisfied. Strengthen and slightly increase the laser power.
L = nλ / 2 (1)
In formula (1), n is an integer. This phenomenon can be sufficiently observed even if the scattered light from the measurement object 104 is very weak, because the apparent reflectance in the resonator 101 of the semiconductor laser increases, causing an amplification effect.

半導体レーザは、注入電流の大きさに応じて周波数の異なるレーザ光を放射するので、発振周波数を変調する際に、外部変調器を必要とせず、注入電流によって直接変調が可能である。図26は、半導体レーザの発振波長をある一定の割合で変化させたときの発振波長とフォトダイオード103の出力波形との関係を示す図である。式(1)に示したL=nλ/2を満足したときに、戻り光と共振器101内のレーザ光の位相差が0°(同位相)になって、戻り光と共振器101内のレーザ光とが最も強め合い、L=nλ/2+λ/4のときに、位相差が180°(逆位相)になって、戻り光と共振器101内のレーザ光とが最も弱め合う。そのため、半導体レーザの発振波長を変化させていくと、レーザ出力が強くなるところと弱くなるところとが交互に繰り返し現れ、このときのレーザ出力を共振器101に設けられたフォトダイオード103で検出すると、図26に示すように一定周期の階段状の波形が得られる。このような波形は一般的には干渉縞と呼ばれる。   Since the semiconductor laser emits laser beams having different frequencies according to the magnitude of the injection current, an external modulator is not required when modulating the oscillation frequency, and direct modulation is possible by the injection current. FIG. 26 is a diagram showing the relationship between the oscillation wavelength and the output waveform of the photodiode 103 when the oscillation wavelength of the semiconductor laser is changed at a certain rate. When L = nλ / 2 shown in Expression (1) is satisfied, the phase difference between the return light and the laser light in the resonator 101 becomes 0 ° (same phase), and the return light and the resonator 101 When L = nλ / 2 + λ / 4, the phase difference is 180 ° (opposite phase), and the return light and the laser light in the resonator 101 are the weakest. Therefore, when the oscillation wavelength of the semiconductor laser is changed, a place where the laser output becomes strong and a place where the laser output becomes weak appear alternately, and the laser output at this time is detected by the photodiode 103 provided in the resonator 101. As shown in FIG. 26, a stepped waveform having a constant period is obtained. Such a waveform is generally called an interference fringe.

この階段状の波形、すなわち干渉縞の1つ1つをモードポップパルス(以下、MHP)と呼ぶ。MHPはモードホッピング現象とは異なる現象である。例えば、測定対象104までの距離がL1のとき、MHPの数が10個であったとすれば、半分の距離L2では、MHPの数は5個になる。すなわち、ある一定時間において半導体レーザの発振波長を変化させた場合、測定距離に比例してMHPの数は変わる。したがって、MHPをフォトダイオード103で検出し、MHPの周波数を測定すれば、容易に距離計測が可能となる。   Each stepped waveform, that is, each interference fringe is called a mode pop pulse (hereinafter referred to as MHP). MHP is a phenomenon different from the mode hopping phenomenon. For example, if the number of MHPs is 10 when the distance to the measurement object 104 is L1, the number of MHPs is 5 at half the distance L2. That is, when the oscillation wavelength of the semiconductor laser is changed for a certain time, the number of MHPs changes in proportion to the measurement distance. Therefore, if the MHP is detected by the photodiode 103 and the frequency of the MHP is measured, the distance can be easily measured.

上田正,山田諄,紫藤進,「半導体レーザの自己結合効果を利用した距離計」,1994年度電気関係学会東海支部連合大会講演論文集,1994年Tadashi Ueda, Satoshi Yamada, Susumu Shito, “Distance Meter Using Self-Coupling Effect of Semiconductor Laser”, Proceedings of the 1994 Tokai Branch Joint Conference of Electrical Engineering Society, 1994 山田諄,紫藤進,津田紀生,上田正,「半導体レーザの自己結合効果を利用した小型距離計に関する研究」,愛知工業大学研究報告,第31号B,p.35−42,1996年Satoshi Yamada, Susumu Shito, Norio Tsuda, Tadashi Ueda, “Study on a small rangefinder using the self-coupling effect of a semiconductor laser”, Aichi Institute of Technology research report, No. 31 B, p. 35-42, 1996 Guido Giuliani,Michele Norgia,Silvano Donati and Thierry Bosch,「Laser diode self-mixing technique for sensing applications」,JOURNAL OF OPTICS A:PURE AND APPLIED OPTICS,p.283−294,2002年Guido Giuliani, Michele Norgia, Silvano Donati and Thierry Bosch, “Laser diode self-mixing technique for sensing applications”, JOURNAL OF OPTICS A: PURE AND APPLIED OPTICS, p. 283-294, 2002

自己結合型を含む従来の干渉型の距離計では、計数装置を用いてMHPの数を測定するか、あるいはFFT(Fast Fourier Transform)を用いてMHPの周波数を測定することにより、測定対象との距離を求めるようにしている。
しかし、FFTを用いる距離計では、レーザの発振波長の変化が線形でない場合、FFTで算出されるピーク周波数と本来求められるべきMHPの平均周波数とに差が生じ、測定した距離に誤差が生じるという問題点があった。
In a conventional interference type distance meter including a self-coupling type, the number of MHPs is measured using a counting device, or the frequency of the MHP is measured using FFT (Fast Fourier Transform). I try to find the distance.
However, in a distance meter using FFT, if the change in the oscillation wavelength of the laser is not linear, there is a difference between the peak frequency calculated by FFT and the average frequency of MHP that should be originally obtained, resulting in an error in the measured distance. There was a problem.

また、計数装置を用いる距離計では、例えば外乱光などのノイズをMHPとして数えたり、信号の歯抜けのために数えられないMHPがあったりして、計数装置で数えるMHPの数に誤差が生じ、測定した距離に誤差が生じるという問題点があった。
なお、このような計数誤差は距離計に限らず、他の計数装置でも同様に発生することがある。
In a distance meter using a counting device, for example, noise such as ambient light is counted as MHP, or there are MHPs that cannot be counted due to signal tooth loss, resulting in an error in the number of MHPs counted by the counting device. There is a problem that an error occurs in the measured distance.
Such a counting error is not limited to the distance meter and may occur in other counting devices as well.

本発明は、上記課題を解決するためになされたもので、計数誤差を補正することができる計数装置および計数方法、MHPの計数誤差を補正して距離の測定精度を向上させることができる距離計および距離計測方法を提供することを目的とする。   The present invention has been made in order to solve the above-described problems. A counting device and a counting method capable of correcting a counting error, and a distance meter capable of improving a distance measurement accuracy by correcting a counting error of an MHP. It is another object of the present invention to provide a distance measuring method.

本発明は、特定の物理量と線形の関係を有し前記物理量が一定の場合は略単一周波数となる信号の数を数える計数装置において、一定の計数期間における入力信号の数を数える計数手段と、前記計数期間中の前記入力信号の周期を信号が入力される度に測定する周期測定手段と、この周期測定手段の測定結果から前記計数期間中の信号周期の度数分布を作成する度数分布作成手段と、前記度数分布をFFTにより解析して、前記度数分布の極大値の間隔に相当する周波数に極大値が出現する周波数スペクトルを求める周期度数分布周波数解析手段と、この周期度数分布周波数解析手段が求めた周波数スペクトルに出現した極大値から前記入力信号の周期の分布の代表値を求める代表値導出手段と、前記度数分布から、前記代表値の第1の所定数倍以下である階級の度数の総和Nsと、前記代表値の第2の所定数倍以上である階級の度数の総和Nwとを求め、これらの度数NsとNwに基づいて前記計数手段の計数結果を補正する補正値算出手段とを有するものである。
また、本発明の計数装置の1構成例において、前記代表値導出手段は、前記周期度数分布周波数解析手段が求めた周波数スペクトルにおいて強度が最大の極大値を前記入力信号の周期の分布の代表値を示すものと見なし、前記代表値を求めるものである。
また、本発明の計数装置の1構成例において、前記代表値導出手段は、前記周期度数分布周波数解析手段が求めた周波数スペクトルにおいて周波数が最も高い極大値を前記入力信号の周期の分布の代表値を示すものと見なし、前記代表値を求めるものである。
また、本発明の計数装置の1構成例において、前記補正値算出手段は、前記計数手段の計数結果をNとしたとき、補正後の計数結果N’を、N’=N+Nw−Nsにより求めるものである。
The present invention relates to a counting device for counting the number of signals having a linear relationship with a specific physical quantity and having a substantially single frequency when the physical quantity is constant, and counting means for counting the number of input signals in a constant counting period. A period measuring unit that measures the period of the input signal during the counting period every time a signal is input, and a frequency distribution generator that generates a frequency distribution of the signal period during the counting period from the measurement result of the period measuring unit means and analyzes the FFT the frequency distribution, a frequency analysis means of the periodic frequency distribution to determine the frequency spectrum maximum value appears in the frequency corresponding to the spacing of the maxima of the frequency distribution, the frequency of the periodic frequency distribution from the maximum value analyzing means appeared to the frequency spectrum obtained representative value deriving means for obtaining a representative value of a distribution of periods of the input signal, from the frequency distribution, a first of said representative value The sum Ns of the frequencies of the class that is less than or equal to a constant multiple and the sum Nw of the frequencies of the class that are greater than or equal to the second predetermined number times the representative value are obtained, and the counting means counts based on these frequencies Ns and Nw. Correction value calculating means for correcting the result.
Further, in one configuration example of the counting device of the present invention, the representative value deriving unit displays a local maximum value having the maximum intensity in the frequency spectrum obtained by the frequency analysis unit of the periodic frequency distribution as a representative of the period distribution of the input signal. It is assumed that the value is shown, and the representative value is obtained.
Further, in one configuration example of the counting device according to the present invention, the representative value deriving means represents a maximum value having the highest frequency in the frequency spectrum obtained by the frequency analysis means of the periodic frequency distribution as a representative of the period distribution of the input signal. It is assumed that the value is shown, and the representative value is obtained.
Further, in one configuration example of the counting device of the present invention, the correction value calculating means obtains a corrected counting result N ′ by N ′ = N + Nw−Ns, where N is the counting result of the counting means. It is.

また、本発明の距離計は、測定対象にレーザ光を放射する半導体レーザと、発振波長が連続的に単調増加する期間を少なくとも含む第1の発振期間と発振波長が連続的に単調減少する期間を少なくとも含む第2の発振期間とが交互に存在するように前記半導体レーザを動作させるレーザドライバと、前記半導体レーザから放射されたレーザ光と前記測定対象からの戻り光との干渉光を電気信号に変換する受光器と、この受光器の出力信号に含まれる、前記半導体レーザから放射されたレーザ光と前記測定対象からの戻り光とによって生じる干渉波形の数を数える計数手段と、前記干渉波形の数を数える計数期間中の前記干渉波形の周期を干渉波形が入力される度に測定する周期測定手段と、この周期測定手段の測定結果から前記計数期間中の干渉波形の周期の度数分布を作成する度数分布作成手段と、前記度数分布をFFTにより解析して、前記度数分布の極大値の間隔に相当する周波数に極大値が出現する周波数スペクトルを求める周期度数分布周波数解析手段と、この周期度数分布周波数解析手段が求めた周波数スペクトルに出現した極大値から前記干渉波形の周期の分布の代表値を求める代表値導出手段と、前記度数分布から、前記代表値の第1の所定数倍以下である階級の度数の総和Nsと、前記代表値の第2の所定数倍以上である階級の度数の総和Nwとを求め、これらの度数NsとNwに基づいて前記計数手段の計数結果を補正する補正値算出手段と、この補正値算出手段で補正された計数結果から前記測定対象との距離を求める演算手段とを有するものである。 The distance meter of the present invention includes a semiconductor laser that emits laser light to a measurement target, a first oscillation period that includes at least a period in which the oscillation wavelength continuously increases monotonously, and a period in which the oscillation wavelength continuously decreases monotonously. A laser driver that operates the semiconductor laser so that second oscillation periods including at least the second oscillation period alternately exist, and an interference signal between the laser light emitted from the semiconductor laser and the return light from the measurement object A light receiving device that converts the interference waveform generated by the laser light emitted from the semiconductor laser and the return light from the measurement object, and the interference waveform included in the output signal of the light receiving device. A period measuring means for measuring the period of the interference waveform during the counting period for counting the number of times each time the interference waveform is input, and a measurement result of the period measuring means from the measurement result of the period measuring means. A frequency distribution generating means for generating a frequency distribution of the periods of Wataru waveform, the frequency distribution by analyzing the FFT and cycle frequency for determining the frequency spectrum maxima in frequency corresponding to the spacing of the maxima of the frequency distribution appears A frequency analysis means for distribution, a representative value deriving means for obtaining a representative value of the period distribution of the interference waveform from a maximum value that appears in the frequency spectrum obtained by the frequency analysis means for the periodic frequency distribution , and the frequency distribution, A sum Ns of the frequencies of the class that is less than or equal to the first predetermined number times the representative value and a sum Nw of the frequencies of the class that are the second predetermined number of times or more of the representative value are obtained, and these frequencies Ns and Nw are obtained. A correction value calculating means for correcting the counting result of the counting means based on the calculation result, and a calculating means for calculating a distance from the measurement object from the counting result corrected by the correction value calculating means.

また、本発明の距離計は、測定対象にレーザ光を放射する半導体レーザと、発振波長が連続的に単調増加する期間を少なくとも含む第1の発振期間と発振波長が連続的に単調減少する期間を少なくとも含む第2の発振期間とが交互に存在するように前記半導体レーザを動作させるレーザドライバと、前記半導体レーザの光出力を電気信号に変換する受光器と、この受光器の出力信号に含まれる、前記半導体レーザから放射されたレーザ光と前記測定対象からの戻り光との自己結合効果によって生じる干渉波形の数を数える計数手段と、前記干渉波形の数を数える計数期間中の前記干渉波形の周期を干渉波形が入力される度に測定する周期測定手段と、この周期測定手段の測定結果から前記計数期間中の干渉波形の周期の度数分布を作成する度数分布作成手段と、前記度数分布をFFTにより解析して、前記度数分布の極大値の間隔に相当する周波数に極大値が出現する周波数スペクトルを求める周期度数分布周波数解析手段と、この周期度数分布周波数解析手段が求めた周波数スペクトルに出現した極大値から前記干渉波形の周期の分布の代表値を求める代表値導出手段と、前記度数分布から、前記代表値の第1の所定数倍以下である階級の度数の総和Nsと、前記代表値の第2の所定数倍以上である階級の度数の総和Nwとを求め、これらの度数NsとNwに基づいて前記計数手段の計数結果を補正する補正値算出手段と、この補正値算出手段で補正された計数結果から前記測定対象との距離を求める演算手段とを有するものである。
また、本発明の距離計の1構成例において、前記代表値導出手段は、前記周期度数分布周波数解析手段が求めた周波数スペクトルにおいて強度が最大の極大値を前記干渉波形の周期の分布の代表値を示すものと見なし、前記代表値を求めるものである。
また、本発明の距離計の1構成例において、前記代表値導出手段は、前記周期度数分布周波数解析手段が求めた周波数スペクトルにおいて周波数が最も高い極大値を前記干渉波形の周期の分布の代表値を示すものと見なし、前記代表値を求めるものである。
また、本発明の距離計の1構成例において、前記補正値算出手段は、前記計数手段の計数結果をNとしたとき、補正後の計数結果N’を、N’=N+Nw−Nsにより求めるものである。
The distance meter of the present invention includes a semiconductor laser that emits laser light to a measurement target, a first oscillation period that includes at least a period in which the oscillation wavelength continuously increases monotonously, and a period in which the oscillation wavelength continuously decreases monotonously. Included in the output signal of the laser receiver, the laser driver that operates the semiconductor laser so that the second oscillation period including at least the second oscillation period alternately exists, the optical receiver that converts the optical output of the semiconductor laser into an electrical signal, and Counting means for counting the number of interference waveforms generated by the self-coupling effect between the laser light emitted from the semiconductor laser and the return light from the measurement object, and the interference waveform during the counting period for counting the number of the interference waveforms A period measuring means for measuring the period of the interference waveform, and a frequency distribution of the period of the interference waveform during the counting period from the measurement result of the period measuring means. A distribution creation means, and analyzed by FFT the frequency distribution, a frequency analysis means of the periodic frequency distribution to determine the frequency spectrum maximum value appears in the frequency corresponding to the spacing of the maxima of the frequency distribution, the cycle frequency distribution Representative value deriving means for obtaining a representative value of the distribution of the period of the interference waveform from the maximum value appearing in the frequency spectrum obtained by the frequency analysis means, and from the frequency distribution at a first predetermined number times or less of the representative value. A total sum Ns of frequencies of a certain class and a total frequency Nw of classes that are equal to or more than the second predetermined number times the representative value are obtained, and the counting result of the counting means is corrected based on these frequencies Ns and Nw. It has a correction value calculation means and a calculation means for obtaining a distance from the measurement object from the count result corrected by the correction value calculation means.
Further, in one configuration example of the distance meter of the present invention, the representative value deriving unit displays a maximum value having the maximum intensity in the frequency spectrum obtained by the frequency analysis unit of the periodic frequency distribution as a representative of the period distribution of the interference waveform. It is assumed that the value is shown, and the representative value is obtained.
Further, in one configuration example of the distance meter of the present invention, the representative value deriving means obtains a maximum value having the highest frequency in the frequency spectrum obtained by the frequency analysis means of the periodic frequency distribution as a representative of the period distribution of the interference waveform. It is assumed that the value is shown, and the representative value is obtained.
Further, in one configuration example of the distance meter of the present invention, the correction value calculation means obtains the corrected count result N ′ by N ′ = N + Nw−Ns, where N is the counting result of the counting means. It is.

また、本発明の計数方法は、一定の計数期間における入力信号の数を数える計数手順と、前記計数期間中の前記入力信号の周期を信号が入力される度に測定する周期測定手順と、この周期測定手順の測定結果から前記計数期間中の信号周期の度数分布を作成する度数分布作成手順と、前記度数分布をFFTにより解析して、前記度数分布の極大値の間隔に相当する周波数に極大値が出現する周波数スペクトルを求める周期度数分布周波数解析手順と、この周期度数分布周波数解析手順で求めた周波数スペクトルに出現した極大値から前記入力信号の周期の分布の代表値を求める代表値導出手順と、前記度数分布から、前記代表値の第1の所定数倍以下である階級の度数の総和Nsと、前記代表値の第2の所定数倍以上である階級の度数の総和Nwとを求め、これらの度数NsとNwに基づいて前記計数手順の計数結果を補正する補正値算出手順とを備えるものである。 The counting method of the present invention includes a counting procedure for counting the number of input signals in a certain counting period, a period measuring procedure for measuring the period of the input signal in the counting period every time a signal is input, A frequency distribution creating procedure for creating a frequency distribution of the signal period during the counting period from the measurement result of the period measuring procedure, and the frequency distribution is analyzed by FFT to maximize the frequency corresponding to the interval between the maximum values of the frequency distribution. a frequency analyzing procedure of the periodic frequency distribution to determine the frequency spectrum value appears, the representative value determining a representative value of a distribution of periods of the input signal from the maximum value that appeared frequency spectrum obtained by the frequency analysis procedure of this cycle frequency distribution From the derivation procedure and the frequency distribution, the sum Ns of the frequencies of the class that is not more than a first predetermined number of times of the representative value and the frequencies of the class that are not less than the second predetermined number of times of the representative value It obtains a sum Nw, in which and a correction value calculation procedure for correcting the counting result of the counting procedure on the basis of the frequencies Ns and Nw.

また、本発明の距離計測方法は、発振波長が連続的に単調増加する期間を少なくとも含む第1の発振期間と発振波長が連続的に単調減少する期間を少なくとも含む第2の発振期間とが交互に存在するように前記半導体レーザを動作させる発振手順と、前記半導体レーザから放射されたレーザ光と前記測定対象からの戻り光との干渉光を電気信号に変換する受光器の出力信号に含まれる、前記半導体レーザから放射されたレーザ光と前記測定対象からの戻り光とによって生じる干渉波形の数を数える計数手順と、前記干渉波形の数を数える計数期間中の前記干渉波形の周期を干渉波形が入力される度に測定する周期測定手順と、この周期測定手順の測定結果から前記計数期間中の干渉波形の周期の度数分布を作成する度数分布作成手順と、前記度数分布をFFTにより解析して、前記度数分布の極大値の間隔に相当する周波数に極大値が出現する周波数スペクトルを求める周期度数分布周波数解析手順と、この周期度数分布周波数解析手順で求めた周波数スペクトルに出現した極大値から前記干渉波形の周期の分布の代表値を求める代表値導出手順と、前記度数分布から、前記代表値の第1の所定数倍以下である階級の度数の総和Nsと、前記代表値の第2の所定数倍以上である階級の度数の総和Nwとを求め、これらの度数NsとNwに基づいて前記計数手順の計数結果を補正する補正値算出手順と、この補正値算出手順で補正された計数結果から前記測定対象との距離を求める演算手順とを備えるものである。 In the distance measuring method of the present invention, the first oscillation period including at least a period in which the oscillation wavelength continuously increases monotonically and the second oscillation period including at least a period in which the oscillation wavelength continuously decreases monotonically are alternated. Included in the oscillation procedure for operating the semiconductor laser to exist in the light source and the output signal of the light receiver that converts the interference light between the laser light emitted from the semiconductor laser and the return light from the measurement object into an electrical signal A counting procedure for counting the number of interference waveforms caused by the laser light emitted from the semiconductor laser and the return light from the measurement object, and the period of the interference waveform during the counting period for counting the number of interference waveforms as an interference waveform A period measurement procedure for measuring the frequency of each input, a frequency distribution creation procedure for creating a frequency distribution of the period of the interference waveform during the counting period from the measurement result of the period measurement procedure, And analyzed by FFT number distribution, and the frequency analyzing procedure of the periodic frequency distribution to determine the frequency spectrum maximum at a frequency appears corresponding to the spacing of the maxima of the frequency distribution, determined by the frequency analysis procedure of this cycle frequency distribution A representative value derivation procedure for obtaining a representative value of the distribution of the period of the interference waveform from the maximum value appearing in the frequency spectrum, and a sum of the frequencies of the class that is not more than a first predetermined number times the representative value from the frequency distribution. A correction value calculation procedure for determining Ns and a total sum Nw of the frequencies of the class that is a second predetermined number times or more of the representative value, and correcting the counting result of the counting procedure based on these frequencies Ns and Nw; And a calculation procedure for obtaining a distance from the measurement target from the count result corrected by the correction value calculation procedure.

また、本発明の距離計測方法は、発振波長が連続的に単調増加する期間を少なくとも含む第1の発振期間と発振波長が連続的に単調減少する期間を少なくとも含む第2の発振期間とが交互に存在するように前記半導体レーザを動作させる発振手順と、前記半導体レーザの光出力を電気信号に変換する受光器の出力信号に含まれる、前記半導体レーザから放射されたレーザ光と前記測定対象からの戻り光との自己結合効果によって生じる干渉波形の数を数える計数手順と、前記干渉波形の数を数える計数期間中の前記干渉波形の周期を干渉波形が入力される度に測定する周期測定手順と、この周期測定手順の測定結果から前記計数期間中の干渉波形の周期の度数分布を作成する度数分布作成手順と、前記度数分布をFFTにより解析して、前記度数分布の極大値の間隔に相当する周波数に極大値が出現する周波数スペクトルを求める周期度数分布周波数解析手順と、この周期度数分布周波数解析手順で求めた周波数スペクトルに出現した極大値から前記干渉波形の周期の分布の代表値を求める代表値導出手順と、前記度数分布から、前記代表値の第1の所定数倍以下である階級の度数の総和Nsと、前記代表値の第2の所定数倍以上である階級の度数の総和Nwとを求め、これらの度数NsとNwに基づいて前記計数手順の計数結果を補正する補正値算出手順と、この補正値算出手順で補正された計数結果から前記測定対象との距離を求める演算手順とを備えるものである。 In the distance measuring method of the present invention, the first oscillation period including at least a period in which the oscillation wavelength continuously increases monotonically and the second oscillation period including at least a period in which the oscillation wavelength continuously decreases monotonically are alternated. An oscillation procedure for operating the semiconductor laser to be present in the laser beam, and a laser beam emitted from the semiconductor laser and an object to be measured included in an output signal of a light receiver that converts the optical output of the semiconductor laser into an electrical signal. Counting procedure for counting the number of interference waveforms caused by the self-coupling effect with the return light, and a period measuring procedure for measuring the period of the interference waveform during the counting period for counting the number of interference waveforms each time the interference waveform is input When the frequency distribution generating procedure for creating a frequency distribution of the periods of the interference waveforms during the counting interval from the measurement result of the cycle measuring procedure, and analyzed by FFT the frequency distribution, the The number distribution and frequency analysis procedure of the periodic frequency distribution to determine the frequency spectrum maximum value appears in the frequency corresponding to the spacing of the maxima of the local maximum value appearing in the frequency spectrum obtained by the frequency analysis procedure of this cycle frequency distribution A representative value derivation procedure for obtaining a representative value of the distribution of the period of the interference waveform, a total sum Ns of the frequencies of the class that is not more than a first predetermined number times the representative value, and a second value of the representative value from the frequency distribution. A correction value calculation procedure for obtaining a total sum Nw of class frequencies that is a predetermined number of times or more and correcting the counting result of the counting procedure based on these frequencies Ns and Nw, and a count corrected by this correction value calculation procedure And a calculation procedure for obtaining a distance from the measurement object from the result.

本発明によれば、計数期間中の入力信号の周期を測定し、この測定結果から計数期間中の信号周期の度数分布を作成し、度数分布の周波数成分の解析結果から入力信号の周期の代表値と見なすことができる値を求めるので、例えば入力信号のS/N比が低下している場合であっても、入力信号の真の周期を求めることができる。そして、本発明では、度数分布から、代表値の第1の所定数倍以下である階級の度数の総和Nsと、代表値の第2の所定数倍以上である階級の度数の総和Nwとを求め、これらの度数NsとNwに基づいて計数手段の計数結果を補正することにより、計数時の欠落や過剰な計数の影響を除去して、計数装置の計数誤差を補正することができる。   According to the present invention, the period of the input signal during the counting period is measured, the frequency distribution of the signal period during the counting period is created from the measurement result, and the input signal period is representative from the analysis result of the frequency component of the frequency distribution. Since a value that can be regarded as a value is obtained, the true period of the input signal can be obtained even when, for example, the S / N ratio of the input signal is reduced. In the present invention, from the frequency distribution, the sum Ns of the frequencies of the class that is less than or equal to the first predetermined number of times of the representative value and the sum of the frequencies of the class Nw that is greater than or equal to the second predetermined number of times of the representative value are obtained. By obtaining and correcting the counting result of the counting means based on these frequencies Ns and Nw, the counting error of the counting device can be corrected by removing the influence of missing or excessive counting at the time of counting.

また、本発明では、計数期間中の干渉波形の周期を測定し、この測定結果から計数期間中の干渉波形の周期の度数分布を作成し、度数分布の周波数成分の解析結果から干渉波形の周期の代表値と見なすことができる値を求めるので、例えば遠距離測定時に信号のS/N比が低下している場合であっても、干渉波形の真の周期を求めることができる。そして、本発明では、度数分布から、代表値の第1の所定数倍以下である階級の度数の総和Nsと、代表値の第2の所定数倍以上である階級の度数の総和Nwとを求め、これらの度数NsとNwに基づいて計数手段の計数結果を補正することにより、計数時の欠落や過剰な計数の影響を除去して、干渉波形の計数誤差を補正することができるので、計数手段を用いて干渉波形の数を測定し測定対象との距離を求める距離計において、距離の測定精度を向上させることができる。   In the present invention, the period of the interference waveform during the counting period is measured, a frequency distribution of the period of the interference waveform during the counting period is created from the measurement result, and the period of the interference waveform is calculated from the analysis result of the frequency component of the frequency distribution. Since a value that can be regarded as a representative value of the signal is obtained, the true period of the interference waveform can be obtained even when, for example, the signal-to-noise ratio of the signal is reduced during long-distance measurement. In the present invention, from the frequency distribution, the sum Ns of the frequencies of the class that is less than or equal to the first predetermined number of times of the representative value and the sum of the frequencies of the class Nw that is greater than or equal to the second predetermined number of times of the representative value are obtained. By finding and correcting the counting result of the counting means based on these frequencies Ns and Nw, it is possible to eliminate the influence of missing or excessive counting at the time of counting, and to correct the counting error of the interference waveform. In a distance meter that measures the number of interference waveforms using a counting means and obtains the distance to the measurement object, the distance measurement accuracy can be improved.

[第1の実施の形態]
以下、本発明の第1の実施の形態について図面を参照して説明する。図1は本発明の第1の実施の形態となる距離計の構成を示すブロック図である。図1の距離計は、測定対象10にレーザ光を放射する半導体レーザ1と、半導体レーザ1の光出力を電気信号に変換する受光器であるフォトダイオード2と、半導体レーザ1からの光を集光して測定対象10に照射すると共に、測定対象10からの戻り光を集光して半導体レーザ1に入射させるレンズ3と、半導体レーザ1に発振波長が連続的に増加する第1の発振期間と発振波長が連続的に減少する第2の発振期間とを交互に繰り返させるレーザドライバ4と、フォトダイオード2の出力電流を電圧に変換して増幅する電流−電圧変換増幅器5と、電流−電圧変換増幅器5の出力電圧から搬送波を除去するフィルタ回路6と、フィルタ回路6の出力電圧に含まれるMHPの数を数える計数装置7と、MHPの数から測定対象10との距離を算出する演算装置8と、演算装置8の算出結果を表示する表示装置9とを有する。
[First Embodiment]
Hereinafter, a first embodiment of the present invention will be described with reference to the drawings. FIG. 1 is a block diagram showing the configuration of a distance meter according to the first embodiment of the present invention. The distance meter of FIG. 1 collects light from a semiconductor laser 1 that emits laser light to a measurement object 10, a photodiode 2 that is a light receiver that converts the optical output of the semiconductor laser 1 into an electrical signal, and light from the semiconductor laser 1. A light 3 that irradiates the measurement target 10 with light, collects return light from the measurement target 10 and makes it incident on the semiconductor laser 1, and a first oscillation period in which the oscillation wavelength continuously increases in the semiconductor laser 1 And a laser driver 4 that alternately repeats the second oscillation period in which the oscillation wavelength continuously decreases, a current-voltage conversion amplifier 5 that converts and amplifies the output current of the photodiode 2 into a voltage, and a current-voltage The filter circuit 6 that removes the carrier wave from the output voltage of the conversion amplifier 5, the counter 7 that counts the number of MHPs included in the output voltage of the filter circuit 6, and the distance to the measurement object 10 from the number of MHPs Having an arithmetic unit 8 for output, and a display device 9 for displaying the calculation result of the arithmetic unit 8.

以下、説明容易にするために、半導体レーザ1には、モードホッピング現象を持たない型(VCSEL型、DFBレーザ型)のものが用いられているものと想定する。   Hereinafter, for ease of explanation, it is assumed that a semiconductor laser 1 of a type that does not have a mode hopping phenomenon (VCSEL type, DFB laser type) is used.

例えば、レーザドライバ4は、時間に関して一定の変化率で増減を繰り返す三角波駆動電流を注入電流として半導体レーザ1に供給する。これにより、半導体レーザ1は、注入電流の大きさに比例して発振波長が一定の変化率で連続的に増加する第1の発振期間と発振波長が一定の変化率で連続的に減少する第2の発振期間とを交互に繰り返すように駆動される。   For example, the laser driver 4 supplies a triangular wave drive current that repeatedly increases and decreases at a constant change rate with respect to time to the semiconductor laser 1 as an injection current. As a result, the semiconductor laser 1 has a first oscillation period in which the oscillation wavelength continuously increases at a constant change rate in proportion to the magnitude of the injection current, and a first oscillation period in which the oscillation wavelength continuously decreases at a constant change rate. It is driven to alternately repeat the two oscillation periods.

図2は、半導体レーザ1の発振波長の時間変化を示す図である。図2において、t−1はt−1番目の発振期間、tはt番目の発振期間、t+1はt+1番目の発振期間、t+2はt+2番目の発振期間、t+3はt+3番目の発振期間、t+4はt+4番目の発振期間、λaは各期間における発振波長の最小値、λbは各期間における発振波長の最大値、Tは三角波の周期である。本実施の形態では、発振波長の最大値λb及び発振波長の最小値λaはそれぞれ常に一定になされており、それらの差λb−λaも常に一定になされている。   FIG. 2 is a diagram showing the change over time of the oscillation wavelength of the semiconductor laser 1. In FIG. 2, t-1 is the t-1th oscillation period, t is the tth oscillation period, t + 1 is the t + 1th oscillation period, t + 2 is the t + 2nd oscillation period, t + 3 is the t + 3rd oscillation period, and t + 4 is The t + 4th oscillation period, λa is the minimum value of the oscillation wavelength in each period, λb is the maximum value of the oscillation wavelength in each period, and T is the period of the triangular wave. In the present embodiment, the maximum value λb of the oscillation wavelength and the minimum value λa of the oscillation wavelength are always constant, and the difference λb−λa is also always constant.

半導体レーザ1から出射したレーザ光は、レンズ3によって集光され、測定対象10に入射する。測定対象10で反射された光は、レンズ3によって集光され、半導体レーザ1に入射する。ただし、レンズ3による集光は必須ではない。フォトダイオード2は、半導体レーザ1の内部又はその近傍に配置され、半導体レーザ1の光出力を電流に変換する。電流−電圧変換増幅器5は、フォトダイオード2の出力電流を電圧に変換して増幅する。   Laser light emitted from the semiconductor laser 1 is collected by the lens 3 and enters the measurement object 10. The light reflected by the measurement object 10 is collected by the lens 3 and enters the semiconductor laser 1. However, condensing by the lens 3 is not essential. The photodiode 2 is disposed in the semiconductor laser 1 or in the vicinity thereof, and converts the optical output of the semiconductor laser 1 into a current. The current-voltage conversion amplifier 5 converts the output current of the photodiode 2 into a voltage and amplifies it.

フィルタ回路6は、変調波から重畳信号を抽出する機能を有するものである。図3(A)は電流−電圧変換増幅器5の出力電圧波形を模式的に示す図、図3(B)はフィルタ回路6の出力電圧波形を模式的に示す図である。これらは、フォトダイオード2の出力である図3(A)の波形(変調波)から、図2の半導体レーザ1の発振波形(搬送波)を除去して、図3(B)のMHP波形(重畳波)を抽出する過程を表している。   The filter circuit 6 has a function of extracting a superimposed signal from the modulated wave. 3A is a diagram schematically showing an output voltage waveform of the current-voltage conversion amplifier 5, and FIG. 3B is a diagram schematically showing an output voltage waveform of the filter circuit 6. As shown in FIG. These are obtained by removing the oscillation waveform (carrier wave) of the semiconductor laser 1 of FIG. 2 from the waveform (modulated wave) of FIG. 3A, which is the output of the photodiode 2, and the MHP waveform (superposition) of FIG. (Wave) is extracted.

計数装置7は、フィルタ回路6の出力電圧に含まれるMHPの数を第1の発振期間t−1,t+1,t+3と第2の発振期間t,t+2,t+4の各々について数える。図4は計数装置7の構成の1例を示すブロック図である。計数装置7は、判定部71と、論理積演算部(AND)72と、カウンタ73と、計数結果補正部74と、記憶部75とから構成される。電流−電圧変換増幅器5とフィルタ回路6と計数装置7の判定部71とAND72とカウンタ73とは、計数手段を構成している。   The counting device 7 counts the number of MHPs included in the output voltage of the filter circuit 6 for each of the first oscillation periods t-1, t + 1, t + 3 and the second oscillation periods t, t + 2, t + 4. FIG. 4 is a block diagram showing an example of the configuration of the counting device 7. The counting device 7 includes a determination unit 71, an AND operation unit (AND) 72, a counter 73, a counting result correction unit 74, and a storage unit 75. The current-voltage conversion amplifier 5, the filter circuit 6, the determination unit 71 of the counting device 7, the AND 72, and the counter 73 constitute counting means.

図5は計数結果補正部74の構成の1例を示すブロック図である。計数結果補正部74は、周期測定部740と、度数分布作成部741と、度数分布周波数解析部742と、代表値導出部743と、補正値算出部744とから構成される。   FIG. 5 is a block diagram showing an example of the configuration of the counting result correction unit 74. The counting result correction unit 74 includes a period measurement unit 740, a frequency distribution creation unit 741, a frequency distribution frequency analysis unit 742, a representative value derivation unit 743, and a correction value calculation unit 744.

図6(A)〜図6(F)は計数装置7の動作を説明するための図であり、図6(A)はフィルタ回路6の出力電圧の波形、すなわちMHPの波形を模式的に示す図、図6(B)は図6(A)に対応する判定部71の出力を示す図、図6(C)は計数装置7に入力されるゲート信号GSを示す図、図6(D)は図6(B)に対応するカウンタ73の計数結果を示す図、図6(E)は計数装置7に入力されるクロック信号CLKを示す図、図6(F)は図6(B)に対応する周期測定部740の測定結果を示す図である。   6A to 6F are diagrams for explaining the operation of the counting device 7. FIG. 6A schematically shows the waveform of the output voltage of the filter circuit 6, that is, the waveform of MHP. FIG. 6B is a diagram showing the output of the determination unit 71 corresponding to FIG. 6A, FIG. 6C is a diagram showing the gate signal GS input to the counting device 7, and FIG. Is a diagram showing the counting result of the counter 73 corresponding to FIG. 6 (B), FIG. 6 (E) is a diagram showing the clock signal CLK input to the counting device 7, and FIG. 6 (F) is FIG. 6 (B). It is a figure which shows the measurement result of the corresponding period measurement part 740. FIG.

まず、計数装置7の判定部71は、図6(A)に示すフィルタ回路6の出力電圧がハイレベル(H)かローレベル(L)かを判定して、図6(B)のような判定結果を出力する。このとき、判定部71は、フィルタ回路6の出力電圧が上昇してしきい値TH1以上になったときにハイレベルと判定し、フィルタ回路6の出力電圧が下降してしきい値TH2(TH2<TH1)以下になったときにローレベルと判定することにより、フィルタ回路6の出力を2値化する。   First, the determination unit 71 of the counting device 7 determines whether the output voltage of the filter circuit 6 shown in FIG. 6A is high level (H) or low level (L), as shown in FIG. Output the judgment result. At this time, the determination unit 71 determines that the output voltage of the filter circuit 6 is high level when the output voltage becomes equal to or higher than the threshold value TH1, and the output voltage of the filter circuit 6 decreases and the threshold value TH2 (TH2 <TH1) The output of the filter circuit 6 is binarized by determining the low level when it becomes below.

AND72は、判定部71の出力と図6(C)のようなゲート信号GSとの論理積演算の結果を出力し、カウンタ73は、AND72の出力の立ち上がりをカウントする(図6(D))。ここで、ゲート信号GSは、計数期間(本実施の形態では第1の発振期間と第2の発振期間の各々)の先頭で立ち上がり、計数期間の終わりで立ち下がる信号である。したがって、カウンタ73は、計数期間中のAND72の出力の立ち上がりエッジの数(すなわち、MHPの立ち上がりエッジの数)を数えることになる。   The AND 72 outputs the logical product operation result of the output of the determination unit 71 and the gate signal GS as shown in FIG. 6C, and the counter 73 counts the rise of the output of the AND 72 (FIG. 6D). . Here, the gate signal GS is a signal that rises at the beginning of the counting period (in this embodiment, each of the first oscillation period and the second oscillation period) and falls at the end of the counting period. Therefore, the counter 73 counts the number of rising edges of the output of the AND 72 during the counting period (that is, the number of rising edges of MHP).

一方、計数結果補正部74の周期測定部740は、計数期間中のAND72の出力の立ち上がりエッジの周期(すなわち、MHPの周期)を立ち上がりエッジが発生する度に測定する。このとき、周期測定部740は、図6(E)に示すクロック信号CLKの周期を1単位としてMHPの周期を測定する。図6(F)の例では、周期測定部740は、MHPの周期としてTα,Tβ,Tγを順次測定している。図6(E)、図6(F)から明らかなように、周期Tα,Tβ,Tγの大きさは、それぞれ5クロック、4クロック、2クロックである。クロック信号CLKの周波数は、MHPの取り得る最高周波数に対して十分に高いものとする。
記憶部75は、カウンタ73の計数結果と周期測定部740の測定結果を記憶する。
On the other hand, the period measuring unit 740 of the counting result correcting unit 74 measures the period of the rising edge of the output of the AND 72 during the counting period (that is, the MHP period) every time a rising edge occurs. At this time, the period measuring unit 740 measures the MHP period with the period of the clock signal CLK shown in FIG. In the example of FIG. 6F, the period measurement unit 740 sequentially measures Tα, Tβ, and Tγ as the MHP period. As is apparent from FIGS. 6E and 6F, the periods Tα, Tβ, and Tγ are 5 clocks, 4 clocks, and 2 clocks, respectively. The frequency of the clock signal CLK is assumed to be sufficiently higher than the highest frequency that the MHP can take.
The storage unit 75 stores the count result of the counter 73 and the measurement result of the period measurement unit 740.

ゲート信号GSが立ち下がり、計数期間が終了した後、計数結果補正部74の度数分布作成部741は、記憶部75に記憶された測定結果から計数期間中のMHPの周期の度数分布を作成する。
続いて、計数結果補正部74の度数分布周波数解析部742は、度数分布作成部741が作成した度数分布の周波数成分を解析する。この周波数解析の手法としては、具体的にはFFT(Fast Fourier Transform)を用いればよい。
After the gate signal GS falls and the counting period ends, the frequency distribution creation unit 741 of the counting result correction unit 74 creates a frequency distribution of the MHP period during the counting period from the measurement result stored in the storage unit 75. .
Subsequently, the frequency distribution frequency analysis unit 742 of the counting result correction unit 74 analyzes the frequency component of the frequency distribution created by the frequency distribution creation unit 741. Specifically, FFT (Fast Fourier Transform) may be used as the frequency analysis method.

度数分布周波数解析部742と後述する代表値導出部743の目的は、MHPの真の周期と見なせる値を求めることである。
通常の測定条件の場合、真の信号周期の近傍に信号周期分布の代表値(最頻値もしくは中央値)が現れる。図7は通常の測定条件におけるMHPの周期の度数分布の1例を示している。図7の例では、最頻値MOと中央値MEは、MHPの真の周期30[samplings]とおおむね一致している。
The purpose of the frequency distribution frequency analyzing unit 742 and the representative value deriving unit 743 described later is to obtain a value that can be regarded as a true period of the MHP.
Under normal measurement conditions, a representative value (mode or median) of the signal period distribution appears in the vicinity of the true signal period. FIG. 7 shows an example of the frequency distribution of the MHP cycle under normal measurement conditions. In the example of FIG. 7, the mode value MO and the median value ME are approximately the same as the true period 30 [samplings] of the MHP.

一方、例えば遠距離測定時に信号のS/N比(Signal to Noise Ratio)が低下すると、MHPの波形に大量の欠落が生じ、MHPの周期の度数分布の代表値が真の周期からかけ離れた値となる。図8はS/N比低下時のMHPの周期の度数分布の1例を示している。図8の例では、最頻値MOと中央値MEは、MHPの真の周期9[samplings]とは一致しない。   On the other hand, for example, when the signal to noise ratio (Signal to Noise Ratio) is lowered during long-distance measurement, a large amount of loss occurs in the MHP waveform, and the representative value of the frequency distribution of the MHP period is far from the true period. It becomes. FIG. 8 shows an example of the frequency distribution of the MHP period when the S / N ratio is lowered. In the example of FIG. 8, the mode value MO and the median value ME do not coincide with the true period 9 [samplings] of the MHP.

ただし、MHPの欠落は、MHPの真の周期の整数倍近傍の周期を持つため、MHPの周期の度数分布の極大値の間隔は真の周期とおおむね等しい値になる。図8の例によると、度数分布の極大値の間隔は、MHPの真の周期9[samplings]とおおむね一致していることが分かる。   However, since the lack of MHP has a period in the vicinity of an integral multiple of the true period of MHP, the interval between the maximum values of the frequency distribution of the MHP period is approximately equal to the true period. According to the example of FIG. 8, it can be seen that the interval between the maximum values of the frequency distribution roughly matches the true period 9 [samplings] of the MHP.

そこで、本実施の形態では、MHPの周期の度数分布の周波数成分を解析することにより、MHPの真の周期と見なせる値を求める。度数分布周波数解析部742が度数分布を周波数解析すると、図9に示すように度数分布の極大値の間隔に相当する周波数に極大値が出現する。
計数結果補正部74の代表値導出部743は、度数分布周波数解析部742の周波数解析結果において強度が最大の信号をMHPの周期の代表値(最頻値もしくは中央値)を示すものと見なし、MHPの周期の代表値T0を求める。
Therefore, in this embodiment, by analyzing the frequency component of the frequency distribution of the MHP cycle, a value that can be regarded as the true cycle of the MHP is obtained. When the frequency distribution frequency analysis unit 742 performs frequency analysis on the frequency distribution, a maximum value appears at a frequency corresponding to the interval between the maximum values of the frequency distribution as shown in FIG.
The representative value deriving unit 743 of the counting result correcting unit 74 regards the signal having the maximum intensity in the frequency analysis result of the frequency distribution frequency analyzing unit 742 as indicating the representative value (mode or median) of the MHP cycle, A representative value T0 of the MHP cycle is obtained.

次に、計数結果補正部74の補正値算出部744は、度数分布作成部741が作成した度数分布から、周期の代表値T0の0.5倍以下である階級の度数の総和Nsと、周期の代表値T0の1.5倍以上である階級の度数の総和Nwとを求め、カウンタ73の計数結果を次式のように補正する。
N’=N+Nw−Ns ・・・(2)
式(2)において、Nはカウンタ73の計数結果であるMHPの数、N’は補正後の計数結果である。
Next, the correction value calculation unit 744 of the counting result correction unit 74 calculates the total frequency Ns of the frequencies that are 0.5 times or less of the representative value T0 of the cycle from the frequency distribution created by the frequency distribution creation unit 741, and the cycle. And a total sum Nw of class frequencies that is 1.5 times or more of the representative value T0 is corrected, and the count result of the counter 73 is corrected as follows.
N ′ = N + Nw−Ns (2)
In Equation (2), N is the number of MHPs that are the counting result of the counter 73, and N ′ is the corrected counting result.

図10にMHPの周期の度数分布の1例を示す。図10において、Tsは周期の代表値T0の0.5倍の階級値、Twは代表値T0の1.5倍の階級値である。図10における階級が、MHPの周期の代表値であることは言うまでもない。なお、図10では記載を簡略化するため、代表値T0とTsとの間、及び代表値T0とTwとの間の度数分布を省略している。   FIG. 10 shows an example of the frequency distribution of the MHP cycle. In FIG. 10, Ts is a class value 0.5 times the representative value T0 of the period, and Tw is a class value 1.5 times the representative value T0. Needless to say, the class in FIG. 10 is a representative value of the MHP cycle. In order to simplify the description in FIG. 10, the frequency distribution between the representative values T0 and Ts and between the representative values T0 and Tw is omitted.

図11はカウンタ73の計数結果の補正原理を説明するための図であり、図11(A)はフィルタ回路6の出力電圧の波形、すなわちMHPの波形を模式的に示す図、図11(B)は図11(A)に対応するカウンタ73の計数結果を示す図である。
本来、MHPの周期は測定対象10との距離によって異なるが、測定対象10との距離が不変であれば、MHPは同じ周期で出現する。しかし、ノイズのために、MHPの波形には欠落が生じたり、信号として数えるべきでない波形が生じたりして、MHPの数に誤差が生じる。
FIG. 11 is a diagram for explaining the principle of correcting the counting result of the counter 73. FIG. 11A is a diagram schematically showing the waveform of the output voltage of the filter circuit 6, that is, the waveform of MHP. () Is a diagram showing the counting result of the counter 73 corresponding to FIG.
Originally, the period of MHP differs depending on the distance to the measurement object 10, but if the distance to the measurement object 10 is unchanged, the MHP appears in the same period. However, due to noise, the MHP waveform may be missing or a waveform that should not be counted as a signal may be generated, resulting in an error in the number of MHPs.

信号の欠落が生じると、欠落が生じた箇所でのMHPの周期Twは、本来の周期のおよそ2倍になる。つまり、MHPの周期が代表値T0のおよそ2倍以上の場合には、信号に欠落が生じていると判断できる。そこで、周期Tw以上の階級の度数の総和Nwを信号が欠落した回数と見なし、このNwをカウンタ73の計数結果Nに加算することで、信号の欠落を補正することができる。   When signal loss occurs, the MHP cycle Tw at the location where the loss occurs is approximately twice the original cycle. That is, when the MHP cycle is approximately twice or more than the representative value T0, it can be determined that the signal is missing. Accordingly, the total frequency Nw of the classes having the period Tw or more is regarded as the number of missing signals, and the missing signal can be corrected by adding this Nw to the counting result N of the counter 73.

また、ノイズをカウントした箇所でのMHPの周期Tsは、本来の周期のおよそ0.5倍になる。つまり、MHPの周期が代表値のおよそ0.5倍以下の場合には、信号を過剰に数えていると判断できる。そこで、周期Ts以下の階級の度数の総和Nsを信号を過剰に数えた回数と見なし、このNsをカウンタ73の計数結果Nから減算することで、誤って数えたノイズを補正することができる。   Further, the MHP cycle Ts at the location where noise is counted is approximately 0.5 times the original cycle. That is, when the MHP cycle is approximately 0.5 times or less of the representative value, it can be determined that the signals are excessively counted. Accordingly, the sum Ns of the frequencies of the class having a period of Ts or less is regarded as the number of times that the signal is excessively counted, and the Ns is subtracted from the count result N of the counter 73, thereby correcting the erroneously counted noise.

以上が、式(2)に示した計数結果の補正原理である。なお、本実施の形態では、Tsを周期の代表値T0の0.5倍の値とし、Twを代表値T0の2倍の値とせずに、1.5倍の値としているが、1.5倍とした理由については後述する。   The above is the correction principle of the counting result shown in Expression (2). In this embodiment, Ts is set to a value that is 0.5 times the representative value T0 of the cycle, and Tw is set to a value that is 1.5 times instead of a value that is twice the representative value T0. The reason for the five times will be described later.

補正値算出部744は、式(2)により計算した補正後の計数結果N’を演算装置8に出力する。計数装置7は、以上のような処理を第1の発振期間t−1,t+1,t+3毎及び第2の発振期間t,t+2,t+4毎に行う。   The correction value calculation unit 744 outputs the corrected count result N ′ calculated by Expression (2) to the arithmetic device 8. The counting device 7 performs the above processing every first oscillation period t-1, t + 1, t + 3 and every second oscillation period t, t + 2, t + 4.

次に、演算装置8は、計数装置7によって計測されたMHPの数N’に基づいて測定対象10との距離を求める。一定期間におけるMHPの数は測定距離に比例する。そこで、一定の計数期間(本実施の形態では第1の発振期間と第2の発振期間の各々)におけるMHPの数と距離との関係を予め求めて演算装置8のデータベース(不図示)に登録しておけば、演算装置8は、計数装置7によって計測されたMHPの数N’に対応する距離の値をデータベースから取得することにより、測定対象10との距離を求めることができる。   Next, the arithmetic device 8 obtains the distance from the measurement object 10 based on the number N ′ of MHPs measured by the counting device 7. The number of MHPs in a certain period is proportional to the measurement distance. Therefore, the relationship between the number of MHPs and the distance in a certain counting period (in this embodiment, each of the first oscillation period and the second oscillation period) is obtained in advance and registered in a database (not shown) of the arithmetic unit 8. Then, the arithmetic device 8 can obtain the distance to the measurement object 10 by obtaining the distance value corresponding to the number N ′ of MHPs measured by the counting device 7 from the database.

あるいは、計数期間におけるMHPの数と距離との関係を示す数式を予め求めて設定しておけば、演算装置8は、計数装置7によって計測されたMHPの数N’を数式に代入することにより、測定対象10との距離を算出することができる。演算装置8は、以上のような処理を第1の発振期間t−1,t+1,t+3毎及び第2の発振期間t,t+2,t+4毎に行う。
表示装置9は、演算装置8によって算出された測定対象10との距離(変位)をリアルタイムで表示する。
Alternatively, if a mathematical expression indicating the relationship between the number of MHPs and the distance in the counting period is obtained and set in advance, the arithmetic device 8 substitutes the number N ′ of MHPs measured by the counting device 7 into the mathematical formula. The distance to the measurement object 10 can be calculated. The arithmetic unit 8 performs the above-described processing every first oscillation period t-1, t + 1, t + 3 and every second oscillation period t, t + 2, t + 4.
The display device 9 displays the distance (displacement) from the measurement object 10 calculated by the arithmetic device 8 in real time.

以上のように、本実施の形態では、計数期間中のMHPの周期を測定し、この測定結果から計数期間中のMHPの周期の度数分布を作成し、度数分布の周波数成分の解析結果からMHPの周期の代表値と見なすことができる値を求めるので、例えば遠距離測定時に信号のS/N比が低下している場合であっても、MHPの真の周期を求めることができる。そして、本実施の形態では、度数分布から、代表値の0.5倍以下である階級の度数の総和Nsと、代表値の1.5倍以上である階級の度数の総和Nwとを求め、これらの度数NsとNwに基づいてカウンタの計数結果を補正することにより、MHPの計数誤差を補正することができるので、距離の測定精度を向上させることができる。   As described above, in the present embodiment, the MHP cycle during the counting period is measured, a frequency distribution of the MHP cycle during the counting period is created from the measurement result, and the MHP is analyzed from the frequency component analysis result of the frequency distribution. Since a value that can be regarded as a representative value of the period is obtained, the true period of the MHP can be obtained even when the S / N ratio of the signal is reduced, for example, during long-distance measurement. Then, in the present embodiment, from the frequency distribution, a total frequency Ns of class frequencies that is 0.5 times or less of the representative value and a total frequency Nw of class frequencies that are 1.5 times or more of the representative value are obtained, By correcting the counting result of the counter based on these frequencies Ns and Nw, the counting error of MHP can be corrected, so that the distance measurement accuracy can be improved.

なお、本実施の形態における計数装置7と演算装置8は、例えばCPU、記憶装置およびインタフェースを備えたコンピュータとこれらのハードウェア資源を制御するプログラムによって実現することができる。このようなコンピュータを動作させるためのプログラムは、フレキシブルディスク、CD−ROM、DVD−ROM、メモリカードなどの記録媒体に記録された状態で提供される。CPUは、読み込んだプログラムを記憶装置に書き込み、このプログラムに従って本実施の形態で説明した処理を実行する。   Note that the counting device 7 and the arithmetic device 8 in the present embodiment can be realized by, for example, a computer including a CPU, a storage device, and an interface, and a program for controlling these hardware resources. A program for operating such a computer is provided in a state of being recorded on a recording medium such as a flexible disk, a CD-ROM, a DVD-ROM, or a memory card. The CPU writes the read program into the storage device, and executes the processing described in this embodiment in accordance with this program.

次に、MHPの基準周期として周期の度数分布の代表値を用いる理由、及び度数Nwを求める際の周期のしきい値を代表値の1.5倍とする理由について説明する。
最初に、誤ってノイズを数えてしまったために、MHPの周期が2分割された場合の計数結果の補正について説明する。半導体レーザの発振波長変化が線形である場合、MHPの周期Tは計測期間TcをMHPの数Nで除算したT0を中心にして正規分布する(図12)。
Next, the reason why the representative value of the frequency distribution of the cycle is used as the MHP reference cycle and the reason why the threshold value of the cycle when determining the frequency Nw is 1.5 times the representative value will be described.
First, correction of the counting result when the MHP cycle is divided into two because the noise is erroneously counted will be described. When the oscillation wavelength change of the semiconductor laser is linear, the period T of MHP is normally distributed around T0 obtained by dividing the measurement period Tc by the number N of MHPs (FIG. 12).

次に、ノイズによって2分割されたMHPの周期を考える。ノイズを過剰に数えた結果として2分割されたMHPの周期は、ランダムな割合で2分割されるが、分割される前の周期がT0を中心とした正規分布であるために、0.5T0に対して対称な度数分布になる(図13のa)。   Next, consider the period of MHP divided into two by noise. The period of MHP divided into two as a result of excessive noise counting is divided into two at a random rate, but since the period before the division is a normal distribution centered on T0, the period is 0.5T0. On the other hand, the frequency distribution is symmetrical (a in FIG. 13).

このノイズを含むMHPの周期の度数分布について、MHPのn%がノイズによって周期が2分割されたと仮定したとき、MHPの周期の平均値及び中央値を算出する。
全ての周期の和は常に計測期間Tcであり、変化はないが、MHPのn%がノイズによって周期が2分割されると、度数の積分値は(1+n[%])Nになるため、MHPの周期の平均値は(1/(1+n[%]))T0になる。
Assuming that the frequency distribution of the MHP cycle including noise is such that n% of the MHP is divided into two by the noise, the average value and the median value of the MHP cycle are calculated.
The sum of all the periods is always the measurement period Tc, and there is no change. However, when n% of MHP is divided into two periods by noise, the integral value of the frequency becomes (1 + n [%]) N. The average value of the period is (1 / (1 + n [%])) T0.

一方、ノイズの分布で正規分布に重なったところを無視した場合、2分割されたノイズの累積度数は中央値とT0との間の階級に含まれる度数の2倍になるため、MHPの周期の中央値は図14のbの面積がaの面積の2倍になる位置になる。   On the other hand, when ignoring the noise distribution that overlaps the normal distribution, the cumulative frequency of the noise divided into two is twice the frequency included in the class between the median and T0. The median is a position where the area of b in FIG. 14 is twice the area of a.

マイクロソフト社のソフトウェアであるExcel(登録商標)に、正規分布の平均値からασ間の両側値の内部割合が「(1−(1−NORMSDIST(α))*2)*100[%]」で表現できるNORMSDIST()という関数があり、この関数を利用すると、MHPの周期の中央値を以下の式で表すことができる。
(1−(1−NORMSDIST((中央値−T0)/σ))*2)
*(100−n)/2=n[%] ・・・(3)
Excel (registered trademark), which is Microsoft software, has an internal ratio of two-sided values between the mean value of α distribution and ασ as “(1- (1-NORMDIST (α)) * 2) * 100 [%]”. There is a function called NORMSDIST () that can be expressed. By using this function, the median value of the MHP cycle can be expressed by the following equation.
(1- (1-NORMDIST ((median-T0) / σ)) * 2)
* (100-n) / 2 = n [%] (3)

以上をもとに、標準偏差σを0.02T0とし、MHPの10%がノイズによって周期が2分割されたとしたときのMHPの周期の平均値T0’および中央値T0’を算出すると、以下のようになる。
T0’=(1/(1+0.1))T0=0.91T0 ・・・(4)
T0’=0.995T0 ・・・(5)
なお、ここでは平均値、中央値共にT0’で表すものとする。カウンタ値(度数の積分値)は、1.1Nとなり、カウント誤差は10%となる。
Based on the above, when the standard deviation σ is 0.02T0, and the average value T0 ′ and median value T0 ′ of the MHP period when 10% of MHP is divided into two periods by noise, the following values are calculated. It becomes like this.
T0 ′ = (1 / (1 + 0.1)) T0 = 0.91T0 (4)
T0 '= 0.995T0 (5)
Here, both the average value and the median value are represented by T0 ′. The counter value (frequency integrated value) is 1.1 N, and the count error is 10%.

ここで、ある周期TaのMHPが2分割された後の2つの周期T1,T2(T1≧T2とする)のとり得る期間の確率を考える。ノイズはランダムに生じると仮定すると、図15に示すようにT2は0<T2≦Ta/2の値を同じ確率でとり得る。T1も同様にT/2≦T1<Taの値を同じ確率でとり得る。図15におけるT1の取り得る確率分布の面積とT2の取り得る確率分布の面積は共に1である。   Here, consider the probability of a period that can be taken by two periods T1 and T2 (T1 ≧ T2) after an MHP of a certain period Ta is divided into two. Assuming that noise occurs randomly, T2 can take a value of 0 <T2 ≦ Ta / 2 with the same probability as shown in FIG. Similarly, T1 can take the value of T / 2 ≦ T1 <Ta with the same probability. In FIG. 15, the area of the probability distribution that T1 can take and the area of the probability distribution that T2 can take are both 1.

周期TaはT0を中心とした正規分布をしているので、Taを集合としてとらえると、T2のとり得る確率の度数分布は、平均値が0.5T0、標準偏差0.5σの正規分布の累積度数分布と同じ形状になる。
また、図16に示すように、T1のとり得る確率の度数分布は平均値が0.5T0、標準偏差0.5σの正規分布の累積度数分布と平均値がT0、標準偏差σの正規分布の累積度数分布を重ねたような形状になる。ここで、T1、T2それぞれの数は、周期が2分割されたMHPの数n[%]・Nに等しい。
Since the period Ta has a normal distribution centered on T0, if Ta is regarded as a set, the frequency distribution of the probability that T2 can take is the accumulation of normal distributions with an average value of 0.5T0 and a standard deviation of 0.5σ. It has the same shape as the frequency distribution.
Further, as shown in FIG. 16, the frequency distribution of the probability that T1 can take is a cumulative frequency distribution of a normal distribution with an average value of 0.5T0 and a standard deviation of 0.5σ, and a normal distribution with an average value of T0 and a standard deviation of σ. The shape is such that the cumulative frequency distribution is superimposed. Here, the number of each of T1 and T2 is equal to the number n [%] · N of MHPs whose period is divided into two.

ノイズによって周期が2分割されたMHPの数n[%]・Nを数えることができれば、以下の式を用いてMHPの数Nを導出することができる。
N=N’−n[%]・N ・・・(6)
図17に示すように、Tb以下の周期を持つMHPの数Nsが2分割されたMHPの数n[%]・Nと等しくなるようにTbを設定することができれば、Tb以下の周期を持つMHPの数Nsを数えることで、周期が2分割されたMHPの数n[%]・Nを間接的に数えることができる。
If the number n [%] · N of MHPs whose period is divided into two by noise can be counted, the number N of MHPs can be derived using the following equation.
N = N′−n [%] · N (6)
As shown in FIG. 17, if Tb can be set so that the number Ns of MHPs having a period equal to or less than Tb is equal to the number n [%] · N of MHPs divided into two, the period is equal to or less than Tb. By counting the number Ns of MHPs, the number n [%] · N of MHPs whose periods are divided into two can be indirectly counted.

図17において、Tb以上の周期を持つMHPの周期T2の度数(図17のc)とTb未満の周期を持つMHPの周期T1の度数(図17のd)が同じになるとき、Tb以下の周期を持つMHPの数は、T2の数、つまり周期が2分割されたMHPの数Ns(=n[%]・N)と等しくなる。つまり、MHPの数Nは以下の式で表すことができる。
N=N’−n[%]・N=N’−Ns ・・・(7)
T1およびT2の度数形状は、0.5Taで対称の形状であるため、0.5Taをしきい値にして判断すると、周期が2分割されたMHPの度数Ns(=n[%]・N)を正確に数えることができる。
In FIG. 17, when the frequency of the MHP cycle T2 having a cycle of Tb or more (c in FIG. 17) and the frequency of the MHP cycle T1 having a cycle of less than Tb (d in FIG. 17) are equal to or less than Tb The number of MHPs having a period is equal to the number of T2, that is, the number Ns (= n [%] · N) of MHPs whose periods are divided into two. That is, the number N of MHPs can be expressed by the following formula.
N = N′−n [%] · N = N′−Ns (7)
Since the frequency shapes of T1 and T2 are symmetric shapes with 0.5 Ta, if the determination is made with 0.5 Ta as a threshold value, the frequency Ns (= n [%] · N) of the MHP whose period is divided into two Can be counted accurately.

次に、0.5T0以下の周期を持つMHPの数を数えることで、周期が2分割されたMHPの数n[%]・Nの数を間接的に数えることができるが、ノイズを含むMHPの周期の度数分布(図13)からは、T0を算出することができない。MHPの母集団が図13の度数分布のように最頻値(モード)がT0と等しくなるほど理想的でかつ母数が大きければ、最頻値をT0’として用いることができる。   Next, by counting the number of MHPs having a period of 0.5T0 or less, the number of MHPs n [%] · N whose period is divided into two can be indirectly counted. T0 cannot be calculated from the frequency distribution of the period (FIG. 13). If the MHP population is ideal and the parameter is large enough that the mode (mode) is equal to T0 as in the frequency distribution of FIG. 13, the mode can be used as T0 '.

ここでは、平均値又は中央値T0’を用いたMHPの数n[%]・Nの計数について記載する。T0’=y・T0で表し、T0の代わりにT0’を代入してNsを求めると、周期が2分割されたMHPの数として判断する0.5T0’よりも小さな周期の度数Ns’は、y・n[%]・Nになる(図18)。   Here, the counting of the number n [%] · N of MHPs using the average value or the median value T0 ′ will be described. T0 ′ = y · T0 and substituting T0 ′ for T0 to obtain Ns, the frequency Ns ′ having a period smaller than 0.5T0 ′, which is determined as the number of MHPs whose periods are divided into two, is y · n [%] · N (FIG. 18).

平均値又は中央値T0’を用いた場合、補正後のカウント値Ntは以下のように表される。
Nt=N’−Ns’=(1+n[%])N−yn[%]N
=(1+(1−y)n[%])N=N+(1−y)n[%]N ・・・(8)
なお、補正後の誤差である(1−y)n[%]Nは、図19のeの部分の度数である。
When the average value or the median value T0 ′ is used, the corrected count value Nt is expressed as follows.
Nt = N′−Ns ′ = (1 + n [%]) N−yn [%] N
= (1+ (1-y) n [%]) N = N + (1-y) n [%] N (8)
Note that (1-y) n [%] N, which is an error after correction, is the frequency of the portion e in FIG.

ここで、平均値又は中央値T0’を用いたカウンタ73の計数結果の補正例について説明する。
標準偏差をσ=0.02T0とし、MHPの10%がノイズによって周期が2分割されたとすると(計数結果は10%の誤差)、MHPの周期の平均値T0’は0.91T0、中央値T0’は0.9949T0であるから、平均値T0’を用いる場合のyは0.91、中央値T0’を用いる場合のyは0.9949であり、補正後の計数結果N’は以下のように算出される。
N’=(1+0.1(1−0.91))N=1.009N ・・・(9)
N’=(1+0.1(1−0.995))N=1.0005N ・・・(10)
Here, an example of correcting the counting result of the counter 73 using the average value or the median value T0 ′ will be described.
Assuming that the standard deviation is σ = 0.02T0 and 10% of the MHP is divided into two periods by noise (the count result is an error of 10%), the average value T0 ′ of the MHP period is 0.91T0, and the median value T0 Since 'is 0.9949T0, y when using the average value T0' is 0.91, and y when using the median value T0 'is 0.9949, and the corrected count result N' is as follows: Is calculated.
N ′ = (1 + 0.1 (1−0.91)) N = 1.09N (9)
N ′ = (1 + 0.1 (1−0.995)) N = 1.0005N (10)

式(9)は平均値T0’を用いた場合の補正後の計数結果N’を示し、式(10)は中央値T0’を用いた場合の補正後の計数結果N’を示している。平均値T0’を用いた場合の計数結果N’の誤差は0.9%であり、中央値T0’を用いた場合の計数結果N’の誤差は0.05%である。   Equation (9) shows the corrected count result N 'when the average value T0' is used, and Equation (10) shows the corrected count result N 'when the median value T0' is used. The error of the count result N ′ when the average value T0 ′ is used is 0.9%, and the error of the count result N ′ when the median value T0 ′ is used is 0.05%.

次に、標準偏差をσ=0.05T0とし、MHPの20%がノイズによって周期が2分割されたとすると(計数結果は20%の誤差)、MHPの周期の平均値T0’は0.83T0、中央値T0’は0.9682T0であるから、平均値T0’を用いる場合のyは0.83、中央値T0’を用いる場合のyは0.968であり、補正後の計数結果N’は以下のように算出される。
N’=(1+0.2(1−0.83))N=1.034N ・・・(11)
N’=(1+0.2(1−0.968))N=1.0064N ・・・(12)
Next, assuming that the standard deviation is σ = 0.05T0 and the period of the MHP is divided into two by noise (the counting result is an error of 20%), the average value T0 ′ of the MHP period is 0.83T0, Since the median value T0 ′ is 0.9682T0, y when using the average value T0 ′ is 0.83, and y when using the median value T0 ′ is 0.968, and the corrected count result N ′ is It is calculated as follows.
N ′ = (1 + 0.2 (1−0.83)) N = 1.304N (11)
N ′ = (1 + 0.2 (1−0.968)) N = 1.0064N (12)

式(11)は平均値T0’を用いた場合の補正後の計数結果N’を示し、式(12)は中央値T0’を用いた場合の補正後の計数結果N’を示している。平均値T0’を用いた場合の計数結果N’の誤差は3.4%であり、中央値T0’を用いた場合の計数結果N’の誤差は0.64%である。
以上のことから、MHPの周期の代表値を使用して計数結果Nを補正すれば、補正後の計数結果N’の誤差を小さくできることが分かる。
Equation (11) shows the corrected count result N ′ when the average value T0 ′ is used, and Equation (12) shows the corrected count result N ′ when the median value T0 ′ is used. The error of the count result N ′ when the average value T0 ′ is used is 3.4%, and the error of the count result N ′ when the median value T0 ′ is used is 0.64%.
From the above, it can be seen that if the count result N is corrected using the representative value of the MHP cycle, the error in the corrected count result N ′ can be reduced.

次に、MHPの波形に欠落が生じた場合の計数結果の補正について説明する。MHPの強度が小さいために計数時に欠落が生じた場合のMHPの周期は、本来のMHPの周期がT0を中心とした正規分布であるために、平均値が2T0、標準偏差2σの正規分布(図20のf)になる。m[%]のMHPが欠落したとすると、この欠落によって周期が2倍になったMHPの周期の度数はNw(=m[%]・N)である。また、計数時の欠落によって減少した後のおおよそT0の周期の度数は、図20に示すgであり、図20のhに示す度数の減少分は2Nw(=2m[%])である。したがって、計数時にMHPの欠落が生じなかった場合の本来のMHPの数N’は以下の式で表すことができる。
N’=N+m[%]=N+Nw ・・・(13)
Next, correction of the counting result when a loss occurs in the MHP waveform will be described. The period of MHP when a loss occurs during counting because the intensity of MHP is small is a normal distribution centered on T0 of the original MHP period, and therefore a normal distribution with an average value of 2T0 and a standard deviation of 2σ ( It becomes f) of FIG. Assuming that m [%] MHP is missing, the frequency of the MHP period whose period is doubled due to this missing is Nw (= m [%] · N). Further, the frequency of the period of approximately T0 after being decreased due to omission at the time of counting is g shown in FIG. 20, and the decrease in the frequency shown in h of FIG. 20 is 2Nw (= 2m [%]). Therefore, the original number N ′ of MHPs when no MHP is lost during counting can be expressed by the following equation.
N ′ = N + m [%] = N + Nw (13)

次に、計数結果を補正するためのNwを数える際の周期のしきい値について考える。ここで、計数時の欠落によって周期が2倍になったMHPの周期の度数Nwのうちノイズによってp[%]が2分割された場合を仮定する。欠落したMHPのうち2分割されたMHPの周期の度数は、Nw’(=m・p[%]・N)である。再度2分割されたMHPの周期の度数分布は、図21のようになる。Nwとみなす周期のしきい値を1.5T0にすると、周期が0.5T0以下のMHPの周期の度数は0.5Nw’(=0.5p[%]・Nw)、周期が0.5T0から1.5T0までのMHPの周期の度数はNw’(=p[%]・Nw)、周期が1.5T0以上のMHPの周期の度数は0.5Nw’(=0.5p[%]・Nw)となる。   Next, let us consider the threshold of the period when counting Nw for correcting the counting result. Here, it is assumed that p [%] is divided into two by noise out of the frequency Nw of the MHP period whose period is doubled due to omission at the time of counting. Of the missing MHPs, the frequency of the MHP period divided into two is Nw ′ (= m · p [%] · N). The frequency distribution of the period of the MHP divided into two again is as shown in FIG. When the threshold value of the period regarded as Nw is 1.5T0, the frequency of the MHP period with a period of 0.5T0 or less is 0.5Nw ′ (= 0.5p [%] · Nw), and the period is from 0.5T0. The frequency of the MHP period up to 1.5T0 is Nw ′ (= p [%] · Nw), and the frequency of the MHP period of 1.5T0 or more is 0.5 Nw ′ (= 0.5 p [%] · Nw). )

よって、全てのMHPの周期の度数分布は図22のようになり、Nsのしきい値を0.5T0、Nwのしきい値を1.5T0にすると、計数結果Nは以下の式で表すことができる。
N=(N’−2Nw)+(Nw−Nw’)+2Nw’=N’−Nw+Nw’
・・・(14)
Therefore, the frequency distribution of all MHP cycles is as shown in FIG. 22. When the threshold value of Ns is 0.5T0 and the threshold value of Nw is 1.5T0, the counting result N is expressed by the following equation. Can do.
N = (N′−2Nw) + (Nw−Nw ′) + 2Nw ′ = N′−Nw + Nw ′
(14)

式(14)より補正された結果は以下のようになり、計数時にMHPの欠落が生じなかった場合の本来のMHPの数N’が算出されることが分かる。
N−0.5Nw’+(0.5Nw’+(Nw−Nw’))
=(N−Nw+Nw’)+(0.5Nw’+(Nw−Nw’))
=N’ ・・・(15)
The corrected result from the equation (14) is as follows, and it can be seen that the original number N ′ of MHPs is calculated when no MHP is lost during the counting.
N−0.5Nw ′ + (0.5Nw ′ + (Nw−Nw ′))
= (N-Nw + Nw ') + (0.5Nw' + (Nw-Nw '))
= N '(15)

以上のことから、度数Nwを求める際の周期のしきい値を代表値の1.5倍とすれば、計数結果Nを補正できることが分かる。なお、ノイズによってMHPの周期が2分割された場合と同様に、T0の代わりに代表値を用いて補正するため、同様の誤差が生じる。   From the above, it can be seen that the counting result N can be corrected by setting the threshold value of the cycle for obtaining the frequency Nw to 1.5 times the representative value. Similar to the case where the period of MHP is divided into two by noise, correction is performed using a representative value instead of T0, so that the same error occurs.

以上の説明では、ノイズを過剰に数えた結果MHPの周期が2分割された場合と計数時の欠落によってMHPの周期が2倍になった場合を別々に説明したが、これらは独立して生じるため、これらの場合を1つの度数分布に表現すると、図23のようになる。Nsのしきい値を0.5T0、Nwのしきい値を1.5T0にすると、計数結果Nは以下の式で表すことができる。
N=(N’−2Nw−Ns)+(Nw−Nw’)+2Nw’+2Ns
=N’−Nw+Nw’+Ns ・・・(16)
In the above description, the case where the MHP cycle is divided into two as a result of excessive noise counting and the case where the MHP cycle is doubled due to omission at the time of counting have been described separately, but these occur independently. Therefore, these cases are expressed as one frequency distribution as shown in FIG. When the threshold value of Ns is 0.5T0 and the threshold value of Nw is 1.5T0, the count result N can be expressed by the following equation.
N = (N′−2Nw−Ns) + (Nw−Nw ′) + 2Nw ′ + 2Ns
= N'-Nw + Nw '+ Ns (16)

式(16)より補正された結果は以下のようになり、計数時に欠落や過剰な計数が生じなかった場合の本来のMHPの数N’が算出されることが分かる。
N−{0.5Nw’+Ns}+{0.5Nw’+(Nw−Nw’)}
={N−Nw+Nw’+Ns}−{0.5Nw’+Ns}
+{0.5Nw’+(Nw−Nw’)}
=N’ ・・・(17)
The corrected result from the equation (16) is as follows, and it can be seen that the original number N ′ of MHPs when there is no omission or excessive counting at the time of counting is calculated.
N− {0.5Nw ′ + Ns} + {0.5Nw ′ + (Nw−Nw ′)}
= {N−Nw + Nw ′ + Ns} − {0.5Nw ′ + Ns}
+ {0.5 Nw ′ + (Nw−Nw ′)}
= N '(17)

[第2の実施の形態]
次に、本発明の第2の実施の形態について説明する。第1の実施の形態では、度数分布周波数解析部742の周波数解析結果において強度が最大の信号をMHPの周期の代表値を示すものと見なした。しかし、ノイズが極めて多い環境では、図24に示すようにMHPの周期の度数分布の極大値の間隔に相当する周波数1/T0ではなく、度数分布の極大値の2倍の間隔に相当する周波数1/(2T0)の方が強度が最大になる場合がある。
[Second Embodiment]
Next, a second embodiment of the present invention will be described. In the first embodiment, the signal having the maximum intensity in the frequency analysis result of the frequency distribution frequency analysis unit 742 is regarded as indicating the representative value of the MHP cycle. However, in an extremely noisy environment, as shown in FIG. 24, not the frequency 1 / T0 corresponding to the maximum value interval of the frequency distribution of the MHP cycle, but the frequency corresponding to twice the maximum value of the frequency distribution. In some cases, 1 / (2T0) has the maximum intensity.

そこで、ノイズが極めて多いと考えられる環境では、計数結果補正部74の代表値導出部743は、度数分布周波数解析部742の周波数解析結果において複数の極大値のうち周波数が最も高い信号をMHPの周期の代表値を示すものと見なし、MHPの周期の代表値T0を求めるようにしてもよい。   Therefore, in an environment that is considered to be extremely noisy, the representative value deriving unit 743 of the counting result correcting unit 74 outputs a signal having the highest frequency among a plurality of maximum values in the frequency analysis result of the frequency distribution frequency analyzing unit 742 to the MHP. The representative value T0 of the MHP period may be obtained by regarding the representative value of the period.

なお、第1、第2の実施の形態では、本発明の計数装置を距離計に適用した場合について説明したが、これに限るものではなく、本発明の計数装置は他の分野にも適用することができる。本発明の計数装置が有効な場合は、計数の対象となる信号が特定の物理量(第1、第2の実施の形態の場合は距離)と線形の関係を有し、物理量が一定の場合は信号が略単一周波数となる場合である。
また、信号が単一周波数でなくても、計数対象となる特定の物理量が計数期間と比較して十分低い周波数で、例えば1/10以下の周波数で振動している対象物の速度のように周期分布の広がりが小さい場合も略単一周波数として本発明の計数装置は有効である。
In the first and second embodiments, the case where the counting device of the present invention is applied to a distance meter has been described. However, the present invention is not limited to this, and the counting device of the present invention is also applied to other fields. be able to. When the counting device of the present invention is effective, the signal to be counted has a linear relationship with a specific physical quantity (distance in the case of the first and second embodiments), and the physical quantity is constant. This is the case where the signal has a substantially single frequency.
Even if the signal is not a single frequency, the specific physical quantity to be counted is sufficiently low compared to the counting period, for example, the speed of an object vibrating at a frequency of 1/10 or less. Even when the spread of the periodic distribution is small, the counting device of the present invention is effective as a substantially single frequency.

また、第1、第2の実施の形態では、MHPの欠落の補正については、1個の欠落によってMHPの周期が本来の周期のおよそ2倍になった場合について説明しているが、連続して2個以上の欠落が生じた場合にも本発明を適用することができる。MHPが連続して2個欠落した場合、中央値の3倍の周期のMHPは3個のMHPが1つになったものだと考えられる。この場合は、周期の中央値のおよそ3倍以上である階級の度数を求めて、この度数を2倍すれば、MHPの欠落を補正することができる。このような考え方を一般化すると、式(2)の代わりに次式を用いればよい。
N’=N+Nw1+Nw2+Nw3+・・・・−Ns ・・・(18)
Nw1は周期の中央値の1.5倍以上である階級の度数の総和、Nw2は周期の中央値の2.5倍以上である階級の度数の総和、Nw3は周期の中央値のおよそ3.5倍以上である階級の度数の総和である。
また、第1、第2の実施の形態では、自己結合型の干渉計に本発明を適用する場合について説明したが、自己結合型以外の干渉計に本発明を適用することもできる。
In the first and second embodiments, the correction of the missing MHP has been described for the case where the MHP period is approximately twice the original period due to one missing. Thus, the present invention can be applied even when two or more missing parts occur. If two MHPs are missing in succession, the MHP with a period three times the median is considered to be a combination of three MHPs. In this case, if the frequency of the class that is approximately three times or more the median value of the period is obtained and this frequency is doubled, the missing MHP can be corrected. If such a way of thinking is generalized, the following equation may be used instead of equation (2).
N ′ = N + Nw1 + Nw2 + Nw3 +... -Ns (18)
Nw1 is the sum of the frequencies of the class that is 1.5 times or more of the median of the period, Nw2 is the sum of the frequencies of the class that is 2.5 or more of the median of the period, and Nw3 is about 3 of the median of the period. It is the sum of the frequency of the class that is 5 times or more.
In the first and second embodiments, the case where the present invention is applied to a self-coupled interferometer has been described. However, the present invention can also be applied to interferometers other than the self-coupled interferometer.

本発明は、信号の数を数える計数装置や、計数装置を用いて干渉波形の数を測定し測定対象との距離を求める干渉型の距離計に適用することができる。   The present invention can be applied to a counting device that counts the number of signals and an interference-type distance meter that measures the number of interference waveforms by using the counting device and obtains a distance from a measurement object.

本発明の第1の実施の形態となる距離計の構成を示すブロック図である。It is a block diagram which shows the structure of the distance meter used as the 1st Embodiment of this invention. 本発明の第1の実施の形態における半導体レーザの発振波長の時間変化の1例を示す図である。It is a figure which shows one example of the time change of the oscillation wavelength of the semiconductor laser in the 1st Embodiment of this invention. 本発明の第1の実施の形態における電流−電圧変換増幅器の出力電圧波形及びフィルタ回路の出力電圧波形を模式的に示す図である。It is a figure which shows typically the output voltage waveform of the current-voltage conversion amplifier in the 1st Embodiment of this invention, and the output voltage waveform of a filter circuit. 本発明の第1の実施の形態における計数装置の構成の1例を示すブロック図である。It is a block diagram which shows one example of a structure of the counting device in the 1st Embodiment of this invention. 図4の計数装置における計数結果補正部の構成の1例を示すブロック図である。It is a block diagram which shows one example of a structure of the count result correction | amendment part in the counting device of FIG. 図4の計数装置の動作を説明するための図である。It is a figure for demonstrating operation | movement of the counting apparatus of FIG. 通常の測定条件におけるモードホップパルスの周期の度数分布の1例を示す図である。It is a figure which shows an example of the frequency distribution of the period of the mode hop pulse in normal measurement conditions. S/N比低下時のモードホップパルスの周期の度数分布の1例を示す図である。It is a figure which shows an example of the frequency distribution of the period of the mode hop pulse at the time of S / N ratio fall. モードホップパルスの周期の度数分布を周波数解析した概略結果を示す周波数スペクトル図である。It is a frequency spectrum figure which shows the general | schematic result which frequency-analyzed the frequency distribution of the period of the mode hop pulse. 本発明の第1の実施の形態におけるモードホップパルスの周期の度数分布の1例を示す図である。It is a figure which shows one example of the frequency distribution of the period of the mode hop pulse in the 1st Embodiment of this invention. 本発明の第1の実施の形態におけるカウンタの計数結果の補正原理を説明するための図である。It is a figure for demonstrating the correction principle of the count result of the counter in the 1st Embodiment of this invention. モードホップパルスの周期の度数分布を示す図である。It is a figure which shows frequency distribution of the period of a mode hop pulse. ノイズを含むモードホップパルスの周期の度数分布を示す図である。It is a figure which shows frequency distribution of the period of the mode hop pulse containing noise. ノイズを含むモードホップパルスの周期の中央値を示す図である。It is a figure which shows the median of the period of the mode hop pulse containing a noise. 周期が2分割されたモードホップパルスの周期の確率分布を示す図である。It is a figure which shows the probability distribution of the period of the mode hop pulse by which the period was divided into two. 周期が2分割されたモードホップパルスの周期の度数分布を示す図である。It is a figure which shows frequency distribution of the period of the mode hop pulse by which the period was divided into two. 周期が2分割されたモードホップパルスの周期の度数分布を示す図である。It is a figure which shows frequency distribution of the period of the mode hop pulse by which the period was divided into two. 周期が2分割されたモードホップパルスの周期の度数分布を示す図である。It is a figure which shows frequency distribution of the period of the mode hop pulse by which the period was divided into two. カウンタ値補正後の誤差を示す図である。It is a figure which shows the error after counter value correction | amendment. 2倍の周期になったモードホップパルスの周期の度数分布を示す図である。It is a figure which shows the frequency distribution of the period of the mode hop pulse which became a 2 times period. 計数時に欠落したモードホップパルスのうち2分割されたモードホップパルスの周期の度数分布を示す図である。It is a figure which shows frequency distribution of the period of the mode hop pulse divided into two among the mode hop pulses lost at the time of counting. 計数時に欠落したモードホップパルスのうち2分割されたモードホップパルスの周期の度数分布を示す図である。It is a figure which shows frequency distribution of the period of the mode hop pulse divided into two among the mode hop pulses lost at the time of counting. 計数時に欠落と過剰な計数が同時に発生した場合のモードホップパルスの周期の度数分布を示す図である。It is a figure which shows the frequency distribution of the period of a mode hop pulse when a loss and an excessive count generate | occur | produce simultaneously at the time of a count. モードホップパルスの周期の度数分布を周波数解析した他の概略結果を示す周波数スペクトル図である。It is a frequency spectrum figure which shows the other general result which frequency-analyzed the frequency distribution of the period of the mode hop pulse. 従来のレーザ計測器における半導体レーザの複合共振器モデルを示す図である。It is a figure which shows the compound resonator model of the semiconductor laser in the conventional laser measuring device. 半導体レーザの発振波長と内蔵フォトダイオードの出力波形との関係を示す図である。It is a figure which shows the relationship between the oscillation wavelength of a semiconductor laser, and the output waveform of a built-in photodiode.

符号の説明Explanation of symbols

1…半導体レーザ、2…フォトダイオード、3…レンズ、4…レーザドライバ、5…電流−電圧変換増幅器、6…フィルタ回路、7…計数装置、8…演算装置、9…表示装置、10…測定対象、71…判定部、72…論理積演算部、73…カウンタ、74…計数結果補正部、75…記憶部、740…周期測定部、741…度数分布作成部、742…度数分布周波数解析部、743…代表値導出部、744…補正値算出部。   DESCRIPTION OF SYMBOLS 1 ... Semiconductor laser, 2 ... Photodiode, 3 ... Lens, 4 ... Laser driver, 5 ... Current-voltage conversion amplifier, 6 ... Filter circuit, 7 ... Counting device, 8 ... Arithmetic unit, 9 ... Display device, 10 ... Measurement Object 71, determination unit 72, logical product operation unit 73, counter, 74, counting result correction unit, 75, storage unit, 740, period measurement unit, 741, frequency distribution creation unit, 742, frequency distribution frequency analysis unit 743 ... representative value deriving unit, 744 ... correction value calculating unit.

Claims (18)

特定の物理量と線形の関係を有し前記物理量が一定の場合は略単一周波数となる信号の数を数える計数装置において、
一定の計数期間における入力信号の数を数える計数手段と、
前記計数期間中の前記入力信号の周期を信号が入力される度に測定する周期測定手段と、
この周期測定手段の測定結果から前記計数期間中の信号周期の度数分布を作成する度数分布作成手段と、
前記度数分布をFFTにより解析して、前記度数分布の極大値の間隔に相当する周波数に極大値が出現する周波数スペクトルを求める周期度数分布周波数解析手段と、
この周期度数分布周波数解析手段が求めた周波数スペクトルに出現した極大値から前記入力信号の周期の分布の代表値を求める代表値導出手段と、
前記度数分布から、前記代表値の第1の所定数倍以下である階級の度数の総和Nsと、前記代表値の第2の所定数倍以上である階級の度数の総和Nwとを求め、これらの度数NsとNwに基づいて前記計数手段の計数結果を補正する補正値算出手段とを有することを特徴とする計数装置。
In a counting device that counts the number of signals having a linear relationship with a specific physical quantity and having a substantially single frequency when the physical quantity is constant,
Counting means for counting the number of input signals in a certain counting period;
Period measuring means for measuring the period of the input signal during the counting period each time a signal is input;
Frequency distribution creating means for creating a frequency distribution of the signal period during the counting period from the measurement result of the period measuring means,
And analyzed by FFT the frequency distribution, a frequency analysis means of the periodic frequency distribution to determine the frequency spectrum maximum value appears in the frequency corresponding to the spacing of the maxima of the frequency distribution,
Representative value deriving means for obtaining a representative value of the distribution of the period of the input signal from the maximum value appearing in the frequency spectrum obtained by the frequency analysis means of the periodic frequency distribution;
From the frequency distribution, the sum Ns of the frequencies of the class that is less than or equal to the first predetermined number of times of the representative value and the sum Nw of the frequencies of the class that are greater than or equal to the second predetermined number of times of the representative value are obtained, and And a correction value calculating means for correcting the counting result of the counting means based on the frequencies Ns and Nw.
請求項1記載の計数装置において、
前記代表値導出手段は、前記周期度数分布周波数解析手段が求めた周波数スペクトルにおいて強度が最大の極大値を前記入力信号の周期の分布の代表値を示すものと見なし、前記代表値を求めることを特徴とする計数装置。
The counting device according to claim 1, wherein
The representative value deriving means regards the maximum value having the maximum intensity in the frequency spectrum obtained by the frequency analysis means of the periodic frequency distribution as a representative value of the period distribution of the input signal, and obtains the representative value. A counting device characterized by.
請求項1記載の計数装置において、
前記代表値導出手段は、前記周期度数分布周波数解析手段が求めた周波数スペクトルにおいて周波数が最も高い極大値を前記入力信号の周期の分布の代表値を示すものと見なし、前記代表値を求めることを特徴とする計数装置。
The counting device according to claim 1, wherein
The representative value deriving means regards the maximum value having the highest frequency in the frequency spectrum obtained by the frequency analysis means of the periodic frequency distribution as the representative value of the period distribution of the input signal, and obtains the representative value. A counting device characterized by.
請求項1乃至3のいずれか1項に記載の計数装置において、
前記補正値算出手段は、前記計数手段の計数結果をNとしたとき、補正後の計数結果N’を、N’=N+Nw−Nsにより求めることを特徴とする計数装置。
The counting device according to any one of claims 1 to 3,
The said correction value calculation means calculates | requires the count result N 'after correction | amendment by N' = N + Nw-Ns when the count result of the said count means is set to N. The counting apparatus characterized by the above-mentioned.
測定対象にレーザ光を放射する半導体レーザと、
発振波長が連続的に単調増加する期間を少なくとも含む第1の発振期間と発振波長が連続的に単調減少する期間を少なくとも含む第2の発振期間とが交互に存在するように前記半導体レーザを動作させるレーザドライバと、
前記半導体レーザから放射されたレーザ光と前記測定対象からの戻り光との干渉光を電気信号に変換する受光器と、
この受光器の出力信号に含まれる、前記半導体レーザから放射されたレーザ光と前記測定対象からの戻り光とによって生じる干渉波形の数を数える計数手段と、
前記干渉波形の数を数える計数期間中の前記干渉波形の周期を干渉波形が入力される度に測定する周期測定手段と、
この周期測定手段の測定結果から前記計数期間中の干渉波形の周期の度数分布を作成する度数分布作成手段と、
前記度数分布をFFTにより解析して、前記度数分布の極大値の間隔に相当する周波数に極大値が出現する周波数スペクトルを求める周期度数分布周波数解析手段と、
この周期度数分布周波数解析手段が求めた周波数スペクトルに出現した極大値から前記干渉波形の周期の分布の代表値を求める代表値導出手段と、
前記度数分布から、前記代表値の第1の所定数倍以下である階級の度数の総和Nsと、前記代表値の第2の所定数倍以上である階級の度数の総和Nwとを求め、これらの度数NsとNwに基づいて前記計数手段の計数結果を補正する補正値算出手段と、
この補正値算出手段で補正された計数結果から前記測定対象との距離を求める演算手段とを有することを特徴とする距離計。
A semiconductor laser that emits laser light to the object to be measured;
The semiconductor laser is operated so that a first oscillation period including at least a period in which the oscillation wavelength continuously increases monotonically and a second oscillation period including at least a period in which the oscillation wavelength continuously decreases monotonously exist. A laser driver,
A light receiver that converts interference light between laser light emitted from the semiconductor laser and return light from the measurement object into an electrical signal;
Counting means for counting the number of interference waveforms generated by the laser light emitted from the semiconductor laser and the return light from the measurement object included in the output signal of the light receiver;
A period measuring means for measuring a period of the interference waveform during a counting period for counting the number of the interference waveforms every time the interference waveform is input;
Frequency distribution creating means for creating a frequency distribution of the period of the interference waveform during the counting period from the measurement result of the period measuring means,
And analyzed by FFT the frequency distribution, a frequency analysis means of the periodic frequency distribution to determine the frequency spectrum maximum value appears in the frequency corresponding to the spacing of the maxima of the frequency distribution,
Representative value deriving means for obtaining a representative value of the period distribution of the interference waveform from the maximum value appearing in the frequency spectrum obtained by the frequency analysis means of the periodic frequency distribution;
From the frequency distribution, the sum Ns of the frequencies of the class that is less than or equal to the first predetermined number of times of the representative value and the sum Nw of the frequencies of the class that are greater than or equal to the second predetermined number of times of the representative value are obtained, and Correction value calculating means for correcting the counting result of the counting means based on the frequencies Ns and Nw of
A distance meter comprising: a calculation unit that obtains a distance from the measurement object from a count result corrected by the correction value calculation unit.
測定対象にレーザ光を放射する半導体レーザと、
発振波長が連続的に単調増加する期間を少なくとも含む第1の発振期間と発振波長が連続的に単調減少する期間を少なくとも含む第2の発振期間とが交互に存在するように前記半導体レーザを動作させるレーザドライバと、
前記半導体レーザの光出力を電気信号に変換する受光器と、
この受光器の出力信号に含まれる、前記半導体レーザから放射されたレーザ光と前記測定対象からの戻り光との自己結合効果によって生じる干渉波形の数を数える計数手段と、
前記干渉波形の数を数える計数期間中の前記干渉波形の周期を干渉波形が入力される度に測定する周期測定手段と、
この周期測定手段の測定結果から前記計数期間中の干渉波形の周期の度数分布を作成する度数分布作成手段と、
前記度数分布をFFTにより解析して、前記度数分布の極大値の間隔に相当する周波数に極大値が出現する周波数スペクトルを求める周期度数分布周波数解析手段と、
この周期度数分布周波数解析手段が求めた周波数スペクトルに出現した極大値から前記干渉波形の周期の分布の代表値を求める代表値導出手段と、
前記度数分布から、前記代表値の第1の所定数倍以下である階級の度数の総和Nsと、前記代表値の第2の所定数倍以上である階級の度数の総和Nwとを求め、これらの度数NsとNwに基づいて前記計数手段の計数結果を補正する補正値算出手段と、
この補正値算出手段で補正された計数結果から前記測定対象との距離を求める演算手段とを有することを特徴とする距離計。
A semiconductor laser that emits laser light to the object to be measured;
The semiconductor laser is operated so that a first oscillation period including at least a period in which the oscillation wavelength continuously increases monotonically and a second oscillation period including at least a period in which the oscillation wavelength continuously decreases monotonously exist. A laser driver,
A light receiver for converting the optical output of the semiconductor laser into an electrical signal;
Counting means for counting the number of interference waveforms generated by the self-coupling effect between the laser light emitted from the semiconductor laser and the return light from the measurement object, included in the output signal of the light receiver,
A period measuring means for measuring a period of the interference waveform during a counting period for counting the number of the interference waveforms every time the interference waveform is input;
Frequency distribution creating means for creating a frequency distribution of the period of the interference waveform during the counting period from the measurement result of the period measuring means,
And analyzed by FFT the frequency distribution, a frequency analysis means of the periodic frequency distribution to determine the frequency spectrum maximum value appears in the frequency corresponding to the spacing of the maxima of the frequency distribution,
Representative value deriving means for obtaining a representative value of the period distribution of the interference waveform from the maximum value appearing in the frequency spectrum obtained by the frequency analysis means of the periodic frequency distribution;
From the frequency distribution, the sum Ns of the frequencies of the class that is less than or equal to the first predetermined number of times of the representative value and the sum Nw of the frequencies of the class that are greater than or equal to the second predetermined number of times of the representative value are obtained, and Correction value calculating means for correcting the counting result of the counting means based on the frequencies Ns and Nw of
A distance meter comprising: a calculation unit that obtains a distance from the measurement object from a count result corrected by the correction value calculation unit.
請求項5又は6記載の距離計において、
前記代表値導出手段は、前記周期度数分布周波数解析手段が求めた周波数スペクトルにおいて強度が最大の極大値を前記干渉波形の周期の分布の代表値を示すものと見なし、前記代表値を求めることを特徴とする距離計。
The distance meter according to claim 5 or 6,
The representative value deriving means regards the maximum value having the maximum intensity in the frequency spectrum obtained by the frequency analysis means of the periodic frequency distribution as a representative value of the period distribution of the interference waveform, and obtains the representative value. Distance meter characterized by.
請求項5又は6記載の距離計において、
前記代表値導出手段は、前記周期度数分布周波数解析手段が求めた周波数スペクトルにおいて周波数が最も高い極大値を前記干渉波形の周期の分布の代表値を示すものと見なし、前記代表値を求めることを特徴とする距離計。
The distance meter according to claim 5 or 6,
The representative value deriving means regards the maximum value having the highest frequency in the frequency spectrum obtained by the frequency analysis means of the periodic frequency distribution as a representative value of the period distribution of the interference waveform, and obtains the representative value. Distance meter characterized by.
請求項6乃至8のいずれか1項に記載の距離計において、
前記補正値算出手段は、前記計数手段の計数結果をNとしたとき、補正後の計数結果N’を、N’=N+Nw−Nsにより求めることを特徴とする距離計。
The distance meter according to any one of claims 6 to 8,
The correction value calculating means obtains the corrected count result N ′ by N ′ = N + Nw−Ns, where N is the count result of the counting means.
特定の物理量と線形の関係を有し前記物理量が一定の場合は略単一周波数となる信号の数を数える計数方法において、
一定の計数期間における入力信号の数を数える計数手順と、
前記計数期間中の前記入力信号の周期を信号が入力される度に測定する周期測定手順と、
この周期測定手順の測定結果から前記計数期間中の信号周期の度数分布を作成する度数分布作成手順と、
前記度数分布をFFTにより解析して、前記度数分布の極大値の間隔に相当する周波数に極大値が出現する周波数スペクトルを求める周期度数分布周波数解析手順と、
この周期度数分布周波数解析手順で求めた周波数スペクトルに出現した極大値から前記入力信号の周期の分布の代表値を求める代表値導出手順と、
前記度数分布から、前記代表値の第1の所定数倍以下である階級の度数の総和Nsと、前記代表値の第2の所定数倍以上である階級の度数の総和Nwとを求め、これらの度数NsとNwに基づいて前記計数手順の計数結果を補正する補正値算出手順とを備えることを特徴とする計数方法。
In a counting method for counting the number of signals having a linear relationship with a specific physical quantity and having a substantially single frequency when the physical quantity is constant,
A counting procedure for counting the number of input signals in a certain counting period;
A period measurement procedure for measuring the period of the input signal during the counting period each time a signal is input;
A frequency distribution creating procedure for creating a frequency distribution of the signal period during the counting period from the measurement result of the period measuring procedure;
Analyzing the frequency distribution by FFT to obtain a frequency analysis procedure of a periodic frequency distribution to obtain a frequency spectrum in which a maximum value appears at a frequency corresponding to the interval between the maximum values of the frequency distribution;
A representative value derivation procedure for obtaining a representative value of the period distribution of the input signal from the maximum value that appears in the frequency spectrum obtained by the frequency analysis procedure of the periodic frequency distribution;
From the frequency distribution, the sum Ns of the frequencies of the class that is less than or equal to the first predetermined number of times of the representative value and the sum Nw of the frequencies of the class that are greater than or equal to the second predetermined number of times of the representative value are obtained, and And a correction value calculation procedure for correcting the counting result of the counting procedure based on the frequencies Ns and Nw.
請求項10記載の計数方法において、
前記代表値導出手順は、前記周期度数分布周波数解析手順で求めた周波数スペクトルにおいて強度が最大の極大値を前記入力信号の周期の分布の代表値を示すものと見なし、前記代表値を求めることを特徴とする計数方法。
The counting method according to claim 10, wherein
In the representative value deriving procedure, the maximum value having the maximum intensity in the frequency spectrum obtained by the frequency analysis procedure of the periodic frequency distribution is regarded as indicating the representative value of the period distribution of the input signal, and the representative value is obtained. A counting method characterized by the above.
請求項10記載の計数方法において、
前記代表値導出手順は、前記周期度数分布周波数解析手順で求めた周波数スペクトルにおいて周波数が最も高い極大値を前記入力信号の周期の分布の代表値を示すものと見なし、前記代表値を求めることを特徴とする計数方法。
The counting method according to claim 10, wherein
In the representative value derivation procedure, the maximum value having the highest frequency in the frequency spectrum obtained by the frequency analysis procedure of the periodic frequency distribution is regarded as indicating the representative value of the period distribution of the input signal, and the representative value is obtained. A counting method characterized by the above.
請求項10乃至12のいずれか1項に記載の計数方法において、
前記補正値算出手順は、前記計数手順の計数結果をNとしたとき、補正後の計数結果N’を、N’=N+Nw−Nsにより求めることを特徴とする計数方法。
The counting method according to any one of claims 10 to 12,
In the correction value calculating procedure, when the counting result of the counting procedure is N, the corrected counting result N ′ is obtained by N ′ = N + Nw−Ns.
半導体レーザを用いて測定対象にレーザ光を放射する距離計測方法において、
発振波長が連続的に単調増加する期間を少なくとも含む第1の発振期間と発振波長が連続的に単調減少する期間を少なくとも含む第2の発振期間とが交互に存在するように前記半導体レーザを動作させる発振手順と、
前記半導体レーザから放射されたレーザ光と前記測定対象からの戻り光との干渉光を電気信号に変換する受光器の出力信号に含まれる、前記半導体レーザから放射されたレーザ光と前記測定対象からの戻り光とによって生じる干渉波形の数を数える計数手順と、
前記干渉波形の数を数える計数期間中の前記干渉波形の周期を干渉波形が入力される度に測定する周期測定手順と、
この周期測定手順の測定結果から前記計数期間中の干渉波形の周期の度数分布を作成する度数分布作成手順と、
前記度数分布をFFTにより解析して、前記度数分布の極大値の間隔に相当する周波数に極大値が出現する周波数スペクトルを求める周期度数分布周波数解析手順と、
この周期度数分布周波数解析手順で求めた周波数スペクトルに出現した極大値から前記干渉波形の周期の分布の代表値を求める代表値導出手順と、
前記度数分布から、前記代表値の第1の所定数倍以下である階級の度数の総和Nsと、前記代表値の第2の所定数倍以上である階級の度数の総和Nwとを求め、これらの度数NsとNwに基づいて前記計数手順の計数結果を補正する補正値算出手順と、
この補正値算出手順で補正された計数結果から前記測定対象との距離を求める演算手順とを備えることを特徴とする距離計測方法。
In a distance measurement method for emitting laser light to a measurement object using a semiconductor laser,
The semiconductor laser is operated so that a first oscillation period including at least a period in which the oscillation wavelength continuously increases monotonically and a second oscillation period including at least a period in which the oscillation wavelength continuously decreases monotonously exist. Oscillation procedure
The laser light emitted from the semiconductor laser and the measurement object included in the output signal of the light receiver that converts the interference light between the laser light emitted from the semiconductor laser and the return light from the measurement object into an electrical signal. A counting procedure for counting the number of interference waveforms caused by the return light of
A period measurement procedure for measuring a period of the interference waveform during a counting period for counting the number of the interference waveforms every time the interference waveform is input;
A frequency distribution creating procedure for creating a frequency distribution of the period of the interference waveform during the counting period from the measurement result of this cycle measuring procedure,
Analyzing the frequency distribution by FFT to obtain a frequency analysis procedure of a periodic frequency distribution to obtain a frequency spectrum in which a maximum value appears at a frequency corresponding to the interval between the maximum values of the frequency distribution;
A representative value derivation procedure for obtaining a representative value of the period distribution of the interference waveform from a local maximum value appearing in the frequency spectrum obtained by the frequency analysis procedure of the periodic frequency distribution;
From the frequency distribution, the sum Ns of the frequencies of the class that is less than or equal to the first predetermined number of times of the representative value and the sum Nw of the frequencies of the class that are greater than or equal to the second predetermined number of times of the representative value are obtained, and A correction value calculation procedure for correcting the counting result of the counting procedure based on the frequencies Ns and Nw of
A distance measurement method comprising: a calculation procedure for obtaining a distance to the measurement object from a count result corrected by the correction value calculation procedure.
半導体レーザを用いて測定対象にレーザ光を放射する距離計測方法において、
発振波長が連続的に単調増加する期間を少なくとも含む第1の発振期間と発振波長が連続的に単調減少する期間を少なくとも含む第2の発振期間とが交互に存在するように前記半導体レーザを動作させる発振手順と、
前記半導体レーザの光出力を電気信号に変換する受光器の出力信号に含まれる、前記半導体レーザから放射されたレーザ光と前記測定対象からの戻り光との自己結合効果によって生じる干渉波形の数を数える計数手順と、
前記干渉波形の数を数える計数期間中の前記干渉波形の周期を干渉波形が入力される度に測定する周期測定手順と、
この周期測定手順の測定結果から前記計数期間中の干渉波形の周期の度数分布を作成する度数分布作成手順と、
前記度数分布をFFTにより解析して、前記度数分布の極大値の間隔に相当する周波数に極大値が出現する周波数スペクトルを求める周期度数分布周波数解析手順と、
この周期度数分布周波数解析手順で求めた周波数スペクトルに出現した極大値から前記干渉波形の周期の分布の代表値を求める代表値導出手順と、
前記度数分布から、前記代表値の第1の所定数倍以下である階級の度数の総和Nsと、前記代表値の第2の所定数倍以上である階級の度数の総和Nwとを求め、これらの度数NsとNwに基づいて前記計数手順の計数結果を補正する補正値算出手順と、
この補正値算出手順で補正された計数結果から前記測定対象との距離を求める演算手順とを備えることを特徴とする距離計測方法。
In a distance measurement method for emitting laser light to a measurement object using a semiconductor laser,
The semiconductor laser is operated so that a first oscillation period including at least a period in which the oscillation wavelength continuously increases monotonically and a second oscillation period including at least a period in which the oscillation wavelength continuously decreases monotonously exist. Oscillation procedure
The number of interference waveforms generated by the self-coupling effect between the laser light emitted from the semiconductor laser and the return light from the measurement object, included in the output signal of the light receiver that converts the optical output of the semiconductor laser into an electrical signal. Counting procedure;
A period measurement procedure for measuring a period of the interference waveform during a counting period for counting the number of the interference waveforms every time the interference waveform is input;
A frequency distribution creating procedure for creating a frequency distribution of the period of the interference waveform during the counting period from the measurement result of this cycle measuring procedure,
Analyzing the frequency distribution by FFT to obtain a frequency analysis procedure of a periodic frequency distribution to obtain a frequency spectrum in which a maximum value appears at a frequency corresponding to the interval between the maximum values of the frequency distribution;
A representative value derivation procedure for obtaining a representative value of the period distribution of the interference waveform from a local maximum value appearing in the frequency spectrum obtained by the frequency analysis procedure of the periodic frequency distribution;
From the frequency distribution, the sum Ns of the frequencies of the class that is less than or equal to the first predetermined number of times of the representative value and the sum Nw of the frequencies of the class that are greater than or equal to the second predetermined number of times of the representative value are obtained, and A correction value calculation procedure for correcting the counting result of the counting procedure based on the frequencies Ns and Nw of
A distance measurement method comprising: a calculation procedure for obtaining a distance to the measurement object from a count result corrected by the correction value calculation procedure.
請求項14又は15記載の距離計測方法において、
前記代表値導出手順は、前記周期度数分布周波数解析手順で求めた周波数スペクトルにおいて強度が最大の極大値を前記干渉波形の周期の分布の代表値を示すものと見なし、前記代表値を求めることを特徴とする距離計測方法。
The distance measuring method according to claim 14 or 15,
In the representative value deriving procedure, the maximum value having the maximum intensity in the frequency spectrum obtained by the frequency analysis procedure of the periodic frequency distribution is regarded as indicating the representative value of the period distribution of the interference waveform, and the representative value is obtained. A distance measurement method characterized by
請求項14又は15記載の距離計測方法において、
前記代表値導出手順は、前記周期度数分布周波数解析手順で求めた周波数スペクトルにおいて周波数が最も高い極大値を前記干渉波形の周期の分布の代表値を示すものと見なし、前記代表値を求めることを特徴とする距離計測方法。
The distance measuring method according to claim 14 or 15,
In the representative value derivation procedure, the maximum value having the highest frequency in the frequency spectrum obtained by the frequency analysis procedure of the periodic frequency distribution is regarded as a representative value of the period distribution of the interference waveform, and the representative value is obtained. A distance measurement method characterized by
請求項14乃至17のいずれか1項に記載の距離計測方法において、
前記補正値算出手順は、前記計数手順の計数結果をNとしたとき、補正後の計数結果N’を、N’=N+Nw−Nsにより求めることを特徴とする距離計測方法。
The distance measuring method according to any one of claims 14 to 17,
In the correction value calculation procedure, when the counting result of the counting procedure is N, the corrected counting result N ′ is obtained by N ′ = N + Nw−Ns.
JP2007259549A 2007-10-03 2007-10-03 Counting device, distance meter, counting method, and distance measuring method Expired - Fee Related JP5545913B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2007259549A JP5545913B2 (en) 2007-10-03 2007-10-03 Counting device, distance meter, counting method, and distance measuring method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2007259549A JP5545913B2 (en) 2007-10-03 2007-10-03 Counting device, distance meter, counting method, and distance measuring method

Publications (2)

Publication Number Publication Date
JP2009085910A JP2009085910A (en) 2009-04-23
JP5545913B2 true JP5545913B2 (en) 2014-07-09

Family

ID=40659501

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2007259549A Expired - Fee Related JP5545913B2 (en) 2007-10-03 2007-10-03 Counting device, distance meter, counting method, and distance measuring method

Country Status (1)

Country Link
JP (1) JP5545913B2 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5798668B2 (en) * 2014-05-19 2015-10-21 アズビル株式会社 Counting device, physical quantity sensor, counting method and physical quantity measuring method
JP6018662B2 (en) * 2015-03-31 2016-11-02 アズビル株式会社 Counting device, physical quantity sensor, counting method and physical quantity measuring method

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5915411B2 (en) * 1975-11-10 1984-04-09 株式会社山武 Keisu Cairo
JPS55166344A (en) * 1979-06-14 1980-12-25 Toshiba Corp Pulse fall-out correcting system by scanning of pulse counter
JPS62101123A (en) * 1985-10-28 1987-05-11 Mitsubishi Electric Corp Pulse type revolution counter
JP2885807B2 (en) * 1988-10-21 1999-04-26 東京航空計器株式会社 Distance detection device
JPH08274627A (en) * 1995-03-31 1996-10-18 Omron Corp Counter
JPH10186039A (en) * 1996-12-20 1998-07-14 Shimadzu Corp Nuclear medical imaging device
JP2980576B2 (en) * 1997-09-12 1999-11-22 株式会社東芝 Physical random number generating apparatus and method, and physical random number recording medium
JP2006322932A (en) * 2006-04-25 2006-11-30 Yamatake Corp Signal detection circuit
JP5545915B2 (en) * 2007-01-25 2014-07-09 アズビル株式会社 Counting device, distance meter, counting method, and distance measuring method

Also Published As

Publication number Publication date
JP2009085910A (en) 2009-04-23

Similar Documents

Publication Publication Date Title
JP5545915B2 (en) Counting device, distance meter, counting method, and distance measuring method
KR100945209B1 (en) Counting Device, Distance Meter, Counting Method, And Distance Measuring Method
JP5545916B2 (en) Physical quantity sensor and physical quantity measuring method
US7649631B2 (en) Distance/speed meter and distance/speed measuring method
JP2006313080A (en) Distance/speed meter and method for distance/speed measurement
JP2009014701A (en) Distance/speed meter, and distance/speed measuring method
JP5545914B2 (en) Physical quantity sensor and physical quantity measuring method
JP5461762B2 (en) Distance meter and distance measuring method
JP5184785B2 (en) Distance meter and distance measuring method
JP5545913B2 (en) Counting device, distance meter, counting method, and distance measuring method
JP5596915B2 (en) Physical quantity sensor and physical quantity measuring method
JP5461761B2 (en) Distance / speed meter and distance / speed measurement method
JP5530068B2 (en) Distance / speed meter and distance / speed measurement method
JP5421568B2 (en) Physical quantity sensor and physical quantity measuring method
JP5612240B2 (en) Distance meter and distance measuring method
JP2008175547A (en) Range finder and distance measuring method
JP5096123B2 (en) Distance / speed meter and distance / speed measurement method
TWI359943B (en)
WO2011111180A1 (en) Physical quantity sensor and physical quantity measuring method
JP5663130B2 (en) Reflective photoelectric switch and object detection method
EP2251707A1 (en) Reflective photoelectric switch and object detection method
JP5484661B2 (en) Physical quantity sensor and physical quantity measuring method
JP2010160138A (en) Physical quantity sensor and physical quantity measurement method
JP2009063346A (en) Physical quantity sensor and physical quantity measuring method

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20100826

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20130514

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20130716

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20140507

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20140512

R150 Certificate of patent or registration of utility model

Ref document number: 5545913

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees