CN102192707A - Physical quantity transducer and physical quantity measuring method - Google Patents
Physical quantity transducer and physical quantity measuring method Download PDFInfo
- Publication number
- CN102192707A CN102192707A CN2010101314864A CN201010131486A CN102192707A CN 102192707 A CN102192707 A CN 102192707A CN 2010101314864 A CN2010101314864 A CN 2010101314864A CN 201010131486 A CN201010131486 A CN 201010131486A CN 102192707 A CN102192707 A CN 102192707A
- Authority
- CN
- China
- Prior art keywords
- period
- unit
- cycle
- oscillation
- distance
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims description 19
- 230000010355 oscillation Effects 0.000 claims abstract description 125
- 238000004364 calculation method Methods 0.000 claims abstract description 92
- 238000006073 displacement reaction Methods 0.000 claims abstract description 92
- 239000004065 semiconductor Substances 0.000 claims abstract description 86
- 238000005259 measurement Methods 0.000 claims abstract description 65
- 238000000605 extraction Methods 0.000 claims abstract description 43
- 230000000694 effects Effects 0.000 claims abstract description 8
- 238000012937 correction Methods 0.000 claims description 30
- 238000005070 sampling Methods 0.000 claims description 19
- 238000001514 detection method Methods 0.000 claims description 17
- 230000007423 decrease Effects 0.000 claims description 9
- 238000000691 measurement method Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 abstract 3
- 238000010586 diagram Methods 0.000 description 28
- 230000014509 gene expression Effects 0.000 description 15
- 238000006243 chemical reaction Methods 0.000 description 12
- 238000012545 processing Methods 0.000 description 9
- 238000009826 distribution Methods 0.000 description 8
- 238000003860 storage Methods 0.000 description 8
- 230000004304 visual acuity Effects 0.000 description 5
- 230000003287 optical effect Effects 0.000 description 4
- 230000002123 temporal effect Effects 0.000 description 4
- 238000002347 injection Methods 0.000 description 3
- 239000007924 injection Substances 0.000 description 3
- 230000003321 amplification Effects 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 230000000630 rising effect Effects 0.000 description 2
- 230000003213 activating effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000035559 beat frequency Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
- 230000010356 wave oscillation Effects 0.000 description 1
Images
Landscapes
- Optical Radar Systems And Details Thereof (AREA)
- Instruments For Measurement Of Length By Optical Means (AREA)
Abstract
Description
技术领域technical field
本发明涉及物理量传感器和物理量计测方法,该物理量传感器根据半导体激光器放射的激光与物体的返回光的自混合效应产生的干涉的信息,计测物体的位移、速度。The present invention relates to a physical quantity sensor and a method for measuring physical quantity. The physical quantity sensor measures the displacement and velocity of an object based on interference information generated by the self-mixing effect of laser light emitted from a semiconductor laser and returned light from the object.
背景技术Background technique
FMCW(Frequency Modulated Continuous Wave)雷达和驻波雷达、自混合型激光传感器等的采用干涉原理进行位移(速度)测定的方法中,在基于拍频和干涉条纹的频率计算测定对象的位移和速度时,一般采用FFT(Fast Fourier Transform)等的信号处理或干涉条纹的计数处理。但是存在着这样的问题,为了通过FFT实现高分辨能力,需要很长的取样时间和较高的取样周期的数据,因此需要非常多的处理时间。又,在干涉条纹的计数处理中,为了测定未满半波长的位移,需要使传感器进行物理振动,或对干涉条纹的振幅进行解析,这样则存在着只能计测作为测定对象的周期运动的振动的问题以及、干涉条纹的计数处理耗费时间的问题。FMCW (Frequency Modulated Continuous Wave) radar, standing wave radar, self-mixing laser sensor and other methods of measuring displacement (velocity) using the principle of interference, when calculating the displacement and velocity of the measurement object based on the beat frequency and the frequency of interference fringes Generally, signal processing such as FFT (Fast Fourier Transform) or counting processing of interference fringes is adopted. However, there is a problem that in order to realize high-resolution capability by FFT, data with a long sampling time and a relatively high sampling period is required, and thus a very large processing time is required. In addition, in the counting process of interference fringes, in order to measure the displacement of less than half the wavelength, it is necessary to physically vibrate the sensor or analyze the amplitude of the interference fringes, so there is a possibility that only the periodic motion that is the measurement object can be measured. The problem of vibration and the problem that the counting process of interference fringes takes time.
一方面,发明者提出了利用半导体激光的自混合效应的波长调制型的激光计测器(参照专利文献1)。该激光计测器的结构如图20所示。图20的激光计测器具有:对物体210放射激光的半导体激光器201;将半导体激光器201的光输出变换为电信号的光电二极管202;透镜203,该透镜对来自半导体激光器201的光进行集光并将该光照射到物体210,同时对物体210的返回光进行激光使其入射到半导体激光器201;使得半导体激光器201在振荡波长连续增加的第一振动期间和振动波长连续减少的第二振荡期间交替反复的激光驱动器204;将光电二极管202的输出电流变换为电压并放大的电流-电压变换放大器205;将电流-电压变换放大器205的输出电压二次微分的信号提取电路206;对信号提取电路206的输出电压中包含的MHP的个数进行计数的计数电路207;计算与物体210的距离和物体210的速度的计算装置208;显示计算装置208的计算结果的显示装置209。On the other hand, the inventors proposed a wavelength-modulated laser measurement device utilizing the self-mixing effect of semiconductor lasers (see Patent Document 1). The structure of this laser measuring device is shown in FIG. 20 . The laser measuring device in FIG. 20 has: a semiconductor laser 201 that radiates laser light to an object 210; a photodiode 202 that converts the light output of the semiconductor laser 201 into an electrical signal; and a lens 203 that collects light from the semiconductor laser 201. And the light is irradiated to the object 210, and the return light of the object 210 is lased to make it incident on the semiconductor laser 201; so that the semiconductor laser 201 is in the first oscillation period in which the oscillation wavelength is continuously increased and the second oscillation period in which the oscillation wavelength is continuously reduced Alternately repeated laser driver 204; the output current of the photodiode 202 is converted into a voltage and amplified current-voltage conversion amplifier 205; the output voltage of the current-voltage conversion amplifier 205 is secondarily differentiated. The signal extraction circuit 206; the signal extraction circuit The counting circuit 207 for counting the number of MHP contained in the output voltage of 206; the calculation device 208 for calculating the distance to the object 210 and the speed of the object 210; and the display device 209 for displaying the calculation result of the calculation device 208.
激光驱动器204将随着时间以与规定的变化率反复增减的三角波驱动电流作为注入电流提供给半导体激光器201。由此,对半导体激光器201进行驱动,使得振荡波长以规定的变化率连续增加的第一振动期间和振荡波长以规定的变化率连续减少的第二振动期间交替反复。图21是显示半导体激光器201的振荡波长的时间变化的示意图。图21中,P1为第一振荡期间,P2为第二振荡期间,λa为各期间振荡波长的最小值、λb为各期间振荡波长的最大值、Tt为三角波的周期。The laser driver 204 supplies the semiconductor laser 201 with a triangular wave drive current that repeatedly increases and decreases at a predetermined rate of change over time as an injection current. Thus, the semiconductor laser 201 is driven so that the first oscillation period in which the oscillation wavelength continuously increases at a predetermined rate of change and the second oscillation period in which the oscillation wavelength continuously decreases at a predetermined rate of change are alternately repeated. FIG. 21 is a schematic diagram showing temporal changes in the oscillation wavelength of the semiconductor laser 201 . In Fig. 21, P1 is the first oscillation period, P2 is the second 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 Tt is the period of the triangle wave.
从半导体激光器201出射的激光由透镜203集光,入射到物体210。由物体210反射的光通过透镜203集光入射到半导体激光器201。光电二极管202将半导体激光器201的光输出变换为电流。电流-电压变换放大器205将光电二极管202的输出电流变换为电压并放大,信号提取电路206对电流-电压变换放大器205的输出电压进行二次微分。计数电路207针对第一振荡期间P1和第二振荡期间P2分别对包含在信号提取电路206的输出电压中的模跳脉冲(MHP)的个数进行计数。计算装置208根据半导体激光器201的最小振荡波长λa、最大振荡波长λb、第一振荡期间P1中的MHP的个数和第二振荡期间P2中的MHP的个数,计算出与物体210的距离和物体210的速度。根据这样的自混合型的激光计测器,可以进行半导体激光器201的半波长程度的分辨能力的位移计测,和与半导体激光器201的波长调制量成反比的分辨能力的距离计测。Laser light emitted from the semiconductor laser 201 is collected by the lens 203 and enters the object 210 . The light reflected by the object 210 is collected by the lens 203 and enters the semiconductor laser 201 . The photodiode 202 converts the light output of the semiconductor laser 201 into an electric current. The current-voltage conversion amplifier 205 converts the output current of the photodiode 202 into a voltage and amplifies it, and the signal extraction circuit 206 secondarily differentiates the output voltage of the current-voltage conversion amplifier 205 . The counting circuit 207 counts the number of mode-hop pulses (MHP) included in the output voltage of the signal extraction circuit 206 for the first oscillation period P1 and the second oscillation period P2 respectively. Calculation device 208 calculates the distance and The velocity of object 210 . According to such a self-mixing type laser measuring device, it is possible to measure the displacement of the resolving power of the semiconductor laser 201 at approximately half wavelength and the distance measurement of the resolving power inversely proportional to the wavelength modulation amount of the semiconductor laser 201 .
专利文献1日本特开2006-313080号公报
发明内容Contents of the invention
发明所要解决的问题The problem to be solved by the invention
根据自结合型激光计测器,和从前的FMCW雷达和驻波雷达、自混合型激光传感器等相比,可以以高分辨能力计测测定对象的位移和速度。但是,自结合型激光计测器由于和FFT一样其计算位移和速度需要一定程度的计测时间(在专利文献1的例子中,是指半导体激光器的振荡波长调制的载波的半周期),因此存在着在速度变化快的测定对象的计测中产生计测误差的问题。又,由于信号处理中需要对MHP的个数进行计数,因此还存在着难以实现半导体激光器的未满半波长的分辨能力的问题。Compared with conventional FMCW radar, standing wave radar, self-mixing laser sensor, etc., the self-integrated laser measuring device can measure the displacement and velocity of the measuring object with high resolution. However, since the self-integrated laser measuring device requires a certain amount of measurement time (in the example of
本发明为了解决上述课题,旨在提供一种能够以高的分辨能力计测物体的位移和速度,并能缩短计测所需的时间的物理量传感器和物理量计测方法。In order to solve the above problems, the present invention aims to provide a physical quantity sensor and a physical quantity measurement method capable of measuring the displacement and velocity of an object with high resolution and shortening the time required for the measurement.
解决问题的手段means of solving problems
本发明提供一种物理量传感器,包括:对测定对象发射激光的半导体激光器;振荡波长调制单元,启动所述半导体激光器动作,使振荡波长连续单调增加的第一振动期间和振动波长连续单调减少的第二振荡期间中的至少一种反复出现;检测单元,其检测包含干涉波形的电信号,所述干涉波形由从所述半导体激光器发射的激光和所述测定对象的返回光的自混合效应产生;信号提取单元,其在每次输入所述干涉波形时对包含在该检测单元的输出信号中的所述干涉波形的周期进行计测;周期修正单元,其通过将该信号提取单元的计测结果和基准周期进行比较对所述计测结果进行修正;计算单元,其基于由该周期修正单元修正后的各个周期,计算出所述测定对象的位移和速度中的至少一项。The present invention provides a physical quantity sensor, comprising: a semiconductor laser that emits laser light to an object to be measured; an oscillation wavelength modulation unit that activates the operation of the semiconductor laser to make the first oscillation period in which the oscillation wavelength continuously monotonically increases and the second oscillation period in which the oscillation wavelength continuously monotonically decreases at least one of the two oscillation periods occurs repeatedly; a detection unit that detects an electrical signal including an interference waveform generated by a self-mixing effect of laser light emitted from the semiconductor laser and returned light from the measurement object; a signal extraction unit that measures the period of the interference waveform included in the output signal of the detection unit each time the interference waveform is input; a period correction unit that uses the measurement result of the signal extraction unit The measurement result is corrected by comparing it with a reference period, and a calculating unit calculates at least one of a displacement and a velocity of the measurement object based on each period corrected by the period correcting unit.
又,在本发明的物理量传感器的一构成例中,所述计算单元,根据计测所述干涉波形的周期的采样时钟的频率、所述基准周期、所述半导体激光的平均波长、由所述周期修正单元修正后的周期相对于所述基准周期的变化量,计算出所述测定对象的位移和速度中的至少一项。Furthermore, in a configuration example of the physical quantity sensor according to the present invention, the calculating means is based on the frequency of the sampling clock for measuring the period of the interference waveform, the reference period, the average wavelength of the semiconductor laser, and the The period correcting unit calculates at least one of the displacement and the velocity of the measurement object from the amount of change in the period after correction with respect to the reference period.
又,在本发明的物理量传感器的一构成例中,所述周期修正单元,在由所述信号提取单元计测到的干涉波形的周期小于所述基准周期的规定数k倍时,其中k是小于1的正值,则将该干涉波形的周期和之后所计测得到的干涉波形的周期合并后的周期作为修正后的干涉波形的周期,将组合周期得到的波形作为一个波形;在由所述信号提取单元计测到的干涉波形的周期为所述基准周期的(m-k)倍以上且小于所述基准周期的(m+k)倍时,其中m是2以上的自然数,将该干涉波形的周期m等分后得到的周期分别作为修正后的周期,修正后的周期的波形具有m个。Also, in a configuration example of the physical quantity sensor of the present invention, the period correcting unit is configured to be smaller than a predetermined number k times of the reference period when the period of the interference waveform measured by the signal extracting unit, where k is If the positive value is less than 1, then the period after combining the period of the interference waveform and the period of the measured interference waveform is taken as the period of the corrected interference waveform, and the waveform obtained by combining the period is regarded as a waveform; When the period of the interference waveform measured by the signal extraction unit is more than (m-k) times the reference period and less than (m+k) times the reference period, wherein m is a natural number greater than 2, the interference waveform Periods obtained by dividing the period of m into equal parts are respectively used as corrected periods, and there are m waveforms of the corrected periods.
又,在本发明的物理量传感器的一构成例中,所述规定数k为0.5。In addition, in one configuration example of the physical quantity sensor of the present invention, the predetermined number k is 0.5.
又,在本发明的物理量传感器的一构成例中,所述周期修正单元,将所述测定对象静止时的所述干涉波形的周期或者所述即将修正之前所计测到的规定数的干涉波形的周期的平均值作为所述基准周期。Furthermore, in a configuration example of the physical quantity sensor of the present invention, the period correcting means adjusts the period of the interference waveform when the measurement object is stationary or a predetermined number of interference waveforms measured immediately before the correction. The average value of the cycle is used as the reference cycle.
又,在本发明的物理量传感器的一构成例中,还包括:计数单元,其对包含在所述检测单元的输出信号中的所述干涉波形的个数,在所述第一振荡期间和所述第二振荡期间分别进行计数;距离计算单元,其根据由该计数单元对干涉波形的个数进行计数的期间的最小振荡波长、最大振荡波长、以及所述计数单元的计数结果来计算与所述测定对象之间的距离;周期计算单元,其根据该距离计算单元所计算得到的距离求得所述干涉波形的周期,所述周期修正单元,将所述周期计算单元求得的周期作为所述基准周期。In addition, in one configuration example of the physical quantity sensor of the present invention, it further includes: a counting unit that counts the number of the interference waveforms included in the output signal of the detection unit during the first oscillation period and the Counting is carried out respectively during the second oscillation period; the distance calculation unit calculates the minimum oscillation wavelength, the maximum oscillation wavelength, and the counting result of the counting unit according to the counting unit during which the number of interference waveforms is counted. The distance between the measurement objects; the period calculation unit, which obtains the period of the interference waveform according to the distance calculated by the distance calculation unit, and the period correction unit uses the period obtained by the period calculation unit as the base cycle.
又,在本发明的物理量传感器的一构成例中,还包括:计数单元,其对包含在所述检测单元的输出信号中的所述干涉波形的个数,在所述第一振荡期间和所述第二振荡期间分别进行计数;距离比例个数计算单元,其通过计算所述干涉波形的个数的平均值,求得与所述半导体激光器和所述测定对象的平均距离成比例的干涉波形的个数即距离比例个数;周期计算单元,其根据所述距离比例个数计算所述干涉波形的周期,所述周期修正单元将所述周期计算单元所求得的周期作为所述基准周期。In addition, in one configuration example of the physical quantity sensor of the present invention, it further includes: a counting unit that counts the number of the interference waveforms included in the output signal of the detection unit during the first oscillation period and the Counting is performed during the second oscillation period; the distance proportional number calculation unit calculates the average value of the number of the interference waveforms to obtain an interference waveform proportional to the average distance between the semiconductor laser and the measurement object The number is the distance ratio number; the period calculation unit calculates the period of the interference waveform according to the distance ratio number, and the period correction unit uses the period obtained by the period calculation unit as the reference period .
本发明提供一种物理量计测方法,包括:振荡步骤,启动所述半导体激光器,使振荡波长连续单调增加的第一振荡期间和振动波长连续单调减少的第二振荡期间中的至少一种反复出现;检测步骤,检测包含干涉波形的电信号,所述干涉波形由从所述半导体激光器发射的激光和测定对象的返回光的自混合效应产生;信号提取步骤,在每次输入干涉波形时对包含于在所述检测步骤得到的输出信号中的所述干涉波形的周期进行计测;周期修正步骤,通过将该信号提取步骤的计测结果和基准周期进行比较,对所述计测结果进行修正;计算步骤,基于该周期修正步骤修正后的各个周期,计算出所述测定对象的位移和速度中的至少一项。The present invention provides a method for measuring a physical quantity, comprising: an oscillating step of activating the semiconductor laser so that at least one of the first oscillation period in which the oscillation wavelength continuously monotonically increases and the second oscillation period in which the oscillation wavelength continuously monotonically decreases occurs repeatedly ; The detection step detects an electrical signal containing an interference waveform generated by the self-mixing effect of the laser light emitted from the semiconductor laser and the return light of the measurement object; The period of the interference waveform in the output signal obtained in the detection step is measured; the period correction step is to correct the measurement result by comparing the measurement result of the signal extraction step with a reference period a calculating step of calculating at least one of the displacement and velocity of the measuring object based on each period corrected by the period correcting step;
发明效果Invention effect
根据本发明,通过基于各干涉波形的周期进行计算,能够以比从前高的分辨能力测定对象的位移和速度。又,相比于在从前的自混合型激光计测器中,需要耗费载波的半周期的计测时间,在本发明中,由于可以根据各个干涉波形的周期求得测定对象的位移和速度,因此大幅度缩短了计测所需要的时间,并且还可以与速度变化快的测定对象相对应。进一步的,在本发明中,通过将信号提取单元的计测结果和基准周期进行比较,可以修正干涉波形的周期误差,因此可提高位移和速度的计测精度。According to the present invention, by performing calculation based on the period of each interference waveform, the displacement and velocity of the object can be measured with higher resolution than before. In addition, compared with the conventional self-mixing type laser measuring instrument, which takes half a cycle of the carrier wave to measure, in the present invention, the displacement and velocity of the measurement object can be obtained from the cycle of each interference waveform, Therefore, the time required for measurement is greatly shortened, and it is also possible to cope with measurement objects whose speed changes rapidly. Furthermore, in the present invention, by comparing the measurement result of the signal extraction unit with the reference period, the period error of the interference waveform can be corrected, so the measurement accuracy of displacement and velocity can be improved.
附图说明Description of drawings
图1是显示本发明的第一实施方式涉及的物理量传感器的结构的框图。FIG. 1 is a block diagram showing the configuration of a physical quantity sensor according to a first embodiment of the present invention.
图2是示意性显示本发明第一实施方式涉及的电流-电压变换放大部的输出电压波形和滤波部的输出电压波形的波形图。2 is a waveform diagram schematically showing an output voltage waveform of a current-voltage conversion amplifying unit and an output voltage waveform of a filter unit according to the first embodiment of the present invention.
图3是用于说明模跳脉冲的图。Fig. 3 is a diagram for explaining mode-hopping pulses.
图4是显示半导体激光器的振荡波长和光电二极管的输出波形的关系的图。FIG. 4 is a graph showing the relationship between the oscillation wavelength of a semiconductor laser and the output waveform of a photodiode.
图5是显示本发明的第一实施方式的信号提取部的构成例的框图。5 is a block diagram showing a configuration example of a signal extraction unit according to the first embodiment of the present invention.
图6是用于说明本发明第一实施方式的信号提取部的动作的图。FIG. 6 is a diagram for explaining the operation of the signal extraction unit according to the first embodiment of the present invention.
图7是显示本发明的第一实施方式的计算部的构成例的框图。7 is a block diagram showing a configuration example of a calculation unit according to the first embodiment of the present invention.
图8是用于说明本发明的第一实施方式的周期修正部的动作的图。FIG. 8 is a diagram for explaining the operation of the period correcting unit according to the first embodiment of the present invention.
图9是用于说明本发明的第一实施方式的信号提取部的计测结果的修正原理的图。FIG. 9 is a diagram for explaining the principle of correction of the measurement result of the signal extraction unit according to the first embodiment of the present invention.
图10是显示模跳脉冲的周期的度数分布的图。Figure 10 is a graph showing the degree distribution of the period of a mode skip pulse.
图11是显示成为两倍周期的模跳脉冲周期的度数分布的图。Fig. 11 is a graph showing the degree distribution of the mode-hop pulse period which becomes twice the period.
图12是显示计数时欠缺的模跳脉冲中被两等分的模跳脉冲的周期的度数分布的图。Fig. 12 is a graph showing degree distribution of the period of the mode skip pulse divided into two equal parts among the mode skip pulses lacking in counting.
图13是显示计数时欠缺的模跳脉冲中被两等分的模跳脉冲的周期的度数分布的图。Fig. 13 is a graph showing the degree distribution of the period of the mode skip pulse divided into two equal parts among the mode skip pulses lacking in counting.
图14是显示本发明的第二实施方式中的计算部的构成例的框图。FIG. 14 is a block diagram showing a configuration example of a calculation unit in a second embodiment of the present invention.
图15是显示本发明的第三实施方式中的计算部的构成例的框图。15 is a block diagram showing a configuration example of a calculation unit in a third embodiment of the present invention.
图16是显示本发明的第三实施方式中的计数部的计数结果的时间变化的一例的图。FIG. 16 is a diagram showing an example of time changes in count results of the count unit in the third embodiment of the present invention.
图17是显示本发明第三实施方式中的计数部的计数结果的时间变化的其他实例的图。Fig. 17 is a diagram showing another example of the time change of the counting result of the counting unit in the third embodiment of the present invention.
图18是显示本发明的第四实施方式中的半导体激光器的振荡波长的时间变化的其他实例的图。FIG. 18 is a graph showing another example of the temporal change of the oscillation wavelength of the semiconductor laser in the fourth embodiment of the present invention.
图19是显示本发明第五实施方式涉及的物理量传感器的构成的框图。FIG. 19 is a block diagram showing the configuration of a physical quantity sensor according to a fifth embodiment of the present invention.
图20是显示现有的激光计测器的构成的框图。FIG. 20 is a block diagram showing the configuration of a conventional laser measuring device.
图21是显示图20的激光计测器中的半导体激光器的振荡波长的时间变化的1个实例的图。FIG. 21 is a graph showing an example of temporal changes in the oscillation wavelength of the semiconductor laser in the laser measuring device of FIG. 20 .
具体实施方式Detailed ways
第一实施方式first embodiment
下面参考附图对本发明的实施方式进行说明。图1是显示本发明第一实施方式涉及的物理量传感器的结构的框图。Embodiments of the present invention will be described below with reference to the drawings. FIG. 1 is a block diagram showing the configuration of a physical quantity sensor according to a first embodiment of the present invention.
图1的物理量传感器具有:对作为测定对象的物体10发射激光的半导体激光器1;将半导体激光器1的光输出变换为电信号的光电二极管2;透镜3,该透镜3对半导体激光器1发出的光进行集光并发射,同时对物体10的返回光进行集光并使该返回光入射到半导体激光器1;驱动半导体激光器1的作为振荡波长调制单元的激光驱动器4;将光电二极管2的输出电流变换为电压并放大的电流-电压变换放大部5;从电流-电压变换放大部5的输出电压中除去载波的滤波部6;对包含在滤波部6的输出电压中的作为自混合信号的模跳脉冲(下面称为MHP)的周期进行计测的信号提取部7;基于信号提取部7计测到的各周期计算物体10的位移和速度的计算部8;显示计算部8的计算结果的显示部9。The physical quantity sensor of FIG. 1 has: a
光电二极管2和电流-电压变换放大部5构成检测单元。下面为了便于说明,假设半导体激光器1采用不具有模跳现象的类型(VCSEL型、DFB激光型)。The
激光驱动器4将随着时间以一定的变化率反复增减三角波驱动电流作为注入电流提供给半导体激光器1。这样,半导体激光器1被驱动为,与注入电流大小成正比地使振荡波长以一定的变化率连续增加的第一振荡期间P1和振荡波长以一定的变化率连续减少的第二振荡期间P2交替重复。此时,半导体激光器1的振荡波长的时间变化如图21所示。在本实施方式中,半导体激光器1的振荡波长的变化速度必须是一定的。The
从半导体激光器1出射的激光通过透镜3集光并入射到物体10。由物体10反射的光通过透镜3集光,入射到半导体激光器1。但是也不是必须通过透镜3进行集光。光电二极管2设置在半导体激光器1的内部或者其附近,将半导体激光器1的光输出变换为电流。电流-电压变换放大部5将光电二极管2的输出电流变换为电压并放大。The laser light emitted from the
滤波部6具有从调制波中抽取重叠信号的功能。图2(A)示意性地显示电流-电压变换放大部5的输出电压波形,图2(B)示意性显示滤波部6的输出电压波形。这些图显示了从相当于光电二极管2的输出的图2(A)的波形(调制波)中,除去图2的半导体激光器1的振荡波形(载波)后,提取如图2(B)所示的MHP波形(干涉波形)的过程。The
接着,信号提取部7在每次MHP发生时计测包含在滤波部6的输出电压中的MHP的周期。这里对作为自混合信号的MHP进行说明。如图3所示,设镜面层1013到物体10的距离为L,激光的振荡波长为λ,满足下列的共振条件时,物体10的返回光和半导体激光器1的光共振器内的激光相互增强,使得激光输出稍有增加。Next, the
L=qλ/2…(1)L=qλ/2...(1)
在式子(1)中,q为整数。该现象,即使在来自物体10的散射光非常微弱时,通过半导体激光器1的共振器内体现的反射率的增加,产生放大作用,从而能充分地进行观测。In the formula (1), q is an integer. Even when the scattered light from the
图4是显示使得半导体激光器1的振荡波长以某个一定的比例变化时的振荡波长和光电二极管2的输出波形之间的关系的图。当满足式子(1)显示的L=qλ/2时,返回光和光共振器内的激光的相位差为0°(同相位),返回光和光共振器内的激光为最大相互增强,当L=qλ/2+λ/4时,相位差为180°(逆相位),返回光和光共振器内的激光为最大相互减弱。因此,如果使得半导体激光器1的振荡波长变化,激光输出交替重复出现增强和减弱的现象,通过发光二极管2检测此时的激光输出的话,如图4所示可得到一定周期的阶梯状的波形。这样的波形一般称为干涉条纹。该阶梯状的波形即各个干涉条纹的为MHP。如前所述,一定时间之内使得半导体激光器1的振荡波长变化时,MHP的个数与测定距离成正比地变化。FIG. 4 is a graph showing the relationship between the oscillation wavelength and the output waveform of the
图5是显示信号提取部7的构成例的框图。信号提取部7包括二值化部70和周期测定部71。FIG. 5 is a block diagram showing a configuration example of the
图6(A)~图6(D)是用于对信号提取部7的动作进行说明的图,图6(A)示意性显示滤波部6的输出电压的波形,即MHP的波形,图6(B)显示对应于图6(A)的二值化部70的输出,图6(C)是显示输入到信号提取部7的采样时钟CLK的图,图6(D)是显示对应于图6(B)的周期测定部71的测定结果的图。6(A) to 6(D) are diagrams for explaining the operation of the
首先,信号提取部7的二值化部70判定图6(A)所示的滤波部6的输出电压为高电平(H)还是低电平(L),并输出如图6(B)那样的判定结果。此时,二值化部70在滤波部6的输出电压上升达到了阈值TH1以上时判定为高电平,当滤波部6的输出电压下降达到了阈值TH2(TH2<TH1)时判定为低电平,从而使得滤波部6的输出二值化。First, the
周期测定部71在每次上升边缘发生时测定二值化部70的输出的上升边缘的周期(即,MHP的周期)。此时,周期测定部71以图6(C)所示的采样时钟CLK的周期为一个单位测定MHP的周期。在图6(D)的示例中,周期测定部71,依次测定Tα、Tβ、Tγ作为MHP的周期。根据图6(C)、图6(D)可知,周期Tα、Tβ、Tγ的大小分别为5[samplings]、4[samplings]、2[samplings]。采样时钟CLK的频率相对于MHP的能得到的最高频率为足够高。The
接着,计算部8基于信号提取部7的计测结果,根据各个MHP的周期的变化计算物体10的位移和速度。图7是显示计算部8的构成例的框图。计算部8包括存储部80、周期修正部81、物理量计算部82。Next, the
存储部80存储信号提取部7的计测结果。周期修正部81将物体10静止时的MHP的周期、计算得到的距离中的MHP的周期、或本次修正之前所计测到的规定数量的MHP的周期的移动平均值中的任意一个作为基准周期T0,通过比较信号提取部7本次的计测结果和基准周期T0,对信号提取部7的计测结果进行修正。图8(A)~图8(F)是用于说明周期修正部81的动作的图。The
当如图8(A)所示由信号提取部7计测得到的MHP的周期T小于0.5T0时,周期修正部81如图8(B)所示将MHP的周期T和下次计测得到的MHP的周期Tnext合成后的周期作为修正后的MHP的周期T’。When the period T of the MHP measured by the
又,当如图8(C)所示由信号提取部7计测得到的MHP的周期T在1.5T0以上且小于2.5T0时,周期修正部81如图8(D)所示将MHP的周期T两等分后得到的周期分别作为修正后的周期T1’、T2’。Also, when the period T of MHP measured by the
又,当如图8(E)所示由信号提取部7所计测到的MHP的周期T在2.5T0以上且小于3.5T0的时候,周期修正部81如图8(F)所示将MHP的周期T三等分后得到的周期分别作为修正后的周期T1’、T2’、T3’。在3.5T0以上的情况下也是同样的。即周期修正部81,在由信号提取部7所计测得到的MHP的周期T在(m-0.5)T0以上且小于(m+0.5)T0的情况下(m为2以上的自然数),将MHP的周期T进行m等分后得到的周期分别作为修正后的周期。周期修正部81在每次信号提取部7输出计测结果时进行如上的修正处理。Also, when the period T of the MHP measured by the
图9是说明信号提取部7的计测结果的修正原理的图,其示意性地显示了滤波部6的输出电压波形,即MHP的波形。FIG. 9 is a diagram illustrating the principle of correction of the measurement result of the
原本,MHP的周期根据与物体10之间的距离不同而不同,但是如果和物体10之间的距离不变,则MHP以同样的周期出现。然而,由于噪音的缘故,MHP的波形会产生欠缺,或产生不应当作为信号的波形,从而导致MHP的周期产生误差。Originally, the cycle of the MHP varies depending on the distance from the
如果产生信号的欠缺,则产生欠缺的地方的MHP的周期Tw是原本周期的大致两倍。即,MHP的周期为基准周期T0的大致两倍以上时,可判断为信号产生了欠缺。这样,即可通过将周期Tw两等分,来修正信号的欠缺。When a signal gap occurs, the cycle Tw of the MHP at the point where the gap occurs is approximately twice the original cycle. That is, when the period of the MHP is approximately twice or more than the reference period T0, it can be determined that the signal is missing. In this way, by dividing the period Tw into two equal parts, it is possible to correct the lack of a signal.
又,在计算了噪音的地方的MHP的周期Ts为原本周期的大致0.5倍。即,MHP的周期小于基准周期T0的大致0.5倍时,可判断为信号计数过剩。这样,通过使得周期Ts和下一次计测得到的周期Tnext相加,对错误计数了的噪声进行修正。Also, the period Ts of the MHP at the place where the noise is calculated is approximately 0.5 times the original period. That is, when the period of the MHP is less than approximately 0.5 times the reference period T0, it can be determined that the signal count is excessive. In this way, the erroneously counted noise is corrected by adding the period Ts to the period Tnext measured next time.
以上,是信号提取部7的计测结果的修正原理。在本实施方式中,使得用于确定视为计算了噪音的周期Ts的阈值为基准周期T0的0.5倍的值,使得用于确定视为产生了信号的欠缺的周期Tw的阈值不为基准周期T0的2倍的值,而是1.5倍的值,确定1.5倍的理由将在之后进行论述。The above is the principle of correction of the measurement result of the
接着,物理量计算部82根据由周期修正部81修正得到的MHP的各个周期相对于基准周期T0的变化计算物体10的位移和速度。设采样时钟的频率为fad[Hz]、基准周期为T0[samplings]、半导体激光器1的振荡平均波长为λ[m]、修正后的MHP的周期由基准周期T0延长n[samplings]时,该修正后的MHP的周期中物体10的位移D[m]如下式所示。Next, the physical
D=n×λ/(2×T0)…(2)D=n×λ/(2×T0)…(2)
修正后MHP的周期由基准周期T0缩短n[samplings]时,使得式子(2)的周期变化量n的符号为负即可。半导体激光器1的振动波长增加的第一振荡期间P1中,位移D为正时,物体10的移动方向是远离半导体激光器1的方向,当位移D为负时,物体10的移动方向为接近半导体激光器1的方向。另,在振荡波长减少的第二振荡期间P2中,当位移D为正,物体10的移动方向为接近半导体激光器1的方向,当位移D为负时,物体10的移动方向为远离半导体激光器1的方向。When the period of the modified MHP is shortened by n[samplings] from the reference period T0, the sign of the period change n in the formula (2) can be negative. In the first oscillation period P1 in which the vibration wavelength of the
又,由于修正后MHP的周期为(T0+n)/fad,该修正后的MHP的周期中物体10的速度V[m/s]为下式所示。Also, since the period of the corrected MHP is (T0+n)/fad, the velocity V [m/s] of the
V=n×λ/(2×T0)×fad/(T0+n) …(3)V=n×λ/(2×T0)×fad/(T0+n) …(3)
物理量计算部82可根据式子(2)计算物体10的位移D,根据式子(3)计算物体10的速度V。例如,设采样时钟的频率fad为16[MHz]、基准周期T0为160[samplings]、半导体激光器1的平均波长为850[nm],若修正后的MHP的周期由基准周期T0增长1[samplings],则可计算得到该MHP的周期中物体10的位移D为5.31[nm]、速度V为1.05[mm/s]。物理量计算部82对修正后的MHP的各周期进行上述那样的计算处理。The physical
显示部9显示计算部8的计算结果。The
这里,令半导体激光器1的振荡波长调制的载波(三角波)的每一半周期的和与物体10的距离相关联的MHP的个数为N1。在将物体10的最大速度的绝对值改为了载波每一周期的位移时则为λ/2×Na,每一载波半周期的MHP的个数为N1±Na。载波每一周期的位移以λ/2×Nb的速度移动时,每一载波半周期的MHP的个数为N1+Nb,这样即可观察与该个数对应的MHP的周期。为了求得物体10的位移D和速度V,可根据各个MHP的周期逆算每一载波半周期的MHP的个数,并根据该MHP的个数计算物体10的位移D和速度V。上述的式子(2)、(3)是基于这样的导出原理。Here, let N1 be the number of MHPs associated with the distance from the
在专利文献1所揭示的自混合型激光计测器中,物体的位移和速度的分辨能力大致为半导体激光器的半波长λ/2。与之相对的,本发明实施方式中,由于位移D和速度V的分辨能力为λ/2×n/T0,因此可实现不到波长λ/2的分辨能力,从而实现以比现有技术更高的分辨能力进行计测。In the self-mixing laser measuring device disclosed in
如上所述,在本实施方式中,物体10的位移D和速度可以以比现有技术高的分辨能力进行计测。又,在专利文献1所揭示的自混合型激光计测器中,需要耗费载波的半周期的计测时间,而在本实施方式中,由于可根据各个MHP的周期求得物体10的位移D和速度V,因此计测所需要的时间能够大幅度缩短,从而可以适用于速度变化快的物体10。进一步的,在本实施方式中,由于可修正MHP的周期的误差,从而能够提高位移D和速度V的计测精度。As described above, in the present embodiment, the displacement D and velocity of the
又,由于即使在物体10静止时,各个MHP的周期在正态分布上也存在偏差,因此对计算得到的位移进行移动平均等处理即可。Also, even when the
又,在本实施方式中,对物体10的位移和速度进行计测,但是也可仅对其中一项进行计测。Also, in this embodiment, the displacement and the velocity of the
接着,对将用于确定视为信号产生了欠缺的周期Tw的阈值设为基准周期T0的1.5倍的理由进行说明。半导体激光器1的振荡波长的变化为线性变化时,MHP的周期以基准周期T0为中心正态分布(图10)。Next, the reason why the threshold value for specifying the period Tw in which the signal is considered to be missing is set to 1.5 times the reference period T0 will be described. When the oscillation wavelength of the
这里,考虑MHP的波形产生了欠缺的情况。由于原本MHP的周期以T0为中心的正态分布,所以因为MHP的强度小而在计测时产生了欠缺时的MHP的周期为平均值2T0、标准偏差2σ的正态分布(图11中的f)。当欠缺了j[%]的MHP时,在第一振荡期间P1或第二振荡期间P2的任一个中,信号提取部7对MHP的个数进行计数,其结果,MHP的个数为N时,由于该欠缺周期为2倍的MHP的周期的度数为Nw(=j[%]·N)。又,由于计测时的欠缺减少后的大致为T0的周期的度数为如图11所示的g,图11的h所示的度数的减少量为2Nw(=2j[%])。因此,在第一振荡期间P1和第二振荡期间P2中的任一个中,未产生MHP欠缺时的原本的MHP的个数N’能够以以下的式子来表示。Here, consider a case where the waveform of the MHP is lacking. Originally, the cycle of the MHP is normally distributed around T0, so the cycle of the MHP when there is a gap in the measurement due to the low intensity of the MHP is a normal distribution with an average value of 2T0 and a standard deviation of 2σ (Fig. 11 f). When the MHP of j [%] is missing, the
N’=N+j[%]=N+Nw …(4)N'=N+j[%]=N+Nw ...(4)
接着,考虑用于修正MHP的周期的计测结果的阈值。此处,假设的是由于计测时的欠缺周期为2倍的MHP的周期的度数Nw中p[%]由噪音一分为二的情况。欠缺了的MHP中被一分为二的MHP的周期的度数为Nw’(=j·p[%]·N)。再度一分为二的MHP的周期的度数分布如图12所示。如果将视为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, a threshold value for correcting the measurement result of the cycle of the MHP is considered. Here, a case is assumed in which p[%] is divided into two by noise in the frequency Nw of the cycle of the MHP which is doubled due to the missing cycle at the time of measurement. Among the missing MHPs, the degree of cycle of the MHP divided into two is Nw' (=j·p[%]·N). The degree distribution of the cycle of the MHP divided again in two is shown in FIG. 12 . If the threshold value of the period regarded as Nw is set to 1.5T0, the degree of the period of the MHP with a period of 0.5T0 or less is 0.5Nw' (=0.5p[%]·Nw), and the degree of the MHP with a period of 0.5T0 to 1.5T0 The degree of cycle is Nw' (=p[%]·Nw), and the degree of cycle of MHP having a cycle of 1.5T0 or more is 0.5Nw'(=0.5p[%]·Nw).
因此,所有MHP的周期的度数分布如图13所示,设与上述Ts对应的周期的度数Ns的阈值为0.5T0,与上述Tw对应的周期的度数Nw的阈值为1.5T0时,计数结果N能够以以下式子表示。Therefore, the degree distribution of all MHP cycles is shown in Figure 13. Assuming that the threshold value of the cycle degree Ns corresponding to the above Ts is 0.5T0, and the threshold value of the cycle degree Nw corresponding to the above Tw is 1.5T0, the counting result N It can be represented by the following formula.
N=(N’-2Nw)+(Nw-Nw’)+2Nw’=N’-Nw+Nw’N=(N'-2Nw)+(Nw-Nw')+2Nw'=N'-Nw+Nw'
…(5)...(5)
根据式子(5),修正后的结果如以下所示,计数时未产生MHP欠缺的情况下原本的MHP的个数N’可计算得到。According to the formula (5), the corrected result is as follows, and the original number N' of MHP can be calculated when there is no shortage of MHP during counting.
N-0.5Nw’+(0.5Nw’+(Nw-Nw’))N-0.5Nw’+(0.5Nw’+(Nw-Nw’))
=(N-Nw+Nw’)+(0.5Nw’+(Nw-Nw’))=(N-Nw+Nw’)+(0.5Nw’+(Nw-Nw’))
=N’ …(6)=N' ...(6)
根据以上,如果使得计算度数Nw时的周期的阈值为基准周期T0的1.5倍,则可以修正计数结果N。如果三角波的每个半周期的采样时钟数为M,则MHP的周期T和计数结果N之间具有T=M/N的关系,由于M为一定值,可知,用于确定视为产生信号欠缺的周期Tw的阈值和计数结果N的情况相同,可以为基准周期T0的1.5倍。Based on the above, if the threshold value of the cycle when calculating the degree Nw is set to be 1.5 times the reference cycle T0, the counting result N can be corrected. If the number of sampling clocks in each half cycle of the triangular wave is M, then there is a relationship of T=M/N between the period T of the MHP and the counting result N. Since M is a certain value, it can be seen that it is used to determine the lack of signal The threshold value of the cycle Tw is the same as that of the counting result N, and may be 1.5 times the reference cycle T0.
又,在本实施方式中,基准周期T0设为物体10静止状态下的MHP的周期,但是不限于此,计算部8也可将修正之前所计测得到的规定个数的MHP的周期的移动平均值作为基准周期T0。根据这种方法,即使在物体10无法静止的情况下,也可求得基准周期T0。Also, in the present embodiment, the reference period T0 is set as the period of the MHP in the stationary state of the
第二实施方式second embodiment
接着,对本发明的第二实施方式进行说明。图14是显示第二实施方式涉及的计算部8的构成例的框图。计算部8包括存储部80、周期修正部81、物理量计算部82、计数部83、距离计算部84、周期计算部85。物理量传感器整体构成和第一实施方式相同,但是半导体激光器1的振荡波长的变化速度为一定,且振荡波长的最大值λb和振荡波长的最小值λa分别为一定,其差λb-λa也需要为一定。Next, a second embodiment of the present invention will be described. FIG. 14 is a block diagram showing a configuration example of the
计数部83针对第一振荡期间P1和第二振荡期间P2分别对包含在滤波部6的输出中的MHP的个数进行计数。计数部83可以利用由逻辑门构成的计数器,也可以利用FFT(Fast FourierTransform)对MHP的频率(即,单位时间内MHP的个数)进行计测。The
接着,距离计算部84基于半导体激光器1的最小振荡波长λa和最大振荡波长λb和计数部83所计数得到的MHP的个数,计算与物体10的距离。在本实施方式中,物体10的状态可以是满足规定条件的微小位移状态,或是比微小位移状态移动大的位移状态中的任何一种。振荡期间P1和振动期间P2的每个期间的物体10的平均位移设为V时,微小位移状态是指满足(λb-λa)/λb>V/Lb的状态(但,Lb是指时刻t时的距离),位移状态是指满足(λb-λa)/λb≤V/Lb的状态。Next, the
首先,距离计算部84如以下式子计算现时刻t的距离的候补值Lα(t)、Lβ(t)和速度的候补值Vα(t)、Vβ(t)。First, the
Lα(t)=λa×λb×(MHP(t-1)+MHP(t))Lα(t)=λa×λb×(MHP(t-1)+MHP(t))
/{4×(λb-λa)} …(7)/{4×(λb-λa)} ... (7)
Lβ(t)=λa×λb×(|MHP(t-1)-MHP(t)|)/{4×(λb-λa)} …(8)Lβ(t)=λa×λb×(|MHP(t-1)-MHP(t)|)/{4×(λb-λa)} …(8)
Vα(t)=(MHP(t-1)-MHP(t))×λb/4 …(9)Vα(t)=(MHP(t-1)-MHP(t))×λb/4 …(9)
Vβ(t)=(MHP(t-1)+MHP(t))×λb/4 …(10)Vβ(t)=(MHP(t-1)+MHP(t))×λb/4 …(10)
在式子(7)~(10)中,MHP(t)是现时刻t所计算得到的MHP的个数,MHP(t-1)是指在MHP(t)的前一次所计算得到的MHP的个数。例如,如果MHP(t)为第一振荡期间P1的计数结果的话,MHP(t-1)为第二振荡期间P2的计数结果,反之,如果MHP(t)为第二振荡期间P2的计数结果的话,则MHP(t-1)为第一振荡期间P1的计数结果。In formulas (7) to (10), MHP(t) is the number of MHPs calculated at the current time t, and MHP(t-1) refers to the MHPs calculated at the previous time of MHP(t). the number of . For example, if MHP(t) is the counting result of the first oscillation period P1, MHP(t-1) is the counting result of the second oscillation period P2, on the contrary, if MHP(t) is the counting result of the second oscillation period P2 If , then MHP(t-1) is the counting result of the first oscillation period P1.
候补值Lα(t)、Vα(t)是假设物体10为微小位移状态时所计算得到的值,候补值Lβ(t)、Vβ(t)是假设物体10为位移状态时所计算得到的值。距离计算部84在通过计数部83对MHP的个数进行测定的各个时刻(各振荡期间)进行式子(7)~(10)的计算。The candidate values Lα(t) and Vα(t) are values calculated assuming that the
接着,距离计算部84对微小位移状态和位移状态,分别通过下列式子计算现时刻t的距离的候补值与之前时刻的距离的候补值之间的差即历史位移。又,在式子(11)、式子(12)中,现时刻t的前一次所计算得到的距离的候补值为Lα(t-1)、Lβ(t-1)。Next, the
Vcalα(t)=Lα(t)-Lα(t-1) …(11)Vcalα(t)=Lα(t)-Lα(t-1) ...(11)
Vcalβ(t)=Lβ(t)-Lβ(t-1) …(12)Vcalβ(t)=Lβ(t)-Lβ(t-1) ...(12)
历史位移Vcalα(t)为假设物体10为微小位移状态时计算得到的值,历史位移Vcalβ(t)是假设物体10为位移状态时计算得到的值。距离计算部84在由计数部83测定MHP的个数的各个时刻进行式子(11)~式子(12)的计算。又,式子(9)~式子(12)中,将物体10靠近本实施方式的物理量传感器的方向确定为正速度,将远离的方向确定为负速度。The historical displacement Vcalα(t) is a value calculated assuming that the
接着,距离计算部84采用式子(7)~式子(12)的计算结果,判定物体10的状态。Next, the
如专利文献1所记载的那样,在假设物体10在微小位移状态计算得到的历史位移Vcalα(t)的符号为一定,且假设物体10在微小位移状态计算得到的速度的候补值Vα(t)和历史位移Vcalα(t)的绝对值的平均值相等时,距离计算部84判定为物体10以微小位移状态做等速度运动。As described in
又,如专利文献1记载的那样,在假设物体10在位移状态计算得到的历史位移Vcalβ(t)的符号为一定,且假设物体10在位移状态计算得到的速度的候补值Vβ(t)和历史位移Vcalβ(t)的绝对值的平均值相等时,距离计算部84判定为物体10以位移状态做等速度运动。Also, as described in
又,如专利文献1记载的那样,在假设物体10在微小位移状态计算得到的历史位移Vcalα(t)的符号在测定MHP的个数的各时刻反转,且假设物体10在微小位移状态计算得到的速度的候补值Vα(t)和历史位移Vcalα(t)的绝对值的平均值不一致时,距离计算部84判定为物体10以微小位移状态做等速度运动以外的运动。Also, as described in
又,如果着眼于速度的候补值Vβ(t),Vβ(t)的绝对值为常数,该值和半导体激光器1的波长变化量(λb-λa)/λb相等。此时,在假设物体10在位移状态而计算得到的速度的候补值Vβ(t)的绝对值和波长变化率相等,且假设物体10在微小位移状态而计算得到的速度的候补值Vα(t)和历史位移Vcalα(t)的绝对值的平均值不一致时,距离计算部84判定为物体10以微小位移状态做等速度运动以外的运动。Also, focusing on the candidate velocity value Vβ(t), the absolute value of Vβ(t) is constant, and this value is equal to the wavelength variation (λb−λa)/λb of the
又,如专利文献1所记载的,在假设物体10在位移状态而计算得到的历史位移Vcalβ(t)的符号在每个测定MHP的个数的时刻反转,且假设物体10在位移状态而计算得到的速度的候补值Vβ(t)和历史位移Vcalβ(t)的绝对值的平均值不一致时,距离计算部84判定为物体10以位移状态做等速度运动以外的运动。Also, as described in
又,如果着眼于速度的候补值Vα(t),Vα(t)的绝对值为常数,该值和半导体激光器1的波长变化率(λb-λa)/λb相等。从而距离计算部84,在假设物体10在微小位移状态而计算得到的速度的候补值V α(t)的绝对值与波长变化率相等,且假设物体10在位移状态而计算得到的速度的候补值Vβ(t)和历史位移Vcalβ(t)的绝对值的平均值不一致时,可判定物体10以位移状态做等速度运动以外的运动。Also, focusing on the candidate velocity value Vα(t), the absolute value of Vα(t) is constant, and this value is equal to the wavelength change rate (λb−λa)/λb of the
距离计算部84根据上述判定结果确定与物体10之间的距离。即,距离计算部84在判定为物体10以微小位移状态做等速度运动时,将距离的候补值Lα(t)作为与物体10的距离,当判定为物体以位移状态做等速度运动时,以距离的候补值Lβ(t)作为与物体10的距离。The
又,距离计算部84,在判定为物体10以微小位移状态做等速度运动以外的运动时,以距离的候补值Lα(t)作为与物体10的距离。然而,实际的距离为距离的候补值Lα(t)的平均值。又,距离计算部84,在判定为物体以位移状态做等速度运动以外的运动时,以距离的候补值Lβ(t)作为与物体10的距离。但是,实际的距离是距离的候补值Lβ(t)的平均值。Further, the
接着,周期计算部85根据距离计算部84所计算得到的距离求得MHP的周期。MHP的频率和测定距离成正比,MHP的周期和测定距离成反比。这里,预先求得MHP的周期和距离的关系并存储在周期计算部85的数据库(未图示)中的话,周期计算部85可通过从数据库中取得与距离计算部84所计算得到的距离相对应的MHP的周期,来求得MHP的周期。或者,如果预先求得并设定表示MHP的周期和距离之间关系的式子,周期计算部85可通过将距离计算部84计算得到的距离带入数式,求得MHP的周期。Next, the
周期修正部81将周期计算部85求得的周期作为基准周期T0,如第一实施方式中说明的那样对信号提取部7的计测结果进行修正即可。物理量计算部82的动作和第一实施方式的相同。在本实施方式中,即使在无法使得物体10静止的情况下,也可求得基准周期T0。The
第三实施方式third embodiment
接着,对本发明的第三实施方式进行说明。图15是本发明第三实施方式涉及的计算部8的构成例的框图。计算部8具有:存储部80;周期修正部81;物理量计算部82;对包括在滤波部6的输出电压中的MHP的个数进行计数的计数部86;存储计数部86的计数结果等的存储部87;通过计算计数部86的计数结果的平均值,求得与半导体激光器1和物体10之间的平均距离成正比的MHP的个数(下面称为距离比例个数)NL的距离比例个数计算部88;符号赋予部89,其根据计数部86的前一次的计数结果与使用比该计数结果更之前的计数结果计算得到的距离比例个数NL的两倍数的大小关系,对计数部86的最新计数结果赋予正负符号;根据距离比例个数NL计算MHP的周期的周期计算部90。物理量传感器的整体结构可与第一实施方式的相同。Next, a third embodiment of the present invention will be described. FIG. 15 is a block diagram of a configuration example of the
计数部86针对第一振荡期间P1和第二振荡期间P2分别对包含于滤波部6的输出中的MHP的个数进行计数。计数部86可采用由逻辑门构成的计数器,也可以是采用FFT对MHP的频率(即单位时间内MHP的个数)进行计数的装置。计数部86的计数结果存储在存储部87中。The
距离比例个数计算部88根据计数部86的计数结果求得距离比例个数NL。图16是对距离比例个数计算部88的动作进行说明的图,其显示了计数部86的计数结果的时间变化。在图16中,Nu为第一振荡期间P1的计数结果,Nd是第二振荡期间P2的计数结果。物体10的距离变化率比半导体激光器1的振荡波长变化率小,物体10做简谐振动时,计数结果Nu的时间变化和计数结果Nd的时间变化如图16所示为相互的相位差为180度的正弦波形。专利文献1中,将此时物体19的状态作为微小位移状态。The distance-proportional
由图21可知,由于第一振荡期间P1和第二振荡期间P2交替出现,计数结果Nu和计数结果Nd也交替出现。计数结果Nu、Nd为,距离比例个数NL和与物体的位移成正比的MHP的个数(下面称为位移比例个数)NV的和或是差。距离比例个数NL相当于图16所示的正弦波形的平均值。又,计数结果Nu或Nd与距离比例个数NL的差,相当于位移比例个数NV。It can be seen from FIG. 21 that since the first oscillation period P1 and the second oscillation period P2 appear alternately, the counting result Nu and the counting result Nd also appear alternately. The counting results Nu and Nd are the sum or difference of the distance proportional number NL and the number of MHPs proportional to the displacement of the object (hereinafter referred to as the displacement proportional number) NV. The distance proportional number NL corresponds to the average value of the sinusoidal waveform shown in FIG. 16 . Also, the difference between the count result Nu or Nd and the distance proportional number NL corresponds to the displacement proportional number NV.
距离比例个数计算部88,通过如以下式子那样对到现时刻t的前两次为止所计测到的偶数次计数结果的平均值进行计算,来计算距离比例个数NL。The distance-proportional
NL={N(t-2)+N(t-3)}/2 …(13)NL={N(t-2)+N(t-3)}/2 ...(13)
式子(13)中,N(t-2)表示现时刻t的前两次所计测得到的MHP的个数N,N(t-3)表示现时刻t的前三次所计测得到的MHP的个数N。若现时刻t的计数结果N(t)为第一振动期间P1的计数结果Nu,则前两次的计数结果N(t-2)也是第一振荡期间P1的计数结果Nu,前三次的计数结果N(t-3)为第二振荡期间P2的计数结果Nd。相反地,如果现时刻t的计数结果N(t)为第二振动期间P2的计数结果Nd,则前两次的计数结果N(t-2)也是第二振荡期间P2的计数结果Nd,前三次的计数结果N(t-3)为第一振荡期间P1的计数结果Nu。In formula (13), N(t-2) represents the number N of MHPs measured two times before the current time t, and N(t-3) represents the number of MHPs measured three times before the current time t. The number N of MHPs. If the counting result N(t) at the current moment t is the counting result Nu of the first oscillation period P1, then the counting result N(t-2) of the first two times is also the counting result Nu of the first oscillation period P1, and the counting result of the first three times is also the counting result Nu of the first oscillation period P1. The result N(t-3) is the counting result Nd of the second oscillation period P2. Conversely, if the counting result N(t) at the current moment t is the counting result Nd of the second oscillation period P2, then the counting results N(t-2) of the previous two times are also the counting results Nd of the second oscillation period P2, and the previous The counting result N(t-3) of three times is the counting result Nu of the first oscillation period P1.
式子(13)是以两次计数结果求距离比例个数NL的情形的式子,在使用2m(m为正整数)次的计数结果的情况下,距离比例个数计算部88如下式那样计算距离比例个数NL。Equation (13) is an expression for the case of calculating the distance proportional number NL from two counting results. When using 2m (m is a positive integer) counting results, the distance proportional
NL={N(t-2m-1)+N(t-2m)+…+N(t-2)}/2mNL={N(t-2m-1)+N(t-2m)+…+N(t-2)}/2m
…(14)...(14)
但是,式子(13)、式子(14)是用于与物体10之间的距离和物体10的速度的计测开始初期所使用的式子,从中途采用使用后述的带符号的计数结果的下列式子代替式子(13),从而进行距离比例个数NL的计算。However, expressions (13) and (14) are expressions used at the beginning of the measurement of the distance to the
NL={N’(t-2)+N’(t-3)}/2 …(15)NL={N’(t-2)+N’(t-3)}/2 …(15)
N’(t-2)是对前两次的计数结果N(t-2)进行后述的符号赋予处理之后的带符号的计数结果,N’(t-3)是对前三次的计数结果N(t-3)进行后述的符号赋予处理之后的带符号的计数结果。式子(15)在现时刻t的计数结果N(t)成为从MHP的个数计测开始第七次的计数结果之后使用。N'(t-2) is the signed counting result after performing the sign assignment process described later on the previous two counting results N(t-2), and N'(t-3) is the counting result of the first three times N(t-3) is the signed count result after performing the sign assignment process described later. Expression (15) is used after the count result N(t) at the current time t becomes the seventh count result from the measurement of the number of MHPs.
又,计测开始初期采用式子(14)的情况下,中途开始不采用式子(14),而采用使用带符号的计数结果的下列式子计算距离比例个数NL。Also, when the formula (14) is used at the beginning of the measurement, the distance ratio number NL is calculated using the following formula using the signed count result without using the formula (14).
NL={N’(t-2m-1)+N’(t-2m)+…+N’(t-2)}/2mNL={N'(t-2m-1)+N'(t-2m)+...+N'(t-2)}/2m
…(16)...(16)
式子(16)在现时刻t的计数结果N(t)成为从MHP的个数计测开始第(2m×2+3)次的计数结果之后使用。Expression (16) is used after the count result N(t) at the current time t becomes the (2m×2+3)-th count result from the start of the counting of the number of MHPs.
距离比例个数NL存储在存储部87。距离比例个数计算部88在各个由计数部86测定MHP的个数的时刻(各振荡周期),进行如上的距离比例个数NL的计算处理。The distance ratio number NL is stored in the
另外,用于距离比例个数NL的计算的计数结果足够多时,也可以以奇数次的计数结果来计算距离比例个数NL。In addition, when there are enough counting results for calculating the distance proportional number NL, the distance proportional number NL can also be calculated with an odd number of counting results.
接着,符号赋予部89根据现时刻t的前一次所计测得到的计数结果N(t-1)和距离比例个数NL的2倍数2NL的大小关系,对计数部86的计数结果N(t)赋予正负符号。符号赋予部89具体执行下列式子。Next, the
如果N(t-1)≥2NL那么N’(t)→-N(t) …(17)If N(t-1)≥2NL then N’(t)→-N(t) …(17)
如果N(t-1)<2NL那么N’(t)→+N(t) …(18)If N(t-1)<2NL then N’(t)→+N(t) …(18)
图17是用于说明符号赋予部9的动作,并显示计数部86的计数结果的时间变化的图。物体10的距离变化率比半导体激光器1的振荡波长变化率大时,计数结果Nu的时间变化为图17的170所示的负侧的波形折返到正侧的形状,同样计数结果Nd的时间变化为图17的171所示的负侧的波形折返到正侧的形状。专利文献1中,产生该计数结果的折返的部分的物体10的状态为位移状态。另一方面,不产生计数结果的折返的部分的物体10的状态为所述的微小位移状态。FIG. 17 is a diagram for explaining the operation of the
为了求得包含位移状态的振动中物体10的物理量,判定物体10为位移状态还是微小位移状态,当物体10为位移状态时,向正侧折返的计数结果如图17的170、171所示的轨迹那样需要修正。式子(17)、式子(18)是用于判定物体10是位移状态还是微小位移状态的式子。在图17中产生计数结果折返的位移状态下,N(t-1)≥2NL成立。从而,如式子(17)所示的,当N(t-1)≥2NL成立时,对计数部86的现时刻t的计数结果N(t)赋予负符号,作为带符号的计数结果N’(t)。In order to obtain the physical quantity of the
另一方面,在图16和图17中未产生计数结果折返的微小位移状态下,N(t-1)<2NL成立。从而,如式子(18)所示,当N(t-1)<2NL成立时,对计数部86的现时刻t的计数结果N(t)为赋予了正符号,作为带符号的计数结果N’(t)。On the other hand, N(t-1)<2NL is established in the minute displacement state in which the count result does not wrap around in FIG. 16 and FIG. 17 . Therefore, as shown in the formula (18), when N(t-1)<2NL holds true, the counting result N(t) of the current time t of the
带符号的计数结果N’(t)存储在存储部87。符号赋予部89在每个由计数部86测定MHP的个数的时刻(各振荡周期),进行上述那样的符号赋予处理。The signed count result N'(t) is stored in the
又,式子(17)的成立条件可以是N(t-1)>2NL,式子(18)的成立条件可以是N(t-1)≤2NL。Also, the establishment condition of formula (17) may be N(t-1)>2NL, and the establishment condition of formula (18) may be N(t-1)≤2NL.
接着,周期计算部90根据距离比例个数NL如下式所示计算MHP的周期T。Next, the
T=C/(2×f×NL) …(19)T=C/(2×f×NL) …(19)
此处,f为三角波的频率,C为光速。Here, f is the frequency of the triangle wave, and C is the speed of light.
周期修正部81以周期计算部90计算得到的周期作为基准周期T0,如第一实施方式所说明的那样对信号提取部7的计测结果进行修正。物理量计算部82的动作和第一实施方式相同。在本实施方式中,即使在无法使物体10静止的情况下,也可求得基准周期T0。The
第四实施方式Fourth Embodiment
接着,对本发明的第四实施方式进行说明。第一~第三实施方式中,使半导体激光器1振荡为三角波状,但是不限于此,也可在第三实施方式中如图18所示使得半导体激光器1振荡为锯齿波状。即,在本实施方式中,只要驱动半导体激光器1使得第一振荡期间P1或第二振荡期间P2中的任一个反复出现即可。但是,在第二实施方式中,需要使半导体激光器1振荡为三角波状。Next, a fourth embodiment of the present invention will be described. In the first to third embodiments, the
即使在如本实施方式一样使半导体激光器1振荡为锯齿波状的情况下,半导体激光器1的振荡波长的变化速度也需要是一定的。第一振荡期间P1或第二振动期间P2中的动作和三角波振荡的情况相同。也可如图18所示,仅第一振荡期间P1重复出现的锯齿波状的振荡的情况下,重复进行第一振荡期间P1的处理,当然,仅第二振荡期间P2重复出现的锯齿波状的振荡的情况下,也可重复进行第二振荡期间P2的处理。Even when the
第五实施方式Fifth Embodiment
接着,对本发明的第五实施方式进行说明。第一~第四实施方式中,虽然采用光电二极管2和电流-电压变换放大部5作为检测包含MHP波形的电信号的检测单元,但也可不采用光电二极管提取MHP波形。图19是显示本发明的第五实施方式涉及的物理量传感器的构成的框图,和图1相同的结构标注同样的符号。本实施方式的物理量传感器,采用电压检测部12作为检测单元,来代替第一实施方式中的光电二极管2和电流-电压变换放大部5。Next, a fifth embodiment of the present invention will be described. In the first to fourth embodiments, although the
电压检测部12对半导体激光器1的端子间电压,即阳极-阴极间电压进行检测并放大。由半导体激光器1发射出的激光和物体10的返回光产生干涉时,半导体激光器1的端子间电压出现MHP波形。从而,可从半导体激光器1的端子间电压提取MHP波形。滤波部6从电压检测部12的输出电压去除载波。物理量传感器的其他结构和第一实施方式的相同。The
这样,在本实施方式中,可不使用光电二极管提取MHP波形,和第一实施方式比较能够减少物理量传感器的部件数量,从而降低物理量传感器的制造成本。又,在本实施方式中,由于不采用光电二极管,从而可去除干扰光的影响。In this way, in this embodiment, the MHP waveform can be extracted without using a photodiode, and the number of components of the physical quantity sensor can be reduced compared with the first embodiment, thereby reducing the manufacturing cost of the physical quantity sensor. Also, in this embodiment, since no photodiode is used, the influence of disturbance light can be eliminated.
又,在第一~第五实施方式中至少信号提取部7和计算部8可以由,例如,具有CPU、存储器和接口的计算机和对这些硬件资源进行控制的程序来实现。CPU根据存储在存储器中的程序进行第一~第五实施方式中所说明的处理。Also, in the first to fifth embodiments, at least the
Claims (14)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101314864A CN102192707B (en) | 2010-03-04 | 2010-03-04 | Physical quantity sensor and physical quantity measurement method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101314864A CN102192707B (en) | 2010-03-04 | 2010-03-04 | Physical quantity sensor and physical quantity measurement method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102192707A true CN102192707A (en) | 2011-09-21 |
CN102192707B CN102192707B (en) | 2013-11-27 |
Family
ID=44601284
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010101314864A Expired - Fee Related CN102192707B (en) | 2010-03-04 | 2010-03-04 | Physical quantity sensor and physical quantity measurement method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102192707B (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102792183A (en) * | 2010-03-10 | 2012-11-21 | 阿自倍尔株式会社 | Physical quantity sensor and physical quantity measuring method |
US8996326B2 (en) | 2009-06-29 | 2015-03-31 | Azbil Corporation | Counting device, physical quantity sensor, counting method, and physical quantity measuring method |
CN105759072A (en) * | 2014-12-02 | 2016-07-13 | 财团法人工业技术研究院 | Optical anemometry system |
CN106643543A (en) * | 2016-12-30 | 2017-05-10 | 东北石油大学 | Method for detecting deformation of surface of magnetic disk storage |
CN111637906A (en) * | 2020-05-11 | 2020-09-08 | 清华大学深圳国际研究生院 | Fiber grating demodulation device and method based on self-differencing coherence |
CN115079145A (en) * | 2022-08-18 | 2022-09-20 | 深圳煜炜光学科技有限公司 | Method and device for improving anti-interference capability of laser radar |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0753804A1 (en) * | 1995-07-13 | 1997-01-15 | RENISHAW plc | Laser measuring apparatus |
CN2310325Y (en) * | 1997-07-09 | 1999-03-10 | 张应宁 | Slope auxiliary time-measuring precision pulse laser speed-measuring instrument |
CN102116775A (en) * | 2010-01-05 | 2011-07-06 | 株式会社山武 | Velocity measuring device and method |
-
2010
- 2010-03-04 CN CN2010101314864A patent/CN102192707B/en not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0753804A1 (en) * | 1995-07-13 | 1997-01-15 | RENISHAW plc | Laser measuring apparatus |
CN2310325Y (en) * | 1997-07-09 | 1999-03-10 | 张应宁 | Slope auxiliary time-measuring precision pulse laser speed-measuring instrument |
CN102116775A (en) * | 2010-01-05 | 2011-07-06 | 株式会社山武 | Velocity measuring device and method |
Non-Patent Citations (1)
Title |
---|
杨新刚: "激光测速仪在秒流量厚度控制系统中的应用", 《酒钢科技》, no. 02, 28 February 2011 (2011-02-28), pages 75 - 78 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8996326B2 (en) | 2009-06-29 | 2015-03-31 | Azbil Corporation | Counting device, physical quantity sensor, counting method, and physical quantity measuring method |
CN102792183A (en) * | 2010-03-10 | 2012-11-21 | 阿自倍尔株式会社 | Physical quantity sensor and physical quantity measuring method |
US8982336B2 (en) | 2010-03-10 | 2015-03-17 | Azbil Corporation | Physical quantity sensor and physical quantity measuring method |
CN105759072A (en) * | 2014-12-02 | 2016-07-13 | 财团法人工业技术研究院 | Optical anemometry system |
CN105759072B (en) * | 2014-12-02 | 2020-04-14 | 财团法人工业技术研究院 | Optical wind measurement system |
CN106643543A (en) * | 2016-12-30 | 2017-05-10 | 东北石油大学 | Method for detecting deformation of surface of magnetic disk storage |
CN111637906A (en) * | 2020-05-11 | 2020-09-08 | 清华大学深圳国际研究生院 | Fiber grating demodulation device and method based on self-differencing coherence |
CN115079145A (en) * | 2022-08-18 | 2022-09-20 | 深圳煜炜光学科技有限公司 | Method and device for improving anti-interference capability of laser radar |
Also Published As
Publication number | Publication date |
---|---|
CN102192707B (en) | 2013-11-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102792183B (en) | Physical quantity sensor and physical quantity measuring method | |
KR100945209B1 (en) | Counting device, rangefinder, counting method and distance measuring method | |
US7961302B2 (en) | Physical quantity sensor and physical quantity measuring method | |
JP5545915B2 (en) | Counting device, distance meter, counting method, and distance measuring method | |
US7649631B2 (en) | Distance/speed meter and distance/speed measuring method | |
CN102192707A (en) | Physical quantity transducer and physical quantity measuring method | |
JP2009014701A (en) | Distance / speed meter and distance / speed measurement method | |
US8537341B2 (en) | Physical quantity sensor and physical quantity measuring method | |
JP5596915B2 (en) | Physical quantity sensor and physical quantity measuring method | |
JP5081778B2 (en) | Vibration amplitude measuring apparatus and vibration amplitude measuring method | |
JP5421568B2 (en) | Physical quantity sensor and physical quantity measuring method | |
JP5541774B2 (en) | Physical quantity sensor and physical quantity measuring method | |
JP5545913B2 (en) | Counting device, distance meter, counting method, and distance measuring method | |
JP5596917B2 (en) | Physical quantity sensor and physical quantity measuring method | |
JP5612240B2 (en) | Distance meter and distance measuring method | |
JP5421577B2 (en) | Physical quantity sensor and physical quantity measuring method | |
WO2011111180A1 (en) | Physical quantity sensor and physical quantity measuring method | |
JP5426345B2 (en) | Vibration amplitude measuring apparatus and vibration amplitude measuring method | |
TWI359943B (en) | ||
JP2010117144A (en) | Physical quantity sensor and method for measuring physical quantity | |
JP5484661B2 (en) | Physical quantity sensor and physical quantity measuring method | |
JP2011106851A (en) | Physical quantity sensor and method for measuring physical quantity | |
EP2357490A1 (en) | Physical quantity sensor and physical quantity measuring method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C53 | Correction of patent of invention or patent application | ||
CB02 | Change of applicant information |
Address after: Japan Tokyo Marunouchi Chiyoda 2 chome 7 No. 3 Applicant after: Azbil Corporation Address before: Japan Tokyo Marunouchi Chiyoda 2 chome 7 No. 3 Applicant before: Yamatake Corp. |
|
COR | Change of bibliographic data |
Free format text: CORRECT: APPLICANT; FROM: YAMATAKE K. K. TO: AZBIL CORPORATION |
|
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20131127 Termination date: 20180304 |