JP7074311B2 - Optical distance measuring device and measuring method - Google Patents

Optical distance measuring device and measuring method Download PDF

Info

Publication number
JP7074311B2
JP7074311B2 JP2017165940A JP2017165940A JP7074311B2 JP 7074311 B2 JP7074311 B2 JP 7074311B2 JP 2017165940 A JP2017165940 A JP 2017165940A JP 2017165940 A JP2017165940 A JP 2017165940A JP 7074311 B2 JP7074311 B2 JP 7074311B2
Authority
JP
Japan
Prior art keywords
frequency
beat
distance
beat frequency
optical
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.)
Active
Application number
JP2017165940A
Other languages
Japanese (ja)
Other versions
JP2019045200A (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.)
National Institute of Advanced Industrial Science and Technology AIST
Original Assignee
National Institute of Advanced Industrial Science and Technology AIST
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 National Institute of Advanced Industrial Science and Technology AIST filed Critical National Institute of Advanced Industrial Science and Technology AIST
Priority to JP2017165940A priority Critical patent/JP7074311B2/en
Publication of JP2019045200A publication Critical patent/JP2019045200A/en
Application granted granted Critical
Publication of JP7074311B2 publication Critical patent/JP7074311B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Description

本発明は、自動車や自律ロボット等に用いられる環境認識センサ等に適する光学的距離測定装置、及び測定方法に関する。 The present invention relates to an optical distance measuring device suitable for an environment recognition sensor or the like used in an automobile, an autonomous robot, or the like, and a measuring method.

自動車や自律ロボットの環境認識センサや、建設・土木現場における形状計測などへの応用を目的として、ライダ(LiDAR:Laser Imaging Detection and Ranging)の開発が進展している。すでに実用化されているTOF(Time of Flight)方式のライダは、対象物に光パルスを照射して、反射して戻ってくるまでの時間から距離を測定し、照射する光パルスを空間的に走査して、3次元距離データを生成するものである。 The development of lidar (LiDAR: Laser Imaging Detection and Ranking) is progressing for the purpose of application to environment recognition sensors for automobiles and autonomous robots, and shape measurement at construction and civil engineering sites. The TOF (Time of Flyght) type rider, which has already been put into practical use, irradiates an object with an optical pulse, measures the distance from the time until it is reflected and returns, and spatially irradiates the optical pulse. It scans and generates three-dimensional distance data.

TOF方式ライダでは、直接検波により対象物からの反射光を検出する。一方、コヒーレント検波を用いるFMCW(Frequency Modulated Continuous Wave)方式は、より高感度の反射光検出が可能であり、対象物までの距離に加えて、ドップラーシフトから運動速度も測定できる特徴を有している。ミリ波領域のFMCWレーダは、車載用の衝突防止センサとして実用化されている。光波領域でFMCWライダを実現できれば、空間分解能の格段の向上が期待できる。現状のFMCWライダは、装置の構成が複雑であり、高コヒーレンスのレーザ光源が要求されるため、応用分野は限定されている。 The TOF rider detects the reflected light from the object by direct detection. On the other hand, the FMCW (Frequency Modified Continuous Wave) method using coherent detection has the feature that it can detect reflected light with higher sensitivity and can measure the motion speed from the Doppler shift in addition to the distance to the object. There is. The FMCW radar in the millimeter wave region has been put into practical use as an in-vehicle collision prevention sensor. If the FMCW rider can be realized in the light wave region, it can be expected that the spatial resolution will be significantly improved. The current FMCW rider has a complicated device configuration and requires a high coherence laser light source, so that the field of application is limited.

図14(a)(b)(c)は、FMCWライダの原理を説明する図である。図14(a)のFMCWライダ装置は、変調信号発生器1と、注入電流源2と、半導体レーザ3と、ビームスプリッタ5a、5bと、光サーキュレータ6と、反射鏡7と、光検出器12を備える。次に、光源として半導体レーザ3を用いた場合について、FMCWライダの動作原理を説明する。三角波を発生する変調信号発生器1の出力を、注入電流源2に入力し、半導体レーザ3の注入電流を変調する。半導体レーザ3の出力光を2分し、一方をプローブ光9、他方を参照光8とする。参照光とは、プローブ光と位相が同期し、光学的遅延の基準となる光をいう。ビームスプリッタ5a、5bと、光サーキュレータ6と、反射鏡7と、対象物11により、干渉計を構成する。プローブ光9を対象物11に照射し、対象物11からの反射光10と参照光8を合波し、光検出器12に入力してビート信号を検出する。図14(b)に半導体レーザ3の光周波数の波形を示す。半導体レーザ3の光周波数は、変調信号発生器1から出力される三角波に対応して、アップ、ダウンチャープを周期的に繰り返す。参照光8と反射光10との間には時間差が発生し、光検出器12の出力には、時間差に比例した周波数を有するビート信号が発生する。図14(c)に、プローブ光の反射光10と参照光の光周波数、およびビート信号波形を示す。図14(c)の上部に、参照光(実線)と反射光(点線)の光周波数を、ビート信号と対応させて示す。したがって、ビート信号の周波数(以下、ビート周波数という。)を測定すれば、時間差、すなわち対象物11までの距離を求めることができる。アップ、ダウンチャープの両方を使用するのは、対象物11が運動している場合に、ドップラーシフトから運動速度を求めるためである。ここでは半導体レーザを用いた場合について説明したが、周波数変調機能を有するレーザであれば、同様の測定を行うことができる。 14 (a), (b) and 14 (c) are diagrams illustrating the principle of the FMCW rider. The FMCW rider device of FIG. 14A includes a modulation signal generator 1, an injection current source 2, a semiconductor laser 3, beam splitters 5a and 5b, an optical circulator 6, a reflector 7, and a photodetector 12. To prepare for. Next, the operating principle of the FMCW rider will be described with respect to the case where the semiconductor laser 3 is used as the light source. The output of the modulation signal generator 1 that generates a triangular wave is input to the injection current source 2 to modulate the injection current of the semiconductor laser 3. The output light of the semiconductor laser 3 is divided into two, one of which is the probe light 9 and the other of which is the reference light 8. The reference light is light that is in phase with the probe light and serves as a reference for optical delay. A beam splitter 5a and 5b, an optical circulator 6, a reflector 7, and an object 11 constitute an interferometer. The probe light 9 is irradiated to the object 11, the reflected light 10 from the object 11 and the reference light 8 are combined and input to the photodetector 12 to detect the beat signal. FIG. 14B shows the waveform of the optical frequency of the semiconductor laser 3. The optical frequency of the semiconductor laser 3 periodically repeats up and down chirps corresponding to the triangular wave output from the modulation signal generator 1. A time difference occurs between the reference light 8 and the reflected light 10, and a beat signal having a frequency proportional to the time difference is generated at the output of the photodetector 12. FIG. 14C shows the reflected light 10 of the probe light, the optical frequency of the reference light, and the beat signal waveform. The optical frequencies of the reference light (solid line) and the reflected light (dotted line) are shown in the upper part of FIG. 14 (c) in correspondence with the beat signal. Therefore, if the frequency of the beat signal (hereinafter referred to as the beat frequency) is measured, the time difference, that is, the distance to the object 11 can be obtained. Both the up and down chirps are used in order to obtain the motion velocity from the Doppler shift when the object 11 is in motion. Here, the case where a semiconductor laser is used has been described, but the same measurement can be performed if the laser has a frequency modulation function.

FMCWライダにより測定されるビート周波数fFMCWは、次式により表すことができる。 The beat frequency f FMCW measured by the FMCW rider can be expressed by the following equation.

Figure 0007074311000001
Figure 0007074311000001

ここで、Δνはチャープ帯域幅、T=1/fは変調周期、fは変調周波数、Lは対象物11までの距離、cは光速度である。(数1)式において、Δν/(T/2)は単位時間当たりの周波数変化、すなわち、チャープ率を表す。ビート周波数から距離を算出するためには、チャープ率を事前に求めておく必要がある。 Here, Δν is the chirp bandwidth, T m = 1 / fm is the modulation period, fm is the modulation frequency, L is the distance to the object 11, and c is the optical velocity. In the equation (Equation 1), Δν / (T m / 2) represents the frequency change per unit time, that is, the chirp rate. In order to calculate the distance from the beat frequency, it is necessary to obtain the chirp rate in advance.

FMCWライダにおける距離分解能δLは、次式により表すことができる。 The distance resolution δL in the FMCW rider can be expressed by the following equation.

Figure 0007074311000002
Figure 0007074311000002

(数2)式における分解能の意味は、近接する2つの反射点を分離して検出する能力である。反射点が1つの場合は、さらに高い精度で距離を測定することができる。距離分解能δLはチャープ帯域幅Δνに反比例するので、高い分解能を得るためには、大きなチャープ帯域幅が必要である。例えば、分解能10cm、1cmを得るのに必要なチャープ帯域幅は、それぞれ1.5GHz、15GHzである。反射点が1つの場合においても、精度はチャープ帯域幅に反比例する。 The meaning of the resolution in the equation (Equation 2) is the ability to detect two adjacent reflection points separately. When there is one reflection point, the distance can be measured with higher accuracy. Since the distance resolution δL is inversely proportional to the chirp bandwidth Δν, a large chirp bandwidth is required to obtain high resolution. For example, the chirp bandwidths required to obtain a resolution of 10 cm and 1 cm are 1.5 GHz and 15 GHz, respectively. Even with one reflection point, the accuracy is inversely proportional to the chirp bandwidth.

チャープ帯域幅に加えて、チャープの直線性もFMCWライダの性能を決定する重要な特性である。(数1)式のビート周波数は、光周波数が時間に比例して増加(アップチャープ)、または減少(ダウンチャープ)することを前提にしている。光周波数が時間に対して非線形に変化する場合は、一定値であるはずのビート周波数が変化して、距離算出の誤差要因となる。 In addition to the chirp bandwidth, the chirp linearity is also an important characteristic that determines the performance of the FMCW rider. The beat frequency of the equation (Equation 1) is based on the premise that the optical frequency increases (up chirp) or decreases (down chirp) in proportion to time. When the optical frequency changes non-linearly with time, the beat frequency, which should be a constant value, changes, which causes an error in the distance calculation.

アイセーフ波長域で動作する半導体レーザは、注入電流変調により周波数を直接変調できることから、小型で低価格のFMCWライダ用の光源として期待されている。ところが、半導体レーザの周波数変調は熱効果に起因し、周波数応答特性が平坦ではないため、非線形チャープが顕著に現れることが知られている。三角波により周波数変調した場合は、非線形チャープに起因して、変調信号に含まれない周波数成分が現れることが報告されている(非特許文献1参照)。 A semiconductor laser operating in the eye-safe wavelength range is expected as a light source for a compact and low-cost FMCW rider because the frequency can be directly modulated by injection current modulation. However, it is known that the frequency modulation of the semiconductor laser is caused by the thermal effect and the frequency response characteristic is not flat, so that the non-linear chirp appears remarkably. It has been reported that when the frequency is modulated by a triangular wave, a frequency component not included in the modulated signal appears due to the nonlinear chirp (see Non-Patent Document 1).

FMCWライダにおいて、このような非線形チャープの影響を抑圧または低減する手法は、2つに大別できる。一つは、半導体レーザの変調を制御して、所望の線形チャープを得る方法である。もう一つは、検出したビート信号を処理して、非線形チャープの影響を除去する方法である。 In the FMCW rider, the method of suppressing or reducing the influence of such a non-linear chirp can be roughly divided into two. One is to control the modulation of the semiconductor laser to obtain the desired linear chirp. The other is to process the detected beat signal to eliminate the influence of the non-linear chirp.

レーザの周波数変調を光学的に検出して、基準信号との誤差をレーザに負帰還して制御する装置、及び方法が、次のように報告されている(特許文献1乃至3参照)。レーザの周波数変調を検出するためのホモダイン、またはヘテロダイン干渉計を用意し、周波数変調光を入力して、干渉計出力からビート周波数を検出する。ビート周波数は、本来一定値になるはずであるが、非線形チャープの影響により変動する。検出したビート周波数を基準信号源の周波数と比較して誤差信号を生成し、レーザの注入電流を負帰還制御して、非線形チャープを抑圧する。 A device and a method for optically detecting the frequency modulation of a laser and controlling an error with a reference signal by negatively feeding back to the laser have been reported as follows (see Patent Documents 1 to 3). A homodyne or heterodyne interferometer for detecting the frequency modulation of the laser is prepared, frequency-modulated light is input, and the beat frequency is detected from the interferometer output. The beat frequency should be a constant value, but it fluctuates due to the influence of the non-linear chirp. The detected beat frequency is compared with the frequency of the reference signal source to generate an error signal, and the injection current of the laser is controlled by negative feedback to suppress the non-linear chirp.

周波数変調信号、またはレーザ出力をモニタして、変調信号発生器を制御するとともに、検出したビート信号を補正して、距離を算出する方法が次のように報告されている(特許文献4参照)。レーザ出力光の位相を数学的にモデル化し、モニタ結果からモデルに含まれるパラメータを推定し、制御と信号処理を行い、距離を算出する。 A method of monitoring a frequency-modulated signal or a laser output to control a modulated signal generator and correcting a detected beat signal to calculate a distance has been reported as follows (see Patent Document 4). .. The phase of the laser output light is mathematically modeled, the parameters included in the model are estimated from the monitor results, control and signal processing are performed, and the distance is calculated.

ミリ波FMCWレーダにおいて、光学的にミリ波周波数を検出し、信号処理により、非線形チャープの影響を抑圧する装置が、以下のように報告されている(特許文献5参照)。送出するミリ波信号を光信号に変換し、ホモダイン干渉計によりビート信号を検出して、パルス信号に変換する。パルス信号には非線形チャープの情報が含まれており、このパルス信号をクロックとして、対象物からの反射により生成されるビート信号をAD変換することにより、非線形チャープの影響を抑圧することができる。 In the millimeter wave FMCW radar, a device that optically detects the millimeter wave frequency and suppresses the influence of the non-linear chirp by signal processing has been reported as follows (see Patent Document 5). The millimeter wave signal to be transmitted is converted into an optical signal, the beat signal is detected by a homodyne interferometer, and the pulse signal is converted. The pulse signal contains information on the non-linear chirp, and the influence of the non-linear chirp can be suppressed by AD-converting the beat signal generated by the reflection from the object using this pulse signal as a clock.

特開2000-111312号公報Japanese Unexamined Patent Publication No. 2000-11132 米国特許出願公開第2010/0085992号明細書U.S. Patent Application Publication No. 2010/0085992 米国特許出願公開第2012/0106579号明細書U.S. Patent Application Publication No. 2012/0106579 米国特許出願公開第2009/0135403号明細書U.S. Patent Application Publication No. 2009/0135403 特表2008-514910号公報Japanese Patent Publication No. 2008-514910

H.Tsuchida、“Waveform measurement technique for phase/frequency-modulated lights based on heterodyne interferometry、”Optics Express、vol.25,no.5、pp.4793-4799(2017年3月)H. Tsuchida, "Waveform measurement technology for phase / frequency-modulated lights based on heterodyne interferometry," Optics Express. 25, no. 5, pp. 4793-4799 (March 2017)

現実のFMCWライダでは、光周波数変化が時間に対して非線形に変化することが原因で、距離算出の誤差要因となる問題がある。 In an actual FMCW rider, there is a problem that the change in optical frequency changes non-linearly with time, which causes an error in distance calculation.

前述した、非線形チャープの影響を抑圧または低減する方法(特許文献1乃至3参照)では、レーザの負帰還制御により非線形チャープを抑圧するため、誤差信号を実時間で生成する必要がある。このため、ホモダイン、またはヘテロダイン干渉計を備える必要がある。また、変調信号発生器の制御と検出信号の補正を行う方法(特許文献4参照)では、干渉計などの光学装置が必要である。また、特許文献5はミリ波レーダ装置に関する技術であるが、光領域のFMCWライダにも適用できる。しかし、AD変換器のクロックを実時間で生成する必要があるため、ホモダイン干渉計が必要である。 In the above-mentioned method of suppressing or reducing the influence of the non-linear chirp (see Patent Documents 1 to 3), since the non-linear chirp is suppressed by the negative feedback control of the laser, it is necessary to generate an error signal in real time. Therefore, it is necessary to provide a homodyne or heterodyne interferometer. Further, in the method of controlling the modulation signal generator and correcting the detection signal (see Patent Document 4), an optical device such as an interferometer is required. Further, although Patent Document 5 is a technique relating to a millimeter wave radar device, it can also be applied to an FMCW rider in an optical region. However, since it is necessary to generate the clock of the AD converter in real time, a homodyne interferometer is required.

このように、従来技術においては、距離測定用の干渉計とは別に、周波数変調を検出して制御するための装置が必要であり、装置の構成が複雑になる。光学的距離測定装置として、装置構成が複雑化せず、小型で低価格のライダシステムが実現できれば、車載用の衝突防止、歩行者検知センサなど、民生分野への展開が期待できる。 As described above, in the prior art, a device for detecting and controlling frequency modulation is required separately from the interferometer for distance measurement, and the configuration of the device becomes complicated. If a compact and low-cost rider system can be realized as an optical distance measuring device without complicating the device configuration, it can be expected to be applied to the consumer field such as in-vehicle collision prevention and pedestrian detection sensor.

本発明は、FMCWライダにおける上述の問題を解決しようとするものであり、干渉計などの付加的な装置を用いることなく、レーザの非線形チャープの影響を除去して、正確な距離測定を可能とする、光学的距離測定装置、及び測定方法を提供することを目的とする。 The present invention is intended to solve the above-mentioned problems in the FMCW rider, and enables accurate distance measurement by eliminating the influence of the non-linear charp of the laser without using an additional device such as an interferometer. It is an object of the present invention to provide an optical distance measuring device and a measuring method.

本発明は、前記目的を達成するために、以下の特徴を有するものである。 The present invention has the following features in order to achieve the above object.

(1) 周波数変調したレーザと、光検出器と、前記レーザの出力光を2分し、一方を参照光、他方をプローブ光とし、前記プローブ光を対象物に照射して、前記対象物からの反射光と、前記参照光を合波して前記光検出器に入射する干渉計と、前記光検出器から発生するビート信号の周波数を検出するビート周波数検出装置と、前記ビート周波数検出装置により検出したビート周波数の平均値を求め、該平均値を用いて、事前に校正した距離と平均ビート周波数との関係を基に、前記対象物までの距離を算出する演算処理装置とを、備えることを特徴とする光学的距離測定装置。
(2) 前記ビート周波数検出装置は、前記ビート信号の同相成分I(t)と直交成分Q(t)を求めて、逆正接tan-1{Q(t)/I(t)}から前記ビート信号の瞬時位相を計算し、該瞬時位相の時間微分から前記ビート周波数を求めることを特徴とする、(1)記載の光学的距離測定装置。
(3) 前記ビート周波数検出装置は、前記ビート信号をデジタル信号に変換するAD変換器と、該デジタル信号のゼロクロス点を挟む2つのサンプル値から、内挿によりゼロクロス点の時刻を求め、該ゼロクロス点の時刻から前記ビート周波数を求める演算処理部とを備えることを特徴とする、(1)記載の光学的距離測定装置。
(4) 前記レーザの周波数変調波形は三角波であることを特徴とする(1)乃至(3)のいずれか1項記載の光学的距離測定装置。
(5) 周波数変調したレーザの出力光を2分し、一方を参照光、他方をプローブ光とし、前記プローブ光を対象物に照射して、前記対象物からの反射光と、前記参照光を合波して光検出器に入射し、該光検出器から発生するビート信号の周波数を検出し、前記ビート信号の周波数の平均値を求め、該平均値を用いて、事前に校正した距離と平均ビート周波数との関係を基に、前記対象物までの距離を算出することを特徴とする光学的距離測定方法。
(6) 前記ビート信号周波数の検出は、前記ビート信号の同相成分I(t)と直交成分Q(t)を求めて、逆正接tan-1{Q(t)/I(t)}から前記ビート信号の瞬時位相を計算し、該瞬時位相の時間微分から前記ビート周波数を求めることを特徴とする、(5)記載の光学的距離測定方法。
(7) 前記ビート信号周波数の検出は、前記ビート信号をデジタル信号に変換し、該デジタル信号のゼロクロス点を挟む2つのサンプル値から、内挿によりゼロクロス点の時刻を求め、該ゼロクロス点の時刻から前記ビート周波数を求めることを特徴とする(5)記載の光学的距離測定方法。
(8) 前記レーザの周波数変調波形は三角波であることを特徴とする(5)乃至(7)のいずれか1項記載の光学的距離測定方法。
(1) The frequency-modulated laser, the light detector, and the output light of the laser are divided into two, one of which is the reference light and the other of which is the probe light. The reflected light of the above, the interferometer that combines the reference light and is incident on the optical detector, the beat frequency detection device that detects the frequency of the beat signal generated from the optical detector, and the beat frequency detection device. Provided with an arithmetic processing device that obtains the average value of the detected beat frequencies and calculates the distance to the object based on the relationship between the distance calibrated in advance and the average beat frequency using the average value. An optical distance measuring device characterized by.
(2) The beat frequency detection device obtains the in-phase component I (t) and the orthogonal component Q (t) of the beat signal, and obtains the beat from the inverse positive tangent tan -1 {Q (t) / I (t)}. The optical distance measuring apparatus according to (1), wherein the instantaneous phase of a signal is calculated and the beat frequency is obtained from the time derivative of the instantaneous phase.
(3) The beat frequency detection device obtains the time of the zero cross point by intercalation from an AD converter that converts the beat signal into a digital signal and two sample values that sandwich the zero cross point of the digital signal, and the zero cross. The optical distance measuring apparatus according to (1), comprising an arithmetic processing unit for obtaining the beat frequency from the time of a point.
(4) The optical distance measuring device according to any one of (1) to (3), wherein the frequency modulation waveform of the laser is a triangular wave.
(5) The output light of the frequency-modulated laser is divided into two, one of which is the reference light and the other of which is the probe light, and the probe light is irradiated to the object to obtain the reflected light from the object and the reference light. The combined wave is incident on the optical detector, the frequency of the beat signal generated from the optical detector is detected, the average value of the frequencies of the beat signal is obtained, and the average value is used to obtain the distance calibrated in advance. An optical distance measuring method characterized in that the distance to the object is calculated based on the relationship with the average beat frequency.
(6) In the detection of the beat signal frequency, the in-phase component I (t) and the orthogonal component Q (t) of the beat signal are obtained, and the reverse positive tangent tan -1 {Q (t) / I (t)} is used. The optical distance measuring method according to (5), wherein the instantaneous phase of the beat signal is calculated and the beat frequency is obtained from the time derivative of the instantaneous phase.
(7) In the detection of the beat signal frequency, the beat signal is converted into a digital signal, the time of the zero cross point is obtained by interpolation from two sample values sandwiching the zero cross point of the digital signal, and the time of the zero cross point is obtained. The optical distance measuring method according to (5), wherein the beat frequency is obtained from the above.
(8) The optical distance measuring method according to any one of (5) to (7), wherein the frequency modulation waveform of the laser is a triangular wave.

本発明の光学的距離測定装置および測定方法においては、FMCWライダにより生成されるビート信号周波数の平均値を測定し、事前に校正した距離と平均周波数との関係を基にして距離を算出する。このため、レーザの周波数変調をモニタする光学系と、周波数変調を制御する電子回路が不要である。その結果、装置構成を格段に簡素化することができて、小型化と低価格化を実現できる。 In the optical distance measuring device and measuring method of the present invention, the average value of the beat signal frequencies generated by the FMCW rider is measured, and the distance is calculated based on the relationship between the pre-calibrated distance and the average frequency. Therefore, an optical system for monitoring the frequency modulation of the laser and an electronic circuit for controlling the frequency modulation are unnecessary. As a result, the device configuration can be remarkably simplified, and miniaturization and cost reduction can be realized.

また、直接周波数変調により、非線形チャープが顕著に発生する半導体レーザを用いる場合は、信号処理だけで非線形チャープの影響を抑制できるため、より小型化で高精度の装置が実現できる。 Further, when a semiconductor laser in which a non-linear chirp is remarkably generated by direct frequency modulation is used, the influence of the non-linear chirp can be suppressed only by signal processing, so that a smaller and more accurate device can be realized.

本発明に係る光学的距離測定装置を説明する図である。It is a figure explaining the optical distance measuring apparatus which concerns on this invention. 理想的な三角波について、理論計算の結果を表す図である。(a)に光周波数とビート周波数の波形を、(b)に平均ビート周波数と相対誤差を示す。It is a figure which shows the result of the theoretical calculation about an ideal triangular wave. (A) shows the waveforms of the optical frequency and the beat frequency, and (b) shows the average beat frequency and the relative error. 非線形チャープを有する三角波について、理論計算の結果を表す図である。(a)に光周波数とビート周波数の波形を、(b)に平均ビート周波数と相対誤差を示す。It is a figure which shows the result of the theoretical calculation about the triangular wave which has a nonlinear chirp. (A) shows the waveforms of the optical frequency and the beat frequency, and (b) shows the average beat frequency and the relative error. 第1の実施の形態の光学的距離測定装置を構成するビート周波数検出装置を説明する図である。It is a figure explaining the beat frequency detection apparatus which constitutes the optical distance measuring apparatus of 1st Embodiment. 第1の実施の形態の光学的距離測定装置と、距離に対する平均ビート周波数の変化率を校正する装置を説明する図である。It is a figure explaining the optical distance measuring apparatus of 1st Embodiment and the apparatus which calibrates the rate of change of the average beat frequency with respect to a distance. 第1の実施の形態に係る、手動式の可変光遅延線を用いた場合のビート周波数の波形を表す図である。(b)は(a)の一部(点線で囲んだ部分)の拡大図である。It is a figure which shows the waveform of the beat frequency when the manual variable optical delay line is used which concerns on 1st Embodiment. (B) is an enlarged view of a part (a part surrounded by a dotted line) of (a). 第1の実施の形態に係る、手動式の可変光遅延線を用いた場合の距離測定値と、測定誤差を表す図である。It is a figure which shows the distance measurement value and the measurement error at the time of using the manual variable optical delay line which concerns on 1st Embodiment. 第1の実施の形態に係る、光ファイバパッチコードを光遅延線として用いた場合の校正結果を表す図である。(a)にビート周波数の波形、(b)にファイバ長と平均ビート周波数との関係を示す。It is a figure which shows the calibration result when the optical fiber patch cord which concerns on 1st Embodiment is used as an optical delay line. (A) shows the waveform of the beat frequency, and (b) shows the relationship between the fiber length and the average beat frequency. 第2の実施の形態の光学的距離測定装置を構成するビート周波数検出装置を説明する図である。(a)にビート周波数検出装置、(b)に、光検出器から出力されるビート信号の波形、(c)に(b)の一部拡大図を示す。It is a figure explaining the beat frequency detection apparatus which constitutes the optical distance measuring apparatus of 2nd Embodiment. (A) shows a beat frequency detector, (b) shows a waveform of a beat signal output from a photodetector, and (c) shows a partially enlarged view of (b). 第2の実施の形態の光学的距離測定装置と、距離に対する平均ビート周波数の変化率を校正する装置を説明する図である。It is a figure explaining the optical distance measuring apparatus of 2nd Embodiment, and the apparatus which calibrates the rate of change of the average beat frequency with respect to the distance. 第2の実施の形態に係る、手動式の可変光遅延線を用いた場合の波形を表す図である。(a)にビート信号の波形、(b)にビート周波数の波形を示す。It is a figure which shows the waveform at the time of using the manual variable optical delay line which concerns on 2nd Embodiment. (A) shows the waveform of the beat signal, and (b) shows the waveform of the beat frequency. 第2の実施の形態に係る、手動式の可変光遅延線を用いた場合の距離測定値と、測定誤差を表す図である。It is a figure which shows the distance measurement value and the measurement error at the time of using the manual variable optical delay line which concerns on 2nd Embodiment. 第2の実施の形態に係る、光ファイバパッチコードを光遅延線として用いた場合の校正結果を表す図である。(a)にビート周波数の波形、(b)にファイバ長と平均ビート周波数との関係を示す。It is a figure which shows the calibration result when the optical fiber patch cord which concerns on 2nd Embodiment is used as an optical delay line. (A) shows the waveform of the beat frequency, and (b) shows the relationship between the fiber length and the average beat frequency. 従来技術におけるFMCWライダの原理を説明する図である。(a)にFMCWライダ装置、(b)に半導体レーザの光周波数の波形、(c)に、反射光と参照光の光周波数、およびビート信号波形を示す。It is a figure explaining the principle of FMCW rider in the prior art. (A) shows the FMCW rider device, (b) shows the waveform of the optical frequency of the semiconductor laser, and (c) shows the optical frequency of the reflected light and the reference light, and the beat signal waveform.

以下、本発明の実施の形態を詳細に説明する。本発明では、周波数変調したレーザと、干渉計と、光検出器と、ビート周波数検出装置と、演算処理装置とを具備する装置を用いて、対象物までの距離を測定する。 Hereinafter, embodiments of the present invention will be described in detail. In the present invention, the distance to an object is measured by using a device including a frequency-modulated laser, an interferometer, a photodetector, a beat frequency detection device, and an arithmetic processing device.

図1は、本発明の本実施の形態に係る光学的距離測定装置の基本構成を説明する図である。図1の光学的距離測定装置は、変調信号発生器1と、注入電流源2と、半導体レーザ3と、干渉計4と、光検出器12と、ビート周波数検出装置13と、信号処理装置14とを備える。図1では、周波数変調したレーザとして、直接変調半導体レーザを例に示している。干渉計4は、ビームスプリッタ5aと光サーキュレータ6とビームスプリッタ5bと反射鏡7と、対象物11とで、主に構成する。図1に示すように、変調信号発生器1の出力を、注入電流源2を介して、半導体レーザ3に入力し、出力光の周波数を変調する。半導体レーザ3の出力光をビームスプリッタ5aにより2分し、一方を参照光8、他方をプローブ光9とする。プローブ光9を光サーキュレータ6を介して、対象物11に照射する。対象物11からの反射光を、光サーキュレータ6を介して、ビームスプリッタ5bに導き、該反射光10と参照光8とを合波して、光検出器12により受光する。周波数変調された参照光8と反射光10との間には、対象物11までの距離に応じた時間差が存在するため、周波数差が生じる。光検出器12の出力には、周波数差に対応したビート信号が発生する。光検出器12の出力をビート周波数検出装置13に入力して、ビート周波数を検出した後、演算処理装置14に入力する。演算処理装置14には、半導体レーザ3を用いて、事前に校正した距離と平均ビート周波数との関係が記録されている。演算処理装置14において、検出した前記ビート周波数の平均値を計算して、事前に校正した距離と平均ビート周波数の関係を基にして、対象物11までの距離を計算する。 FIG. 1 is a diagram illustrating a basic configuration of an optical distance measuring device according to an embodiment of the present invention. The optical distance measuring device of FIG. 1 includes a modulation signal generator 1, an injection current source 2, a semiconductor laser 3, an interferometer 4, an optical detector 12, a beat frequency detection device 13, and a signal processing device 14. And prepare. In FIG. 1, a directly modulated semiconductor laser is shown as an example of a frequency-modulated laser. The interferometer 4 is mainly composed of a beam splitter 5a, an optical circulator 6, a beam splitter 5b, a reflector 7, and an object 11. As shown in FIG. 1, the output of the modulation signal generator 1 is input to the semiconductor laser 3 via the injection current source 2 to modulate the frequency of the output light. The output light of the semiconductor laser 3 is split into two by a beam splitter 5a, one of which is the reference light 8 and the other of which is the probe light 9. The probe light 9 is irradiated to the object 11 via the optical circulator 6. The reflected light from the object 11 is guided to the beam splitter 5b via the optical circulator 6, and the reflected light 10 and the reference light 8 are combined and received by the photodetector 12. Since there is a time difference between the frequency-modulated reference light 8 and the reflected light 10 according to the distance to the object 11, a frequency difference occurs. A beat signal corresponding to the frequency difference is generated at the output of the photodetector 12. The output of the photodetector 12 is input to the beat frequency detection device 13, the beat frequency is detected, and then the output is input to the arithmetic processing device 14. The arithmetic processing apparatus 14 records the relationship between the pre-calibrated distance and the average beat frequency using the semiconductor laser 3. The arithmetic processing apparatus 14 calculates the average value of the detected beat frequencies, and calculates the distance to the object 11 based on the relationship between the pre-calibrated distance and the average beat frequency.

数式を用いて、以下詳しく説明する。図1において参照光8の電場E(t)は次式により表すことができる。 The following will be described in detail using mathematical formulas. In FIG. 1, the electric field ER (t) of the reference light 8 can be expressed by the following equation.

Figure 0007074311000003
Figure 0007074311000003

一方、反射光10の電場E(t)は次式により表すことができる。 On the other hand, the electric field ES (t) of the reflected light 10 can be expressed by the following equation.

Figure 0007074311000004
Figure 0007074311000004

ここで、AとAは電場の振幅、νは中心周波数、φ(t)は半導体レーザ3の周波数変調に起因する瞬時位相、τは対象物11までの光の往復時間を表す。対象物11までの距離をLとすれば、光の往復時間τは、光速度cを用いて、τ=2L/cにより表すことができる。φ(t)は、半導体レーザ3の光周波数ν(t)と次式により関係づけられる。 Here, AR and AS are the amplitude of the electric field, ν 0 is the center frequency, φ (t) is the instantaneous phase caused by the frequency modulation of the semiconductor laser 3, and τ d is the round-trip time of light to the object 11. .. Assuming that the distance to the object 11 is L, the round-trip time τ d of light can be expressed by τ d = 2 L / c using the speed of light c. φ (t) is related to the optical frequency ν (t) of the semiconductor laser 3 by the following equation.

Figure 0007074311000005
Figure 0007074311000005

FMCWライダにおいては、周波数変調信号として、三角波または鋸波を用いることが一般的であるが、後述するように、距離と平均ビート周波数との関係を事前に校正しておけば、これらの波形に限定されるものではない。 In the FMCW rider, it is common to use a triangular wave or a sawtooth wave as the frequency modulation signal, but as will be described later, if the relationship between the distance and the average beat frequency is calibrated in advance, these waveforms can be obtained. Not limited.

光検出器12から出力されるビート信号V(t)は、次式により表すことができる。 The beat signal V B (t) output from the photodetector 12 can be expressed by the following equation.

Figure 0007074311000006
Figure 0007074311000006

ここで、ηは光検出器12の感度などにより決まる定数である。ビート周波数検出装置13において、(数6)式により表されるビート信号から、次式により表されるビート周波数fB(t)を検出する。 Here, η is a constant determined by the sensitivity of the photodetector 12 and the like. The beat frequency detection device 13 detects the beat frequency f B (t) represented by the following equation from the beat signal represented by the equation (Equation 6).

Figure 0007074311000007
Figure 0007074311000007

次に、演算処理装置14において、(数7)式により表されるビート周波数fB(t)から、次式により表される平均ビート周波数fB,avg(t,t)を計算する。 Next, in the arithmetic processing apparatus 14, the average beat frequencies f B , avg (t 1 , t 2 ) expressed by the following equation are calculated from the beat frequency f B (t) expressed by the equation (Equation 7). ..

Figure 0007074311000008
Figure 0007074311000008

平均ビート周波数fB,avg(t,t)は、ビート周波数fB(t)を有限区間[t,t]に渡り積分した後、有限区間の時間t-tにより除した値である。ビート周波数をN個のサンプル値fB(t)(t=1,2,・・,N)として測定した場合は、平均ビート周波数は、N個のサンプル値の総和をサンプル数Nにより除した値になる。平均化の区間[t,t]は任意であるが、ドップラーシフトから対象物11の運動速度を求める場合は、周波数変調のアップ、ダウンチャープに対応した半周期が望ましい。運動速度を求める必要が無い場合は、変調周期の整数倍にしてもよい。平均化の区間が短いほど、測定時間を短縮できる。 The average beat frequency f B , avg (t 1 , t 2 ) is obtained by integrating the beat frequency f B (t) over a finite interval [t 1 , t 2 ] and then dividing by the time t 2 − t 1 in the finite interval. It is the value that was set. When the beat frequency is measured as N sample values f B ( tk ) ( tk = 1, 2, ..., N), the average beat frequency is the sum of the N sample values based on the number of samples N. It becomes the value divided. The averaging interval [t 1 , t 2 ] is arbitrary, but when the motion velocity of the object 11 is obtained from the Doppler shift, a half cycle corresponding to the up and down chirps of the frequency modulation is desirable. If it is not necessary to obtain the motion velocity, it may be an integral multiple of the modulation period. The shorter the averaging interval, the shorter the measurement time.

演算処理装置14には、事前に校正した距離Lと平均ビート周波数fB,avg(t,t)との関係が記録又は記憶されている。両者の関係は、距離に対する平均周波数のデータを、表の形式で記録することも可能であるが、解析的に表現できる形式で記録しておくことが望ましい。例えば、次式のように、平均ビート周波数fB,avg(t,t)が距離Lに比例する場合を考える。 The arithmetic processing apparatus 14 records or stores the relationship between the pre-calibrated distance L and the average beat frequencies f B , avg (t 1 , t 2 ). Regarding the relationship between the two, it is possible to record the data of the average frequency with respect to the distance in the format of a table, but it is desirable to record it in a format that can be expressed analytically. For example, consider the case where the average beat frequencies f B , avg (t 1 , t 2 ) are proportional to the distance L as shown in the following equation.

Figure 0007074311000009
Figure 0007074311000009

ここで、γは距離に対する平均ビート周波数の変化率を表す。γの値が既知であれば、測定した平均ビート周波数から、次式を用いて距離L’を算出できる。 Here, γ represents the rate of change of the average beat frequency with respect to the distance. If the value of γ is known, the distance L'can be calculated from the measured average beat frequency using the following equation.

Figure 0007074311000010
Figure 0007074311000010

このように、単一のパラメータγのみで距離を算出できるので、演算処理に要する時間を短縮できる。平均ビート周波数を計算する区間[t,t]は、事前に行う校正と同一に設定する必要がある。 In this way, since the distance can be calculated using only a single parameter γ, the time required for the arithmetic processing can be shortened. The interval [t 1 , t 2 ] for calculating the average beat frequency must be set to be the same as the calibration performed in advance.

(数10)式により算出した距離L’と、実際の距離Lとの相対誤差εは、次式により表すことができる。 The relative error ε between the distance L'calculated by the equation (Equation 10) and the actual distance L can be expressed by the following equation.

Figure 0007074311000011
Figure 0007074311000011

以下、平均ビート周波数fB,avg(t,t)と距離Lの関係について、具体的な波形を用いて、理論的に計算した結果を説明する。最初に、光周波数ν(t)として、理想的な三角波を考える。理想的な三角波は、次式のフーリエ級数展開で表すことができる。 Hereinafter, the results of theoretically calculating the relationship between the average beat frequencies f B , avg (t 1 , t 2 ) and the distance L will be described using a specific waveform. First, consider an ideal triangular wave as the optical frequency ν (t). The ideal triangle wave can be expressed by the Fourier series expansion of the following equation.

Figure 0007074311000012
Figure 0007074311000012

ここで、Δνはチャープ帯域幅、fは変調周波数である。(数12)式が示すように、理想的な三角波は、基本波と奇数次の高調波により展開できる。 Here, Δν is the chirp bandwidth and fm is the modulation frequency. As shown by Eq. (Equation 12), an ideal triangular wave can be developed by a fundamental wave and odd-order harmonics.

図2は、理想的な三角波について、理論的に計算した結果を表す図である。図2(a)は、(数12)式の光周波数(破線)と、(数12)式を(数7)式に代入して計算したビート周波数(実線)を表す。変調周波数はf=5kHz、チャープ帯域幅はΔν=15GHzである。Lは対象物までの距離を表し、ビート周波数は、距離5、50、500、5000mに対する計算結果である。 FIG. 2 is a diagram showing the results of theoretical calculations for an ideal triangular wave. FIG. 2A shows the optical frequency (broken line) of the equation (Equation 12) and the beat frequency (solid line) calculated by substituting the equation (Equation 12) into the equation (Equation 7). The modulation frequency is fm = 5 kHz, and the chirp bandwidth is Δν = 15 GHz. L represents the distance to the object, and the beat frequency is the calculation result for the distances of 5, 50, 500, and 5000 m.

図2(a)のビート周波数が示すように、三角波の線形チャープ領域では、ビート周波数は一定値になる。一方、三角波の頂点では、アップチャープとダウンチャープが入れ替わるため、ビート周波数がスパイク状に変化する。FMCWライダにおいては、ビート周波数が一定の領域を抽出してスペクトル解析を行い、距離を算出する。距離の増大とともに、ビート周波数が一定の領域が減少し、スペクトル解析に利用するデータが減少する。 As shown by the beat frequency in FIG. 2A, the beat frequency becomes a constant value in the linear chirp region of the triangular wave. On the other hand, at the apex of the triangular wave, the up chirp and the down chirp are switched, so that the beat frequency changes like a spike. In the FMCW rider, a region having a constant beat frequency is extracted, spectrum analysis is performed, and a distance is calculated. As the distance increases, the region where the beat frequency is constant decreases, and the data used for spectrum analysis decreases.

図2(b)は、理想的な三角波について、(数8)式を用いて計算した平均ビート周波数(実線)と、(数11)式を用いて計算した相対誤差(一点鎖線)を距離の関数としてプロットしたものである。点線は(数1)式により表されるビート周波数fFMCWである。変調周波数はf=5kHz、チャープ帯域幅はΔν=15GHz、(数8)式における積分範囲は[0、200μs]である。 FIG. 2B shows the distance between the average beat frequency (solid line) calculated using Eq. (Equation 8) and the relative error (single point chain line) calculated using Eq. (Equation 11) for an ideal triangular wave. It is plotted as a function. The dotted line is the beat frequency f FMCW expressed by the equation (Equation 1). The modulation frequency is fm = 5 kHz, the chirp bandwidth is Δν = 15 GHz, and the integration range in the equation (Equation 8) is [ 0,200 μs].

図2(b)の平均ビート周波数は、距離に比例して増加し、(数1)式のビート周波数fFMCWとほぼ一致するが、距離1km以上の領域では、飽和する傾向にある。これはビート周波数fFMCWが、図2(a)のビート周波数が一定の領域のみを考慮しているのに対して、平均ビート周波数は、スパイク状の変化を含めて計算しているためである。 The average beat frequency in FIG. 2B increases in proportion to the distance and is almost the same as the beat frequency f FMCW of the equation (Equation 1), but tends to be saturated in the region of the distance of 1 km or more. This is because the beat frequency f FMCW considers only the region where the beat frequency in FIG. 2 (a) is constant, whereas the average beat frequency is calculated including spike-like changes. ..

距離1km以下の領域では、図2(b)の平均ビート周波数は、次式により距離と関係づけることができる。 In the region where the distance is 1 km or less, the average beat frequency in FIG. 2B can be related to the distance by the following equation.

Figure 0007074311000013
Figure 0007074311000013

(数1)式のビート周波数fFMCWは、距離に対して1.000MHz/mの割合で増加するが、図2(b)の平均ビート周波数は、距離に対して0.9978MHz/mの割合で変化する。このように変化の割合がわずかに小さくなるのは、スパイク状のビート周波数変化を含めているためである。 The beat frequency f FMCW of the equation (Equation 1) increases at a rate of 1.000 MHz / m with respect to the distance, while the average beat frequency of FIG. 2 (b) increases at a rate of 0.9978 MHz / m with respect to the distance. It changes with. The reason why the rate of change is slightly small is that the spike-like beat frequency change is included.

図2(b)の相対誤差は、距離とともに増大する傾向にあり、距離1km以上では、測定精度が急激に劣化する。実用的な測定精度(5%程度)や、レーザのコヒーレンスを考慮すると、周波数変化率γから距離を算出する方式は、距離1km以下の測定に適用できると考えられる。 The relative error in FIG. 2B tends to increase with a distance, and the measurement accuracy deteriorates sharply at a distance of 1 km or more. Considering practical measurement accuracy (about 5%) and laser coherence, it is considered that the method of calculating the distance from the frequency change rate γ can be applied to the measurement of the distance of 1 km or less.

次に、光周波数ν(t)として、非線形チャープを有する三角波を考える。非線形チャープを有する三角波は、次式のフーリエ級数展開で表すことができる。 Next, consider a triangular wave having a nonlinear chirp as the optical frequency ν (t). A triangular wave with a non-linear chirp can be represented by the Fourier series expansion of the following equation.

Figure 0007074311000014
Figure 0007074311000014

理想的な三角波のフーリエ級数展開は、基本波と奇数次高調波からなるが、非線形チャープを有する場合は、(数14)式が示すように、非線形に起因する偶数次高調波が現れるとともに、フーリエ係数の振幅C、Dと、位相ψ、θ、の値が変化する。 The Fourier series expansion of an ideal triangular wave consists of a fundamental wave and odd-order harmonics, but when it has a non-linear charp, as shown by Eq. (Equation 14), even-order harmonics due to non-linearity appear, and at the same time, The values of the amplitudes C k and D k of the Fourier coefficient and the phases ψ k and θ k change.

図3は、非線形チャープを有する三角波について、理論的に計算した結果を表す図である。図3(a)は、(数14)式の光周波数(破線)と、(数14)式を(数7)式に代入して計算したビート周波数(実線)を表す。変調周波数はf=5kHz、チャープ帯域幅はΔν=15GHzである。Lは対象物までの距離を表し、ビート周波数は、距離5、50、500、5000mに対する計算結果である。 FIG. 3 is a diagram showing the result of theoretical calculation for a triangular wave having a non-linear chirp. FIG. 3A shows the optical frequency (broken line) of the equation (Equation 14) and the beat frequency (solid line) calculated by substituting the equation (Equation 14) into the equation (Equation 7). The modulation frequency is fm = 5 kHz, and the chirp bandwidth is Δν = 15 GHz. L represents the distance to the object, and the beat frequency is the calculation result for the distances of 5, 50, 500, and 5000 m.

図3(a)のビート周波数が示すように、三角波の頂点付近では、ビート周波数がスパイク状に変化する。一方、それ以外の領域では、非線形チャープの影響を受けて、ビート周波数は緩やかに変化する。また、(数14)式における偶数次高調波の存在により、アップチャープ領域よりも、ダウンチャープ領域の周波数変化が大きくなる。非線形チャープの影響により、ビート周波数から距離を正確に算出することは困難である。例えば、対象物11までの距離が5mの場合、(数1)式のビート周波数fFMCWが5.00MHzであるのに対して、ビート周波数は1.05MHzから6.77MHzの範囲で変化する。この周波数変化は、距離1.05mから6.77mに相当する。 As shown by the beat frequency in FIG. 3A, the beat frequency changes like a spike near the apex of the triangular wave. On the other hand, in the other regions, the beat frequency changes slowly under the influence of the nonlinear chirp. Further, due to the existence of the even-order harmonics in the equation (Equation 14), the frequency change in the down chirp region becomes larger than that in the up chirp region. Due to the influence of non-linear chirps, it is difficult to accurately calculate the distance from the beat frequency. For example, when the distance to the object 11 is 5 m, the beat frequency f FMCW of the equation (Equation 1) is 5.00 MHz, whereas the beat frequency changes in the range of 1.05 MHz to 6.77 MHz. This frequency change corresponds to a distance of 1.05 m to 6.77 m.

図3(b)は、非線形チャープを有する三角波について、(数8)式を用いて計算した平均ビート周波数(実線)と、(数11)式を用いて計算した相対誤差(一点鎖線)を距離の関数としてプロットしたものである。点線は(数1)式により表されるビート周波数fFMCWである。変調周波数はf=5kHz、チャープ帯域幅はΔν=15GHz、(数8)式における積分範囲は[0、200μs]である。 FIG. 3B shows the distance between the average beat frequency (solid line) calculated using Eq. (Equation 8) and the relative error (single point chain line) calculated using Eq. (Equation 11) for a triangular wave having a non-linear chirp. It is plotted as a function of. The dotted line is the beat frequency f FMCW expressed by the equation (Equation 1). The modulation frequency is fm = 5 kHz, the chirp bandwidth is Δν = 15 GHz, and the integration range in the equation (Equation 8) is [ 0,200 μs].

理想的な三角波の場合と同様に、図3(b)の平均ビート周波数は、距離に比例して増加し、(数1)式のビート周波数fFMCWとほぼ一致するが、距離1km以上の領域では、飽和する傾向にある。 As in the case of the ideal triangular wave, the average beat frequency in FIG. 3 (b) increases in proportion to the distance and is almost the same as the beat frequency f FMCW in the equation (Equation 1), but in a region of 1 km or more in the distance. Then, it tends to be saturated.

距離1km以下の領域では、図3(b)の平均ビート周波数は、次式により距離と関係づけることができる。 In the region where the distance is 1 km or less, the average beat frequency in FIG. 3B can be related to the distance by the following equation.

Figure 0007074311000015
Figure 0007074311000015

(数1)式のビート周波数fFMCWは、距離に対して1.000MHz/mの割合で増加する。一方、図3(b)の平均ビート周波数は、距離に対して1.0018MHz/mの割合で変化する。スパイク状のビート周波数変化を含めているのに関わらず、変化の割合がわずかに大きくなるのは、非線形チャープによるビート周波数変化が大きいためである。 The beat frequency f FMCW of the equation (Equation 1) increases at a rate of 1.000 MHz / m with respect to the distance. On the other hand, the average beat frequency in FIG. 3B changes at a rate of 1.0018 MHz / m with respect to the distance. The reason why the rate of change is slightly large even though the spike-shaped beat frequency change is included is that the beat frequency change due to the non-linear chirp is large.

図3(b)の相対誤差は、距離とともに増大する傾向にあり、距離1km以上では測定精度が急激に劣化する。距離1km以下の領域における相対誤差は、理想的な三角波の場合の1/2程度である。これは、スパイク状の周波数変化による平均ビート周波数の減少を、非線形チャープに起因する増加により、打ち消しているためである。非線形チャープを有する三角波を用いた場合でも、周波数変化率γから距離を算出する方式は、1km以下の距離測定に適用できると考えられる。 The relative error in FIG. 3B tends to increase with a distance, and the measurement accuracy deteriorates sharply at a distance of 1 km or more. The relative error in the region with a distance of 1 km or less is about ½ of that in the case of an ideal triangular wave. This is because the decrease in the average beat frequency due to the spike-like frequency change is canceled by the increase due to the nonlinear chirp. Even when a triangular wave having a non-linear chirp is used, it is considered that the method of calculating the distance from the frequency change rate γ can be applied to the distance measurement of 1 km or less.

従来技術であるスペクトル解析による距離算出では、ビート周波数が一定値の領域の信号を抽出するので、距離に応じて利用するデータ領域を判断する必要がある。このため、距離の増大とともにデータ量が減少し、周波数分解能が低下する問題がある。一方、平均ビート周波数から距離を算出する場合は、距離に関わらず、すべてデータを利用するので、利用可能なデータ領域を判断する必要がなく、周波数分解能低下の問題も無い。 In the distance calculation by spectrum analysis, which is a conventional technique, a signal in a region where the beat frequency is a constant value is extracted, so it is necessary to determine the data region to be used according to the distance. Therefore, there is a problem that the amount of data decreases as the distance increases and the frequency resolution decreases. On the other hand, when the distance is calculated from the average beat frequency, all the data is used regardless of the distance, so that it is not necessary to determine the available data area and there is no problem of deterioration of the frequency resolution.

図2と図3の例では、三角波を変調信号として用いた場合を説明したが、事前に校正した距離と平均ビート周波数の関係から、距離を一意的に算出できるのであれば、他の波形の変調信号を用いた場合も、同様に実施できる。 In the examples of FIGS. 2 and 3, a case where a triangular wave is used as a modulation signal has been described, but if the distance can be uniquely calculated from the relationship between the pre-calibrated distance and the average beat frequency, other waveforms can be used. The same can be performed when a modulated signal is used.

(第1の実施の形態)
第1の実施の形態の光学的距離測定装置は、ビート周波数検出装置として、ビート信号の瞬時位相を計算し、該瞬時位相の時間微分から前記ビート周波数を求める装置を用いた場合であり、図4乃至8を参照して、以下説明する。
(First Embodiment)
The optical distance measuring device of the first embodiment is a case where a device for calculating the instantaneous phase of a beat signal and obtaining the beat frequency from the time derivative of the instantaneous phase is used as the beat frequency detecting device. The following will be described with reference to 4 to 8.

図4は、ビート周波数検出装置13において行う処理を表す図である。ビート周波数検出装置13に入力したビート信号を2分し、ミキサ17aとミキサ17bにそれぞれ入力する。ミキサ17aの局発信号として、周波数がfの局部発振器15の出力を、ミキサ17bの局発信号として、局部発振器15の出力を位相シフタ16によりπ/2の位相シフトを与えた後、それぞれ入力する。ミキサ17aの出力と、ミキサ17bの出力を、それぞれローパスフィルタ18aとローパスフィルタ18bに入力し、高周波成分を除去して、ビート信号と局発信号の差周波成分を出力する。ローパスフィルタ18aの出力には、次式で表されるビート信号の同相成分19(I(t))が出力される。 FIG. 4 is a diagram showing a process performed by the beat frequency detection device 13. The beat signal input to the beat frequency detection device 13 is divided into two and input to the mixer 17a and the mixer 17b, respectively. The output of the local oscillator 15 having a frequency of f C as the local oscillator signal of the mixer 17a and the output of the local oscillator 15 as the local oscillator signal of the mixer 17b are given a phase shift of π / 2 by the phase shifter 16, respectively. input. The output of the mixer 17a and the output of the mixer 17b are input to the low-pass filter 18a and the low-pass filter 18b, respectively, the high-frequency component is removed, and the difference frequency component between the beat signal and the local signal is output. The common mode component 19 (I (t)) of the beat signal represented by the following equation is output to the output of the low-pass filter 18a.

Figure 0007074311000016
Figure 0007074311000016

一方、ローパスフィルタ18bの出力には、次式で表されるビート信号の直交成分20(Q(t))が出力される。 On the other hand, the orthogonal component 20 (Q (t)) of the beat signal represented by the following equation is output to the output of the low-pass filter 18b.

Figure 0007074311000017
Figure 0007074311000017

(数16)式と(数17)式におけるΦ(t)は次式で表される。Φ(t)はビート信号の瞬時位相である。 Φ (t) in the equations (Equation 16) and (Equation 17) is expressed by the following equation. Φ (t) is the instantaneous phase of the beat signal.

Figure 0007074311000018
Figure 0007074311000018

演算処理部21において、同相成分I(t)と直交成分Q(t)とからビート周波数f(t)を求める。(数19)(数20)式を用いて具体的に説明する。 In the arithmetic processing unit 21, the beat frequency f B (t) is obtained from the common mode component I (t) and the orthogonal component Q (t). (Equation 19) and (Equation 20) will be specifically described.

まず、演算処理部21において、上述のようにして求めた同相成分I(t)と、直交成分Q(t)と、次式を用いて、(数17)式のΦ(t)を計算する。 First, in the arithmetic processing unit 21, the common mode component I (t) obtained as described above, the orthogonal component Q (t), and the following equation are used to calculate Φ (t) of the equation (Equation 17). ..

Figure 0007074311000019
Figure 0007074311000019

注入電流を変調した半導体レーザでは、周波数だけでなく、光強度も同時に変調される。ビート信号には、光強度変調に起因する成分も含まれるが、(数19)式の処理により除去できる。 In a semiconductor laser in which the injection current is modulated, not only the frequency but also the light intensity is modulated at the same time. The beat signal also includes a component caused by light intensity modulation, which can be removed by the process of equation (Equation 19).

次に、演算処理部21において、±πを越える位相変化を検出するため、(数19)式により求めた位相に対して、アンラッピング処理を行う。アンラッピング処理とは、不連続な位相変化に対して、2πの整数倍の位相を付加して、連続的な位相変化にする処理をいう。アンラッピング処理した位相Φ(t)を時間微分して、局部発振器15の周波数fを加えることにより、次式で表されるビート周波数f(t)を求めることができる。 Next, in order to detect a phase change exceeding ± π, the arithmetic processing unit 21 performs unwrapping processing on the phase obtained by the equation (Equation 19). The unwrapping process is a process of adding a phase that is an integral multiple of 2π to a discontinuous phase change to make a continuous phase change. By time-differentiating the unwrapped phase Φ (t) and adding the frequency f C of the local oscillator 15, the beat frequency f B (t) represented by the following equation can be obtained.

Figure 0007074311000020
Figure 0007074311000020

図5は、第1の実施の形態に係る光学的距離測定装置と、距離と平均ビート周波数の関係を校正する装置を説明する図である。図5の光学的距離測定装置は、変調信号発生器1と、注入電流源2と、半導体レーザ3と、干渉計に対応する光学系と、バランス型光検出器24と、ベクトル信号解析装置26と、演算処理装置14とを備える。図1の干渉計4に対応する光学系は、光ファイバを用いて構成し、光方向性結合器22a、22bと、偏波コントローラ23と、光サーキュレータ6とを備える。変調信号発生器1の出力を、注入電流源2を介して、半導体レーザ3に入力し、出力光の周波数を変調する。半導体レーザ3の出力光を、光方向性結合器22aにより二分し、一方を参照光として、偏波コントローラ23を介して、光方向性結合器22bに導く。もう一方はプローブ光として、光サーキュレータ6を介して、光遅延線25に入力する。光遅延線25の出力端面でフレネル反射を受けた光を、光サーキュレータ6を介して光方向性結合器22bに導き、参照光と合波した後、バランス型光検出器24により受光して、ビート信号として出力する。バランス型光検出器24から出力されるビート信号を、ベクトル信号解析装置26に入力して、デジタル信号への変換(分解能12ビット)と、ビート周波数f(t)の計算を実行する。ここでは、ビート信号をデジタル信号に変換した後、ソフトウェア処理により、ビート周波数f(t)の検出を行うが、同様の処理はアナログ電子回路を用いて行うことも可能である。(数8)式に基づく平均ビート周波数の算出は、演算処理装置14を用いて、オフライン処理により実行する。 FIG. 5 is a diagram illustrating an optical distance measuring device according to the first embodiment and a device for calibrating the relationship between the distance and the average beat frequency. The optical distance measuring device of FIG. 5 includes a modulation signal generator 1, an injection current source 2, a semiconductor laser 3, an optical system corresponding to an interferometer, a balanced optical detector 24, and a vector signal analyzer 26. And the arithmetic processing device 14. The optical system corresponding to the interferometer 4 in FIG. 1 is configured by using an optical fiber, and includes optical directional couplers 22a and 22b, a polarization controller 23, and an optical circulator 6. The output of the modulation signal generator 1 is input to the semiconductor laser 3 via the injection current source 2 to modulate the frequency of the output light. The output light of the semiconductor laser 3 is divided into two by the optical directional coupler 22a, and one of them is guided to the optical directional coupler 22b via the polarization controller 23 as reference light. The other is input as probe light to the optical delay line 25 via the optical circulator 6. The light subjected to Fresnel reflection at the output end face of the optical delay line 25 is guided to the optical directional coupler 22b via the optical circulator 6, combined with the reference light, and then received by the balanced photodetector 24. Output as a beat signal. The beat signal output from the balanced photodetector 24 is input to the vector signal analyzer 26 to convert it into a digital signal (resolution 12 bits) and calculate the beat frequency f B (t). Here, after the beat signal is converted into a digital signal, the beat frequency f B (t) is detected by software processing, but the same processing can also be performed using an analog electronic circuit. The calculation of the average beat frequency based on the equation (Equation 8) is executed by offline processing using the arithmetic processing device 14.

具体的な実施例1を示す。変調信号発生器1の出力として、周波数5kHzの三角波を用い、半導体レーザ3の変調電流振幅は60mApp、出力光のチャープ帯域幅は15GHzである。半導体レーザ3の出力光は非線形チャープを有し、光周波数は(数14)式のフーリエ級数により表すことができる。チャープ帯域幅15GHzは、(数2)式の分解能10cmに対応するが、図5に示した光学的距離測定装置においては、反射点1つであるため、さらに高い分解能で距離を測定することができる。 Specific Example 1 is shown. As the output of the modulation signal generator 1, a triangular wave having a frequency of 5 kHz is used, the modulation current amplitude of the semiconductor laser 3 is 60 mA pp , and the chirp bandwidth of the output light is 15 GHz. The output light of the semiconductor laser 3 has a non-linear chirp, and the optical frequency can be expressed by the Fourier series of the equation (Equation 14). The chap bandwidth of 15 GHz corresponds to the resolution of equation (Equation 2) of 10 cm, but in the optical distance measuring device shown in FIG. 5, since there is only one reflection point, the distance can be measured with a higher resolution. can.

校正の第1の例を示す。最初に、図5の光遅延線25として、光路長を精密に設定できる手動式の可変光遅延線を用いて、周波数変化率の校正を行った。可変光遅延線の最小目盛は0.08mmである。校正においては、光路長を8.0mm単位で変化させて、ビート周波数を測定し、変調の1周期に亘る平均値を算出した。図6(a)は、可変光遅延線の設定が0、及び80.0mmの場合のビート周波数を表す図である。図6(a)のスケールでは、両者はほぼ重なっている。図6(b)は図6(a)の一部(点線で囲んだ部分)の拡大図であり、拡大すると、0mm(実線)と80mm(点線、グレー表示)でわずかな差を生じていることがわかる。 A first example of calibration is shown. First, as the optical delay line 25 in FIG. 5, the frequency change rate was calibrated using a manual variable optical delay line capable of precisely setting the optical path length. The minimum scale of the variable optical delay line is 0.08 mm. In the calibration, the optical path length was changed in units of 8.0 mm, the beat frequency was measured, and the average value over one cycle of modulation was calculated. FIG. 6A is a diagram showing a beat frequency when the variable optical delay line is set to 0 and 80.0 mm. On the scale of FIG. 6 (a), the two are almost overlapped. FIG. 6B is an enlarged view of a part of FIG. 6A (the part surrounded by the dotted line), and when enlarged, there is a slight difference between 0 mm (solid line) and 80 mm (dotted line, grayed out). You can see that.

図7は、可変光遅延線の設定値と、平均ビート周波数の関係を表す図である。黒丸は、平均ビート周波数、白丸は、平均ビート周波数から(数10)式を用いて算出した距離と可変光遅延線の設定値との誤差を表す。可変光遅延線25の遅延が0の場合の平均ビート周波数5.60MHzは、光ファイバ干渉計自体が有する光路差に対応する。平均ビート周波数は、遅延の設定値に比例して増加している。距離に対する平均ビート周波数の変化率は、1.023MHz/mであり、理論値とほぼ一致する。算出した距離の誤差は±1mm以内であり、平均ビート周波数から距離を正確に測定できることを示している。 FIG. 7 is a diagram showing the relationship between the set value of the variable optical delay line and the average beat frequency. The black circles represent the average beat frequency, and the white circles represent the error between the distance calculated from the average beat frequency using the equation (Equation 10) and the set value of the variable optical delay line. The average beat frequency of 5.60 MHz when the delay of the variable optical delay line 25 is 0 corresponds to the optical path difference of the optical fiber interferometer itself. The average beat frequency increases in proportion to the delay setting. The rate of change of the average beat frequency with respect to the distance is 1.023 MHz / m, which is almost the same as the theoretical value. The calculated distance error is within ± 1 mm, indicating that the distance can be measured accurately from the average beat frequency.

次に、校正の第2の例を示す。光遅延線25として、光ファイバパッチコードを利用して、光路長を0から10mまで1m単位で変化させて、周波数変化率の校正を行った。長さ1、2、3、5、10mのパッチコードを組み合わせて使用した。パッチコードは、仕様よりも最大で+5cm程度長く製作されている。図8(a)は、光ファイバ長を変化させた場合のビート周波数の変化を表す図である。光ファイバ長に対応して、ビート周波数が増加していることがわかる。 Next, a second example of calibration is shown. Using an optical fiber patch cord as the optical delay line 25, the optical path length was changed from 0 to 10 m in 1 m units, and the frequency change rate was calibrated. Patch cords with lengths of 1, 2, 3, 5, and 10 m were used in combination. The patch cord is manufactured up to +5 cm longer than the specifications. FIG. 8A is a diagram showing a change in the beat frequency when the optical fiber length is changed. It can be seen that the beat frequency increases corresponding to the optical fiber length.

図8(b)は、光ファイバ長と平均ビート周波数の関係を表す図である。水平軸のファイバ長は、パッチコードの仕様であり、実際のファイバ長とは異なる。ファイバ長が0の場合の周波数2.36MHzは、光ファイバ干渉計自体が有する光路差に対応する。光ファイバの屈折率(1.467)を考慮して計算した距離に対する平均ビート周波数の変化率は、1.0095MHz/mであり、理論値と可変光遅延線で校正した値よりもわずかに大きい。その原因は、光ファイバパッチコードが、仕様よりも長く製作されているためである。図8(b)の結果は、光ファイバ長と平均ビート周波数との比例関係を示しており、1m単位の光路長変化に対しても、平均ビート周波数から距離算出が可能であることを示している。 FIG. 8B is a diagram showing the relationship between the optical fiber length and the average beat frequency. The fiber length on the horizontal axis is a patch cord specification and differs from the actual fiber length. The frequency of 2.36 MHz when the fiber length is 0 corresponds to the optical path difference of the optical fiber interferometer itself. The rate of change of the average beat frequency with respect to the distance calculated in consideration of the refractive index (1.467) of the optical fiber is 1.0095 MHz / m, which is slightly larger than the theoretical value and the value calibrated by the variable optical delay line. .. The reason is that the optical fiber patch cord is manufactured longer than the specification. The result of FIG. 8B shows the proportional relationship between the optical fiber length and the average beat frequency, and shows that the distance can be calculated from the average beat frequency even for a change in the optical path length in units of 1 m. There is.

(第2の実施の形態)
第2の実施の形態の光学的距離測定装置は、ビート周波数検出装置として、ビート信号のゼロクロス点の時刻からビート周波数を求める装置を用いた場合であり、図9乃至13を参照して、以下説明する。
(Second embodiment)
The optical distance measuring device of the second embodiment is a case where a device for obtaining the beat frequency from the time of the zero cross point of the beat signal is used as the beat frequency detecting device, and the following is referred to with reference to FIGS. 9 to 13. explain.

図9は、本実施の形態の装置で用いるビート周波数検出装置及びその動作を説明する図である。図9(a)はビート周波数検出装置の回路の模式図である。図9(a)において、ビート周波数検出装置13は、AD変換器27と演算処理部21により構成され、これらには変調信号発生器1の出力から変調信号も入力される。AD変換器27に、図1の光検出器12から出力されるビート信号が入力される。図9(b)は、図1の光検出器12から出力されるビート信号の波形を表す図である。直流成分を除去したビート信号は、複数のゼロクロス点t、t、・・、t、tN+1、・・を有する。時刻にtおけるビート周波数は近似的に次式により表すことができる。 FIG. 9 is a diagram illustrating a beat frequency detection device used in the device of the present embodiment and its operation. FIG. 9A is a schematic diagram of the circuit of the beat frequency detection device. In FIG. 9A, the beat frequency detection device 13 is composed of an AD converter 27 and an arithmetic processing unit 21, and a modulation signal is also input to these from the output of the modulation signal generator 1. The beat signal output from the photodetector 12 of FIG. 1 is input to the AD converter 27. FIG. 9B is a diagram showing a waveform of a beat signal output from the photodetector 12 of FIG. The beat signal from which the DC component has been removed has a plurality of zero cross points t 1 , t 2 , ..., T N , t N + 1 , .... The beat frequency at t N at the time can be approximately expressed by the following equation.

Figure 0007074311000021
Figure 0007074311000021

したがって、ビート信号に含まれるゼロクロス点の時刻を検出すれば、ビート周波数を求めることができる。AD変換器27のサンプリングは、変調信号発生器1の出力と同期しているが、光検出器12から出力されるビート信号とは同期していない。このため、AD変換器27のサンプリング点は、必ずしもゼロクロス点とは一致しない。 Therefore, the beat frequency can be obtained by detecting the time of the zero cross point included in the beat signal. The sampling of the AD converter 27 is synchronized with the output of the modulation signal generator 1, but is not synchronized with the beat signal output from the photodetector 12. Therefore, the sampling point of the AD converter 27 does not always coincide with the zero cross point.

図9(c)は、図9(b)の一部(点線で囲んだ部分)の拡大図であり、AD変換器27のサンプル値から、ゼロクロス点を求める方法を説明する図である。あるゼロクロス点の時刻をtにより表し、時刻tに近接する2つのサンプル点の時刻と電圧を(tNa、VNa)と、(tNb、VNb)により表す(tNb>Na)。2つのサンプリング点の間にゼロクロス点が存在する場合は、電圧の積VNa・VNbは負になる。したがって、すべてのサンプル値に対して、隣り合う2点の電圧の積を計算することにより、ゼロクロス点の位置を特定することができる。ゼロクロス点の時刻をtは、次式によりサンプル値を内挿して求めることができる。 FIG. 9 (c) is an enlarged view of a part (a portion surrounded by a dotted line) of FIG. 9 (b), and is a diagram for explaining a method of obtaining a zero cross point from the sample value of the AD converter 27. The time of a zero cross point is represented by t N , and the time and voltage of two sample points close to the time t N are represented by (t Na , V Na ) and (t Nb , VN b ) (t Nb> t Na ). ). If there is a zero cross point between the two sampling points, the voltage product V Na · V Nb becomes negative. Therefore, the position of the zero cross point can be specified by calculating the product of the voltages of two adjacent points for all the sample values. The time of the zero cross point can be obtained by interpolating the sample value by the following equation for t N.

Figure 0007074311000022
Figure 0007074311000022

図9(a)の演算処理部21においては、AD変換器27のサンプル値からゼロクロス点の時刻を求め、(数21)式を用いてビート周波数f(t)を計算する。図9(a)の演算処理部21は、図1の演算処理装置14と装置を兼用することもできる。 In the arithmetic processing unit 21 of FIG. 9A, the time of the zero cross point is obtained from the sample value of the AD converter 27, and the beat frequency f B (t N ) is calculated using the equation (Equation 21). The arithmetic processing unit 21 of FIG. 9A can also be used as the arithmetic processing apparatus 14 of FIG.

図10は、本実施の形態に係る光学的距離測定装置と、距離と平均ビート周波数の関係を校正する装置を説明する図である。変調信号発生器1から、バランス型光検出器24に至る装置は、図5の装置と同一である。バランス型光検出器24から出力されるビート信号を、デジタルオシロスコープ28(分解能8ビット)に入力して、デジタル信号に変換する。(数22)式によるゼロクロス点の検出と、(数21)式によるビート周波数f(t)の計算と、(数8)式による平均ビート周波数fB,avg(t,t)の計算は、演算処理装置14を用いて、オフライン処理により実行する。図10の演算処理装置14は、図9(a)の演算処理部21の機能も兼ねている。 FIG. 10 is a diagram illustrating an optical distance measuring device according to the present embodiment and a device for calibrating the relationship between the distance and the average beat frequency. The apparatus from the modulation signal generator 1 to the balanced photodetector 24 is the same as the apparatus of FIG. The beat signal output from the balanced photodetector 24 is input to the digital oscilloscope 28 (resolution 8 bits) and converted into a digital signal. Detection of zero cross point by Eq. (Equation 22), calculation of beat frequency f B (t N ) by Eq. (Equation 21), and average beat frequency f B , avg (t 1 , t 2 ) by Eq. (Equation 8) Is calculated by offline processing using the arithmetic processing device 14. The arithmetic processing apparatus 14 of FIG. 10 also has a function of the arithmetic processing unit 21 of FIG. 9A.

具体的な実施例2を示す。変調信号発生器1の出力として、周波数5kHzの三角波を用い、半導体レーザ3の変調電流振幅は60mApp、出力光のチャープ帯域幅は15GHzである。半導体レーザ3の出力光は非線形チャープを有し、光周波数は(数14)式のフーリエ級数により表すことができる。 Specific Example 2 is shown. As the output of the modulation signal generator 1, a triangular wave having a frequency of 5 kHz is used, the modulation current amplitude of the semiconductor laser 3 is 60 mA pp , and the chirp bandwidth of the output light is 15 GHz. The output light of the semiconductor laser 3 has a non-linear chirp, and the optical frequency can be expressed by the Fourier series of the equation (Equation 14).

校正の第3の例を示す。最初に、図10の光遅延線25として、光路長を精密に設定できる手動式の可変光遅延線を用いて、周波数変化率の校正を行った。図11(a)は、可変光遅延線の設定が0mmの場合のビート信号の波形を表す図である。変調周期の前半では、半導体レーザの注入電流が増加するため、ビート信号の包絡線が増加し、ダウンチャープが発生する。変調周期の後半では、注入電流が減少するため、前半とは逆の過程になる。図11(b)は、図11(a)の波形から求めたビート周波数の波形を表す図である。第1の実施の形態により求めた図6(a)と同様の波形であるが、信号対雑音比は小さくなっている。 A third example of calibration is shown. First, as the optical delay line 25 in FIG. 10, the frequency change rate was calibrated using a manual variable optical delay line capable of precisely setting the optical path length. FIG. 11A is a diagram showing a waveform of a beat signal when the variable optical delay line is set to 0 mm. In the first half of the modulation cycle, the injection current of the semiconductor laser increases, so that the envelope of the beat signal increases and down chirp occurs. In the latter half of the modulation cycle, the injection current decreases, so the process is the reverse of the first half. 11 (b) is a diagram showing a waveform of a beat frequency obtained from the waveform of FIG. 11 (a). The waveform is the same as that of FIG. 6A obtained by the first embodiment, but the signal-to-noise ratio is smaller.

図12は、可変光遅延線の設定値と、平均ビート周波数の関係を表す図である。黒丸は、平均ビート周波数、白丸は、平均ビート周波数から(数10)式を用いて算出した距離と可変光遅延線の設定値との誤差を表す。可変光遅延線25の遅延が0の場合の平均ビート周波数5.87MHzは、光ファイバ干渉計自体が有する光路差に対応する。平均ビート周波数は、遅延の設定値に比例して増加している。距離に対する平均ビート周波数の変化率は、1.0323MHz/mであり、理論値と第1の実施の形態により求めた値よりもわずかに大きい。また、算出した距離の誤差は、最大で6.4mmであり、第1の実施の形態よりも大きい。第1の実施の形態と比較して、平均ビート周波数の値と、誤差が大きい原因は、AD変換器として用いたデジタルオシロスコープ28の分解能(8ビット)である。 FIG. 12 is a diagram showing the relationship between the set value of the variable optical delay line and the average beat frequency. The black circles represent the average beat frequency, and the white circles represent the error between the distance calculated from the average beat frequency using the equation (Equation 10) and the set value of the variable optical delay line. The average beat frequency of 5.87 MHz when the delay of the variable optical delay line 25 is 0 corresponds to the optical path difference of the optical fiber interferometer itself. The average beat frequency increases in proportion to the delay setting. The rate of change of the average beat frequency with respect to the distance is 1.0323 MHz / m, which is slightly larger than the theoretical value and the value obtained by the first embodiment. Further, the error of the calculated distance is 6.4 mm at the maximum, which is larger than that of the first embodiment. Compared with the first embodiment, the value of the average beat frequency and the cause of the large error are the resolution (8 bits) of the digital oscilloscope 28 used as the AD converter.

次に、校正の第4の例を示す。光遅延線25として、光ファイバパッチコードを利用して、光路長を0から10mまで1m単位で変化させて、周波数変化率の校正を行った。図13(a)は、光ファイバ長0、2m、4m、6m、8m、10mに対するビート周波数の波形を表す図である。光ファイバ長に対応して、ビート周波数が増加していることがわかる。 Next, a fourth example of calibration is shown. Using an optical fiber patch cord as the optical delay line 25, the optical path length was changed from 0 to 10 m in 1 m units, and the frequency change rate was calibrated. FIG. 13A is a diagram showing waveforms of beat frequencies for optical fiber lengths of 0, 2 m, 4 m, 6 m, 8 m, and 10 m. It can be seen that the beat frequency increases corresponding to the optical fiber length.

図13(b)は、光ファイバ長と平均ビート周波数の関係を表す図である。ファイバ長が0の場合の周波数2.43MHzは、光ファイバ干渉計自体が有する光路差に対応する。光ファイバの屈折率(1.467)を考慮して計算した距離に対する平均ビート周波数の変化率は1.0651MHz/mであり、理論値と第1の実施の形態により求めた値よりもわずかに大きい。原因は、デジタルオシロスコープ28の分解能(8ビット)である。よって、分解能の最適化により変化率は改善できる。図13(b)は、光ファイバ長と平均ビート周波数との比例関係を示しており、1m単位の光路長変化に対しても、平均ビート周波数から距離算出が可能であることを示している。 FIG. 13B is a diagram showing the relationship between the optical fiber length and the average beat frequency. The frequency 2.43 MHz when the fiber length is 0 corresponds to the optical path difference of the optical fiber interferometer itself. The rate of change of the average beat frequency with respect to the distance calculated in consideration of the refractive index (1.467) of the optical fiber is 1.0651 MHz / m, which is slightly larger than the theoretical value and the value obtained by the first embodiment. big. The cause is the resolution (8 bits) of the digital oscilloscope 28. Therefore, the rate of change can be improved by optimizing the resolution. FIG. 13B shows the proportional relationship between the optical fiber length and the average beat frequency, and shows that the distance can be calculated from the average beat frequency even for a change in the optical path length in units of 1 m.

以上、第1と第2の実施の形態においては、変調信号として、三角波を用いた場合を説明した。前述したように、距離と平均ビート周波数の関係が既知であれば、他の変調波形を用いることができる。 In the first and second embodiments, the case where a triangular wave is used as the modulation signal has been described above. As described above, other modulated waveforms can be used as long as the relationship between the distance and the average beat frequency is known.

また、第1と第2の実施の形態の校正は、自律ロボットに搭載する環境認識センサを想定して、最大測定距離10m、分解能1cmに対応する変調条件で行ったものである。車載用センサとして用いる場合は、最大測定距離200m、分解能10cmの仕様が要求されるが、チャープ帯域や変調周波数を最適化して、同様に実施することが可能である。 Further, the calibration of the first and second embodiments is performed under modulation conditions corresponding to a maximum measurement distance of 10 m and a resolution of 1 cm, assuming an environment recognition sensor mounted on an autonomous robot. When used as an in-vehicle sensor, specifications with a maximum measurement distance of 200 m and a resolution of 10 cm are required, but the chirp band and modulation frequency can be optimized and implemented in the same manner.

また、第1と第2の実施の形態においては、光源として半導体レーザを用いた場合について説明したが、周波数変調機能を有するレーザであれば、同様にして実施できる。 Further, in the first and second embodiments, the case where the semiconductor laser is used as the light source has been described, but any laser having a frequency modulation function can be used in the same manner.

上記実施の形態等で示した例は、発明を理解しやすくするために記載したものであり、この形態に限定されるものではない。 The examples shown in the above embodiments and the like are described for the sake of easy understanding of the invention, and are not limited to this embodiment.

本発明の光学的距離測定装置および方法は、レーザの周波数変調を制御するための付加的な装置が不要であるので、小型でかつ高精度で低価格のFMCWライダシステムとして産業上有用である。自動車、自律ロボットなどの環境認識センサとして利用を含め、民生機器への利用可能性を大とするものである。 Since the optical distance measuring device and method of the present invention do not require an additional device for controlling the frequency modulation of the laser, they are industrially useful as a small, highly accurate and low-cost FMCW rider system. The possibility of using it for consumer devices, including its use as an environment recognition sensor for automobiles and autonomous robots, is great.

1 変調信号発生器
2 注入電流源
3 半導体レーザ
4 干渉計
5a、5b ビームスプリッタ
6 光サーキュレータ
7 反射鏡
8 参照光
9 プローブ光
10 反射光
11 対象物
12 光検出器
13 ビート周波数検出装置
14 演算処理装置
15 局部発振器
16 位相シフタ
17a、17b ミキサ
18a、18b ローパスフィルタ
19 同相成分
20 直交成分
21 演算処理部
22a、22b 光方向性結合器
23 偏波コントローラ
24 バランス型光検出器
25 光遅延線
26 ベクトル信号解析装置
27 AD変換器
28 デジタルオシロスコープ
1 Modulation signal generator 2 Injection current source 3 Semiconductor laser 4 Interferometer 5a, 5b Beam splitter 6 Optical circulator 7 Reflector 8 Reference light 9 Probe light 10 Reflected light 11 Object 12 Optical detector 13 Beat frequency detector 14 Arithmetic processing Device 15 Local oscillator 16 Phase shifter 17a, 17b Mixer 18a, 18b Low-pass filter 19 In-phase component 20 Orthogonal component 21 Arithmetic processing unit 22a, 22b Optical directional coupler 23 Polarization controller 24 Balanced optical detector 25 Optical delay line 26 Vector Signal analyzer 27 AD converter 28 Digital oscilloscope

Claims (4)

周波数変調したレーザと、
光検出器と、
前記レーザの出力光を2分し、一方を参照光、他方をプローブ光とし、前記プローブ光を対象物に照射して、前記対象物からの反射光と、前記参照光を合波して前記光検出器に入射する干渉計と、
前記光検出器から発生するビート信号の周波数であって、時間と共に連続的に変化するビート周波数を検出するビート周波数検出装置と、
前記ビート周波数検出装置により検出した前記ビート周波数を、変調周期の整数倍の区間を含む任意の有限区間にわたって積分した後、該有限区間の時間により除すことにより平均値を求め、該平均値を、事前に計算した前記平均値の、距離に対する変化率を示す、事前に校正した距離と前記平均値との関係に照らして、前記対象物までの距離を算出する演算処理装置とを、備え
前記ビート周波数検出装置は、前記ビート信号の同相成分I(t)と直交成分Q(t)を求めて、逆正接tan -1 {Q(t)/I(t)}から前記ビート信号の瞬時位相を計算し、該瞬時位相の時間微分から前記ビート周波数を求めることを特徴とする、光学的距離測定装置。
With a frequency-modulated laser,
With a photodetector
The output light of the laser is divided into two parts, one of which is the reference light and the other of which is the probe light. The interferometer incident on the optical detector and
A beat frequency detection device that detects a beat frequency that is the frequency of a beat signal generated from the photodetector and that changes continuously with time.
The beat frequency detected by the beat frequency detection device is integrated over an arbitrary finite interval including an interval that is an integral multiple of the modulation period , and then divided by the time of the finite interval to obtain an average value, and the average value is calculated. , An arithmetic processing device that calculates the distance to the object in light of the relationship between the pre-calibrated distance and the average value, which indicates the rate of change of the pre-calculated average value with respect to the distance. , Prepare ,
The beat frequency detection device obtains the in-phase component I (t) and the orthogonal component Q (t) of the beat signal, and instantaneously obtains the beat signal from the inverse positive tangent tan -1 {Q (t) / I (t)}. An optical distance measuring device , characterized in that the phase is calculated and the beat frequency is obtained from the time derivative of the instantaneous phase .
前記レーザの周波数変調波形は三角波であることを特徴とする請求項記載の光学的距離測定装置。 The optical distance measuring device according to claim 1 , wherein the frequency modulation waveform of the laser is a triangular wave. 周波数変調したレーザの出力光を2分し、一方を参照光、他方をプローブ光とし、
前記プローブ光を対象物に照射して、前記対象物からの反射光と、前記参照光を合波して光検出器に入射し、
該光検出器から発生するビート信号の周波数であって、時間と共に連続的に変化するビート周波数を検出し、
検出した前記ビート周波数を、変調周期の整数倍の区間を含む任意の有限区間にわたって積分した後、該有限区間の時間により除すことにより平均値を求め、該平均値を、事前に計算した前記平均値の、距離に対する変化率を示す、事前に校正した距離と前記平均値との関係に照らして、前記対象物までの距離を算出し、
前記ビート信号の周波数の検出は、前記ビート信号の同相成分I(t)と直交成分Q(t)を求めて、逆正接tan -1 {Q(t)/I(t)}から前記ビート信号の瞬時位相を計算し、該瞬時位相の時間微分から前記ビート周波数を求めることを特徴とする、光学的距離測定方法。
The output light of the frequency-modulated laser is divided into two, one is the reference light and the other is the probe light.
By irradiating the object with the probe light, the reflected light from the object and the reference light are combined and incident on the photodetector.
The frequency of the beat signal generated from the photodetector, which is the frequency of the beat signal that continuously changes with time, is detected.
The detected beat frequency was integrated over an arbitrary finite interval including an interval of an integral multiple of the modulation period , and then divided by the time of the finite interval to obtain an average value, and the average value was calculated in advance. The distance to the object is calculated in light of the relationship between the pre-calibrated distance and the average value, which indicates the rate of change of the average value with respect to the distance.
To detect the frequency of the beat signal, the in-phase component I (t) and the orthogonal component Q (t) of the beat signal are obtained, and the beat signal is obtained from the inverse positive tangent tan -1 {Q (t) / I (t)}. An optical distance measuring method, characterized in that the instantaneous phase of the above is calculated and the beat frequency is obtained from the time derivative of the instantaneous phase .
前記レーザの周波数変調波形は三角波であることを特徴とする請求項記載の光学的距離測定方法。 The optical distance measuring method according to claim 3 , wherein the frequency modulation waveform of the laser is a triangular wave.
JP2017165940A 2017-08-30 2017-08-30 Optical distance measuring device and measuring method Active JP7074311B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2017165940A JP7074311B2 (en) 2017-08-30 2017-08-30 Optical distance measuring device and measuring method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2017165940A JP7074311B2 (en) 2017-08-30 2017-08-30 Optical distance measuring device and measuring method

Publications (2)

Publication Number Publication Date
JP2019045200A JP2019045200A (en) 2019-03-22
JP7074311B2 true JP7074311B2 (en) 2022-05-24

Family

ID=65814167

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2017165940A Active JP7074311B2 (en) 2017-08-30 2017-08-30 Optical distance measuring device and measuring method

Country Status (1)

Country Link
JP (1) JP7074311B2 (en)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11353558B2 (en) 2018-12-29 2022-06-07 Gm Cruise Holdings Llc Multiple laser, single resonator lidar
US11079480B2 (en) * 2018-12-29 2021-08-03 Gm Cruise Holdings Llc FMCW lidar with wavelength diversity
JP7107268B2 (en) * 2019-04-03 2022-07-27 日本製鉄株式会社 Velocity measuring method and velocity measuring device
AU2020251989A1 (en) * 2019-03-29 2021-10-21 Aurora Operations, Inc. Switchable coherent pixel array for frequency modulated continuous wave light detection and ranging
WO2020261391A1 (en) * 2019-06-25 2020-12-30 三菱電機株式会社 Signal processing device, signal processing method, and radar device
JP7247058B2 (en) * 2019-08-26 2023-03-28 株式会社ミツトヨ Measuring device and method
US20230052690A1 (en) * 2019-12-25 2023-02-16 National Institute Of Advanced Industrial Science And Technology Optical measurement device and measurement method
KR102275387B1 (en) * 2020-02-25 2021-07-09 한화시스템 주식회사 Apparatus and method for detecting small target
US11650318B2 (en) * 2020-04-06 2023-05-16 Gm Cruise Holdings Llc Characterizing linearity of an optical frequency chirp output by an FMCW laser
US11428785B2 (en) * 2020-06-12 2022-08-30 Ours Technology, Llc Lidar pixel with active polarization control
JPWO2022209367A1 (en) 2021-03-30 2022-10-06
CN116981961A (en) 2021-03-30 2023-10-31 松下知识产权经营株式会社 Device and method for measuring distance and/or speed of object
JPWO2022219911A1 (en) 2021-04-15 2022-10-20

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002196067A (en) 1996-06-03 2002-07-10 Mitsubishi Electric Corp Frequency modulation data determining device
JP2003028753A (en) 2001-05-08 2003-01-29 Agilent Technol Inc System and method for measuring group delay based on zero-crossing
JP2003106944A (en) 2001-08-22 2003-04-09 Agilent Technol Inc System and method for measuring optical part characteristic by use of interferometer
JP2005010130A (en) 2003-06-18 2005-01-13 Kumataka Engineering:Kk Measuring apparatus for distance or the like
JP2008531993A (en) 2005-02-17 2008-08-14 メトリス ユーエスエー インク. Compact fiber optic geometry for inverse chirped FMCW coherent laser radar
JP2014202716A (en) 2013-04-09 2014-10-27 株式会社日立ハイテクノロジーズ Distance measuring device
WO2016132520A1 (en) 2015-02-19 2016-08-25 三菱電機株式会社 Fm-cw radar and fm-cw signal generation method
JP2016526150A (en) 2013-04-24 2016-09-01 ローデ ウント シュワルツ ゲゼルシャフト ミット ベシュレンクテル ハフツング ウント コンパニー コマンディット ゲゼルシャフトRohde & Schwarz GmbH & Co.KG Method and apparatus for measuring broadband measurement signals

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS62251686A (en) * 1986-04-24 1987-11-02 Omron Tateisi Electronics Co Laser range finder
JP2733990B2 (en) * 1988-09-06 1998-03-30 日本電気株式会社 Distance measuring device
JP3583906B2 (en) * 1996-07-15 2004-11-04 伊藤 弘昌 Optical rangefinder
US7280186B1 (en) * 2005-02-10 2007-10-09 Hrl Laboratories, Llc High resolution chirped/AM optical FM laser radar

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002196067A (en) 1996-06-03 2002-07-10 Mitsubishi Electric Corp Frequency modulation data determining device
JP2003028753A (en) 2001-05-08 2003-01-29 Agilent Technol Inc System and method for measuring group delay based on zero-crossing
JP2003106944A (en) 2001-08-22 2003-04-09 Agilent Technol Inc System and method for measuring optical part characteristic by use of interferometer
JP2005010130A (en) 2003-06-18 2005-01-13 Kumataka Engineering:Kk Measuring apparatus for distance or the like
JP2008531993A (en) 2005-02-17 2008-08-14 メトリス ユーエスエー インク. Compact fiber optic geometry for inverse chirped FMCW coherent laser radar
JP2014202716A (en) 2013-04-09 2014-10-27 株式会社日立ハイテクノロジーズ Distance measuring device
JP2016526150A (en) 2013-04-24 2016-09-01 ローデ ウント シュワルツ ゲゼルシャフト ミット ベシュレンクテル ハフツング ウント コンパニー コマンディット ゲゼルシャフトRohde & Schwarz GmbH & Co.KG Method and apparatus for measuring broadband measurement signals
WO2016132520A1 (en) 2015-02-19 2016-08-25 三菱電機株式会社 Fm-cw radar and fm-cw signal generation method

Also Published As

Publication number Publication date
JP2019045200A (en) 2019-03-22

Similar Documents

Publication Publication Date Title
JP7074311B2 (en) Optical distance measuring device and measuring method
JP6806347B2 (en) Optical distance measuring device and measuring method
JP7169642B2 (en) Optical measuring device and measuring method
JP7239975B2 (en) Optical angle modulation measuring device and measuring method
EP3540401B1 (en) High resolution interferometric optical frequency domain reflectometry (ofdr)
Norgia et al. Self-mixing instrument for simultaneous distance and speed measurement
US10739448B2 (en) Laser phase estimation and correction
JP7291385B2 (en) Optical measuring device and measuring method
JP2014202716A (en) Distance measuring device
CA2599128C (en) Phase noise compensation for interferometric absolute distance measuring
CN109031340B (en) Continuous frequency modulation laser radar device for measuring object movement speed
JP2019020143A (en) Optical fiber vibration detection sensor and method therefor
JP7426123B2 (en) Optical measuring device and measuring method
Li et al. High-resolution frequency-modulated continuous-wave LiDAR using multiple laser sources simultaneously scanning
CN109612590B (en) Ultrafast optical wavelength measuring system
EP4336209A1 (en) Detection device, radar, and terminal
US9857159B1 (en) Velocity compensated frequency sweeping interferometer and method of using same
JP7061364B2 (en) Distance measuring device and distance measuring method
JP7380382B2 (en) range finder
US11815399B2 (en) Spectrum measurement method and spectrum measurement device
JP2023135347A (en) Measuring apparatus
JP2024007750A (en) Ranging device and ranging method
JP2001235397A (en) Method for measuring group delay time

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20200522

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20210304

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20210324

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210420

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20210929

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20211029

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20220330

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20220407

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: 20220426

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220502

R150 Certificate of patent or registration of utility model

Ref document number: 7074311

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150