Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP2915435B2 - Information processing apparatus and information processing method thereof - Google Patents
[go: Go Back, main page]

JP2915435B2 - Information processing apparatus and information processing method thereof - Google Patents

Information processing apparatus and information processing method thereof

Info

Publication number
JP2915435B2
JP2915435B2 JP17396189A JP17396189A JP2915435B2 JP 2915435 B2 JP2915435 B2 JP 2915435B2 JP 17396189 A JP17396189 A JP 17396189A JP 17396189 A JP17396189 A JP 17396189A JP 2915435 B2 JP2915435 B2 JP 2915435B2
Authority
JP
Japan
Prior art keywords
arrival time
information processing
phase
sensor
group
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
JP17396189A
Other languages
Japanese (ja)
Other versions
JPH0339610A (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.)
Canon Inc
Original Assignee
Canon Inc
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 Canon Inc filed Critical Canon Inc
Priority to JP17396189A priority Critical patent/JP2915435B2/en
Publication of JPH0339610A publication Critical patent/JPH0339610A/en
Application granted granted Critical
Publication of JP2915435B2 publication Critical patent/JP2915435B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Length Measuring Devices Characterised By Use Of Acoustic Means (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Description

【発明の詳細な説明】 [産業上の利用分野] 本発明は超音波距離計測装置、あるいは超音波の伝播
を利用する距離計測技術を取り入れた各種応用装置一般
の情報処理装置及びその情報処理方法に関するものであ
る。
DETAILED DESCRIPTION OF THE INVENTION [Industrial Application Field] The present invention relates to an ultrasonic distance measuring apparatus, or an information processing apparatus in general, and various information processing apparatuses incorporating a distance measuring technique using the propagation of ultrasonic waves. It is about.

[従来の技術] 超音波距離計測の応用装置として考案されたものの一
つに超音波座標入力装置がある。以下、この装置を例に
とつて本発明の説明を行う。
[Related Art] An ultrasonic coordinate input device is one of devices devised as an application device for ultrasonic distance measurement. Hereinafter, the present invention will be described using this apparatus as an example.

第8図は一般に考案されている超音波座標入力装置の
概略図である。1は座標指示具(以後ペンと呼ぶ)であ
り、内部に圧電素子が組み込まれており、先端より所望
の超音波信号を発信する発振具である。2a〜2cはセンス
であり、ペンから発せられた超音波信号を伝播体を介し
て受信する受信具である。3は伝播体で、超音波の伝播
媒体となり、ガラス,アルミ板等が用いられる。4は防
振材で、センサ2a〜2cにペン1からの波以外の反射波等
が混入しないようにする目的で設定されている。以上の
ような構成をとる超音波座標入力装置において、2個以
上の複数の各センサと指示点の距離より指示点の座標が
計算される。
FIG. 8 is a schematic view of a generally devised ultrasonic coordinate input device. Reference numeral 1 denotes a coordinate pointing device (hereinafter referred to as a pen), which has a piezoelectric element incorporated therein and is an oscillator that emits a desired ultrasonic signal from its tip. 2a to 2c are senses, which are receivers that receive an ultrasonic signal emitted from a pen via a propagator. Reference numeral 3 denotes a propagation body, which is a propagation medium of ultrasonic waves, and is made of glass, an aluminum plate, or the like. Reference numeral 4 denotes a vibration isolator, which is set for the purpose of preventing reflected waves other than the waves from the pen 1 from entering the sensors 2a to 2c. In the ultrasonic coordinate input device having the above-described configuration, the coordinates of the designated point are calculated from the distances between the two or more sensors and the designated point.

センサと指示点の距離算出方法はさまざまなものが考
案されているが、基本的にはペンから発せられた超音波
信号のセンサへの到達時間によつて算出する。群速度と
位相速度とが異なる横波超音波を用いる場合は、受信波
形に対しエンベローブをとつて群到達時間を求めて大ま
かな距離計測等を行い、さらに精度を求められる場合に
は、適当な位置の位相を検知して位相到達時間を求め、
より細かい精度を出している。
Various methods for calculating the distance between the sensor and the designated point have been devised, but basically, the distance is calculated based on the arrival time of the ultrasonic signal emitted from the pen to the sensor. When using transverse ultrasonic waves having different group velocities and phase velocities, the envelope of the received waveform is used to determine the group arrival time and rough distance measurement is performed. Phase arrival time by detecting the phase of
It has finer precision.

次に、群到達時間と位相到達時間との2つの時間情報
からペンとセンサ間距離rを算出する手順を説明する。
まず、第9図で示すような電圧をペン内の圧電素子に印
加した場合、センサの受信信号は第10図で示すようにな
る。これに対してそのピーク位置を群到達位置tgとして
検出する場合、第11図に示すようにセンサ出力信号に対
して初段増幅器5,全波整流器6,ローパスフイルタ7,微分
回路8の各回路を通してから、微分信号のゼロクロスを
コンパレータ9によつて検出し、その時刻を群到達時間
tgとして認識する。
Next, a procedure for calculating the distance r between the pen and the sensor from two pieces of time information of the group arrival time and the phase arrival time will be described.
First, when a voltage as shown in FIG. 9 is applied to the piezoelectric element in the pen, the received signal of the sensor becomes as shown in FIG. On the other hand, when the peak position is detected as the group arrival position tg, as shown in FIG. 11, the circuit of the first-stage amplifier 5, the full-wave rectifier 6, the low-pass filter 7, and the differentiating circuit 8 is applied to the sensor output signal. After that, the zero crossing of the differential signal is detected by the comparator 9 and the time is taken as the group arrival time.
Recognized as t g .

これにより、群速度をυとすれば、 r=υ・tg でr算出する事ができるが、エンベローブをもとに時間
を検出するという方法によるため、信号出力の大きさや
フイルタ特性の影響で、どうしてもある程度のゆらぎΔ
tが発生する。一般には特定の位相ゼロクロスポイント
を検出して時間を決定する方が、よりゆらぎの少ない値
を得ることができる。
Thus, if the group velocity and upsilon g, but can be r calculated by r = υ g · t g, because due to the method of detecting the time based on the envelope of the signal output size and the filter characteristics Inevitably, some fluctuation Δ
t occurs. Generally, detecting a specific phase zero crossing point and determining the time can obtain a value with less fluctuation.

そこで、群到達時間tg決定直後の位相ゼロクロスとい
うように検出ポイントを規定すると、υ≠υである
がために距離rと共に群の中の位相がずれるため、位相
到達時間tpとして第12図に示すような段階状のものが観
測される。この段は位相の検出ポイントの移動を示すも
ので、各段のつなぎ目は信号の周期Tだけ平行移動した
ものとなる。υとυが等しく常に一定の位相検出ポ
イントを観測できる場合には、このような段階はできず
直線aのような位相到達時間tpが得られる。
Therefore, since when defining the detection point as that phase zero cross just after the group arrival time t g determined, the phases of the group with upsilon g ≠ upsilon distance for is a p r deviates, the a phase arrival time t p Stepwise ones as shown in Fig. 12 are observed. This stage indicates the movement of the phase detection point, and the joint between the stages is a parallel movement of the signal period T. If υ g and υ p are equal and a constant phase detection point can always be observed, such a step cannot be performed and a phase arrival time t p as shown by a straight line a is obtained.

従つて、段階状に与えられる位相到達時間tpを元の直
線aに変換すればよい。
Accordance connexion, may be converted to a phase arrival time t p given to stepped to the original straight line a.

つまり、 tpa≒(υg)tg−tof (tof:オフセツト値) であるが、群到達時間tgはゆらぎが大きいので、 tpi=nT≒(υg)tg−tof−tp (nは整数) という性質を利用して、 tpa=tp+T×Int(tpi/T+0.5) =tp+T×Int[{(υg)tg −tof−tp}/T+0.5] としてtpaを求め、結局rを次式より計算する。That is, t pa ≒ (υ g / υ p ) t g −t of (t of : offset value). However, since the group arrival time t g has a large fluctuation, t pi = nT ≒ (υ g / υ g Using the property of t g −t of −t p (n is an integer), t pa = t p + T × Int (t pi /T+0.5)=t p + T × Int [{(υ g / υ p ) t g −t of −t p } /T+0.5], and t pa is obtained. Finally, r is calculated by the following equation.

r=υ・tpa−rof(rof:オフセツト値) =Atp+B×Int{Ctg+Dtp+E+0.5}+F A:υ B:υ・T=λ C:(υg)/T=f(υs) D:−1/T=−f E:−tof/T F:−rof [発明が解決しようとしている課題] しかしながら、以上のような方法でペンとセンサ間距
離rを精度よく算出するには、正確な計算定数(A,B,C,
D,E,F)を用いることが当然大前提となる。ところが、
このような計算定数は超音波の伝播を利用するという方
式の性質上、個々の製品の各種の物理定数,調整具合等
の影響で微妙な変動が起こることが必至であり、このた
めに安定した座標精度を確保した製品を大量生産するこ
とは事実上困難となつている。
r = υ p · t pa -r of (r of: offset value) = At p + B × Int {Ct g + Dt p + E + 0.5} + F A: υ p B: υ p · T = λ C: (υ g / υ p) / T = f (υ s / υ p) D: -1 / T = -f E: -t of / TF: -r of [ problems invention is trying to solve] However, as described above In order to calculate the distance r between the pen and the sensor with high accuracy,
D, E, F) is a major premise. However,
Due to the nature of the method of utilizing the propagation of ultrasonic waves, such calculation constants are inevitably subject to subtle variations due to various physical constants of individual products, adjustment conditions, and the like. It is practically difficult to mass-produce products with coordinate accuracy.

本発明は、前記従来の欠点を除去し、個々の装置の各
種の物理定数,調整具合等の影響による微妙な変動に関
わりなく、超音波の到達時間を基に指示点とセンサ間の
距離を正確に算出する情報処理装置及びその情報処理方
法を提供する。
The present invention eliminates the above-mentioned drawbacks of the related art, and makes it possible to determine the distance between the designated point and the sensor based on the arrival time of the ultrasonic wave, regardless of the subtle variations due to various physical constants and adjustment conditions of each device. An information processing apparatus and an information processing method for accurately calculating the information are provided.

[発明が解決しようとしている手段] この課題を解決するために、本発明の情報処理装置
は、超音波の到達時間を基に指示点とセンサ間の距離を
算出する情報処理装置であって、較正モードを判断する
判断手段と、前記判断手段により較正モードであると判
断された場合、任意に入力される連続的な複数の点での
超音波の到達時間から較正パラメータを算出する算出手
段と、指示点とセンサ間の距離を求めるために予め記憶
されている較正パラメータを、前記算出手段により算出
された較正パラメータに更新する更新手段とを有するこ
とを特徴とする。
[Means to be Solved by the Invention] In order to solve this problem, an information processing apparatus of the present invention is an information processing apparatus that calculates a distance between a designated point and a sensor based on an arrival time of an ultrasonic wave, Judgment means for judging the calibration mode, and calculation means for calculating a calibration parameter from arrival times of ultrasonic waves at a plurality of arbitrarily inputted continuous points when the judgment mode judges that the mode is the calibration mode. And updating means for updating a calibration parameter stored in advance to obtain the distance between the designated point and the sensor to the calibration parameter calculated by the calculating means.

ここで、前記算出手段は、予め定められた位置での群
到達時間tg及び位相到達時間tpを所定の演算式に代入し
て、発振源とセンサ間の距離の算出に用いる計算パラメ
ータの群速度υg,位相速度υp,周波数fを求める。ま
た、前記算出手段は、任意の連続的な複数位置で観測さ
れる群到達時間tg及び位相到達時間tpによる所定の統計
計算によつて、計算パラメータの群速度υg,位相速度υ
p,周波数fを求める。また、前記算出手段は、予め定め
られた適当な位置での群到達時間tg及び位相到達時間tp
を所定の演算式に代入して計算パラメータの位相速度υ
求め、任意の連続的な複数位置で得られる群到達時間
tg及び位相到達時間tpによる所定の統計計算によつて計
算パラメータの群速度υを求める。
Here, the calculating means, by substituting the group arrival time t g and the phase arrival time t p at a predetermined position on a predetermined arithmetic expression, the calculation parameters used to calculate the distance between the oscillation source and the sensor The group velocity υ g , the phase velocity υ p , and the frequency f are obtained. In addition, the calculating means performs a predetermined statistical calculation based on the group arrival time t g and the phase arrival time t p observed at arbitrary continuous positions, and calculates the group velocity { g , phase velocity} of the calculation parameter.
Find p and frequency f. Further, the calculating means calculates the group arrival time t g and the phase arrival time t p at a predetermined appropriate position.
Into a predetermined equation to calculate the phase velocity 計算
p arrival, group arrival time obtained at arbitrary continuous plural positions
Request group velocity upsilon g of calculated pursuant parameter to a predetermined statistical computation by t g and the phase arrival time t p.

又、本発明の情報処理装置は、超音波の到達時間を基
に指示点とセンサ間の距離を算出する情報処理装置であ
って、較正モードを判断する判断手段と、前記判断手段
により較正モードであると判断された場合、予め定めら
れた適当な位置での超音波の到達時間から設定された計
算パラメータを用いてセンサの取り付け位置を計算し
て、正規の位置からのズレ量を算出する算出手段と、指
示点とセンサ間の距離を求めるために予め記憶されてい
るセンサの位置補正量を、前記算出手段により算出され
たズレ量に更新する更新手段とを有することを特徴とす
る。
An information processing apparatus according to the present invention is an information processing apparatus that calculates a distance between an indicated point and a sensor based on an arrival time of an ultrasonic wave, wherein a determination unit that determines a calibration mode, and the calibration mode is determined by the determination unit. When it is determined that the sensor position is determined, the sensor mounting position is calculated using a calculation parameter set from the arrival time of the ultrasonic wave at a predetermined appropriate position, and the deviation amount from the normal position is calculated. It is characterized by comprising calculating means, and updating means for updating the position correction amount of the sensor stored in advance for obtaining the distance between the designated point and the sensor to the displacement amount calculated by the calculating means.

ここで、前記算出手段は、予め定められた適当な位置
での群到達時間tg及び位相到達時間tpを基に初期設定さ
れた計算パラメータの群速度υg,位相速度υp,周波数f
を用いてセンサの取り付け位置を計算して、正規の位置
からのズレ量を算出する。
Here, the calculating means, group velocity upsilon g of the initially set calculation parameters based on the group arrival time t g and the phase arrival time t p at an appropriate position predetermined phase velocity upsilon p, the frequency f
Is used to calculate the sensor mounting position, and the amount of deviation from the normal position is calculated.

又、本発明の情報処理方法は、超音波の到達時間を基
に指示点とセンサ間の距離を算出する情報処理装置での
情報処理方法であって、較正モードを有し、該較正モー
ド時には、任意に入力される連続的な複数の点での超音
波の到達時間から較正パラメータを算出し、指示点とセ
ンサ間の距離を求めるために予め記憶されている較正パ
ラメータを、前記算出された較正パラメータに更新する
ことを特徴とする。
Also, the information processing method of the present invention is an information processing method in an information processing apparatus that calculates a distance between an indicated point and a sensor based on an arrival time of an ultrasonic wave, and has a calibration mode. Calculate the calibration parameters from the arrival times of the ultrasonic waves at a plurality of points that are arbitrarily input, and calculate the calibration parameters stored in advance to determine the distance between the designated point and the sensor. It is characterized by updating to calibration parameters.

ここで、前記算出工程では、予め定められた位置での
群到達時間tg及び位相到達時間tpを所定の演算式に代入
して、発振源とセンサ間の距離の算出に用いる計算パラ
メータの群速度υg,位相速度υp,周波数fを求める。ま
た、前記算出工程では、任意の連続的な複数位置で観測
される群到達時間tg及び位相到達時間tpによる所定の統
計計算によつて、計算パラメータの群速度υg,位相速度
υp,周波数fを求める。また、前記算出工程では、予め
定められた適当な位置での群到達時間tg及び位相到達時
間tpを所定の演算式に代入して計算パラメータの位相速
度υ求め、任意の連続的な複数位置で得られる群到達
時間tg及び位相到達時間tpによる所定の統計計算によつ
て計算パラメータの群速度υを求める。
Here, in the calculating step, by substituting the group arrival time t g and the phase arrival time t p at a predetermined position on a predetermined arithmetic expression, the calculation parameters used to calculate the distance between the oscillation source and the sensor The group velocity υ g , the phase velocity υ p , and the frequency f are obtained. Further, in the calculation step, Yotsute a predetermined statistical computation by the group arrival time observed by any continuous multiple positions t g and the phase arrival time t p, the group velocity upsilon g of calculation parameters, phase velocity upsilon p , Frequency f. Further, in the calculating step calculates the phase velocity upsilon p calculation by substituting the group arrival time t g and the phase arrival time t p of a suitable predetermined positions on a predetermined arithmetic expression parameter, any continuous Request group velocity upsilon g of calculated pursuant parameter to a predetermined statistical computation by the time the group reached obtained at a plurality of positions t g and the phase arrival time t p.

又、本発明の情報処理方法は、超音波の到達時間を基
に指示点とセンサ間の距離を算出する情報処理装置での
情報処理方法であって、較正モードを有し、該較正モー
ド時には、予め定められた適当な位置での超音波の到達
時間から設定された計算パラメータを用いてセンサの取
り付け位置を計算して、正規の位置からのズレ量を算出
し、指示点とセンサ間の距離を求めるために予め記憶さ
れているセンサの位置補正量を、前記算出されたズレ量
に更新することを特徴とする。
Also, the information processing method of the present invention is an information processing method in an information processing apparatus that calculates a distance between an indicated point and a sensor based on an arrival time of an ultrasonic wave, and has a calibration mode. Calculate the mounting position of the sensor using the calculation parameters set from the arrival time of the ultrasonic wave at a predetermined appropriate position, calculate the amount of deviation from the normal position, between the designated point and the sensor It is characterized in that the position correction amount of the sensor stored in advance for obtaining the distance is updated to the calculated shift amount.

ここで、前記算出工程では、予め定められた適当な位
置での群到達時間tg及び位相到達時間tpを基に初期設定
された計算パラメータの群速度υg,位相速度υp,周波数
fを用いてセンサの取り付け位置を計算して、正規の位
置からのズレ量を算出する。
Here, in the calculating step, the group velocity upsilon g of the initially set calculation parameters based on the group arrival time t g and the phase arrival time t p at an appropriate position predetermined phase velocity upsilon p, the frequency f Is used to calculate the sensor mounting position, and the amount of deviation from the normal position is calculated.

[作用] かかる構成において、予め定められた適当な位置もし
くは任意の連続的な複数位置で観測される群到達時間
tg,位相到達時間tpを基に、計算パラメータを新たに設
定して較正できるようにし、個々の装置ごとに個別に最
適な計算パラメータを設定・較正できるようにしたもの
である。又、センサの取り付け位置の経年変化を補償で
きるようにしたものである。
[Operation] In such a configuration, the group arrival time observed at a predetermined appropriate position or any continuous plural positions.
Based on t g and phase arrival time t p , calculation parameters are newly set so that calibration can be performed, and optimum calculation parameters can be individually set and calibrated for each device. Further, it is possible to compensate for the secular change of the mounting position of the sensor.

[実施例] 以下、添付図面を参照して本発明の実施例を詳細に説
明する。
Hereinafter, an embodiment of the present invention will be described in detail with reference to the accompanying drawings.

<第1実施例> 第1図は本発明の第1実施例における較正用指示ポイ
ント及び超音波座標入力装置の構成例を示す図である。
図中、1は振動ペン、2はセンサ、3は振動板、4は防
振材であり、振動ペンはペン制御部10によりその制御を
され、センサ制御部20には第11図に示す各回路が含まれ
ていても良い。100は演算制御用のCPU、101は制御プロ
グラムを格納するROM、102は補助記憶用RAMで、例えば
本例の定数A,B,C,D,Eを記憶する。103は入力座標を表示
画面等に出力する出力部である。
<First Embodiment> FIG. 1 is a diagram showing a configuration example of a calibration instruction point and an ultrasonic coordinate input device according to a first embodiment of the present invention.
In the figure, reference numeral 1 denotes a vibration pen, 2 denotes a sensor, 3 denotes a vibration plate, and 4 denotes an anti-vibration material. The vibration pen is controlled by a pen control unit 10, and a sensor control unit 20 includes each of the components shown in FIG. A circuit may be included. 100 is a CPU for arithmetic control, 101 is a ROM for storing a control program, and 102 is a RAM for auxiliary storage, which stores, for example, constants A, B, C, D, and E in this example. An output unit 103 outputs the input coordinates to a display screen or the like.

前述したようにペンとセンサ間距離rは、 r=A tp+B Int{C tg+D tp+E+0.5}+F で求められるが、本式の第2項{ }内を整数化する項
であり、整数化によって整数値がずれてしまう程の物理
定数の変動が起きない限りは、予め決めた点における整
数値nに変動はないはずである。したがつて、例えばP
2,P3,P4の3点での(tg,tp)がわかれば、既知のnp2,n
p3,np4に対し C tg p2+D tp p2+E=np2 C tg p3+D tp p3+E=np3 C tg p4+D tp p4+E=np4 より、C,D,Eを求めることが可能となる。
Pen and the sensor distance r as described above is obtained by r = A t p + B Int {C t g + D t p + E + 0.5} + F, to integer within the second term of the equation {} term As long as the physical constant does not fluctuate such that the integer value shifts due to the conversion into an integer, the integer value n at a predetermined point should not fluctuate. Thus, for example, P
If (t g , t p ) at the three points 2,2, P3, P4 is known, the known n p2 , n
p3, to n p4 than C t g p2 + D t p p2 + E = n p2 C t g p3 + D t p p3 + E = n p3 C t g p4 + D t p p4 + E = n p4, C, D, and E It is possible to ask.

また、rに関しても同様に既知のnp2,np3,np4,rp2,r
p3,rp4に対し A tp p2+B np2+F=rp2 A tp p3+B np3+F=rp3 A tp p4+B np4+F=rp4 より、A,B,Fを求めることが可能となる。
Similarly, regarding r, the known n p2 , n p3 , n p4 , r p2 , r
p3, r p4 respect than A t p p2 + B n p2 + F = r p2 A t p p3 + B n p3 + F = r p3 A t p p4 + B n p4 + F = r p4, A, B, be determined F It becomes possible.

本実施例を行う上で注意することは、言うまでもなく
各較正用指示ポイントとr,nとが正確に対応しているこ
とである。そのため、較正用指示ポイントの場所は第3
図で示すような位相到達時間tp対rの関係のうちで段差
が生じる近辺からずらした場所に設定する必要がある。
尚、位相到達時間tpの段差は、おおむね ごとに発生する。
It should be noted that, in carrying out the present embodiment, it is needless to say that each calibration instruction point and r, n correspond exactly. Therefore, the location of the calibration instruction point is
It is necessary to set a location which is shifted from the vicinity of a step in among the relationship between the phase arrival time t p versus r as shown in FIG.
Incidentally, the step of phase arrival time t p is approximately Occurs every time.

例として、λ=3mm,ν=2350m/s,υ=1350m/sと
すれば、約7mmごとに位相到達時間tpの段差が生じるこ
とになるので、その近辺をうまくはずしてポイントを設
定する。
Examples, λ = 3mm, ν g = 2350m / s, if υ p = 1350m / s, it means that the step of phase arrival time t p occurs about every 7 mm, a point to remove the vicinity well Set.

以上の動作を実現するためのシステムの流れを第2図
のフローチヤートに示す。まず、何らかのスイツチ動作
によりシステムを較正モードにする。その後、ステツプ
S1で3点を入力して、ステツプS2,S3で計算パラメータ
を算出し、ステツプS4で定数の切り替えを行う。なお、
定数の切り替えは、バツクアツプされるRAM上での書き
変えや、何らかの手段でマイコンが読み込み可能なデイ
ツプ・スイツチ等の設定切り替えなどの手段をとればよ
い。
The flow of the system for realizing the above operation is shown in the flowchart of FIG. First, the system is put into the calibration mode by some switch operation. Then step
Three points are input in S1, calculation parameters are calculated in steps S2 and S3, and constants are switched in step S4. In addition,
The switching of the constants may be performed by means of rewriting on the RAM to be backed up, or by switching the setting of a dip switch which can be read by the microcomputer by some means.

<第2実施例> 次に本発明の第2実施例について説明する。Second Embodiment Next, a second embodiment of the present invention will be described.

r算出式の中でInt関数が正しい値を出していれば絶
対精度を出すものはA,B,Fである。ところで、A,B,Fは位
相速度υ及び周波数fに関与するパラメータであり、
これらはrと位相到達時間tpとの関係より求められる
が、一般に位相ゼロクロス検出により位相到達時間tp
求めると、第3図で示すようにrと位相到達時間tpの関
係においては各階段ごとに直線性がよく、またそれぞれ
の傾きυもほぼ一致する。また、階段ごとの切れ目が
示すT(=1/f)も全体にほぼ一様である。したがつ
て、これらA,B,Fは前記第1実施例で述べたような方法
で決定しても、おおむね真値に近い値が得られ実用上問
題ない。
A, B, and F are those that give absolute accuracy if the Int function gives a correct value in the r calculation formula. By the way, A, B, and F are parameters related to the phase velocity 及 びp and the frequency f,
These are determined from the relationship between r and the phase arrival time t p, but in general when obtaining the phase arrival time t p the phase zero-cross detection, in relation r and a phase arrival time t p, as shown in FIG. 3 each The linearity is good for each step, and the respective slopes そ れ ぞ れp are almost the same. In addition, T (= 1 / f) indicated by a break at each step is substantially uniform throughout. Therefore, even if these A, B, and F are determined by the method described in the first embodiment, a value close to the true value is obtained, and there is no practical problem.

しかしながら、Int関数内のC,D,Eは群速度υに関与
するパラメータであり、一般にr対群到達時間tgの関係
が第3図で示すように観測される以上、前記第1実施例
のようにこのうちの数点をとることによりC,D,Eを決定
することは大変危険である。つまり、第3図のようにr
対群到達時間tgの関係が得られるとき、例えば図中のよ
うな2点よりυを決定することが全く無意味であるこ
とに対応する。そこで、一般にC,D,Eを決定するには、
前記方法に加え何らかの統計処理をすることが必要とな
る。
However, C, D, and E in the Int function are parameters related to the group velocity υ p . Generally, the relationship between r and the group arrival time t g is observed as shown in FIG. It is very dangerous to determine C, D, and E by taking several points as in the example. That is, as shown in FIG.
When the relationship of the pair group arrival time t g is obtained, which corresponds to be determined upsilon p than two points as in FIG example is quite meaningless. Therefore, in general, to determine C, D, and E,
It is necessary to perform some statistical processing in addition to the above method.

そこで、本実施例においては次にような手段をとる。
まず、伝播体上の任意の場所で入力ペンにより第4図で
示すような任意の曲線を適当な速度で描き、(tg,tp
を連続入力する。すると、第5図で示すような(tg,
tp)対rの関係図が得られる。これに対し、開始点Aに
おける整数値の初期値nを適当に決め(例えば0)、そ
の後ペンがセンサから遠ざかる方向で位相到達時間tp
段差が生じた場合は+1、近づく方向では−1をして各
段差ごとに整数値nを対応させる。その結果、第6図で
示すようなデータ群ができる。これをもとにして独立変
数を群到達時間tg,位相到達時間tp,従属変数をnとした n=C tg+D tp+E′ (E′=|E|−Int{|E|+0.5}−nint) なる関係式による多重回帰を行い、最適定数C,D,Eを統
計的に決定する。
Therefore, the following means is taken in this embodiment.
First, an arbitrary curve as shown in FIG. 4 is drawn at an appropriate speed at an arbitrary position on the propagator using an input pen, and (t g , t p )
Is input continuously. Then, as shown in FIG. 5, (t g ,
t p ) vs. r. In contrast, decide appropriately the initial value n of the integer value at the start point A (e.g. 0), then if the pen has occurred stepped phase arrival time t p in a direction away from the sensor +1, in approaching direction -1 To associate an integer value n for each step. As a result, a data group as shown in FIG. 6 is created. Group arrival time independent variables it based on t g, the phase arrival time t p, the dependent variable was n n = C t g + D t p + E '(E' = | E | -Int {| E | +0.5} −n int ) Multiple regression is performed by the relational expression, and the optimal constants C, D, and E are statistically determined.

<第3実施例> 以上説明した実施例のいわば逆の適用例として、セン
サ位置のずれ量検出をあげることができる。前記第1,第
2実施例は、適当な位置を指示したときの(tg,tp)を
もとに距離計算の際の定数であるυgp,fを逆算し、R
AMバツクアツプなりデイツプ・スイツチなりの手段によ
り較正を行うというものであつた。ところが、υgp,
fは伝播体の厚みのばらつき,超音波発信ペンの周波数
特性の影響に大きく依存しており、個別の製品ごとに一
度設定を行えばその後は安定していると考えられ、逆に
使用により経時変化が出やすいのはセンサ位置である。
特に、製品自体が静的な安定した環境になく、頻繁に衝
撃を与えられるような場所にはなおさらである。従つ
て、適当な位置に指示したときの(tg,tp)をもとに初
期設定されたυgp,fを用いてセンサ位置を逆算し、
その結果を距離計算式にフイードバツクする方法をとる
ことにより、センサ位置のずれ補正を行うことができ
る。
<Third Embodiment> As a reverse application example of the above-described embodiment, detection of a shift amount of a sensor position can be mentioned. Said first, second embodiment, and calculated back (t g, t p) to a constant when the distance calculated on the basis υ g, υ p, f when instructed appropriate position, R
The calibration was performed by means of an AM backup or a depth switch. However, υ g , υ p ,
f greatly depends on the variation of the thickness of the propagation body and the influence of the frequency characteristics of the ultrasonic transmission pen, and once it is set for each individual product, it is considered to be stable thereafter. It is the sensor position that is likely to change.
This is especially true in locations where the product itself is not in a static, stable environment and can be impacted frequently. Therefore, the sensor position is calculated backward using υ g , υ p , f which is initially set based on (t g , t p ) when the appropriate position is indicated,
By using a method of feeding back the result to a distance calculation formula, it is possible to correct the displacement of the sensor position.

以下、第7図を使つて説明する。まず、何らかのスイ
ツチ動作によりシステムの動作モードを較正モードにす
ると、ステツプS75からステツプS76に進む。ここで、絶
対座標x軸上のP1,P2点をポイントする。その際のP1
センサS2間,P2とセンサS2間の距離を、予め初期設定さ
れたυgp,fを用いて、ステツプS77で取り込んだ(tg
p1,tp p1),(tg p2,tp p2)をもとにステツプS78で計算
し、その結果をもとにセンサS2の絶対座標のx成分を算
出する。同様に絶対座標y軸上のP6,P7点をポイント
し、その際の(tg P6,tp P6),(tg P7,tp P7)をもとに、
センサS2の絶対座標のy成分を算出する。このようにし
て、全てのセンサの絶対座標を求め、本来あるべき位置
からのずれ量を読みとり、ステツプS79でこの結果を通
常の動作モードにおける座標計算式にフイードバツクす
る。
Hereinafter, description will be made with reference to FIG. First, when the operation mode of the system is set to the calibration mode by some switch operation, the process proceeds from step S75 to step S76. Here, points P 1 and P 2 on the absolute coordinate x axis are pointed. Between P 1 and the sensor S 2 at that time, the distance between P 2 and the sensor S 2, advance initially set upsilon g, using upsilon p, f, taken at step S77 (t g
p1, t p p1), and calculates the (t g p2, t p p2 ) calculated in step S78 on the basis of, x component of the absolute coordinate of the sensor S 2 based on the results. Similarly, point P 6 and P 7 on the absolute coordinate y axis, and based on (t g P6 , t p P6 ) and (t g P7 , t p P7 ) at that time,
Calculating the y component of the absolute coordinate of the sensor S 2. In this way, the absolute coordinates of all the sensors are obtained, the deviation from the original position is read, and the result is fed back to the coordinate calculation formula in the normal operation mode in step S79.

通常モードでは、ステツプS72でペンとセンサ間距離
算出、ステツプS73で指示座標算出の後、ステツプS74で
ステツプS79で設定されたセンサずれ補正量による補正
を行う。
In the normal mode, the distance between the pen and the sensor is calculated in step S72, the designated coordinate is calculated in step S73, and then the correction is performed in step S74 based on the sensor shift correction amount set in step S79.

[発明の効果] 本発明により、個々の装置の各種の物理定数,調整具
合等の影響による微妙な変動に関わりなく、超音波の到
達時間を基に指示点とセンサ間の距離を正確に算出する
情報処理装置及びその情報処理方法を提供できる。
[Effects of the Invention] According to the present invention, the distance between the designated point and the sensor can be accurately calculated based on the arrival time of the ultrasonic wave, regardless of subtle variations due to various physical constants and adjustment conditions of each device. And an information processing method thereof.

すなわち、予め定められた適当な位置もしくは任意の
連続的な複数位置で観測される群到達時間tg,位相到達
時間tpを基に、係数定数を新たに設定して較正できるよ
うなシステム構成をとることにより、個々の製品ごとに
個別に最適な計算定数を設定・較正できるようになり、
製品を構成する各部材の物理定数のバラツキ,取り付け
具合のバラツキまた経時変化などによる最適計算定数の
変動を工場出荷時のみならず、サービスマン,ユーザら
が市場で較正を行うことで簡便に吸収することが可能と
なり、安定した精度を確保した超音波距離計測の応用製
品が実用化する。
That is, predetermined appropriate position or any successive plurality of positions in the observed group arrival time t g, based on the phase arrival time t p, the system configuration can be calibrated by setting the constant coefficient newly By taking, the optimal calculation constant can be set and calibrated individually for each product,
Variations in the physical constants of the components that make up the product, variations in the degree of mounting, and changes in the optimal calculation constants due to aging, etc., are easily absorbed not only at the time of factory shipment but also by servicemen and users calibrating in the market. And practical application of ultrasonic distance measurement application products that ensure stable accuracy.

【図面の簡単な説明】[Brief description of the drawings]

第1図は本実施例較正用指示ポイント及び超音波座標入
力装置の構成を示す図、 第2図は第1実施例における処理手順を示すフローチヤ
ート図、 第3図は距離rと群到達時間tg,位相到達時間tpとの関
係を示す図、 第4図は第2実施例における較正用入力例を示す図、 第5図は第2実施例におけるデータ番号と群到達時間
tg,位相到達時間tpの関係例を示す図、 第6図は第2実施例における較正用データ群を示す図、 第7図は第3実施例における処理手順を示すフローチヤ
ート、 第8図は座標入力装置の概略図、 第9図は指示用ペンの駆動電圧を示す図、 第10図はセンサ出力波形を示す図、 第11図は信号処理回路の構成例を示す図、 第12図は距離計算を説明する図である。 図中、1……座標指示具、2,2a〜2c……センサ、3……
伝播体、4……防振材、5……初段増幅器、6……全波
整流器、7……ローパスフイルタ、8……微分回路、9
……コンパレータ、10……ペン制御部、20……センサ制
御部、100……CPU、101……ROM、102……RAM、103……
出力部である。
FIG. 1 is a diagram showing the configuration of an instruction point for calibration and an ultrasonic coordinate input device in this embodiment. FIG. 2 is a flowchart showing a processing procedure in the first embodiment. FIG. 3 is a distance r and a group arrival time. t g, shows a relationship between the phase arrival time t p, Fig. 4 shows a calibration input example of the second embodiment, Figure 5 is a data number and the group arrival time in the second embodiment
t g, diagram illustrating an example of the relationship between the phase arrival time t p, Fig. 6 is a diagram showing the calibration data group in the second embodiment, Figure 7 is flow chart showing a processing procedure in the third embodiment, eighth FIG. 9 is a schematic diagram of a coordinate input device, FIG. 9 is a diagram showing a driving voltage of the pointing pen, FIG. 10 is a diagram showing a sensor output waveform, FIG. 11 is a diagram showing a configuration example of a signal processing circuit, FIG. The figure illustrates the distance calculation. In the figure, 1 ... coordinate indicating tool, 2,2a ~ 2c ... sensor, 3 ...
Propagator, 4… Vibration isolator, 5… First-stage amplifier, 6… Full-wave rectifier, 7… Low-pass filter, 8… Differentiation circuit, 9
...... Comparator, 10 ... Pen control unit, 20 ... Sensor control unit, 100 ... CPU, 101 ... ROM, 102 ... RAM, 103 ...
Output unit.

───────────────────────────────────────────────────── フロントページの続き (58)調査した分野(Int.Cl.6,DB名) G01S 5/22 G01B 17/00 G06F 3/03 ──────────────────────────────────────────────────続 き Continued on the front page (58) Fields investigated (Int.Cl. 6 , DB name) G01S 5/22 G01B 17/00 G06F 3/03

Claims (12)

(57)【特許請求の範囲】(57) [Claims] 【請求項1】超音波の到達時間を基に指示点とセンサ間
の距離を算出する情報処理装置であって、 較正モードを判断する判断手段と、 前記判断手段により較正モードであると判断された場
合、任意に入力される連続的な複数の点での超音波の到
達時間から較正パラメータを算出する算出手段と、 指示点とセンサ間の距離を求めるために予め記憶されて
いる較正パラメータを、前記算出手段により算出された
較正パラメータに更新する更新手段とを有することを特
徴とする情報処理装置。
1. An information processing apparatus for calculating a distance between an indicated point and a sensor based on an arrival time of an ultrasonic wave, comprising: a determination unit for determining a calibration mode; and the determination unit determines that the mode is the calibration mode. Calculation means for calculating a calibration parameter from arrival times of ultrasonic waves at a plurality of points that are arbitrarily input; and a calibration parameter stored in advance to obtain a distance between the designated point and the sensor. And an updating unit for updating the calibration parameters calculated by the calculating unit.
【請求項2】前記算出手段は、予め定められた位置での
群到達時間tg及び位相到達時間tpを所定の演算式に代入
して、発振源とセンサ間の距離の算出に用いる計算パラ
メータの群速度υg,位相速度υp,周波数fを求めること
を特徴とする請求項1記載の情報処理装置。
Wherein said calculation means substitutes a group arrival time t g and the phase arrival time t p at a predetermined position in the predetermined arithmetic expression, calculation used to calculate the distance between the oscillation source and the sensor 2. The information processing apparatus according to claim 1, wherein a group velocity υ g , a phase velocity υ p , and a frequency f of parameters are obtained.
【請求項3】前記算出手段は、任意の連続的な複数位置
で観測される群到達時間tg及び位相到達時間tpによる所
定の統計計算によつて、計算パラメータの群速度υg,位
相速度υp,周波数fを求めることを特徴とする請求項1
記載の情報処理装置。
Wherein said calculating means Yotsute a predetermined statistical computation by the group arrival time observed by any continuous multiple positions t g and the phase arrival time t p, the group velocity upsilon g of calculation parameters, phase 2. The method according to claim 1, wherein the speed υ p and the frequency f are obtained.
An information processing apparatus according to claim 1.
【請求項4】前記算出手段は、予め定められた適当な位
置での群到達時間tg及び位相到達時間tpを所定の演算式
に代入して計算パラメータの位相速度υ求め、任意の
連続的な複数位置で得られる群到達時間tg及び位相到達
時間tpによる所定の統計計算によつて計算パラメータの
群速度υを求めることを特徴とする請求項1記載の情
報処理装置。
Wherein said calculating means obtains the phase velocity upsilon p of calculation parameters by substituting the group arrival time t g and the phase arrival time t p of a suitable predetermined positions on a predetermined arithmetic expression, any the information processing apparatus according to claim 1, wherein the determination of the group velocity upsilon g continuous calculated pursuant parameter to a predetermined statistical computation by the group arrival time obtained t g and the phase arrival time t p at a plurality of positions.
【請求項5】超音波の到達時間を基に指示点とセンサ間
の距離を算出する情報処理装置であって、 較正モードを判断する判断手段と、 前記判断手段により較正モードであると判断された場
合、予め定められた適当な位置での超音波の到達時間か
ら設定された計算パラメータを用いてセンサの取り付け
位置を計算して、正規の位置からのズレ量を算出する算
出手段と、 指示点とセンサ間の距離を求めるために予め記憶されて
いるセンサの位置補正量を、前記算出手段により算出さ
れたズレ量に更新する更新手段とを有することを特徴と
する情報処理装置。
5. An information processing apparatus for calculating a distance between an indicated point and a sensor based on an arrival time of an ultrasonic wave, comprising: a judging unit for judging a calibration mode; and the judging unit judges the calibration mode. Calculating means for calculating a mounting position of the sensor using a calculation parameter set from an arrival time of the ultrasonic wave at a predetermined appropriate position and calculating a deviation amount from a normal position; An information processing apparatus, comprising: updating means for updating a position correction amount of a sensor stored in advance for obtaining a distance between a point and a sensor to a displacement amount calculated by the calculating means.
【請求項6】前記算出手段は、予め定められた適当な位
置での群到達時間tg及び位相到達時間tpを基に初期設定
された計算パラメータの群速度υg,位相速度υp,周波数
fを用いてセンサの取り付け位置を計算して、正規の位
置からのズレ量を算出することを特徴とする請求項5記
載の情報処理装置。
Wherein said calculating means, group velocity upsilon g of the initially set calculation parameters based on the group arrival time t g and the phase arrival time t p of a suitable predetermined positions, phase velocity upsilon p, 6. The information processing apparatus according to claim 5, wherein the sensor mounting position is calculated using the frequency f, and a shift amount from a normal position is calculated.
【請求項7】超音波の到達時間を基に指示点とセンサ間
の距離を算出する情報処理装置での情報処理方法であっ
て、 較正モードを有し、 該較正モード時には、任意に入力される連続的な複数の
点での超音波の到達時間から較正パラメータを算出し、 指示点とセンサ間の距離を求めるために予め記憶されて
いる較正パラメータを、前記算出された較正パラメータ
に更新することを特徴とする情報処理方法。
7. An information processing method in an information processing apparatus for calculating a distance between a designated point and a sensor based on an arrival time of an ultrasonic wave, the method including a calibration mode. The calibration parameters are calculated from the arrival times of the ultrasonic waves at a plurality of continuous points, and the calibration parameters stored in advance for obtaining the distance between the designated point and the sensor are updated to the calculated calibration parameters. An information processing method, comprising:
【請求項8】前記算出工程では、予め定められた位置で
の群到達時間tg及び位相到達時間tpを所定の演算式に代
入して、発振源とセンサ間の距離の算出に用いる計算パ
ラメータの群速度υg,位相速度υp,周波数fを求めるこ
とを特徴とする請求項7記載の情報処理方法。
The method according to claim 8, wherein the calculating step, by substituting the group arrival time t g and the phase arrival time t p at a predetermined position in the predetermined arithmetic expression, calculation used to calculate the distance between the oscillation source and the sensor 8. The information processing method according to claim 7, wherein a group velocity υ g , a phase velocity υ p , and a frequency f of parameters are obtained.
【請求項9】前記算出工程では、任意の連続的な複数位
置で観測される群到達時間tg及び位相到達時間tpによる
所定の統計計算によつて、計算パラメータの群速度υg,
位相速度υp,周波数fを求めることを特徴とする請求項
7記載の情報処理方法。
9. In the calculating step, the group velocity υ g , the group velocity 計算g , of the calculation parameter is obtained by a predetermined statistical calculation based on the group arrival time t g and the phase arrival time t p observed at arbitrary continuous positions.
8. The information processing method according to claim 7, wherein the phase velocity υ p and the frequency f are obtained.
【請求項10】前記算出工程では、予め定められた適当
な位置での群到達時間tg及び位相到達時間tpを所定の演
算式に代入して計算パラメータの位相速度υ求め、任
意の連続的な複数位置で得られる群到達時間tg及び位相
到達時間tpによる所定の統計計算によつて計算パラメー
タの群速度υを求めることを特徴とする請求項7記載
の情報処理方法。
The method according to claim 10, wherein the calculating step calculates the phase velocity upsilon p of calculation parameters by substituting the group arrival time t g and the phase arrival time t p of a suitable predetermined positions on a predetermined arithmetic expression, any the information processing method according to claim 7, wherein the determination of the group velocity upsilon g continuous calculated pursuant parameter to a predetermined statistical computation by the group arrival time obtained t g and the phase arrival time t p at a plurality of positions.
【請求項11】超音波の到達時間を基に指示点とセンサ
間の距離を算出する情報処理装置での情報処理方法であ
って、 較正モードを有し、 該較正モード時には、予め定められた適当な位置での超
音波の到達時間から設定された計算パラメータを用いて
センサの取り付け位置を計算して、正規の位置からのズ
レ量を算出し、 指示点とセンサ間の距離を求めるために予め記憶されて
いるセンサの位置補正量を、前記算出されたズレ量に更
新することを特徴とする情報処理方法。
11. An information processing method in an information processing apparatus for calculating a distance between an indicated point and a sensor based on an arrival time of an ultrasonic wave, comprising a calibration mode, wherein a predetermined mode is set in the calibration mode. Using the calculation parameters set from the arrival time of the ultrasonic wave at the appropriate position, calculate the sensor mounting position, calculate the amount of deviation from the regular position, and obtain the distance between the indicated point and the sensor An information processing method comprising: updating a position correction amount of a sensor stored in advance to the calculated shift amount.
【請求項12】前記算出工程では、予め定められた適当
な位置での群到達時間tg及び位相到達時間tpを基に初期
設定された計算パラメータの群速度υg,位相速度υp,周
波数fを用いてセンサの取り付け位置を計算して、正規
の位置からのズレ量を算出することを特徴とする請求項
11記載の情報処理方法。
The method according to claim 12, wherein the calculating step, the group velocity upsilon g of the initially set calculation parameters based on the group arrival time t g and the phase arrival time t p of a suitable predetermined positions, phase velocity upsilon p, 2. The method according to claim 1, further comprising calculating a mounting position of the sensor using the frequency f, and calculating an amount of deviation from a normal position.
11. The information processing method according to 11.
JP17396189A 1989-07-07 1989-07-07 Information processing apparatus and information processing method thereof Expired - Fee Related JP2915435B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP17396189A JP2915435B2 (en) 1989-07-07 1989-07-07 Information processing apparatus and information processing method thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP17396189A JP2915435B2 (en) 1989-07-07 1989-07-07 Information processing apparatus and information processing method thereof

Publications (2)

Publication Number Publication Date
JPH0339610A JPH0339610A (en) 1991-02-20
JP2915435B2 true JP2915435B2 (en) 1999-07-05

Family

ID=15970253

Family Applications (1)

Application Number Title Priority Date Filing Date
JP17396189A Expired - Fee Related JP2915435B2 (en) 1989-07-07 1989-07-07 Information processing apparatus and information processing method thereof

Country Status (1)

Country Link
JP (1) JP2915435B2 (en)

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS62281023A (en) * 1986-05-30 1987-12-05 Canon Inc coordinate input device
JPS63100531A (en) * 1986-10-17 1988-05-02 Canon Inc input device

Also Published As

Publication number Publication date
JPH0339610A (en) 1991-02-20

Similar Documents

Publication Publication Date Title
EP0467296B1 (en) Coordinates input apparatus
US5239138A (en) Coordinate input apparatus
KR100259693B1 (en) A coordinate input device and a control method therefor
EP0584695A3 (en) Coordinate input apparatus and method, and oscillation generating apparatus
JP2915435B2 (en) Information processing apparatus and information processing method thereof
JP3115307B2 (en) Coordinate input device and coordinate detection method
JP3244895B2 (en) Distance measuring device, coordinate input device using the same, distance measuring method, and coordinate input method using the same
JP3010398B2 (en) Steel sheet shape measurement method
JPH03243840A (en) Temperature compensation method of piezoelectric vibrator pressure sensor
JP3401374B2 (en) Sound velocity correction method for underwater position measurement system
RU2221991C1 (en) Method establishing positions of mobile ground objects and facility for its realization
JP3077989B2 (en) Ultrasonic coordinate input method and device
US12460961B2 (en) Surface acoustic wave scale
JP2558104B2 (en) Method and device for adjusting sensor gain
JPH109852A (en) Straightness measurement method
JP3034271B2 (en) Ultrasonic distance measuring method and apparatus
JPH08313248A (en) Apparatus and method for probe-type measurement of surface roughness
JPH0318917A (en) Coordinate input device
JP2953285B2 (en) Pointing device
JPH03148040A (en) Viscosity measuring apparatus
JP2537537B2 (en) Coordinate input device
SU1270582A1 (en) Method of measuring velocity and attenuation of acoustic waves in material
JPH03192421A (en) Ultrasonic distance measuring instrument
SU643843A1 (en) Device for determining dynamic characteristics of oscillating systems
JPS61138131A (en) Crystal thermometer

Legal Events

Date Code Title Description
LAPS Cancellation because of no payment of annual fees