Apparatus and method for compensating for noise in signals

Abstract

Signal compensation apparatus and method compensates for and removes errors due to error factors included in source signals x and y approximated by the sine function and cosine function, respectively, from a measurement apparatus such as a resolver or an encoder. Compensation is achieved by approximating the source signals x and y as expressions which include the terms due to the error factors, detecting the maximum values of x, x+y, etc., and determining coefficients corresponding to the error factors on the basis of the relationships between the maximum values and the coefficients, and removing the errors from the signals provided by the measurement apparatus.

Description

BACKGROUND OF THE INVENTION

1. Field of the Invention

This invention relates to a signal compensator which automatically deletes noise components of signals approximated by sin .theta. and cos .theta., to produce sin .theta. and cos .theta. signals with low distortion. In particular, this invention can be applied to high-resolution measurement of position or angle for resolver and encoder signals through arc-tangent operation, to reduce measurement error.

2. Background Art

Resolvers and encoders are widely used to measure rotational angles (.theta.). Resolvers output one or more signal cycles approximated by sin .theta. and cos .theta. for each revolution, and these are processed (by arc-tangent operation) to produce a signal equivalent to .theta.. Encoders output several dozen to several thousand signal cycles approximated by sin .theta. and cos .theta. for each revolution, and in most cases the wave numbers of these signals are used to count the amount of rotation. Where higher resolution is required, arc-tangent operation is used to yield .theta. just as for resolvers.

To determine .theta. through arc-tangent operation with high precision, distortion-free sin .theta. and cos .theta. signals are essential. But, in fact, a variety of distortion components are included in the signals. Further, changes in temperature-related and time-related device characteristic can affect the amount of distortion present.

These problems can be avoided to some extent through the use of high-stability devices, temperature compensation for electronic circuitry, and precise adjustment. These measures, however, increase the cost of the equipment and still leave doubt as to reliability.

The inventor has proposed a magnetoresistance element (U.S. Pat. No. 4,843,365) providing low-distortion sin .theta. and cos .theta. signals through the use of a magnetic encoder. But even the use of this element can not prevent the obtained signal from degenerating through amplifier zero drift, gain fluctuation, and non-linearity. In addition, improper device positioning and scale pitch error can prevent realization of the full potential of the device.

SUMMARY OF THE INVENTION

Accordingly, it is an object of this invention to provide a signal compensator which automatically deletes noise components generated in measurement instrumentation by sine and cosine functions signal zero point and amplitude fluctuations, and cross-talk, to produce low-distortion sine and cosine signals.

It is another object of this invention to provide a signal compensator which can measure position or angle with high resolution to reduce measurement error, and which makes possible an inexpensive high-precision measurement system.

In order to attain the above objects, this invention provides a signal compensator for eliminating error due to at least one of amplitude difference, zero point variation, cross-talk, second-order harmonics and third-order harmonics of two signals x and y approximated by sine function and cosine function, respectively, said signal compensator comprising: a first means for approximating the two signals x and y by the following expressions, and detecting respective maximum values of d.sub.0, d.sub.1, d.sub.2, d.sub.3, d.sub.4, d.sub.5, d.sub.6, and d.sub.7, of x, x+y, y, -x+y, -x, -x-y, -y and x-y:

x=G(cos .theta.+g cos .theta.+k sin .theta.-b.sub.x cos 2.theta.+h cos 3.theta.)+z.sub.x y=G(sin .theta.-g sin .theta.+k cos .theta.+b.sub.y cos 2.theta.-h sin 3.theta.)+z.sub.y where G, g, k, b.sub.x, b.sub.y, h, z.sub.x, and z.sub.y are respective coefficients corresponding to the size of respective error factors consisting of average amplitude, amplitude difference between the two signals x and y, cross-talk, second-order harmonics included in the signal x, second-order harmonics included in the signal y, third-order harmonics included in both signals x and y at the same amount respectively, zero point of the signal x and zero point of the signal y; a second means for determining at least one of the coefficients using the following expressions providing the relationships between the values detected by the first means and the respective sizes of the error factors: ##EQU1## and; a third means for compensating the signals x and y, or a signal or signals derived from the signals x and y, using the coefficients determined by the second means.

This invention improves measurement precision by automatically determining the value of the distortion components in the signal processing step and subtracting the distortion components from the signals, and by predicting measurement error due to distortion components and compensating the result. These operations are handled digitally with a microprocessor, providing high reliability and stability with simple adjustment.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a block diagram showing the configuration of one example of the signal compensator according to this invention;

FIGS. 2 and 3, 3(a)-(c) are flow charts showing the data flow and the processing flow, respectively, in the signal compensator in FIG. 1;

FIGS. 4(a) and 4(b) are graphs showing the change of error after power is turned on, for one rotation in the positive direction;

FIG. 5 is a graph showing the size of error after compensation, compared with that before compensation.

FIG. 6 is an illustrative circuit diagram for the signal compensator shown in FIG. 1;

FIG. 7 is an illustrative diagram for the error factor computing circuit of FIG. 6;

FIG. 8 is an illustrative diagram for the compensating circuit of FIG. 6; and

FIG. 9 is an illustrative diagram for the arc-tangent operation circuit of FIG. 6.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

First, the principle of this invention will be explained more fully. The signal compensator of this invention is connected to the output of a measuring apparatus, such as a resolver, an encoder, and so on, through amplifiers and analog-to-digital (A/D) converters for signals x and y. The input signals x and y of this signal compensator can be approximated by the following expressions:

x=G(cos .theta.+g cos .theta.+k sin .theta.-b.sub.x cos 2.theta.+h cos 3.theta.)+z.sub.x (1) y=G(sin .theta.-g sin .theta.+k cos .theta.+b.sub.y cos 2.theta.-h sin 3.theta.)+z.sub.y (2) where G is the average amplitude and does not affect the arc-tangent operation, but the other coefficients g, k, b.sub.x, b.sub.y, h, z.sub.x and z.sub.y all cause errors.

Coefficient g represents the amplitude difference between x and y, and if the gain of their signal amplifiers differs, a non-zero value will result. Amplifier gain is generally determined by the resistance, but this itself contains error which makes precise adjustment difficult. Further, temperature causes change in resistance, and therefore, when the amount of change on the signal x side differs from that on the signal y side, the second terms appear in the expressions (1) and (2).

Coefficient k indicates cross-talk, and is a non-zero value when the phase difference between x and y is not precisely .pi./2. This occurs when the signal lines approach and cause signal mixing between them, or when the encoder scale pitch is different from the value demanded by the detector pitch.

Coefficients b.sub.x and b.sub.y indicate the non-linearity of the amplifiers, and are non-zero values if their non-linearity cannot be neglected. The relation between non-linear amplifier input v.sub.i and output v.sub.0 can be approximated by the second-order function v.sub.o =a v.sub.i -b v.sub.i.sup.2 +c, and when v.sub.i is replaced by sin .theta. and cos .theta. the following expressions result, causing the appearance of the fourth terms in the expressions (1) and (2) above:

a cos .theta.-b cos.sup.2 .theta.+c=a cos .theta.-(b/2)cos 2.theta.+b/2 +c a sin .theta.-b sin.sup.2 .theta.+c=a sin .theta.+(b/2)cos 2.theta.-b/2 +c The amount of non-linearity is dependent on the amplifier itself, and therefore the size of the fourth terms may differ between x and y.

Coefficient h represents the amount of third-order harmonics present, and will be a non-zero value when the source signal (unprocessed signal) from the detector is distorted from a sine wave. Distortion of the signal from the sine wave means harmonic content of a variety of orders and phases. However, for measurement instruments such as resolvers and encoders the detector is generally designed to be symmetric with respect to the signal sign and direction of rotation, so that harmonics of even orders and harmonics of odd orders with phases different from the fundamental wave will be suppressed. In these measurement devices the magnetic field distribution is unified by the air gap to reduce high-order harmonics, leaving third-order harmonics dominant. This means that the distortion of the source signal from the detector can be approximated by the fifth terms in the expressions (1) and (2) above.

The terms z.sub.x and z.sub.y express the zero points of the signals, and will be non-zero values when the zero point fluctuation or adjustment error cannot be ignored. Amplifiers exhibit a phenomenon called drift, where temperature change and other factors can cause the zero point to vary. In many cases the zero point can be adjusted with a variable resistance (trimmer), but precise adjustment is difficult and temperature change will cause variation in resistance as well. In this case, the sixth terms appear in expressions (1) and (2). When the measurement circuit functions from a single power supply, or when A/D converter outputs only positive values, this tenn appears in operation processing, but this is a problem in the definition of the reference point or in the signing system used for integers, and therefore, outside the scope of this invention.

This invention monitors input signals x and y to automatically determine any of the above coefficients that cannot be ignored and compensates accordingly. Where all coefficients are determined the procedure is as described below.

The maximum values for x, x+y, y, -x+y, -x, -x-y, -y and x-y are determined respectively as d.sub.o, d.sub.1, d.sub.2, d.sub.3, d.sub.4, d.sub.5, d.sub.6 and d.sub.7. This may appear to be a random parameter selection, but d.sub.i, corresponds to the maximum value for the general expression x cos i(.pi./4)+y sin i(.pi./4). This is equivalent to .pi./4 radii for the Lissajous' curve obtained when the x and y signals are plotted orthogonally. In this specification, however, where i is an odd number d.sub.i, the value in the above general expression can be treated as times .sqroot.2. This is because the expression itself can be simplified, reducing the operations required. This treatment is employed in the actual operation. The coefficients d.sub.o through d.sub.7 have the following relationships: ##EQU2## And their sums and differences result in the following values s.sub.0 through s.sub.7 :

s.sub.0 =(d.sub.0 -d.sub.4)/2=z.sub.x -Gb.sub.x s.sub.1 =(d.sub.1 -d.sub.5)/2=z.sub.x +x.sub.y s.sub.2 =(d.sub.2 -d.sub.6)/2=z.sub.y -Gb.sub.y s.sub.3 =(d.sub.3 -d.sub.7)/2=-z.sub.x +z.sub.y s.sub.4 =(d.sub.o +d.sub.4)/2=G(1+g+h) ##EQU3## This allows the coefficients to be determined through the following expressions: z.sub.x =(s.sub.1 -s.sub.3)/2 z.sub.y =(s.sub.1 +s.sub.3)/2 Gb.sub.x =z.sub.x -s.sub.0 Gb.sub.y =z.sub.y -s.sub.2 Gg=(s.sub.4 -s.sub.6)/2 Gk=(s.sub.5 -s.sub.7)/2 Gh=s.sub.4 -s.sub.5 +s.sub.6 -s.sub.7)/4 G=(s.sub.4 +s.sub.5 +s.sub.6 +s.sub.7)/4 Once these are known, approximate values for sin .theta. and cos .theta. are determined from the expressions: cos .theta..apprxeq.A=(x-z.sub.x)/G sin .theta..apprxeq.B=(y-z.sub.y)/G making it possible to apply the double angle formula as: cos 2.theta..apprxeq.C=A.sup.2 -B.sup.2 cos 3.theta..apprxeq.D=A(A.sup.2 -3B.sup.2) sin 3.theta..apprxeq.E=B(3A.sup.2 -B.sup.2) and to compensate the source signal through product with Gg, Gk, Gb.sub.x, Gb.sub.y and Gh: G cos .theta..apprxeq.x-G g A-G k B+Gb.sub.x C-Gh D-z.sub.x G cos .theta..apprxeq.y+G g B-G k A-G b.sub.y C+G h E-z.sub.y

The coefficients d.sub.o through d.sub.7 represent the maximum values for x, x+y, y, -x+y, -x, -x-y, -y and x-y. However, since the encoder or resolver shaft cannot be expected to rotate at a fixed speed, simple peak detection methods cannot be used. Since x and y are signals approximating cos .theta. and sin .theta., it is possible to detect these values by sampling x, x+y, etc., where .theta. is i .pi./4.

Of course, there is the possibility that .theta. is exactly equal to i .pi./4, but since cos 2.pi./64=0.995 and cos 2.theta./128=0.999, signal measurement with a 0.1% resolution would allow angles from .+-.2.pi./64 to .+-.2.pi./128 to be used as the maximum values.

Usually, the value .theta. used here is obtained through arc-tangent operation, but the value which has not been subjected to arc-tangent operation may be used in some cases. For example, when .theta. is close to zero, .theta. can be approximated by sin .theta., and when y or its compensated absolute value is small, it is sufficient to sample d.sub.o when x is positive or d.sub.4 when x is negative. In the same way, when the absolute x value is small it is sufficient to sample d.sub.2 or d.sub.6 depending on the sign of y, and d.sub.1 or d.sub.5 depending on the sign of x+y when the absolute value of x-y is small, and d.sub.7 or d.sub.3 depending on the sign of x-y when the absolute value of x+y is small.

While measurement values will contain error due to intensity variation stemming from noise and scaling position, it can be expected that peak value change over time is slight, and therefore, some form of averaging for peak sampling is desirable. In this case a moving average is often used, but this requires storage of a great number of past measurement points, and this is ill-suited to microprocessors with limited storage capacity. If the expression

d.sub.o =d.sub.o +p(x-d.sub.o) is used to update the peak value (given that p is a constant where 0<p<1), it is sufficient to store only the latest peak value, which is extremely advantageous. This method is in fact a weighed average where weighing is reduced exponentially with time.

Use of this method means that as the rotating shaft stops near i .pi./4, the specific peak value would be continuously updated. This would make the averaging to a specific biased data, and generate bad such as causing gain fluctuation to be misinterpreted as zero point variations. To avoid this, it is desirable to perform d.sub.i correction and d.sub.i.+-.4 correction alternately.

Expressions (1) and (2) above take into account all possible factors for error, but in fact in most cases it is sufficient to compensate only for amplitude and zero point variation. In these cases, coefficients k, b.sub.x, b.sub.y, and h can be replaced by 0, as described below.

Where error factors are due only to amplitude and zero point fluctuations, the input signals can be expressed as:

x=g.sub.x cos .theta.+z.sub.x y=g.sub.y sin .theta.+z.sub.y where g.sub.x =G (1+g) and g.sub.y =G (1-g). In this case only d.sub.o, d.sub.2, d.sub.4 and d.sub.6 need be used, with relations: d.sub.o =g.sub.x +z.sub.x d.sub.2 =g.sub.y +z.sub.y d.sub.4 =g.sub.x -z.sub.x d.sub.6 =g.sub.y -z.sub.y This allows the determination of: g.sub.x =(d.sub.o +d.sub.4)/2 g.sub.y =(d.sub.2 -d.sub.6)/2 z.sub.x =(d.sub.o -d.sub.4)/2 z.sub.y =(d.sub.2 -d.sub.6)/2 In this case signal compensation is simplified, and the following expressions can be used: cos .theta.=(x-z.sub.x)/g.sub.x sin .theta.=(y-z.sub.y)/g.sub.y

It is interesting to compare these extremely simple compensation expressions and the complete compensation expressions given above. It can be seen that the third harmonics enlarges the amplitude to (1+h) times, and that the second harmonics changes zero points by -Gb.sub.x and -Gb.sub.y only. This can be taken to mean that where the third harmonics can be ignored G can be determined with the simplified expression, and where the second-order harmonics can be ignored the simplified expression can be used to determine z.sub.x and z.sub.y. Further, because G error is canceled by the arc-tangent operation and because errors in other coefficients are canceled where second terms can be ignored, the same result is obtained regardless of which expression is used to determine G, for all practical purposes.

Where the amplitudes of two input signals are equal and only the cross-talk term is left, in other words where g=b.sub.x =b.sub.y =h=0, the source signals can be given as:

x=G(cos .theta.+k sin .theta.)+z.sub.x y=G(sin .theta.+k cos .theta.)+z.sub.y When the relations for d.sub.1, d.sub.3, d.sub.5 and d.sub.7 are: ##EQU4## this results in the calculation of coefficients as: ##EQU5##

The method described above consists of detection of the amount of the error factors included in the source signal, and compensation by subtracting these from the source signal. If the size of the error factor is known, it is possible to estimate the size of the error when a signal with that error component is arc-tangent processed, allowing compensation of the operation results. This method is described below.

The source signal can be represented in vector notation as follows:

.vertline.r*=.vertline.r+.vertline.g+.vertline.k+.vertline.b+.vertline.h+z Each vector can be expressed as indicated below when the components are given: ##EQU6## Although the angle .theta. between the vector .vertline.r and the X axis is needed, the angle determined through arc-tangent operation using the source signal is the angle .theta.* between .vertline.r* and the X axis. The difference between them can be approximated by the tangential components g.sub.t, k.sub.t, b.sub.t, h.sub.t and z.sub.t of the vectors .vertline.g, .vertline.k, .vertline.b, .vertline.h and z. These tangential components are given by the inner product of each vector and the tangential unit vectors (-sin .theta., cos .theta.) to produce the following expressions: g.sub.t =-g sin 2.theta. k.sub.t =k cos 2.theta. b.sub.t =b.sub.x sin .theta.cos 2.theta.+b.sub.y cos .theta.cos 2.theta. h.sub.t =-h sin 4.theta.z.sub.t =(-z.sub.x sin .theta.+z.sub.y cos .theta.)/G These, then, can be used in the expression: .theta.=.theta.* -g.sub.t -k.sub.t -h.sub.t -z.sub.t to accurately determine the angle .theta.. Note that where the absolute error vector is small compared to 1, the error can be calculated using .theta.* instead of .theta..

Either of the two compensation methods described above may be used individually, or they may be used specifically to match different error factors. In practical application the correction of the effects of third-order harmonics at the source signal requires more calculations, and compensation of angle error may therefore, be more efficient. But direct compensation of other errors at the source signal is generally most efficient. Application requires selection of the optimum method, taking into account the amount of calculation required and remembering that accurate compensation is not possible where error is large because compensation is implemented through approximation.

FIG. 1 shows a block diagram of one example of the signal compensator of this invention, which comprises an error factor eliminating means 1, an arc-tangent operating means 2, a maximum value detecting means 3, a coefficients calculating means 4 and control means 5.

The error factor eliminating means I receives the input signals x and y, approximates the signals in accordance with the above-mentioned expressions, and compensates the signals using the coefficients supplied from the coefficients calculating means 4. The arc-tangent operating means 2 receives the compensated cos .theta. and sin .theta. signals and calculates .theta. by arc-tangent operation. The arc-tangent operation is an established principle used in signal processing wherein the angle of the signal, i.e., .theta., is determined by the familiar mathematical arc-tangent (tan.sup.-1) function. As can be seen from FIG. 2, the angle being calculated is the arc-tangent of the ratio of x and y. Thus, for example, the tan.sup.-1 of 1/2 is 26.56 degrees.

The maximum value detecting means 3 detects the maximum values of x, x+y, etc., using the method mentioned above. The coefficients calculating means 4 calculates the coefficients in accordance with the above-mentioned expressions providing the relationships between the maximum values detected by the maximum value detecting means 3 and the respective size of the error factors. The output of the coefficients calculating means 4 is fed back to the error factor eliminating means 1 for compensation. The control means 5 controls the operations of the means 1 through 4.

FIG. 2 indicates the data flow in the signal compensator of FIG. 1, and FIGS. 3 (a)-(c) indicate the processing flow of a microprocessor handling the signal processing described in this invention, giving C language function names and operators. When the program as shown is started, it initializes the parameters (setting errors to zero, setting amplitude to specified values, setting maximum values to values indicated by amplitude), and then enters a compensation operation loop. The functions in the INITIALIZATION block are shown at the top of FIG. 3(c).

Within the loop the signals x and y are read from A/D converters or other devices (not shown), and the zero point amplitude used to determine cos .theta. and sin .theta.. The double angle formula is then used to determine cos 2.theta.. Next, the values determined above and the coefficients representing the amounts of included error are used to compensate the signals, and .theta. is determined through arc-tangent operation. The arc-tangent operation here is iterated for the angular range of vector (x,y) from -.pi. to .pi., just as with the C language function atan2. This function actually is implemented by case branching and table reference through the signs of both signals and the absolute value size functions. Next, error due to the third harmonics is predicted, compensated for, and the result output as a value from 0 to 2.pi..

Next, the value of .theta. is used to determine if any maximum values need to be compensated. Maximum value compensation is carried out alternately (using a flag such as f04, FIG. 3(b)) for the values where .theta. is different by .pi. from each other. Exponential smoothing is used in compensation to weight the latest value with 1/256 weighing (in this example). After compensation, the alternating compensation flag and the flag indicating that all maximum values have been reset (approximate f bit) are set. Finally, the coefficients are recalculated if all have been updated.

FIG. 3 follows C language expression to make comprehension easier, but Assembler description is actually used to enhance execution speed. Processing flow is simplified in this figure, but in fact actual processing consists of multiple interrupt routines operating together to minimize the time from input of the signals to output of .theta..

FIGS. 4 (a) and (b) show the change in error in an encoder after power is turned on, for one rotation in the positive direction. FIGS. 4(a), (b) indicate, respectively, the error for the first half-rotation and the second half-rotation, at .+-.1/10000 rotations full-scale. In this case, one rotation of the encoder shaft corresponds to 512 cycles of sin .theta. and cos .theta.. In FIGS. 4(a), (b), the arc-tangent error reduces with the rotation of the shaft. FIG. 5 also plots the arc-tangent operation error before and after compensation, with ordinate of .+-.1/10000 rotations and abscissa of .theta. from 0 to 2.pi.. From these figures it is clear that error rapidly decreases with rotation after the power is turned on, demonstrating the effectiveness of the present invention.

FIGS. 6-9 show illustrative circuits which may constitute the present signal compensator, and the flow of signals therein. The term "S/H" denotes a sample and hold circuit, and "Comp" denotes a digital comparator for detecting a value at or near the angle indicated to the right of the respective comparator in FIG. 6. Additional details of the compensating circuit 6, the error factor computing circuit 7 and the arc-tangent operating circuit 8 are shown respectively in FIGS. 7, 8 and 9. The compensating circuit 6 (FIGS. 6 and 8) and the error factor computing circuit 7 (FIGS. 6 and 7) comprise, in part, the error factor eliminating means 1 shown in FIG. 1. In these figures, "Mul" denotes a multiplier; "Div", a divider; "Sq", a square multiplier for calculating the square of the input signal value; and the numerals in the blocks are the numerical values by which the input signal values are multiplied.

In FIG. 7, the error factor computing circuit 7 receives the coefficients d.sub.0 -d.sub.7 from the respective sample-and-hold (S/H) circuits and provides as output the error factor coefficients described above. As shown in FIGS. 6 and 8, the output from the error factor computing circuit 7 is received by the compensating circuit 6, along with the source signals x and y from the measurement device (not shown), and used to compensate the source signals through the product with Gg, Gk, Gb.sub.x, Gb.sub.y and Gh, as described above, to provide as output G cos .theta. and G sin .theta.. This output is received by the arc-tangent operation circuit 8 (FIGS. 6 and 9), an illustrative circuit for the arc-tangent operating means 2 in FIG. 1, and used to provide the value of the angle .theta.. As shown in FIG. 9, the determination of .theta. can be achieved using a table reference wherein the angle-tangent values are stored, for example, on a read-only memory (ROM), a common technique in the art.

It is, of course, possible to embody the invention in specific forms other than described above without departing from the spirit of the present invention. The embodiment described above is merely illustrative and should not be considered restrictive in any way. The scope of the invention is given in the appended claims, rather than the above description, and all variations and equivalents which fall within the range of the claims are intended to be embraced therein.

Claims

1. An apparatus for compensating for measurement errors included in the output from a measurement device for sensing the angular position .theta. (theta) of a rotating member, the output signals x and y from said device being represented as a sine function and a cosine function of the angle .theta., respectively, said apparatus comprising:

compensating means receiving said output signals x and y and generating approximations of said signals as the average amplitude of the input signals and error factors consisting of at least one of amplitude difference between the two signals x and y, cross-talk, second-order harmonics included in the signal x, second-order harmonics included in the signal y, third-order harmonics included in both signals x and y at the same amount, respectively, zero point of the signal x and zero point of the signal y;

first computational means using the approximations from said compensating means for determining the maximum values of x, x+y, y, -x+y, -x, -x-y, -y and x-y and providing output signals of said maximum values as d.sub.0, d.sub.1, d.sub.2, d.sub.3, d.sub.4, d.sub.5, d.sub.6 and d.sub.7, respectively; and

second computational means for receiving the output of said first computational means and determining the magnitude of the respective error factors based upon the input from said first computational means;

said compensating means using said approximation signals and the magnitude of the error factors from said second computational means, removing said error factors from the x and y signals from said measurement device, or a signal or signals derived from the signals x and y, and providing as outputs error-compensated x and y signals as sin .theta. and cos .theta.,

wherein said apparatus monitors the signals x and y, automatically determines the measurement error from the magnitude of the error factors and compensates the signals accordingly by removing said error components to provide output compensated signals of improved precision.

2. An apparatus as claimed in claim 1, wherein said compensating means provides said approximation x and y signals in accordance with the following relationships:

x=G(cos .theta.+g cos .theta.+k sin .theta.-b.sub.x cos 2.theta.+h cos 3.theta.)+z.sub.x

y=G(sin .theta.-g sin .theta.+k cos .theta.+b.sub.y cos 2.theta.-h sin 3.theta.)+z.sub.y

where G, g, k, b.sub.x, b.sub.y, h, z.sub.x, z.sub.y are coefficients corresponding to the size of respective error factors consisting of average amplitude, amplitude difference between the two signals x and y, cross-talk, second-order harmonics included in the signal x, second-order harmonics included in the signal y, third-order harmonics included in both signals x and y at the same amount respectively, zero point of the signal x and zero point of the signal y; and

said second computational means determines the magnitude of the respective error factors in accordance with the following relationships: ##EQU7##

3. An apparatus as claimed in claim 2, wherein the error factors are due to signal amplitude and zero point fluctuation and said compensating means provides said approximation x and y signals in accordance with the following expressions:

x=g.sub.x cos .theta.+z.sub.x

y=g.sub.y sin .theta.+z.sub.y

where g.sub.x =G(1+g) and g.sub.y =G(1-g);

said first computational means determines the maximum values of x, y, -x and -y and provides output signals of said maximum values as d.sub.0, d.sub.2, d.sub.4 and d.sub.6 ; and

said second computational means determines the magnitude of the respective error factors in accordance with the following expressions:

d.sub.0 =g.sub.x +z.sub.x

d.sub.2 =g.sub.y +z.sub.y

d.sub.4 =g.sub.x -z.sub.x

d.sub.6 =g.sub.y -z.sub.y.

4. An apparatus as claimed in claim 2, wherein the amplitudes of the signals from the measurement device are equal and the error factor is due only to cross-talk, said compensating means provides said approximation x and y signals in accordance with the following relationships:

x=G(cos .theta.+k sin .theta.)+z.sub.x

y=G(sin .theta.+k cos .theta.)+z.sub.y

said first computational means determines the maximum values of x+y, -x+y, -x-y and x-y and provides output signals of said maximum values as d.sub.1, d.sub.3, d.sub.5 and d.sub.7 ; and

said second computational means determines the magnitude of the respective error factors in accordance with following expressions: ##EQU8##

5. An apparatus as claimed in claim 2, wherein the signal to be compensated is .theta. derived from the signals x and y through arc-tangent operation.

6. An apparatus as claimed in claim 2, 3 or 4, wherein said first computational means determines the maximum values d.sub.i on the basis of signal values at .theta.=i .pi./4.+-..delta. where .delta. is an angle whose cosine value can be approximated by 1 to adjust for the speed variation of the rotating member.

7. An apparatus as claimed in claim 6, wherein the maximum values are compensated by exponential smoothing.

8. An apparatus as claimed in claim 6, wherein the compensation is carried out alternately for the maximum values wherein .theta. is different from each other by .pi..

9. An apparatus for compensating for measurement errors included in the output from a measurement device for sensing the angular position .theta. (theta) of a rotating member, the output signals x and y from said device being represented as a sine function and a cosine function of the angle .theta., respectively, said apparatus comprising:

(i) compensating means receiving said signals x and y and generating approximations of said signals as the average amplitude of the input signals and the error factors consisting at least of one of amplitude difference between the two signals x and y, cross-talk, second-order harmonics included in the signal x, second-order harmonics included in the signal y, third-order harmonics included in both signals x and y at the same amount respectively, zero point of the signal x and zero point of the signal y;

(ii) first computational means using the approximations from said compensating means for determining the maximum values of x, x+y, y, -x+y, -x, -x-y, -y and x-y and providing output signals of said maximum values as d.sub.0, d.sub.1, d.sub.2, d.sub.3, d.sub.4, d.sub.5, d.sub.6 and d.sub.7, respectively;

(iii) second computational means for receiving the output of said first computational means and determining the magnitude of the respective error factors based upon the input from said first computational means;

(iv) said compensating means using said approximation signals and the magnitude of the error factors from said second computational means, removing said error factors from the x and y signals from said measurement device, or a signal or signals derived from the signals x and y, and providing as outputs error-compensated x and y signals as sin .theta. and cos .theta.;

(v) arc-tangent operating means receiving the compensated signals from said compensating means for determining the value of .theta. by arc-tangent operation; and

(vi) control means for controlling the operations of said compensating means, said first and second computational means, and said arc-tangent operating means,

wherein said apparatus monitors the signals x and y, automatically determines the measurement error from the magnitude of the error factors and compensates the signals accordingly by removing said error components to provide output compensated signals of improved precision.

10. An apparatus as claimed in claim 9, wherein said compensating means provides said approximation x and y signals in accordance with the following relationship:

x=G(cos .theta.+g cos .theta.+k sin .theta.-b.sub.x cos 2.theta.+h cos 3.theta.)+z.sub.x

y=G(sin .theta.-g sin .theta.+k cos .theta.+b.sub.y cos 2.theta.-h sin 3.theta.)+z.sub.y

where G, g, k, b.sub.x, b.sub.y, h, z.sub.x, z.sub.y are coefficients corresponding to the size of respective error factors consisting of average amplitude, amplitude difference between the two signals x and y, cross-talk, second-order harmonics included in the signal x, second-order harmonics included in the signal y, third-order harmonics included in both signals x and y at the same amount respectively, zero point of the signal x and zero point of the signal y; and

said second computational means determines the magnitude of the respective error factors in accordance with the following relationships: ##EQU9##

11. An apparatus as claimed in claim 10, wherein the error factors are due to signal amplitude and zero point fluctuations and said compensating means provides said approximation x and y signals in accordance with the following expressions:

x=g.sub.x cos .theta.+z.sub.x

y=g.sub.y cos .theta.+z.sub.y

where g.sub.x =G(1+g) and g.sub.y =G(1-g);

said first computational means determines the maximum values of x, y, -x and -y and provides output signals of said maximum values as d.sub.0, d.sub.2, d.sub.4 and d.sub.6 ; and

said second computational means determines the magnitude of the respective error factors in accordance with the following equations:

d.sub.0 =g.sub.x +z.sub.x

d.sub.2 =g.sub.y +z.sub.y

d.sub.4 =g.sub.x -z.sub.x

d.sub.6 =g.sub.y -z.sub.y.

12. An apparatus as claimed in claim 10, wherein the amplitude of the signal from the measurement device are equal and the error factor is due only to cross-talk and said compensating means provides said approximation x and y signals in accordance with the following relationships:

x=G(cos .theta.+k sin .theta.)+z.sub.x

y=G(sin .theta.+k cos .theta.)+z.sub.y

said first computational means determines the maximum values of x+y, -x+y, -x-y and x-y and provides output signals of said maximum values as d.sub.1, d.sub.3, d.sub.5 and d.sub.7 ; and

said second computational means determines the magnitude of the respective error factors in accordance with following expressions: ##EQU10##

13. An apparatus as claimed in claim 9, wherein the signal to be compensated is .theta. provided by said arc-tangent operating means.

14. An apparatus as defined in claim 9, 10 or 11, wherein to adjust for the speed variation of the rotating member said first computational means determines the maximum values d.sub.i on the basis of signal values at .theta.=i .pi./4.+-..delta. where .delta. is an angle whose cosine value can be approximated by 1.

15. An apparatus as defined in claim 14, wherein the maximum values are compensated by exponential smoothing.

16. An apparatus as defined in claim 14, wherein the compensation is carried out alternately for the maximum values wherein .theta. is different by .pi. from each other.