Open AccessArticle

Analysis of the Subdivision Errors of Photoelectric Angle Encoders and Improvement of the Tracking Precision of a Telescope Control System

Jiawei Yu

^1,2,3

Qiang Wang

^1,2,

Guozhong Zhou

^1,2,

Dong He

^1,2,

Yunxia Xia

^1,2,

Xiang Liu

^1,2,

Wenyi Lv

^1,2,3 and

Yongmei Huang

^1,2,*

Institute of Optics and Electronics, Chinese Academy of Sciences, No. 1 Guangdian Road, Chengdu 610209, China

Key Laboratory of Optical Engineering, Chinese Academy of Sciences, Chengdu 610209, China

University of Chinese Academy of Sciences, Beijing 100049, China

Author to whom correspondence should be addressed.

Sensors 2018, 18(9), 2998; https://doi.org/10.3390/s18092998

Submission received: 25 July 2018 / Revised: 20 August 2018 / Accepted: 4 September 2018 / Published: 7 September 2018

(This article belongs to the Section Physical Sensors)

Download

Browse Figures

Figure 1
Distribution law of <math display="inline"><semantics> <mrow> <mo>Δ</mo> <msub> <mi>θ</mi> <mi>D</mi> </msub> </mrow> </semantics></math> in one period of <math display="inline"><semantics> <mrow> <msub> <mi>θ</mi> <mi>r</mi> </msub> </mrow> </semantics></math>. "> Figure 2
Distribution law of <math display="inline"><semantics> <mrow> <mo>Δ</mo> <msub> <mi>θ</mi> <mi>M</mi> </msub> </mrow> </semantics></math> in one period of <math display="inline"><semantics> <mrow> <msub> <mi>θ</mi> <mi>r</mi> </msub> </mrow> </semantics></math>. "> Figure 3
Distribution law of <math display="inline"><semantics> <mrow> <mo>Δ</mo> <msub> <mi>θ</mi> <mi>P</mi> </msub> </mrow> </semantics></math> in one period of <math display="inline"><semantics> <mrow> <msub> <mi>θ</mi> <mi>r</mi> </msub> </mrow> </semantics></math>. "> Figure 4
Distribution law of <math display="inline"><semantics> <mrow> <mo>Δ</mo> <msub> <mi>θ</mi> <mi>N</mi> </msub> </mrow> </semantics></math> in one period of <math display="inline"><semantics> <mrow> <msub> <mi>θ</mi> <mi>r</mi> </msub> </mrow> </semantics></math>. "> Figure 5
Simulation of the tracking error with subdivision errors. "> Figure 6
Telescope platform used for the experiment. "> Figure 7
Internal structure of the photoelectric angle encoder. "> Figure 8
Working process of the angle encoder (including subdivision errors). "> Figure 9
Compensation algorithm for the angle encoder. "> Figure 10
Control system based on classical double-loop structure. "> Figure 11
(a) Distribution of the position error with leading speed of 0.25°/s; (b) Distribution of the position error with leading speed of 1°/s. "> Figure 12
(a) Distribution of the compensated position error with a leading speed of 0.25°/s; (b) Distribution of the compensated position error with a leading speed of 1°/s. "> Figure 13
Adding a feed-forward loop into the control system. "> Figure 14
(a) Distribution of the compensated position error after adding the feed-forward loop with a leading speed of 0.25°/s; (b) Distribution of the compensated position error after adding the feed-forward loop with a leading speed of 1°/s. ">

Versions Notes

Abstract

Photoelectric angle encoders, working as position sensors, have a great influence on the accuracy and stability of telescope control systems (TCS). In order to improve the tracking precision of TCS, a method based on subdivision error compensation for photoelectric angle encoders is proposed. First, a mathematical analysis of six types of subdivision errors (DC error, phase error, amplitude error, harmonic error, noise error, and quantization error) is presented, which is different from the previously used analysis based on the Lissajous figure method. In fact, we believe that a mathematical method is more efficient than the figure method for the expression of subdivision errors. Then, the distribution law and period length of each subdivision error are analyzed. Finally, an error compensation algorithm is presented. In a real TCS, the elevation jittering phenomenon occurs, which indicates that compensating for the amplitude error is necessary. A feed-forward loop is then introduced into the TCS, which is position loop- and velocity loop-closed, leading to a decrease of the tracking error by nearly 54.6%, from 2.31” to 1.05”, with a leading speed of 0.25°/s, and by 40.5%, from 3.01” to 1.79”, with a leading speed of 1°/s. This method can realize real-time compensation and improve the ability of TCS without any change of the hardware. In addition, independently of the environment and the kind of control strategy used, this method can also improve the tracking precision presumably because it compensates the measuring error inside the photoelectric angle encoder.

Keywords:

telescope control systems; photoelectric angle encoders; tracking precision; subdivision errors; amplitude error; real-time compensation

1. Introduction

With the development of free space optical communication, a higher tracking precision of telescope control systems (TCS), which usually display arc second-scale accuracy [1,2], is a necessity. The photoelectric angle encoder is a kind of angular sensor, which achieves high resolution by the subdivision of two-way orthogonal moire pattern (sine and cosine) signals from a fine code disk [3]. Obviously, high-precision angle encoders contribute to the high accuracy of TCS. However, the real photoelectric parameters are not well matched with the processing circuit, and the normal subdivision methods are based on standard sine and cosine signals, which causes subdivision errors [4].

Scientists are working to manufacture high-precision and high-resolution photoelectric angle encoders. St. Petersburg State Electrotechnical University developed an encoder error detection system with an accuracy of 0.1” at the speed of up to 10 Round Per Second (RPS), which can realize dynamic error detection [5]. Wan et al. from Changchun Institute of Optics, Fine Mechanical, and Physics, Chinese Academy of Sciences, designed a high-precision photoelectric angle encoder with a dual-reading system, which satisfies the technique requirements of spaceborne equipment [6]. However, all these high-resolution angle measurements depend on large code disks. It is a challenge for scientists to achieve high resolution with small code disks [7].

To compensate the subdivision error of the photoelectric angle encoders, whose code disks are small, Wang et al. proposed an error compensation method based on an improved Back Propagation (BP) neural network, which brought a 15% improvement in accuracy [8]. An experimental–theoretical method was proposed by E. V. Kaplya to identify and compensate static errors in the phase coordinates of the elements of an optical encoder and decreased the error of the encoder by 10.8% [9]. Yu’s team designed a 23-bit angle reference encoder by structure and circuit devising, whose subdivision error was then compensated by a Radical Basis Function (RBF) neural network. The resolution of the angle encoder was 0.15”, and, after compensating, the error of the encoder decreased from 1.3” to 0.6” [10]. An automatic compensation method for sine deviation, proposed by Gao and his cooperators, used particle swam optimization to identify the seven undefined parameters of the signal model and could reduce the subdivision error from 0.923” to 0.316” [11]. A method based on moire triangle subdivision error compensation, which was proposed by Sun, decreased the original error by 1/3 [12].

However, all the studies mentioned above have their shortcomings. For those based on neural or improved neural networks, their training sets are too big and need a long processing time before they can be implemented. Besides, when working in real systems, neural network-based methods need a considerable time to match the data input with the training sets. This is not practical for spaceborne equipment. For some other methods, the compensating strategies are difficult to adapt to changing working sites.

To improve the precision of photoelectric encoders, we designed a compensation method without changing the parameters of the hardware. First, we defined the position error of the encoder and identified the type of subdivision error corresponding to the main error. Then, we compensated the main error through the properties of the angle encoder. Finally, we improved the structure of the control system to further increase its precision. This compensation method required a short time and significantly increased the precision of the system. Most importantly, this method works in real time and can satisfy the requirements of spaceborne equipment.

2. Mathematical Analysis on Subdivision Errors

Ideal grating fringe photoelectric signals are sine and cosine signals with equal amplitude, which can be expressed as follow [11]:

u_{a} = A \sin θ, u_{b} = A \cos θ,

(1)

where

A

is the amplitude of the signal, and

θ

is the theoretic rotary angle of the signal; the angle

θ

can be deduced from:

θ = \arctan \frac{u_{a}}{u_{b}},

(2)

However, some factors influence the precision of subdivision, for example, the precision of the code disk, the precision of the axis, the quality of the photoelectric signals, the performance of the processing circuit, and so on. The real outputs are obtained by the moire fringe technique, which can be expressed as:

u_{a} = A_{0} + A_{1} \sin (θ + φ_{1}) + \sum_{i = 2}^{\infty} A_{i} \sin (i θ + φ_{i a}) + δ_{a},

(3)

u_{b} = B_{0} + B_{1} \sin (θ + φ_{2}) + \sum_{i = 2}^{\infty} B_{i} \sin (i θ + φ_{i b}) + δ_{b},

(4)

where the two signals consist of four parts.

A₀ and B₀ are DC components, which cause the DC subdivision error. A₁ and B₁ are the amplitudes of base waves, which cause the amplitude subdivision error.

\sum_{i = 2}^{\infty} A_{i} \sin (i θ + φ_{i a})

and

\sum_{i = 2}^{\infty} B_{i} \sin (i θ + φ_{i b})

are the sum of higher harmonics, which cause the harmonic subdivision error.

δ_{a}

and

δ_{b}

are the electric noise, which causes the electric noise subdivision error. In addition, the analog-to-digital conversion (A/D) causes the quantization subdivision error. Also, the difference between the phases

φ_{1}

and

φ_{2}

causes the phase subdivision error [8].

Here, we define

Δ θ

as the subdivision error,

θ_{r}

as the real rotary angle, and

θ_{d}

as the measured rotary angle. We have, then, the following equation:

Δ θ = θ_{r} - θ_{d},

(5)

2.1. DC Subdivision Error Analysis

If there is only the DC subdivision error, we can introduce the DC components

a, b

into the expression of the subdivision error

Δ θ

and then simplify it as shown:

Δ θ = θ_{r} - \arctan \frac{r \sin θ_{r} + a}{r \cos θ_{r} + b},

(6)

so : \tan Δ θ = \tan (θ_{r} - \arctan \frac{\sin θ_{r} + \frac{a}{r}}{\cos θ_{r} + \frac{b}{r}}),

(7)

We can simplify the equation above:

\tan Δ θ = \frac{b \sin θ_{r} - a \cos θ_{r}}{r + b \cos θ_{r} + a \sin θ_{r}},

(8)

Here, we define:

\cos φ = \frac{b}{\sqrt{a^{2} + b^{2}}}, \sin φ = \frac{a}{\sqrt{a^{2} + b^{2}}},

(9)

We can use zero shift of the angle (initial angle

α

) to calculate

φ

α = \arctan (\frac{a}{1 + b}),

(10)

Then:

\tan Δ θ = \frac{\sin (θ_{r} - φ)}{\frac{r}{\sqrt{a^{2} + b^{2}}} + \cos (θ_{r} - φ)}, φ = \arctan \frac{a}{b},

(11)

Finally, we can get:

Δ θ = \arctan (\frac{\sin (θ_{r} - φ)}{\frac{r}{\sqrt{a^{2} + b^{2}}} + \cos (θ_{r} - φ)}) = \arctan (\frac{\sin (θ_{d} + Δ θ - φ)}{\frac{r}{\sqrt{a^{2} + b^{2}}} + \cos (θ_{d} + Δ θ - φ)})

(12)

According to the equation above, we can use the measured rotary angle

θ_{d}

and

α

to get the DC subdivision error

Δ θ

and compensate it.

Then, we can determine the change of

Δ θ

when

θ_{r}

ranges from 0 to 2π and obtain the distribution law and period length of the DC subdivision error in Figure 1.

Here, we consider as an example the DC components

a = 0.2 r, b = 0.2 r .

Each type of error has a special distribution law, which can help us to identify the main error. Here, we can see that, in one period of

θ_{r}

Δ θ

is a single period.

2.2. Magnitude Subdivision Error Analysis

If there is only the magnitude subdivision error, we can introduce the magnitude components

d r

into the expression of the subdivision error

Δ θ

and then simplify it as follows:

Δ θ = θ_{r} - \arctan \frac{(r - d r) \sin θ_{r}}{r \cos θ_{r}},

(13)

so : \tan Δ θ = \frac{\tan θ_{r} - \frac{(r - d r) \sin θ_{r}}{r \cos θ_{r}}}{1 + \tan θ_{r} \frac{(r - d r) \sin θ_{r}}{r \cos θ_{r}}},

(14)

We can finally simplify it:

Δ θ = \arctan (\frac{\sin 2 θ_{r}}{\cos 2 θ_{r} + \frac{2 r}{d r} - 1}) = \arctan (\frac{\sin (2 θ_{d} + 2 Δ θ)}{\cos (2 θ_{d} + 2 Δ θ) + \frac{2 r}{d r} - 1}),

(15)

According to the equation above, we can use the measured rotary angle

θ_{d}

and

\frac{d r}{r}

to obtain the magnitude subdivision error

Δ θ

and compensate it.

Then, we can determine the change of

Δ θ

when

θ_{r}

ranges from 0 to 2π and obtain the distribution law and the period length of the magnitude subdivision error in Figure 2.

Here, we consider the magnitude component

d r = 0.2 r

We can see that, in one period of

θ_{r}

Δ θ

is a dual period.

2.3. Phase Subdivision Error Analysis

If there is only the phase subdivision error, we can introduce the phase component

φ_{3}

into the expression of the subdivision error

Δ θ

and then simplify it as follows:

Δ θ = θ_{r} - \arctan \frac{r \sin (θ_{r} + φ_{3})}{r \cos θ_{r}},

(16)

next : \arctan (\frac{r \sin (θ_{r} + φ_{3})}{r \cos θ_{r}}) = θ_{r} - Δ θ = θ_{d},

(17)

so : Δ θ = \arctan (\frac{\tan θ_{d}}{\cos φ_{3}} - \tan φ_{3}) - θ_{d},

(18)

According to the equation above, we can use the measured rotary angle

θ_{d}

and

φ_{3}

to get the phase subdivision error

Δ θ

and compensate it.

Then, we can determine the change of

Δ θ

when

θ_{r}

ranges from 0 to 2π and obtain the distribution law and period length of the phase subdivision error in Figure 3.

Here, we consider the phase components

φ_{3} = 10 °

We can see that, in one period of

θ_{r}

Δ θ

is a dual period.

2.4. Harmonic Subdivision Error Analysis

If only the harmonic subdivision error is present, we can introduce the harmonic component

N θ_{r}

into the expression of the subdivision error

Δ θ

and then simplify it as follows:

Δ θ = θ_{r} - \arctan \frac{A_{1} \sin θ_{r} + \sum_{N = 2}^{\infty} A_{N} \sin (N θ_{r})}{A_{1} \cos θ_{r} + \sum_{N = 2}^{\infty} A_{N} \cos (N θ_{r})},

(19)

next : \tan Δ θ = \frac{\tan θ_{r} - \frac{A_{1} \sin θ_{r} + \sum_{N = 2}^{\infty} A_{N} \sin (N θ_{r})}{A_{1} \cos θ_{r} + \sum_{N = 2}^{\infty} A_{N} \cos (N θ_{r})})}{1 + \tan θ_{r} \frac{A_{1} \sin θ_{r} + \sum_{N = 2}^{\infty} A_{N} \sin (N θ_{r})}{A_{1} \cos θ_{r} + \sum_{N = 2}^{\infty} A_{N} \cos (N θ_{r})})},

(20)

so : Δ θ = \arctan [- \frac{\sum_{N = 2}^{\infty} A_{N} \sin ((N - 1) (θ_{d} + Δ θ))}{A_{1} + \sum_{N = 2}^{\infty} A_{N} \cos ((N - 1) (θ_{d} + Δ θ))}],

(21)

According to Equation (21) above, we can see that the harmonic subdivision error consists of different harmonic components, so it is not as simple as the subdivision errors discussed above.

2.5. Noise Subdivision Error Analysis

If only the harmonic subdivision error is present, we can introduce the noise components

δ_{a}, δ_{b}

into the expression of the subdivision error

Δ θ

and then simplify it as follows:

Δ θ = θ_{r} - \arctan \frac{r \sin θ_{r} + δ_{a}}{r \cos θ_{r} + δ_{b}},

(22)

next : \tan Δ θ = \frac{\tan θ_{r} - \frac{\sin θ_{r} + δ_{a}}{\cos θ_{r} + δ_{b}}}{1 + \tan θ_{r} \frac{\sin θ_{r} + δ_{a}}{\cos θ_{r} + δ_{b}}},

(23)

\begin{array}{l} so : Δ θ = \arctan [\frac{\sqrt{δ_{a}^{2} + δ_{b}^{2}} \sin (θ_{d} + Δ θ - λ)}{r + \sqrt{δ_{a}^{2} + δ_{b}^{2}} \cos (θ_{d} + Δ θ - λ)}], \\ here : λ = \arctan \frac{δ_{a}}{δ_{b}} . \end{array}

(24)

Especially when

δ_{a} = δ_{b} = δ,

Finally : Δ θ = \arctan [\frac{\sqrt{2} δ \sin (θ_{d} + Δ θ - \frac{π}{4})}{r + \sqrt{2} δ \cos (θ_{d} + Δ θ - \frac{π}{4})}],

(25)

According to the equation above, we can use the measured rotary angle

θ_{d}

and

δ_{a}, δ_{b}

to get the subdivision error

Δ θ

and compensate it.

Then, we can determine the change of

Δ θ

when

θ_{r}

ranges from 0 to 2π and obtain the distribution law and period length of the noise subdivision error in Figure 4.

Here, we consider the noise components

δ_{a} = 0.05 r, δ_{b} = 0.08 r .

We can see that, in one period of

θ_{r}

Δ θ

is a single period. In almost all spaceborne equipment, the noise should be a minimal component, which can be neglected.

2.6. Quantization Subdivision Error Analysis

As we know, the digital code

X

of a photoelectric angle encoder is uniformly distributed, which means the probability density function of

X

is as follow [8]:

f (x) = {\begin{cases} \frac{1}{b - a} & a < x < b \\ 0 & o t h e r w i s e \end{cases}} .

(26)

We can obtain the mathematical expectation and variance:

E (x) = \int_{a}^{b} x \frac{1}{b - a} d x = \frac{a + b}{2},

(27)

D (x) = E (x^{2}) - {[E (x)]}^{2} = \frac{{(b - a)}^{2}}{12},

(28)

When the resolution of the angle code is

Q

, the quantization subdivision error can be expressed as [8]:

σ_{x} = \frac{Q}{\sqrt{12}} \approx \frac{Q}{3.46},

(29)

Here, the quantization subdivision error is the root-mean-square (RMS) error [8], which can be neglected in a standard electric system.

After obtaining the main subdivision error and its expression, we should think of how to compensate it.

All the subdivision errors are listed in the following Table 1:

3. Compensation Algorithm and Experiment Set Up

3.1. Simulations

Once completed the above analysis, we set up a simulation. We assumed three kinds of subdivision error in the system to determine the tracking error before and after compensation. In Figure 5, we considered a sine input

x = 0.5 * \sin (0.05 * t) - 10

at the center of minus 10 and added three kinds of subdivision error (amplitude error:

d r = 0.3 r

, DC error:

a = 0.2, b = 0

, and phase error:

φ = 10 °

) in the output to evaluate the tracking error. The transfer function

G (s) = \frac{17.5}{{7.842}^{- 5} s^{2} + 0.02143 s + 1}

and the controller

C_{P I} (s) = 0.1943 * \frac{s + 62.98}{s}

were from reference [13]. We determined that the tracking error was 11.15”. Then, we compensated for the amplitude subdivision error, which is the main error, and obtained a tracking error of 6.83”. This means that, after compensating for the main subdivision error, the tracking error decreased of 38.7%. We were not concern about the type of control strategy, because the tracking error would decrease as a result of the measuring error being compensated. In other words, regardless of the type of control strategy, this method can also improve the tracking precision presumably because it compensates the measuring error inside the photoelectric angle encoder.

3.2. Compensation Algorithm

The experimental telescope platform and photoelectric angle encoder we used are shown in Figure 6 and Figure 7.

To compensate the subdivision error, we should know the structure of the angle encoder and the algorithm to calculate the angle. In Figure 7, we see the code disk of the absolute photoelectric angle encoder, on which there are two coarse-reading heads and two fine-reading heads (indicated in Figure 7 with green and blue lines, respectively). We could read the coarse code and fine code through these reading heads. We removed the processing circuit chip lying on the code disk in order to examine the details of the code disk. The fine code is generated by electronic subdivision, which is processed in the DSP (Digital Signal Processor). The coarse code represents the completed period numbers of subdivision signals A and B the shaft has rotated. The subdivision of A and B is described by the fine code, which includes the subdivision errors [14].

The movement of the telescope shaft is transformed to optical signals through the code disk, on which we could obtain the coarse signals and fine signals (electric signals) by the reading heads. In Figure 8, the processing details are discussed. The coarse signals and fine signals were sent into amplifiers and filters to ensure their purity. Afterwards, the analog signals were converted into digital signals by A/D sampling and were then sent into the DSP. After a series of processes, we could determine the measuring angle. However, subdivision errors were present in the angle we obtained. Therefore, a compensation algorithm was required.

In Figure 9, we see the compensation algorithm. First, according to Section 2, we identified the main kind of subdivision error and calculated the compensation angle. In the following Equations,

Q

is the resolution of the angle encoder,

Q_{C}

is the resolution represented by the coarse code,

Q_{F}

is the resolution represented by the fine code. The total bits of the angle encoder are the combination of the coarse code bits and the fine code bits [14]. These parameters are expressed as follows:

Q = {(360 / 2^{T o t a l B i t s})}^{°},

(30)

Q_{C} = {(360 / 2^{C o a r s e B i t s})}^{°},

(31)

Q_{F} = {(360 / 2^{F i n e B i t s})}^{°},

(32)

T o t a l B i t s = C o a r s e B i t s + F i n e B i t s,

(33)

C o d e A l l = C o a r s e C o d e * 2^{F i n e B i t s} + F i n e C o d e,

(34)

F i n e C o d e C o r r e c t e d = F i n e C o d e + Δ θ * Q_{F},

(35)

C o d e A l l C o r r e c t e d = C o a r s e C o d e * 2^{F i n e B i t s} + F i n e C o d e C o r r e c t e d,

(36)

θ_{d} = C o d e A l l * Q,

(37)

θ_{r} = C o d e A l l C o r r e c t e d * Q,

(38)

The processing steps can be described as follows:

Step 1: Obtain the data of coarse code and fine code.
Step 2: According to the analysis of the subdivision errors, identify the main kind of subdivision error and calculate $Δ θ$ .
Step 3: Calculate the fine code corrected according to Equation (35).
Step 4: Finally, calculate $θ_{r}$ through Equations (36) and (38).

Here, we can see that, if we only determined and calculated the subdivision error

Δ θ

, the real rotary angle

θ_{r}

is easily obtained. However, it is not obvious how to identify the main kind of subdivision error and calculate it.

3.3. Experiment Set Up

In a real telescope system, the elevation jittering cannot be neglected. Compared with the elevation error, the azimuth error was low, i.e., below the expected value (2″). Therefore, we only consider the elevation as the control target. We used a system with the classical double-loop structure (position loop and velocity loop) [15], as shown in Figure 10.

4. Experimental Results

4.1. Results before Adding the Compensation Algorithm

The elevation position error, before we applied our compensation algorithm, is shown in Figure 11a,b. The leading speed of the elevation position errors was 0.25°/s and 1°/s. These values were equivalent to sinusoidal guidance whose amplitude and frequency are 1.99°, 0.02 Hz, and 1.99°, 0.08 Hz.

Here, we recorded the leading speed v, time t, and feedback position

θ_{d}

measured by the photoelectric angle encoder. We expressed the leading position as

θ_{L} = v * t

, so the position error could be expressed as

p o s i t i o n e r r o r = θ_{L} - θ_{d}

. Finally, we considered Error as the RMS (root mean square) of the position error,

E r r o r = {(\sum_{i = 1}^{N} {(p o s i t i o n e r r o r (i))}^{2} / N)}^{1 / 2}

. We could calculate the Error of the two kinds of speed:

E r r o r_{v = 0.25} = {2.31}^{″}

E r r o r_{v = 1} = {3.01}^{″}

. Obviously, the elevation error did not satisfy our expectation. Thus, error compensation was necessary.

4.2. Results after Applying the Compensation Algorithm

From the analysis in Section 2, we know there are six kinds of subdivision errors that decrease the precision of the angle encoder. However, there must be a main subdivision error out of the six kinds of subdivision errors. To ensure real-time compensation, our algorithm has to compensate only for the main subdivision error on the basis of our analysis above.

The angle encoder we used was 23-bit. More precisely, there were 13 bits for the coarse code and 10 bits for the fine code. According to the expression (31), we determined the resolution of the coarse code,

Q_{C} = 360 ° / 2^{13} = 0.0439 ° = {158.04}^{″}

. Together with the leading speed v, we determined the coarse period time

T_{C} = Q_{C} / v

. Then, we recorded the coarse cycles according to the coarse period time

T_{C}

, and the subdivision error cycles according to Figure 11a,b in one second.

For

v = 0.25 ° / s

T_{C} = (0.0439 / 0.25) s = 0.176 s

. From Figure 11a, we could count the periods. We found that, in one second, there were

1 / 0.176 = 5.68

periods for the coarse code and 12 periods for the subdivision error. This means that there were approximately 2.1 periods for the subdivision error in one period of the coarse code.

For

v = 1 ° / s

T_{C} = (0.0439 / 1) s = 0.0439 s

. From Figure 11b, we could count the periods. We found that, in one second, there were

1 / 0.0439 = 22.78

periods for the coarse code and 46 periods for the subdivision error. This means that there were approximately 2.0 periods for the subdivision error in one period of the coarse code.

We observed that the magnitude subdivision error and phase subdivision error are dual periods that are close to the values we calculated. We also found that the distribution of the error was symmetrical with respect to the center of zero, which suggested that the magnitude subdivision error was the main error.

As we reported, the magnitude subdivision error is:

Δ θ = θ_{r} - \arctan \frac{(r - d r) \sin θ_{r}}{r \cos θ_{r}},

so : Δ θ = \arctan (\frac{\sin 2 θ_{r}}{\cos 2 θ_{r} + \frac{2 r}{d r} - 1}) = \arctan (\frac{\sin (2 θ_{d} + 2 Δ θ)}{\cos (2 θ_{d} + 2 Δ θ) + \frac{2 r}{d r} - 1}),

We defined

\frac{d r}{r}

and introduced all the required parameters into the TCS, using the compensation algorithm described in Figure 9. We considered v = 0.25°/s and v = 1°/s after compensation, as indicated below.

As shown in Figure 12a,b, the position errors were clearly smaller after compensation. The RMS of the position error after compensation were

E r r o r_{v = 0.25} = {1.15}^{″}

and

E r r o r_{v = 1} = {1.92}^{″}

. This means that the errors decreased, respectively, of 50.2% and 36.2%. Therefore, the compensation algorithm considerably increased the precision.

4.3. Results after Introducing the Feed-Forward Loop into the TCS

In order to further increase the tracking precision, we added a feed-forward loop into the system (see Figure 13). A feed-forward loop can introduce the object properties into the loop and reduce the error.

After adding the feed-forward loop, the RMS of the position errors were

E r r o r_{v = 0.25} = {1.05}^{″}

and

E r r o r_{v = 1} = {1.79}^{″}

. The feed-forward loop further decreased the previously compensated errors of about 9.5% and 6.8%. Because of the characteristics of the feed-forward loop, if the frequency of the error was lower, the error decrease would be greater [16]. The distribution of the error after adding the feed-forward loop was as follows (Figure 14a,b).

5. Conclusions

In a real TCS, the elevation jittering phenomenon occurs. We considered the elevation as the control target. After analyzing the period of the position error, we found that magnitude subdivision exists, using the analysis in Section 2. After compensation, the tracking error was decreased by nearly 50.2%, from 2.31” to 1.15” with the leading speed of 0.25°/s, and by 36.2%, from 3.01” to 1.92”, with the leading speed of 1°/s. We observed that, when the speed increased, the tracking error increased, which suggested that when the system worked at a lower speed, the compensation algorithm worked better. We then found that adding a feed-forward loop into the system could further decrease the tracking error, even if only by a few percent points (9.5% and 6.8%). The reason of the small decrease was that the frequency of the error was so high that the feed forward-loop could not run perfectly. Therefore, in the future, we will try to improve the functioning of the feed-forward loop.

The compensation algorithm described above is suitable for complex devices and can provide a real-time compensation without any change of the hardware. It is also adaptable to changing environments, regardless of the control strategy used, and can also improve the tracking precision presumably because it compensate the measuring error inside of the photoelectric angle encoder. Moreover, the compensation method described only requires the position of the shaft that can be easily determined by the photoelectric angle encoder; therefore, its cost is limited, which means it works efficiently.

However, the compensation algorithm has its shortcomings, since it cannot compensate all subdivision errors in the same time. In fact, we observed that other subdivision errors were still present after compensation, and that we only compensated the main error. Therefore, although the compensation algorithm worked well, it still requires improvements. We will work on these improvements in the future to enhance the system’s efficiency.

Author Contributions

Conceptualization, J.Y. and Q.W.; Methodology, J.Y. and Q.W.; Software, J.Y., D.H., Y.X., W.L.; Validation, G.Z. and Y.H.; Formal Analysis, J.Y. and Q.W.; Investigation, J.Y. and Q.W.; Resources, D.H., Y.X., X.L.; Data Curation, J.Y. and Q.W.; Writing-Original Draft Preparation, J.Y.; Writing-Review & Editing, J.Y.; Visualization, J.Y., D.H., Y.X., X.L.; Supervision, G.Z., Y.H.; Project Administration, G.Z., Y.H.

Funding

This research received no external funding.

Acknowledgments

Thanks for my tutor and all the colleagues in the Key Laboratory of Optical Engineering, Chinese Academy of Sciences.

Conflicts of Interest

The authors declare no conflict of interest.

References

Jian, Z.; Du, F. The servo control system of KDUST telescope. In Proceedings of the Advances in Optical and Mechanical Technologies for Telescopes and Instrumentation, Montréal, QC, Canada, 18 June 2014; pp. 112–117. [Google Scholar]
Dindar, M.; Helhel, S.; Esenoğlu, H.; Parmaksızoğlu, M. A new software on TUG-T60 autonomous telescope for astronomical transient events. Exp. Astron. 2015, 39, 21–28. [Google Scholar] [CrossRef]
Lopez, J.; Artes, M. A New Methodology for Vibration Error Compensation of Optical Encoders. Sensors 2012, 12, 4918–4933. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Wang, X.-J. Errors and precision analysis of subdivision signals for photoelectric angle encoders. Opt. Precis. Eng. 2012, 20, 379–386. [Google Scholar] [CrossRef]
Filatov, Y.V.; Agapov, M.Y.; Bournachev, M.N.; Loukianov, D.P.; Pavlov, P.A. Laser goniometer systems for dynamic calibration of optical encoders. In Proceedings of the Optical Measurement Systems for Industrial Inspection III, Munich, Germany, 30 May 2003; Volume 5144, pp. 381–390. [Google Scholar]
Wan, Q.H.; Sun, Y.; Wang, S.J.; She, R.; Lu, X.R.; Liang, L.H. Design for spaceborne absolute photoeletric encoder of dual numerical system. Opt. Previs. Eng. 2009, 17, 52–57. [Google Scholar]
Yu, H.; Wan, Q.; Zhao, C.; Liang, L.; Du, Y. A High-resolution Subdivision Algorithm for Photographic Encoders and Its Error Analysis. Acta Opt. Sin. 2017, 37. (In Chinese) [Google Scholar] [CrossRef]
Wang, X.-G.; Cai, T.; Deng, F.; Xu, L.S. Error compensation of photoelectric encoder based on improved BP neural network. In Proceedings of the 24th Chinese Control and Decision Conference (CCDC), Taiyuan, China, 23–25 May 2012; Volume 2074, pp. 3941–3946. [Google Scholar]
Kaplya, E.V. Identification and compensation of static error in the phase coordinates of elements of an optical encoder. Meas. Tech. 2014, 57, 493–500. [Google Scholar] [CrossRef]
Yu, H.; Wan, Q.H.; Wang, S.J.; Lu, X.R.; Du, Y.C. High-precision realtime angle reference in dynamic measurement of photoelectric encoder. Chin. Opt. 2015, 3, 020. [Google Scholar]
Gao, X.; Wan, Q.; Lu, X. Automatic compensation of sine deviation for grating fringe photoelectric signal. Acta Opt. Sin. 2013, 33, 0712001. [Google Scholar] [CrossRef]
Sun, Y. Study on the Interpolation Error Correlation Method for the Minitype spaceflight Photoelectric Encoder. Ph.D. Thesis, University of Chinese Academy of Sciences, Shenzhen, China, 2010. (In Chinese). [Google Scholar]
Xia, Y.-X. Study on Control Techniques of Inertial Stabilization on Moving Platform; Institute of Optics and Electronics. Ph.D. Thesis, Chinese Academy of Sciences, Beijing, China, 2013. (In Chinese). [Google Scholar]
Su, Y.; Wang, Q.; Yan, F.; Liu, X.; Huang, Y. Subdivision error analysis and compensation for photoelectric angle encoder in a telescope control system. Math. Probl. Eng. 2015, 2015, 967034. [Google Scholar] [CrossRef]
Wang, Q.; Cai, H.X.; Huang, Y.M.; Ge, L.; Tang, T.; Su, Y.R.; Liu, X.; Li, J.Y.; He, D.; Du, S.P.; et al. Acceleration feedback control (AFC) enhanced by disturbance observation and compensation (DOC) for high precision tracking in telescope systems. Res. Astron. Astrophys. 2016, 16, 124. [Google Scholar] [CrossRef]
Tang, T. Acceleration Feedback Control for High precision Tracking Control. Ph.D. Thesis, Institute of Optics and Electronics, Chinese Academy of Sciences, Beijing, China, 2009. (In Chinese). [Google Scholar]

Figure 1. Distribution law of

Δ θ_{D}

in one period of

θ_{r}

Figure 1. Distribution law of

Δ θ_{D}

in one period of

θ_{r}

Figure 2. Distribution law of

Δ θ_{M}

in one period of

θ_{r}

Figure 2. Distribution law of

Δ θ_{M}

in one period of

θ_{r}

Figure 3. Distribution law of

Δ θ_{P}

in one period of

θ_{r}

Figure 3. Distribution law of

Δ θ_{P}

in one period of

θ_{r}

Figure 4. Distribution law of

Δ θ_{N}

in one period of

θ_{r}

Figure 4. Distribution law of

Δ θ_{N}

in one period of

θ_{r}

Figure 5. Simulation of the tracking error with subdivision errors.

Figure 6. Telescope platform used for the experiment.

Figure 7. Internal structure of the photoelectric angle encoder.

Figure 8. Working process of the angle encoder (including subdivision errors).

Figure 9. Compensation algorithm for the angle encoder.

Figure 10. Control system based on classical double-loop structure.

Figure 11. (a) Distribution of the position error with leading speed of 0.25°/s; (b) Distribution of the position error with leading speed of 1°/s.

Figure 12. (a) Distribution of the compensated position error with a leading speed of 0.25°/s; (b) Distribution of the compensated position error with a leading speed of 1°/s.

Figure 13. Adding a feed-forward loop into the control system.

Figure 14. (a) Distribution of the compensated position error after adding the feed-forward loop with a leading speed of 0.25°/s; (b) Distribution of the compensated position error after adding the feed-forward loop with a leading speed of 1°/s.

Table 1. Definition of the subdivision errors.

Kind of Errors	Expression of Errors
DC error	$\arctan (\frac{\sin (θ_{d} + Δ θ - φ)}{\frac{r}{\sqrt{a^{2} + b^{2}}} + \cos (θ_{d} + Δ θ - φ)})$
Magnitude error	$\arctan (\frac{\sin (2 θ_{d} + 2 Δ θ)}{\cos (2 θ_{d} + 2 Δ θ) + \frac{2 r}{d r} - 1})$
Phase error	$\arctan (\frac{\tan θ_{d}}{\cos φ_{3}} - \tan φ_{3}) - θ_{d}$
Harmonic error	$\arctan [- \frac{\sum_{N = 2}^{\infty} A_{N} \sin ((N - 1) (θ_{d} + Δ θ))}{A_{1} + \sum_{N = 2}^{\infty} A_{N} \cos ((N - 1) (θ_{d} + Δ θ))}]$
Noise error	$\arctan [\frac{\sqrt{δ_{a}^{2} + δ_{b}^{2}} \sin (θ_{d} + Δ θ - λ)}{r + \sqrt{δ_{a}^{2} + δ_{b}^{2}} \cos (θ_{d} + Δ θ - λ)}]$
Quantization error	×

© 2018 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Yu, J.; Wang, Q.; Zhou, G.; He, D.; Xia, Y.; Liu, X.; Lv, W.; Huang, Y. Analysis of the Subdivision Errors of Photoelectric Angle Encoders and Improvement of the Tracking Precision of a Telescope Control System. Sensors 2018, 18, 2998. https://doi.org/10.3390/s18092998

AMA Style

Yu J, Wang Q, Zhou G, He D, Xia Y, Liu X, Lv W, Huang Y. Analysis of the Subdivision Errors of Photoelectric Angle Encoders and Improvement of the Tracking Precision of a Telescope Control System. Sensors. 2018; 18(9):2998. https://doi.org/10.3390/s18092998

Chicago/Turabian Style

Yu, Jiawei, Qiang Wang, Guozhong Zhou, Dong He, Yunxia Xia, Xiang Liu, Wenyi Lv, and Yongmei Huang. 2018. "Analysis of the Subdivision Errors of Photoelectric Angle Encoders and Improvement of the Tracking Precision of a Telescope Control System" Sensors 18, no. 9: 2998. https://doi.org/10.3390/s18092998

APA Style

Yu, J., Wang, Q., Zhou, G., He, D., Xia, Y., Liu, X., Lv, W., & Huang, Y. (2018). Analysis of the Subdivision Errors of Photoelectric Angle Encoders and Improvement of the Tracking Precision of a Telescope Control System. Sensors, 18(9), 2998. https://doi.org/10.3390/s18092998

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu

Analysis of the Subdivision Errors of Photoelectric Angle Encoders and Improvement of the Tracking Precision of a Telescope Control System

Abstract

1. Introduction

2. Mathematical Analysis on Subdivision Errors

2.1. DC Subdivision Error Analysis

2.2. Magnitude Subdivision Error Analysis

2.3. Phase Subdivision Error Analysis

2.4. Harmonic Subdivision Error Analysis

2.5. Noise Subdivision Error Analysis

2.6. Quantization Subdivision Error Analysis

3. Compensation Algorithm and Experiment Set Up

3.1. Simulations

3.2. Compensation Algorithm

3.3. Experiment Set Up

4. Experimental Results

4.1. Results before Adding the Compensation Algorithm

4.2. Results after Applying the Compensation Algorithm

4.3. Results after Introducing the Feed-Forward Loop into the TCS

5. Conclusions

Author Contributions

Funding

Acknowledgments

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI