JP3092100B2 - Position detector with error correction function - Google Patents
Position detector with error correction functionInfo
- Publication number
- JP3092100B2 JP3092100B2 JP04089391A JP8939192A JP3092100B2 JP 3092100 B2 JP3092100 B2 JP 3092100B2 JP 04089391 A JP04089391 A JP 04089391A JP 8939192 A JP8939192 A JP 8939192A JP 3092100 B2 JP3092100 B2 JP 3092100B2
- Authority
- JP
- Japan
- Prior art keywords
- equation
- sine wave
- phase
- error
- time
- 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
Links
Landscapes
- Length Measuring Devices With Unspecified Measuring Means (AREA)
- Transmission And Conversion Of Sensor Element Output (AREA)
Description
【0001】[0001]
【産業上の利用分野】本発明は、サーボモータの回転角
度の検出等に利用される位置検出装置の改良に関するも
のであり、移動体の位置を検出する際に、その検出信号
である2相交流信号に含まれる各種の誤差成分を検出
し、自動的にこの誤差を補正して位置を検出するように
した誤差補正機能付き位置検出装置に関するものであ
る。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an improvement of a position detecting device used for detecting a rotation angle of a servomotor and the like. The present invention relates to a position detection device with an error correction function that detects various error components included in an AC signal and automatically corrects the error to detect a position.
【0002】[0002]
【従来の技術】従来の位置検出装置の一例を図5に示
す。被検出体の回転体1には透磁材で成るギア2が取付
けられており、ギア2に対向して磁気センサ3,4及び
磁石5が配設されており、磁気センサ3及び4の出力V
x及びVyはそれそれ瞬時値検出部10内のサンプルホ
ールド回路11及び12に入力されると共に、位置上位
桁検出回路20内の比較演算器21及び22に入力され
る。サンプルホールド回路11及び12の出力はそれぞ
れA/D変換器13及び14に入力され、A/D変換さ
れたディジタル瞬時値x=cosθ及びy=sinθは
除算器31に入力され、その除算結果y/xがtan
−1回路32に入力されている。このtan−1回路3
2の出力は、ギア2の1歯に相当する回転角を2π(r
ad)とした場合の1歯内の微小な位置θを内挿検出し
た結果である。比較演算器21及び22の出力はそれぞ
れパルスカウンタ23に入力され、計数された位置の上
位桁データTNは位置検出部30内の加算器33に入力
される。加算器33では上述の1歯内の内挿検出した位
置θと、位置の上位桁データTNとを加算することによ
り、回転体1の回転角度位置PSを得ている。このよう
な構成において、位置検出部30における処理は、マイ
クロプロセッサ等を用いたソフトウェア処理で実現する
ことが考えられる。この場合のソフトウェアの動作は、
図6のフローチャートに示されるようなものである。2. Description of the Related Art An example of a conventional position detecting device is shown in FIG. A gear 2 made of a magnetically permeable material is attached to a rotating body 1 of the detected object, and magnetic sensors 3 and 4 and a magnet 5 are disposed to face the gear 2, and outputs of the magnetic sensors 3 and 4 are provided. V
x and Vy are respectively input to the sample and hold circuits 11 and 12 in the instantaneous value detection unit 10 and also to the comparison arithmetic units 21 and 22 in the position upper digit detection circuit 20. The outputs of the sample and hold circuits 11 and 12 are input to A / D converters 13 and 14, respectively. The digital instantaneous values x = cos θ and y = sin θ that have been A / D converted are input to a divider 31 and the division result y / X is tan
-1 circuit 32. This tan -1 circuit 3
2 indicates that the rotation angle corresponding to one tooth of the gear 2 is 2π (r
This is a result of interpolation detection of a minute position θ in one tooth in the case of (ad). The outputs of the comparators 21 and 22 are input to the pulse counter 23, and the upper digit data TN of the counted position is input to the adder 33 in the position detector 30. The adder 33 obtains the rotation angle position PS of the rotating body 1 by adding the position θ detected by interpolation within one tooth and the upper digit data TN of the position. In such a configuration, the processing in the position detection unit 30 may be realized by software processing using a microprocessor or the like. The operation of the software in this case is
This is as shown in the flowchart of FIG.
【0003】[0003]
【発明が解決しようとする課題】従来の位置検出装置で
は、磁気センサ3及び4の出力Vx及びVyは互いに振
幅が等しく、位相差が正確に90°異なっており、又直
流オフセット成分を0として扱っている。ただし、一般
的には直流オフセットについては、0でなくても予め既
知である一定の値として扱っている例も多く、この場合
でも位置検出部30の内部処理において、既知である一
定の直流オフセットをVx及びVyから減算することに
よって、前述のように直流オフセットが0である場合と
同様に扱っている。ところが実際の位置検出装置におい
ては、磁気センサ3及び4により出力される実際の2相
交流信号について、オフセット又はオフセットの変動が
生じていたり、各振幅が異なっていたり、位相差が90
°からずれていたりするため、その信号は検出されるべ
きV・sinθ,V・cosθとは異なった値として検
出されていた。このような誤差を含んだ検出値に基づい
て演算を行なっていたので、従来は高精度な位置検出が
できないという問題点があった。In the conventional position detecting device, the outputs Vx and Vy of the magnetic sensors 3 and 4 have the same amplitude, the phase difference is exactly 90 °, and the DC offset component is 0. Are dealing. However, in general, the DC offset is often treated as a fixed value that is known in advance even if it is not 0. Even in this case, in the internal processing of the position detection unit 30, the fixed DC offset that is known is also used. Is subtracted from Vx and Vy, so that the DC offset is treated as 0 as described above. However, in an actual position detecting device, an offset or a fluctuation of an offset occurs in an actual two-phase AC signal output from the magnetic sensors 3 and 4, each amplitude is different, and a phase difference is 90%.
°, the signal was detected as a value different from V · sin θ and V · cos θ to be detected. Since the calculation is performed on the basis of the detection value including such an error, there has been a problem that the position cannot be detected with high accuracy conventionally.
【0004】尚、従来そのような場合に、上記各誤差補
正回路を設けることにより対処することもあったが、各
誤差補正回路は可変抵抗等の調整が必要であったため、
人為的な調整ミスが生じたり、回路素子の温度変動等に
よる特性変化や経年変化による特性劣化等により、常に
高精度な位置検出を行なうことができなかった。又、本
出願人は特開平3−170011号において、上記各誤
差を2相交流信号Vx,Vyの複数組の瞬時値から検出
する誤差検出方法を提案しているが、その演算は非常に
複雑であって、実際の位置検出装置においては演算回路
が高価になること、又複数組の瞬時値の選定をどのよう
に行なうかによって所望の検出精度が得られない等の課
題を有していた。Conventionally, in such a case, there has been a case where the above error correction circuits are provided to cope with the problem. However, since each error correction circuit requires adjustment of a variable resistor or the like,
Due to artificial adjustment errors, characteristic changes due to circuit element temperature fluctuations and the like, and characteristic deterioration due to aging, it has not been possible to always perform highly accurate position detection. Further, the applicant of the present invention has proposed in Japanese Patent Application Laid-Open No. Hei 3-170011 an error detection method for detecting each of the above errors from a plurality of sets of instantaneous values of the two-phase AC signals Vx and Vy. However, in an actual position detecting device, there are problems that an arithmetic circuit becomes expensive and that a desired detection accuracy cannot be obtained depending on how to select a plurality of sets of instantaneous values. .
【0005】本発明は上述したような事情から成された
ものであり、本発明の目的は、上記各誤差を簡単かつ安
価な演算処理回路により同時に演算して正確に検出し、
さらに自動的に補正することによって常に正確な位置検
出を行なうことのできる誤差補正機能付き位置検出装置
を提供することにある。The present invention has been made in view of the above circumstances, and an object of the present invention is to simultaneously and accurately detect each of the above errors by a simple and inexpensive arithmetic processing circuit.
It is still another object of the present invention to provide a position detection device with an error correction function that can always perform accurate position detection by automatically correcting the position.
【0006】[0006]
【課題を解決するための手段】本発明は、誤差補正機能
付き位置検出装置に関するもので、本発明の前記目的
は、移動体の移動距離θに応じて、互いに位相の異なる
第1及び第2の2相正弦波信号x及びyを出力する検出
部と;前記第1の正弦波信号をx=1/A・(cosθ
−K・sinθ)+C、前記第2の正弦波信号をy=s
inθ+dと表わしたとき、x=0となる時点における
yの瞬時値y1及びy2と、y=0となる時点における
xの瞬時値x1及びx2と、|x|=|y|となる時点
におけるx又はyの瞬時値x3及びx4、又はy3及び
y4とを記憶する記憶手段と;前記記憶手段に記憶され
ているy1,y2,x1,x2,x3又はy3,x4又
はy4を用いて、前記2相正弦波信号x及びyの振幅の
差を表わすAを数1の演算により求め、前記2相正弦波
信号x及びyの位相の差を表わすKを数2の演算により
求め、前記2相正弦波信号x及びyの直流オフセットを
表わすC及びdを数3及び数4の演算により求める演算
手段と;前記誤差演算手段によって同時に求めたC,
d,A,Kを用いて、前記2相正弦波信号x及びyか
ら、互いに位相が90°異なり振幅が等しく、かつ直流
オフセット成分のない第3及び第4の2相正弦波信号
x′=cosθ及びy′=sinθを求める誤差補正手
段と;前記第3及び第4の2相正弦波信号x′及びy′
からθ=tan−1(y′/x′)の演算に基づいて前
記移動体の移動距離θを求める位置検出部とを具備した
ことにより達成される。そして、前記記憶手段は、前記
移動体の移動速度が所定値以上の場合には前記y1,y
2,x1,x2,x3,x4を更新することなく保持し
ており、前記移動体の移動速度が前記所定値以下の場合
には前記y1,y2,x1,x2,x3,x4をそれぞ
れx=0となる時点、y=0となる時点、|x|=|y
|となる時点において順次更新して記憶するようになっ
ている。SUMMARY OF THE INVENTION The present invention relates to a position detecting device having an error correcting function, and an object of the present invention is to provide a position detecting device having first and second phases different in phase according to a moving distance .theta. A detection unit that outputs two-phase sine wave signals x and y; x = 1 / A · (cos θ
−K · sin θ) + C, and the second sine wave signal is represented by y = s
When expressed as inθ + d, instantaneous values y 1 and y 2 of y at the time when x = 0, instantaneous values x 1 and x 2 of x at the time of y = 0, and | x | = | y | x or y instantaneous value x 3 and x 4 in at the time made, or y 3 and y 4 and storage means for storing a; y 1 stored in the storage means, y 2, x 1, x 2, x 3 or y 3 , x 4 or y 4 to calculate the amplitude of the two-phase sinusoidal signals x and y .
A representing the difference is obtained by the operation of Equation 1, and the two-phase sine wave is obtained.
K representing the phase difference between the signals x and y is calculated by the equation (2).
And the DC offset of the two-phase sinusoidal signals x and y
An operation for finding C and d to be expressed by the operations of Expressions 3 and 4.
Means ; C obtained simultaneously by the error calculating means,
Using d, A, and K, the third and fourth two-phase sine wave signals x '= 90 ° differing in phase from each other, have the same amplitude, and have no DC offset component from the two-phase sine wave signals x and y. error correcting means for obtaining cos θ and y ′ = sin θ; the third and fourth two-phase sine wave signals x ′ and y ′
And a position detection unit for calculating the moving distance θ of the moving body based on the calculation of θ = tan −1 (y ′ / x ′). The storage means stores the y 1 , y when the moving speed of the moving body is equal to or higher than a predetermined value.
2, x 1, x 2, x 3, x 4 holds without updating the said when the moving speed of the moving body is equal to or less than the predetermined value is the y 1, y 2, x 1 , x 2 , X 3 , x 4 at the time when x = 0, y = 0, | x | = | y
At the time of |, the data is sequentially updated and stored.
【0007】[0007]
【作用】本発明においては、移動体の移動距離θに応じ
て検出部より出力される2相の正弦波信号Vx及びVy
は互いの振幅に誤差があり、位相差も90°に対して誤
差があり、かつそれぞれに直流オフセット成分を含んで
いるため、According to the present invention, the two-phase sine wave signals Vx and Vy output from the detector in accordance with the moving distance θ of the moving body.
Since there is an error in the amplitude of each other, the phase difference also has an error with respect to 90 °, and each contains a DC offset component,
【数5】Vx=V1・(cosθ−K・sinθ)+aVx = V 1 · (cos θ−K · sin θ) + a
【数6】Vy=V2・sinθ+b と表わされる。2相正弦波信号Vx及びVyは、位置検
出部においてA/D変換器によってデジタル化された数
値データx及びyとして取扱うが、このデジタル化され
た数値データx及びyは、スケーリングファクタSfを
用いて以下のように規格化して考えることができる。す
なわち、Vy = V 2 · sin θ + b The two-phase sine wave signals Vx and Vy are handled as numerical data x and y digitized by an A / D converter in the position detecting unit. The digitized numerical data x and y are calculated using a scaling factor Sf. Can be standardized as follows. That is,
【数7】 Vx=Sf・x=V1(cosθ−K・sinθ)+aVx = Sf · x = V 1 (cos θ−K · sin θ) + a
【数8】Vy=Sf・y=V2・sinθ+b ここで、スケーリングファクタSfは2相正弦波信号V
x及びVyの振幅、A/D変換器の分解能等によって任
意に仮定できる変数であるので、正弦波Vyを振幅の基
準として考えてSf=V2と仮定すると、数7及び数8
は以下のように書換えることができる。Vy = Sf ・ y = V 2・ sin θ + b where the scaling factor Sf is a two-phase sine wave signal V
The amplitude of x and Vy, since a variable that can assume arbitrarily by the resolution or the like of the A / D converter, when thinking of a sine wave Vy as a reference amplitude assuming Sf = V 2, number 7 and number 8
Can be rewritten as follows:
【0008】[0008]
【数9】(Equation 9)
【数10】y=sinθ+d ただし、A=V2/V1,C=a/V2,d=b/V2 本発明では、2相正弦波信号のxが0となる時点のyの
瞬時値、2相正弦波信号のyが0となる時点のxの瞬時
値、及び|x|=|y|となる時点のx又はyの瞬時値
を用いて、数9及び数10のA,K,C,d、すなわち
2相正弦波信号x及びyの振幅誤差,位相誤差,直流オ
フセットの各値を同時に演算、検出している。そして、
検出された各誤差成分を用いて、上記2相正弦波信号x
及びyからY = sin θ + d where A = V 2 / V 1 , C = a / V 2 , d = b / V 2 In the present invention, the instant of y at the time when x of the two-phase sine wave signal becomes 0 Using the instantaneous value of x at the time when y of the two-phase sine wave signal becomes 0 and the instantaneous value of x or y at the time when | x | = | y | K, C, and d, that is, the values of the amplitude error, phase error, and DC offset of the two-phase sine wave signals x and y are simultaneously calculated and detected. And
Using the detected error components, the two-phase sine wave signal x
And from y
【数11】x′=A(x−C)+K(y−d)=cos
θX '= A (x-C) + K (yd) = cos
θ
【数12】y′=y−d=sinθ という演算によって、互いに振幅が等しく位相差が90
°であって、直流オフセットを含まない2相正弦波信号
x′及びy′を得ることができ、この2相正弦波信号
x′及びy′から[Mathematical formula-see original document] By the operation of y '= yd-sin [theta], the amplitudes are equal to each other and the phase difference is 90.
° and two-phase sine wave signals x 'and y' which do not include a DC offset can be obtained, and from these two-phase sine wave signals x 'and y'
【数13】θ=tan−1(y′/x′)=tan−1
(sinθ/cosθ) なる演算によって、移動体の移動距離θを誤差なく高精
度に検出することができる。## EQU13 ## θ = tan -1 (y '/ x') = tan -1
By the calculation (sin θ / cos θ), the moving distance θ of the moving body can be detected with high accuracy without errors.
【0009】[0009]
【実施例】以下に図面に基づいて本発明の実施例を詳細
に説明する図1は本発明による位置検出装置の一実施例
であり、磁気式回転位置検出装置に応用した例の構成図
である。回転体1に取付けられたギア2、磁気センサ
3,4及び磁石5は従来と同一であり、磁気センサ3,
4の出力は前記数5及び数6で示したように振幅誤差,
位相誤差,直流オフセットなどの誤差を含む2相正弦波
信号である。この2相正弦波信号Vx及びVyはサンプ
ルホールド回路11及び12,A/D変換器13及び1
4によってデジタル化された数値データx及びyに変換
される。なお、x及びyは前述した数9及び数10で表
わされる。この2相正弦波信号x及びyは図2に示すよ
うに、ギア2の1歯分に相当する回転角を1周期、すな
わち2π(rad)とする正弦波であり、理想的な2相
正弦波に対して振幅誤差A,位相誤差K,直流オフセッ
トC及びdを含んでいる。これら各誤差量は、このよう
な位置検出装置において直接的に検出することは不可能
であるので、2相正弦波信号x及びyの瞬時値からこれ
ら各誤差量A,K,c,dを検出することを考える。ま
ず数9及び数10を変形してDETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS An embodiment of the present invention will be described below in detail with reference to the accompanying drawings. FIG. is there. The gear 2, the magnetic sensors 3, 4 and the magnet 5 attached to the rotating body 1 are the same as in the prior art.
The output of 4 is the amplitude error as shown in the above equations (5) and (6),
This is a two-phase sine wave signal including errors such as a phase error and a DC offset. The two-phase sine wave signals Vx and Vy are supplied to sample-hold circuits 11 and 12, and A / D converters 13 and 1, respectively.
4 is converted into digitized numerical data x and y. Note that x and y are represented by the above-described equations 9 and 10. As shown in FIG. 2, the two-phase sine wave signals x and y are sine waves having a rotation angle corresponding to one tooth of the gear 2 for one cycle, that is, 2π (rad). Includes amplitude error A, phase error K, DC offset C and d for the wave. Since these error amounts cannot be directly detected by such a position detecting device, these error amounts A, K, c, and d are calculated from the instantaneous values of the two-phase sine wave signals x and y. Think about detecting. First, transform equations 9 and 10
【数14】sinθ=y−d## EQU14 ## sin θ = yd
【数15】cosθ=A(x−C)+K(y−d) を得る。数14及び数15はsin2θ+cos2θ=
1という関係が成り立つので、[Mathematical formula-see original document] cos [theta] = A (x-C) + K (yd) is obtained. Equations 14 and 15 are sin 2 θ + cos 2 θ =
Since the relationship of 1 holds,
【数16】 (y−d)2+{A(x−C)+K(y−d)2}=1 を得る。この数16を変形すると、(Y−d) 2 + {A (x−C) + K (y−d) 2 } = 1 is obtained. By transforming this equation 16,
【数17】 A2x2+(1+K2)y2+2AKxy−2A(AC+Kd)x −2{(1+K2)d+AKC}y+{(1+K2)d2+A2 C2+2AKCd−1}=0 となる。ここで位相誤差Kについて考えると、一般的な
磁気センサにおいて位相誤差Kはできるだけ0となるよ
うに設計、製造されていることから、Equation 17] A 2 x 2 + (1 + K 2) y 2 + 2AKxy-2A (AC + Kd) x -2 and {(1 + K 2) d + AKC} y + {(1 + K 2) d 2 + A 2 C 2 + 2AKCd-1} = 0 Become. Here, considering the phase error K, since the phase error K is designed and manufactured so as to be as small as possible in a general magnetic sensor,
【数18】 と近似することが可能であり、その結果、数17は更に
以下のように考えることができる。(Equation 18) Can be approximated, and as a result, Equation 17 can be further considered as follows.
【数19】 A2x2+y2+2AKxy−2A(AC+Kd)X −2(d+AKC)y+(d2+A2C2+2AKCd−1} =0 ちなみに、数19は2相正弦波信号x及びyを2次元平
面の直交する座標軸として考えたとき、このxy平面上
における楕円の方程式にほかならない。Equation 19] A 2 x 2 + y 2 + 2AKxy-2A (AC + Kd) X -2 (d + AKC) y + (d 2 + A 2 C 2 + 2AKCd-1} = 0 In this connection, the number 19 two-phase sinusoidal signals x and y When considered as a coordinate axis orthogonal to a two-dimensional plane, this is nothing but the equation of an ellipse on the xy plane.
【0010】さて、数19の各誤差量A,K,C,dを
同時に算出することが本発明の目的であるが、単に2相
正弦波信号x及びyの瞬時値からこれら誤差量A,K,
C,dを算出することは、演算が非常に複雑となり現実
的でない。そこで、本発明においては、次の3つの特定
のタイミングにおける2相正弦波信号x及びyの瞬時値
を用いて演算することを考える。すなわち、3つのタイ
ミングは x=0となる時のyの瞬時値y1及びy2 y=0となる時のxの瞬時値x1及びx2 |x|=|y|となる時のxの瞬時値x3=y3及
びx4=−y4である。なお、タイミングの場合につ
いて、yの瞬時値y3及びy4を用いても以下の演算は
全く同様の結果が得られる。ここで、これらの特定のタ
イミングにおける2相正弦波信号x及びyの瞬時値を得
る方法について、具体的な実施例を説明する。Now, each of the error amounts A, K, C and d in Expression 19 is
It is an object of the present invention to calculate them at the same time , but simply from these instantaneous values of the two-phase sinusoidal signals x and y, these error amounts A, K,
Calculating C and d is not realistic because the calculation becomes very complicated. Therefore, in the present invention, it is considered that the calculation is performed using the instantaneous values of the two-phase sine wave signals x and y at the following three specific timings. That is, three timing is x = 0 and the instantaneous value y 1 of y when made and y 2 y = 0 the instantaneous value of x when the become x 1 and x 2 | x | = | y | become when x instantaneous value x 3 = y 3 and x 4 = a -y 4. Note that for the case of timing, the following operations be used instantaneous value y 3 and y 4 and y is obtained exactly the same results. Here, a specific example of a method for obtaining instantaneous values of the two-phase sine wave signals x and y at these specific timings will be described.
【0011】図1において、瞬時値検出部10から出力
される2相正弦波信号x及びyは、コンパレータ等によ
って構成される瞬時値比較器41と、メモリ素子等によ
って構成される瞬時値保持回路42とから成る瞬時値記
憶部40に入力されている。瞬時値比較器41は、常時
2相正弦波信号x及びyの瞬時値を検出しており、前記
3つの条件〜、すなわちx=0,y=0及び|x|
=|y|となるタイミングで制御信号CSを発生する。
この制御信号CSが瞬時値保持回路42に入力される
と、そのタイミングにおいて入力されている2相正弦波
信号x及びyの瞬時値を取込んで保持する。なお、この
動作によって保持される2相正弦波信号x及びyの瞬時
値は、図2の波形中にy1,y2,x1,x2,x3,
x4で示される値である。このようにして保持記憶され
たy1,y2,x1,x2,x3,x4の値から、前記
各誤差量A,K,C,dを同時に算出するアルゴリズム
を以下に説明する。なお、この演算は図1中の誤差演算
部50において行なわれる。まず数19に条件、すな
わちx=0のときy=y1及びy3を代入すると、次の
2つの式数20及び数21を得る。In FIG. 1, the two-phase sine wave signals x and y output from the instantaneous value detector 10 are converted into an instantaneous value comparator 41 composed of a comparator and the like, and an instantaneous value holding circuit composed of a memory element and the like. 42, and is input to the instantaneous value storage unit 40. The instantaneous value comparator 41 always detects the instantaneous values of the two-phase sine wave signals x and y, and the three conditions or above, ie, x = 0, y = 0 and | x |
The control signal CS is generated at the timing when | = | y |.
When the control signal CS is input to the instantaneous value holding circuit 42, the instantaneous values of the two-phase sine wave signals x and y input at that timing are captured and held. The instantaneous values of the two-phase sine wave signals x and y held by this operation are represented by y 1 , y 2 , x 1 , x 2 , x 3 , and y 3 in the waveform of FIG.
It is a value indicated by x 4. An algorithm for simultaneously calculating the error amounts A, K, C, and d from the values of y 1 , y 2 , x 1 , x 2 , x 3 , and x 4 thus held and stored will be described below. . This calculation is performed in the error calculation unit 50 in FIG. First condition number 19, that is, substituting y = y 1 and y 3 when x = 0, we obtain the following two formulas number 20 and number 21.
【数20】(Equation 20)
【数21】そして、数20と数21の両辺を減算するこ
とによって、[Mathematical formula-see original document] Then, by subtracting both sides of the equations (20) and (21),
【数22】を得る。同様にして数19に条件、すなわ
ちy=0のときx=x1,x2を代入してThe following is obtained. Similarly, substituting x = x 1 , x 2 for the condition in Equation 19, ie, when y = 0,
【数23】(Equation 23)
【数24】を得る。同様に、数23と数24の両辺を減
算し、24 is obtained. Similarly, subtract both sides of Expression 23 and Expression 24,
【数25】を得る。また、数20と数23を用いて次の
数26を得る。Equation 25 is obtained. Further, the following Expression 26 is obtained by using Expression 20 and Expression 23.
【数26】(Equation 26)
【0012】そして、数26の両辺に数22及び数25
を代入することができ、その結果次式の数27を得る。Then, on both sides of Equation 26, Equations 22 and 25
Can be substituted, and as a result, the following Expression 27 is obtained.
【数27】 この数27を変形すると、[Equation 27] By transforming this equation 27,
【数28】−y1y2=A2(−x1x2)-Y 1 y 2 = A 2 (-x 1 x 2 )
【数29】A2=y1y2/x1x2 が得られる。数28より振幅誤差Aを求めることができ
るが、そのためには平方根演算が必要であり、一般的な
マイクロプロセッサを使用することを考えると、更に簡
単な演算であることが望ましい。そこで、振幅誤差Aが
一般的な磁気センサにおいては極めて1に近いというこ
とを利用して、以下のような近似計算を行なう。すなわ
ち、次式が求められるEquation 29] A 2 = y 1 y 2 / x 1 x 2 is obtained. Although the amplitude error A can be obtained from Expression 28, a square root operation is required for that. Considering the use of a general microprocessor, a simpler operation is desirable. Therefore, the following approximate calculation is performed by utilizing that the amplitude error A is extremely close to 1 in a general magnetic sensor. That is, the following equation is obtained
【数30】A=1+α ただし、A = 1 + α where
【数31】 である。この数30を数29に代入すると、(Equation 31) It is. Substituting Equation 30 into Equation 29 gives
【数32】(1+α)2=y1y2/x1x2 [Number 32] (1 + α) 2 = y 1 y 2 / x 1 x 2
【数33】1+2α+α2=y1y2/x1x2 が得られる。そして、数31が成り立つからEquation 33] 1 + 2α + α 2 = y 1 y 2 / x 1 x 2 is obtained. And because equation 31 holds
【数34】 と近似でき、その結果(Equation 34) And the result is
【数35】α=(y1y2−x1x2)/2x1x2 となる。ここで、再びA=1+αに戻すことによって、Α = (y 1 y 2 −x 1 x 2 ) / 2 × 1 x 2 Here, by returning to A = 1 + α again,
【数36】 A=1+α=(x1x2+y1y2)/2x1x2 となる。A = 1 + α = (x 1 x 2 + y 1 y 2 ) / 2x 1 x 2
【0013】次に数19に前記条件、すなわちx3=
y3及びx4=−y4を代入するが、その前に以下の説
明の簡単化のために、数21及び数24をもとに数18
の見直しを行なう。つまり、数19は次の式ように変形
できる。Next, the above condition, that is, x 3 =
y 3 and x 4 = −y 4 are substituted, but before that, for simplification of the following description, Equation 18 is obtained based on Equation 21 and Equation 24.
Review of. That is, Equation 19 can be transformed into the following equation.
【数37】 A2{x2−2(C+Kd/A)x}+y2−2(d+AKC)y +2AKxy+(d2+A2C2+2AKCd−1)=0 そして、この数37に数22及び数25を代入すると、Equation 37] A 2 {x 2 -2 (C + Kd / A) x} + y 2 -2 (d + AKC) y + 2AKxy + (d 2 + A 2 C 2 + 2AKCd-1) = 0 and the number 22 and number this number 37 Substituting 25 gives
【数38】 A2{x2−(x1+x2)x}+y2−(y1+y2)y +2AKxy+(d2+A2C2+2AKCd−1)=0 となる。このようにして得られた数38を数19の代り
として、条件のx3=y3を代入すると、Equation 38] A 2 becomes (y 1 + y 2) y + 2AKxy + (d 2 + A 2 C 2 + 2AKCd-1) = 0 - {x 2 - (x 1 + x 2) x} + y 2. By substituting the condition x 3 = y 3 for the equation 38 obtained in this way instead of the equation 19,
【数39】 が得られ、この数39の変形により、[Equation 39] Is obtained, and by the transformation of Expression 39,
【数40】 が得られる。同様にして数38にx4=−y4を代入す
ると、(Equation 40) Is obtained. Similarly, by substituting x 4 = −y 4 into Equation 38,
【数41】 が得られ、この数41の変形により、[Equation 41] Is obtained, and by the transformation of Expression 41,
【数42】 となる。(Equation 42) Becomes
【0014】そして、数40と数42の両辺を加算する
ことにより、Then, by adding both sides of Expression 40 and Expression 42,
【数43】 となるので、左辺の係数2AKについて考える。Aが極
めて1に近いことを利用して、数30に示すようにA=
1+α(ただし、数31である)とおくと、[Equation 43] Therefore, the coefficient 2AK on the left side is considered. Taking advantage of the fact that A is extremely close to 1, A =
1 + α (however, Equation 31)
【数44】 と近似することができる。そこで、この数44を数43
に代入して整理すると、[Equation 44] Can be approximated. Therefore, this equation 44 is replaced by
Substituting into
【数45】 が得られ、更に数29で求めたA2を代入し、[Equation 45] Is obtained, and further, A 2 obtained by Expression 29 is substituted,
【数46】 が得られる。[Equation 46] Is obtained.
【0015】以上より、数36及び数46によって、2
相正弦波信号x及びyの瞬時値から2相正弦波の誤差成
分のうちA及びKを解くことができた。次に、数22及
び数25を用いて残りのC及びdについて解くと、まず
数21よりFrom the above, according to equations 36 and 46, 2
From the instantaneous values of the phase sine wave signals x and y, A and K of the error components of the two-phase sine wave could be solved. Next, when solving for the remaining C and d using Equations 22 and 25,
【数47】d=(y1+y2)/2−AKC が得られるので、これを数25に代入し、Equation 47] Since d = (y 1 + y 2 ) / 2-AKC is obtained by substituting this into Equation 25,
【数48】 を得、この数48の変形により、[Equation 48] And by the transformation of the expression 48,
【数49】C+K(y1+y2)/2A−K2C=(x
1+x2)/2 を得る。ここで、C + K (y 1 + y 2 ) / 2A−K 2 C = (x
1 + x 2 ) / 2. here,
【数50】 と近似して[Equation 50] Approximate
【数51】 C=(x1+x2)/2−K(y1+y2)/2A を得る。そして、数51を数47に代入してThe following equation is obtained: C = (x 1 + x 2 ) / 2−K (y 1 + y 2 ) / 2A Then, substituting equation 51 into equation 47
【数52】d=(y1+y2)/2−AK{(x1+x
2)/2−K(y1+y2)/2A} を得、再び数50と近似することにより、D = (y 1 + y 2 ) / 2−AK {(x 1 + x
2 ) / 2−K (y 1 + y 2 ) / 2A}, and by approximating again to Equation 50,
【数53】 d=(y1+y2)/2−AK(x1+x2)/2 が得られる。The following equation is obtained: d = (y 1 + y 2 ) / 2−AK (x 1 + x 2 ) / 2
【0016】以上で2相正弦波の誤差成分A,K,C,
dの全てが2相正弦波信号x及びyの瞬時値(x1,x
2,y1,y2,x3,x4)を用いて、数33,数4
3,数47,数49によって同時に演算できた。すなわ
ち、図1中の誤差演算部50における演算処理は上記の
数33,数43,数47,数49を具現化したものであ
り、これらの演算は、マイクロプロセッサ等の演算素子
によって容易に実現できるものである。As described above, the error components A, K, C, and
d is the instantaneous value (x 1 , x 1) of the two-phase sinusoidal signals x and y
2 , y 1 , y 2 , x 3 , x 4 ), and
3, Equation 47 and Equation 49 could be calculated simultaneously . That is, the arithmetic processing in the error arithmetic unit 50 in FIG. 1 embodies the above equations 33, 43, 47 and 49, and these operations are easily realized by an arithmetic element such as a microprocessor. You can do it.
【0017】このようにして誤差演算部50によって演
算された各誤差成分A,K,C,dを用いて、図1中の
誤差補正部60では以下のようなアルゴリズムにより振
幅誤差がなく、位相差が90°であって直流オフセット
のない、理想的な2相正弦波x′及びy′を算出する。
すなわち、前記数11及び数12に示す演算を行なう。
このようにして得られた2相正弦波信号x′及びy′は
誤差のない理想的な2相正弦波であるので、位置検出部
30において除算器31によってUsing the error components A, K, C, and d calculated by the error calculation unit 50 in this manner, the error correction unit 60 in FIG. The ideal two-phase sine waves x ′ and y ′ having a phase difference of 90 ° and no DC offset are calculated.
That is, the operations shown in the above equations 11 and 12 are performed.
Since the two-phase sine wave signals x 'and y' thus obtained are ideal two-phase sine waves having no error, the position detector 30 uses the divider 31 in the position detector 30.
【数54】y′/x′=sinθ/cosθ を求め、tan−1データテーブル32によってY ′ / x ′ = sin θ / cos θ is obtained, and by using the tan −1 data table 32,
【数55】θ=tan−1(y′/x′)=tan−1
(sinθ/cosθ) のように変換することによって、回転体1の回転角θを
正確に検出できる。なお、図1には示していないが、図
5に示したように上位桁検出部20を設けて回転角θと
加算することによって、回転体1の全周にわたる位置検
出を行なうことができる。[Mathematical formula-see original document] θ = tan −1 (y ′ / x ′) = tan −1
By converting as (sin θ / cos θ), the rotation angle θ of the rotating body 1 can be accurately detected. Although not shown in FIG. 1, as shown in FIG. 5, by providing the upper digit detection unit 20 and adding the upper digit detection unit 20 to the rotation angle θ, the position detection over the entire circumference of the rotating body 1 can be performed.
【0018】ところで、以上説明したような誤差演算処
理,誤差補正処理等は、A/D変換器によって2相正弦
波をデジタル化し、数値化されたデータを用いて処理す
ることが一般的である。そして、このようにA/D変換
を行なう場合、A/D変換器の変換速度の制約及び演算
の処理時間の制約から、デジタル数値化された2相正弦
波は連続的な値をとることができない。その結果、図3
に示すように、例えば前記条件においてx=0となる
タイミングが正確に検出できず、そのためyの瞬時値y
1(又はy2)の正しい値を得ることができないという
事態が生じる。この状態を図3を用いて説明する。図中
のVxが0になるタイミングに対して、デジタル数値化
されたデータxが0になったことを検出できるのは図中
のA若しくはBの区間であり、そのため瞬時値保持回路
42が保持するy1,y2には、図3中の△y1,△y
2だけの誤差を含むことになる。このような事態を避け
るために、図4に示すようにサンプリング周期Tに対し
てVx及びVyの周期が大きな時、つまり回転体1の回
転速度が低い場合にのみ、瞬時値保持回路42が2相正
弦波信号x及びyの瞬時値を取込んで保持するようにし
ており、誤差△y1が精度上無視できるほど小さくなる
ようにしている。すなわち、回転体1の回転速度が高い
場合には、瞬時値保持回路42は2相正弦波信号x及び
yの瞬時値を更新せず、回転速度が低い際に取込んだ瞬
時値を保持し続けるようにしている。この結果、回転速
度が高い場合には2相正弦波に含まれる各誤差はリアル
タイムで検出できず、回転数が低い際に検出された値を
そのまま使用して補正することになる。しかしながら、
一般的に位置検出精度が必要となるのは極低速で使用さ
れる場合が多いことから、実用上支障はきたさない。Incidentally, the above-described error calculation processing, error correction processing, and the like are generally performed by digitizing a two-phase sine wave with an A / D converter and using numerical data. . When the A / D conversion is performed as described above, the digital-valued two-phase sine wave may take a continuous value due to the limitation of the conversion speed of the A / D converter and the limitation of the processing time of the operation. Can not. As a result, FIG.
As shown in, for example, the timing at which x = 0 under the above conditions cannot be accurately detected, and therefore, the instantaneous value y of y
A situation arises in which the correct value of 1 (or y 2 ) cannot be obtained. This state will be described with reference to FIG. It is in the section A or B in the figure that it is possible to detect that the digitized data x has become 0 at the timing when Vx becomes 0 in the figure. Y 1 , y 2 in FIG. 3 include △ y 1 , 中 y
It will contain only two errors. In order to avoid such a situation, as shown in FIG. 4, only when the period of Vx and Vy is larger than the sampling period T, that is, when the rotation speed of the rotating body 1 is low, the instantaneous value holding circuit 42 The instantaneous values of the phase sine wave signals x and y are fetched and held, so that the error Δy 1 is so small that it can be ignored in accuracy. That is, when the rotation speed of the rotating body 1 is high, the instantaneous value holding circuit 42 does not update the instantaneous values of the two-phase sine wave signals x and y, but holds the instantaneous value captured when the rotation speed is low. I try to keep going. As a result, when the rotation speed is high, each error included in the two-phase sine wave cannot be detected in real time, and correction is performed using the value detected when the rotation speed is low. However,
Generally, position detection accuracy is required at very low speeds in many cases, so that there is no practical problem.
【0019】以上、本発明の実施例を磁気式回転位置検
出装置について説明したが、本発明の原理は回転型位置
検出装置に限定されず、直線移動型の位置検出装置でも
全く同様に適用できる。また、ギアを利用した磁気式に
限らず、例えば着磁されたドラムを利用した検出器、光
学式の検出器など2相正弦波信号を利用したものであれ
ば同様に適用できる。Although the embodiment of the present invention has been described above with reference to a magnetic rotary position detecting device, the principle of the present invention is not limited to a rotary type position detecting device, and can be applied to a linear moving type position detecting device in the same manner. . The present invention is not limited to the magnetic type using a gear, but may be applied to any type using a two-phase sine wave signal such as a detector using a magnetized drum or an optical type detector.
【0020】[0020]
【発明の効果】以上のように本発明による誤差補正機能
付き位置検出装置によれば、移動体の移動距離に応じて
検出された2相正弦波信号に含まれる振幅誤差,位相誤
差及び直流オフセット誤差を自動的に同時に検出して補
正することにより、これら各誤差の影響を含まない高精
度な位置検出装置を実現できる。また、この誤差の検出
及び補正を行なう演算処理は特別な処理回路を必要とす
ることなく、安価なマイクロプロセッサ等によって実現
できる。さらにまた、2相正弦波信号に上記のような誤
差成分が含まれていても自動的に補正して高精度な位置
検出を実現できることから、従来2相正弦波信号の調整
作業に必要であった多大な時間工数を削減することがで
きると共に、調整作業における人為的な調整不良の影響
がなく、信頼性の高い位置検出装置を実現できる。As described above, according to the position detecting device with an error correcting function according to the present invention, the amplitude error, the phase error and the DC offset included in the two-phase sine wave signal detected according to the moving distance of the moving body. By automatically detecting and correcting the errors at the same time , it is possible to realize a high-accuracy position detection device that does not include the effects of these errors. The arithmetic processing for detecting and correcting the error can be realized by an inexpensive microprocessor or the like without requiring a special processing circuit. Furthermore, even if the two-phase sine wave signal contains the above error component, it can be automatically corrected to realize high-accuracy position detection. A great amount of time and man-hours can be reduced, and a highly reliable position detection device can be realized without the influence of artificial adjustment failure in the adjustment work.
【図1】本発明の一実施例である磁気式回転位置検出装
置の構成図である。FIG. 1 is a configuration diagram of a magnetic rotational position detecting device according to an embodiment of the present invention.
【図2】本発明の原理を説明するための2相正弦波信号
の波形図である。FIG. 2 is a waveform diagram of a two-phase sine wave signal for explaining the principle of the present invention.
【図3】本発明の動作を説明するための図であり、磁気
センサによって検出された2相正弦波信号と、これをA
/D変換器によってデジタル数値化したデータを示す波
形図である。FIG. 3 is a diagram for explaining the operation of the present invention, in which a two-phase sine wave signal detected by a magnetic sensor and this signal is represented by A
FIG. 4 is a waveform diagram showing data converted into digital numerical values by a / D converter.
【図4】図3と同様の波形図であり、回転体の回転速度
が低い場合を示すものである。FIG. 4 is a waveform diagram similar to FIG. 3, showing a case where the rotation speed of a rotating body is low.
【図5】従来技術による磁気式回転位置検出装置の構成
図である。FIG. 5 is a configuration diagram of a magnetic rotational position detecting device according to the related art.
【図6】従来の検出装置をソフトウェア処理によって実
現した場合のフローチャートである。FIG. 6 is a flowchart when a conventional detection device is realized by software processing.
1 回転体 2 ギア 10 瞬時値検出部 30 位置検出部 40 瞬時値記憶部 41 瞬時値比較器 42 瞬時値保持回路 50 誤差演算部 DESCRIPTION OF SYMBOLS 1 Rotating body 2 Gear 10 Instantaneous value detection part 30 Position detection part 40 Instantaneous value storage part 41 Instantaneous value comparator 42 Instantaneous value holding circuit 50 Error calculation part
Claims (2)
相の異なる第1及び第2の2相正弦波信号x及びyを出
力する検出部と;前記第1の正弦波信号をx=1/A・
(cosθ−K・sinθ)+C、前記第2の正弦波信
号をy=sinθ+dと表わしたとき、x=0となる時
点におけるyの瞬時値y1及びy2と、y=0となる時
点におけるxの瞬時値x1及びx2と、|x|=|y|
となる時点におけるx又はyの瞬時値x3及びx4、又
はy3及びy4とを記憶する記憶手段と;前記記憶手段
に記憶されているy1,y2,x1,x2,x3又はy
3,x4又はy4を用いて、前記2相正弦波信号x及び
yの振幅の差を表わすAを数1の演算により求め、前記
2相正弦波信号x及びyの位相の差を表わすKを数2の
演算により求め、前記2相正弦波信号x及びyの直流オ
フセットを表わすc及びdを数3及び数4の演算により
求める演算手段と;前記誤差演算手段によって同時に求
めたC,d,A,Kを用いて、前記2相正弦波信号x及
びyから、互いに位相が90°異なり振幅が等しく、か
つ直流オフセット成分のない第3及び第4の2相正弦波
信号x′=cosθ及びy′=sinθを求める誤差補
正手段と;前記第3及び第4の2相正弦波信号x′及び
y′からθ=tan−1(y′/x′)の演算に基づい
て前記移動体の移動距離θを求める位置検出部とを具備
したことを特徴とする誤差補正機能付き位置検出装置。 【数1】 【数2】 【数3】 【数4】 A detecting unit that outputs first and second two-phase sine wave signals x and y having different phases according to a moving distance θ of a moving object; 1 / A ・
(Cos θ−K · sin θ) + C, when the second sine wave signal is expressed as y = sin θ + d, instantaneous values y 1 and y 2 of y at the time of x = 0 and at the time of y = 0 the instantaneous value x 1 and x 2 of x, | x | = | y |
Become instantaneous value of x or y at time x 3 and x 4, or y 3 and y 4 and storage means for storing a; y 1 stored in the storage means, y 2, x 1, x 2, x 3 or y
3, with x 4 or y 4, determined by calculation the number 1 of the A representative of the difference in amplitude of the two-phase sine wave signals x and y, the
K representing the phase difference between the two-phase sinusoidal signals x and y is expressed by
It is obtained by an operation, and the DC output of the two-phase sine wave signals x and y is obtained.
C and d representing the offset are calculated by the operations of Equations 3 and 4.
Calculating means ; and using C, d, A, and K simultaneously obtained by the error calculating means, the two-phase sine wave signals x and y have phases different from each other by 90 ° and equal amplitudes, Error correction means for obtaining third and fourth two-phase sinusoidal signals x '= cos θ and y' = sin θ without a DC offset component; and the third and fourth two-phase sinusoidal signals x 'and y' And a position detecting unit for calculating the moving distance θ of the moving body based on the calculation of θ = tan −1 (y ′ / x ′). (Equation 1) (Equation 2) (Equation 3) (Equation 4)
所定値以上の場合には前記y1,y2,x1,x2,x
3,x4を更新することなく保持しており、前記移動体
の移動速度が前記所定値以下の場合には前記y1,
y2,x1,x2,x3,x4をそれぞれx=0となる
時点、y=0となる時点、|x|=|y|となる時点に
おいて順次更新して記憶するようになっている請求項1
に記載の誤差補正機能付き位置検出装置。Wherein said storage means, the moving speed of the moving body in the case of more than the predetermined value y 1, y 2, x 1 , x 2, x
3, x 4 holds without updating the said when the moving speed of the moving body is equal to or less than the predetermined value is the y 1,
y 2 , x 1 , x 2 , x 3 , and x 4 are sequentially updated and stored at a time when x = 0, a time when y = 0, and a time when | x | = | y | Claim 1
4. The position detecting device with an error correction function according to 1.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP04089391A JP3092100B2 (en) | 1992-03-13 | 1992-03-13 | Position detector with error correction function |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP04089391A JP3092100B2 (en) | 1992-03-13 | 1992-03-13 | Position detector with error correction function |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH05256638A JPH05256638A (en) | 1993-10-05 |
JP3092100B2 true JP3092100B2 (en) | 2000-09-25 |
Family
ID=13969363
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP04089391A Expired - Fee Related JP3092100B2 (en) | 1992-03-13 | 1992-03-13 | Position detector with error correction function |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3092100B2 (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3893033B2 (en) * | 2001-06-29 | 2007-03-14 | オークマ株式会社 | Position detection device |
JP2006201148A (en) * | 2004-12-22 | 2006-08-03 | Toshiba Mach Co Ltd | Signal processing device, signal processing method, signal processing program, recording medium with the signal processing program stored, speed detection device, and servo mechanism |
CN111336963B (en) * | 2020-03-27 | 2021-05-18 | 重庆市计量质量检测研究院 | Error compensation method and measurement error evaluation method based on biaxial symmetry standard device |
CN111336964B (en) * | 2020-03-27 | 2021-05-18 | 重庆市计量质量检测研究院 | Error compensation method and measurement error evaluation method based on uniaxial symmetry standard device |
-
1992
- 1992-03-13 JP JP04089391A patent/JP3092100B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JPH05256638A (en) | 1993-10-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6188341B1 (en) | Encoder interpolation circuit which corrects an interpolation angle between a received sine-wave encoder signal and a cosine-wave encoder signal | |
JP4138899B2 (en) | Phase difference detection apparatus, position detection system and method for position detection | |
JP6005781B2 (en) | Resolver device | |
JP3772121B2 (en) | Encoder signal processing device | |
JPH06167354A (en) | Interpolation unit for scale | |
US4795954A (en) | Resolver controlling method and apparatus | |
JP3624458B2 (en) | Digital angle detection method | |
JP2014025871A (en) | Encoder output signal correction apparatus | |
JP3092100B2 (en) | Position detector with error correction function | |
JPH05231879A (en) | Correcting method for detecting position | |
US6484120B1 (en) | Position detection data generating method and apparatus based on phase shift principle | |
JP3365063B2 (en) | Rotation angle detection method | |
JPH0658769A (en) | Signal processing method and displacement detector using method thereof | |
JPH0529045B2 (en) | ||
JPS6211113A (en) | Position detector | |
JP3439814B2 (en) | Digital PLL device | |
JP3592432B2 (en) | Amplitude correction method of sine wave signal in encoder | |
JP3349575B2 (en) | Displacement speed detection method and device | |
JP3451851B2 (en) | Rotation speed detector | |
JP2938472B2 (en) | Rotation angle detector | |
JPH1062203A (en) | Position detector | |
JPH03296618A (en) | Correcting method of waveform for interpolation | |
JPH0469079A (en) | Apparatus for sensing motor speed | |
JP3208932B2 (en) | Position measuring device | |
JP2741431B2 (en) | Encoder that detects the movement position from any reference point |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090728 Year of fee payment: 9 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100728 Year of fee payment: 10 |
|
LAPS | Cancellation because of no payment of annual fees |