Research Article  Open Access
Matej Šalamon, Bojan Jarc, "Macromodel of Precise SignalPhase Meter", Modelling and Simulation in Engineering, vol. 2012, Article ID 856497, 7 pages, 2012. https://doi.org/10.1155/2012/856497
Macromodel of Precise SignalPhase Meter
Abstract
This paper presents a precise macromodel of a signalphase meter, which allows continuous phase measurement during simulation. It has been developed as a support tool during the design process of a signalconditioning circuit for incremental position encoders. The development of a signal conditioning circuit requires precise measurements of small signal phases, amplitudes and offsets using the analog/digital circuit simulator. The phase measurement cannot be performed directly with a simulator, therefore an appropriate macromodel is needed for a circuit simulator. The structure of the signalphase meter is based on the conventional signalphase measuring method and is intended for the measuring of a cosinesignal phase with a known frequency. It recommends that the time variations of an input signal’s parameters (amplitude, phase, frequency, and offset voltage) are slow and small as possible. Rapid change of a signal’s parameters decreases the simulation result’s accuracy. A macromodel’s precision mainly depends on the chosen parameters for the macromodel and for the simulation. We show that with the proposed meter’s model, the phase angle can be measured with an accuracy of more than ±0,02%.
1. Introduction
Position encoders are sensors that convert their linear or rotary displacements into appropriate electrical signals and are widely used in many computercontrolled machines. Position signals can be generated by magnetic, capacitive, or optical means, and represent displacement information in absolute or incremental forms [1, 2].
Ideally, the linear optical incremental position encoders produce two 90° shifted analog signals, which enable us to detect the position of an encoder’s head. However, due to encoder’s and optical grating’s mechanical limitations, these signals have unideal sine/cosine shapes, variable amplitudes, offset voltages, and additional phase shifts. The influences of such various unwanted effects are reduced and corrected by a signalconditioning circuit, which is usually an encoder’s integral part, realized as an integrated circuit. Thus, the signalconditioning circuit is an analog/digital circuit, generally developed using an analog/digital circuit simulator. The development of such circuit, always requires very precise signalphase measurements, which cannot be performed directly by a simulator.
This paper deals with the SPICE macromodel of signalphase meter, which enables a precise phase measurement whilst simulation is running. It is based on wellknown practically implemented conventional phasemeasurement techniques, as described in [3–7]. Although structures, similar to our model, can be found in [8], to our knowledge, this type of model has never been developed.
This paper is organized as follows. Section 2 briefly describes the conventional signalphase measurement techniques. A detailed description of the proposed macromodel is given in Section 3. Section 4 comprises examples of signalphase measurement using the proposed model and the SPICE simulator. This model’s accuracy and other properties are also described. The conclusions are given in Section 4.
2. Conventional SignalPhase Measurement Techniques
Conventional phasemeasurement techniques between two signals of the same frequency [3] usually apply a square waveconversion to signals. Phase is then determined by measuring the time difference, either between the zero crossingpoints or between the pulse centers [4–7]. Due to certain disadvantages of conventional phase measurement methods, especially in cases of harmonically distorted signals, many authors suggest signal digitalization and further digital signal processing. The author in [6] presents a theory of highprecision phase measurement and a practical algorithm for phase calculation between two signals. The possibility for its implementation is described in [7]. Similarly, improved ways of phase measuring can also be found in [4, 9, 10].
3. Macromodel Description
The proposed phase meter macromodel consists of basic components (resistors, capacitors, etc.) and analog behavioral functions (differentiators, multipliers, etc.). Their real circuit implementation would be difficult, but usage within a simulator is simple. It is intended for the phasemeasurements of cosine signals with known frequencies and excludes measurement of the sinesignal phase. It could however, also be generalized for such cases.
A generalized expression of macromodel’s input cosine signal can be written using the equation: where : amplitude time waveform, : frequency time waveform, : offset voltage time waveform, and φ(t): signal phase time waveform.
Figure 1 shows a block diagram of the phase meter macromodel. The input signals are the voltage signals (t) and f(t). The signal f(t) in volts is proportional to the inputsignal frequency in hertz. The output signal φ(t) in volts quantifies a measured phase in degrees. The phase is measured between the input (t) and internally generated cosine signal (the latter is not shown in this figure).
The model does not measure the frequency f(t). It is presumed that f(t) is known in advance. This simplifies the model’s structure. It also requires a quantification of the following expected maximal input signal parameters: the frequency , the amplitude , and the phase .
The macromodel was developed modularly with appropriate sub models:(i)peakdetector,(ii)offset voltage corrector and amplitude corrector,(iii)signalphase calculator.
The peak detector generates both the maximal and minimal inputsignal time values. They are required for input signal amplitude normalization performed using the offset voltage corrector and the amplitude corrector. The normalized signal drives the phase calculator, where the phase difference is calculated between the input signal and the internally generated ideal cosine signal. The ideal cosine signal has the amplitude 1 V and the same frequency as the inputsignal (t).
The entire signalphase meter macromodel was constructed in the ICAP (ICAP/4 for Windows 8.1.11) SpiceNet as detailed in the following subsections. The usage of the macromodel requires a schematic form (*.dng file) conversion to the text netlist form (*.cir file) and its inclusion to the SPICE library file (*.lib file) as shown in Figure 2.
3.1. Peak Detector
The main task of the peak detector is the detection and holding back of those extreme values from input signal (t). The current peak value is kept constant until the next peakvalue occurs.
The peak detector submodel works as an analog memory of peak values. Therefore, it could also be used in the other cases [11], where the analog peak values should be temporarily stored.
Figure 3 presents the structure of peak detector submodel. It has two inputs (in, ) and two outputs (peak, ). in is driven by cosine signal (t), from which the minimal and maximal time values are generated. Since the concept of the macromodel requires knowledge of the frequency, there is a second voltage driven input . The input voltage in volts must be proportional to the input signal frequency in hertz. Other than , additional parameters (dv, , ) are required for the peak detector to operate. They are further described in this sub section.
Capacitor C1 is needed for temporal storage of the input signals’ local minimum values, whilst capacitor C2 stores its local maximum values. Appropriate control of switches X1 and X2 ensures, that both capacitors are timely charged to the inputvoltages’ local extremes. The charging period is followed by a holdback period, when the capacitors’ voltages are kept constant. Capacitors C1 and C2 are connected to outputs and peak.
The inputvoltages’ extreme values occur during the time instants when the first derivative of the input signal is equal to zero. A second derivative sign is required in order to distinguish between the local minimum and maximum.
Although the input signal parameters could be time variable values, it was assumed that they are constant, in order to simplify the model. In this case, the first derivative of the input signal (t) can be written using the equation: If it is assumed that , , , , (2) can be simplified to
Switches X1 and X2 are controlled by inputsignal derivative values which have been previously normalized. According to (3), the derivative of the input signal can be normalized using constant , where stands for the expected maximal inputsignal amplitude. The first and second derivatives’ normalizations are performed using the voltagecontrolled sources B4 and B5, respectively. Both sources B4 and B5 are linked to sources B1 and B2, which generate the control voltages for switches X1 and X2. X1, X2 are switched “ON” at a positively control voltage, otherwise they are switched “OFF”.
The “ON” resistances of switches X1 and X2 are inverselyproportional to the values of their control voltages. By choosing a value for the control voltage equal to 1/ and equal to 1 Ω, the “ON” time constant , only depends on capacitance C.
X1 is “ON” when the values of the second and first derivatives are correspondingly positive and negative. Simultaneously, the value of the first derivative must be greater than the negative reference voltage . Vice versa applies for switch X2. X2 is “ON” if the second derivative is negative and the first derivative is positive. Simultaneously, the value of the first derivative must be smaller than the reference voltage , generated by source B3.
The reference voltage is the threshold value of the first derivative, when the switches X1 and X2 must be toggled “ON” to ensure adequate charging time Δt for capacitors C1 and C2. In Δt the capacitors are charged to input voltage extreme values with Δv tolerance.
Since the capacitor’s chargingtime is substantially shorter than the input signal period, it can be presume that during the charging periods, the input signal voltage is constant. In the worstcase scenario, at maximal frequency and maximal amplitude , capacitor C is charging from 0 V to . At the permitted tolerance Δv, the capacitor will be charged to the value (Δv) within the time interval:
Since, within the time Δt, the cosine signal will reach the value: the capacitor’s charging time can be written as
Equalizing (4) and (6), the capacitances C1 and C2 are
Choosing equal to 1 Ω, the capacitances depend solely on the permitted tolerance v, and the expected and .
C1 and C2 start to charge at the time instant: For example, for the time t before extremevalue occurrence. According to (6), (8), and the normalized (3), the reference voltage can be calculated using
After the charging period t, the switches X1 and X2 disconnect the connections between the input and capacitors C1 and C2, correspondingly. Whilst the switches are “OFF”, the capacitors remain charged to the input signals’ extreme values. Since switches X1 and X2 are “OFF”, the resistances are very high (e.g., = 1 GΩ), during the “OFF” period the voltages on the capacitors will not change significantly.
Using the submodel requires setting the parameters: v, , and , which are labeled in Figure 3 as dv, , and , consecutively. These parameters enable the calculation of reference voltage according to (9), and the capacitances C1 and C2 according to (7).
3.2. Offset/Amplitude Corrector
The offset voltage corrector removes the offset voltage from the input cosine signal (t). Assuming a constant offset voltage, it can be written as where stands for the minimal and stands for the maximal input cosine signal amplitudes.
The cosinetime waveform without offset voltage can thus be expressed as
The absolute minimal and maximal peak values are equal and are The signal (t) is normalized with , thus providing a constant equalizer output amplitude of value 1 V. The normalized signal waveform describes (13):
The offset voltage corrector and amplitude corrector submodel is shown in Figure 4. The inputs’ peak and are connected to the peak detector. The output signal is the normalized cosine signal with a constant amplitude value of 1 V and a frequency f.
The offset voltage and amplitude correction are performed by controlled voltage source B17, which generates voltage according to (13). B17 additionally limits the output signal’s amplitude within the range . Limitation is performed for the cases of variable inputsignal amplitude, which will be described later. Should the output voltage exceed the permitted range , the signal phase calculation would not be feasible. The process of amplitude limitation does not affect the signal phase measurements, since they are always performed at zero crossing instant.
3.3. SignalPhase Calculator
The signal phase calculator calculates the phase difference φ(t) between the normalized cosine signal and the internally generated cosine signal of amplitude 1 V. The signals frequencies are equal and known in advance. The ideal cosine signal is generated using a voltagecontrolled oscillator modeled with Laplace transformation and a voltagecontrolled source, as proposed in [12]. Figure 5 shows the signal phase calculator submodel.
The calculation of phase difference is periodic within the period T = 1/f, according to the following: It is performed using the controlled voltage source B9. When the switch X3 is switched ON, the capacitor C3 charges to a voltage value proportional to the phase difference measured over the previous period. When the switch X3 is switched OFF, the precharged value is held constant at the Vphase output. X3 is controlled by the source B7. It ensures that the switch is turned on once per period. This event occurs when a second derivative of input signal (t) is negative and the value of the normalized signal is within the range: . Here, stands for the reference voltage, provided by the source B11, according to (9).
When the signal reaches zero value, X3 switches OFF. The voltage of C3 remains constant. Its value in volts is proportional to the phasedifference value φ(t) in degrees. For the same reason described for the peakdetector submodel, the X3 switch “OFF” resistance has to be high enough. This ensures that the voltage on C3 will not change significantly until the next toggle of the switch X3.
The capacitance of C3 is calculated using an equation similar to (7) where is replaced with the predicted maximum value of the inputsignal phase fimax. In the worst case scenario, C3 must be charged to the voltage value proportionally to the fimax value in degrees, with the permissible phasetolerance dφ. The same values of the permissible tolerances were used for the phase calculator and the peak detector in order to reduce the number of model parameters.
3.4. Macromodel Parameters
An entire macromodel of the precise signalphase meter, suitable for SPICE simulators consists of all previous described submodels. It performs precise cosine signal phase measurements at any frequency, amplitude, and offset voltage on condition that the input signalfrequency time waveform is known in advance. The model’s operation is based on the assumption that the inputsignal parameters are constant. However, its usefulness in situations of time variable parameters was also examined. The signal phase measurement, whilst the signal’s frequency, amplitude, and offset voltage are timevariable, is a very problematic task [13] and in general it has not been fully explored as yet.
Using the macromodel of the signalphase meter requires knowledge of the four model parameters that have an influence on its accuracy and precision. These parameters are(i): expected maximum inputsignal amplitude;(ii): expected maximum inputsignal frequency;(iii): expected maximum inputsignal phase;(iv)dv: admissible tolerance of input signal peak values (equal to the admissible tolerance for the measured signal phase).
The accuracies of signalphase measurements depend on the macromodel’s accuracy, and on the simulation parameter values. They are transient (.TRAN) analysis parameters:(i): time to start recording data;(ii): total analysis time;(iii): datastep time;(iv): maximum timestep.
The most influential simulation parameter is . It must be lower or equal to the capacitors’ C1–C3 chargingtimes. Using (6), the recommended maximum time step can be calculated as where stands for the smallest values:
Any increase of model’s accuracy can be achieved with a decrease in the parameters dv or dφ. According to (15), this decreases which reduces the integration step time. Therefore, more accurate results prolong the simulation duration.
4. PhaseMeasurement Simulation Results
The following section presents the phase measurement simulation results using the proposed macromodel for the SPICE simulator. The focus was on the phase measurements of the inputcosine signal with constant and timevariable parameters. In the further described simulation examples the transient analysis (.TRAN) with the GEAR integration method was used, with a relative error tolerance .
4.1. Phase Measurements of the Input Signal with Constant Parameters
This example focused on the phase measurement of a cosine signal with constant parameters, which were (t) = 1 V, , φ(t) = 0.1°, and . Using these parameters the following macromodel parameters were defined: , = 1 kHz, and = 0.1°. Simulations were performed at three different tolerances dv of values: V, V, and V. Using (15) and (7), the simulation parameters : 2.25 μs, 712 ns and 225 ns, were calculated. The rest of simulation parameters were = 5 ms, = 100 ms, and = 10 μs.
Figure 6 shows the simulation results at three different permissible tolerances dv. Δφ(t) is the phase error between measured and the actual phase of the input signal.
Despite the constant input signal parameters, Δφ is time variable. Its time variability is the consequence of the simulator’s internal integration time step, which varies whilst the transient analysis is running. The integration timestep automatically adapts to the time changes of circuit signals. Therefore, Δφ(t) are not constanttime waveforms.
Table 1 shows the worst case results for all three examples.

From Figure 6 and Table 1, it can be concluded that reducing the parameter dv decreases the phase error Δφ.
4.2. Phase Measurements of Input Signal with TimeVariable Parameters
At the macromodel development stage, it was assumed that the input signal parameters are not time variable. Nevertheless, its suitability for cases of time variable parameters has also been verified.
The chosen intervals of the timevariable input signal parameters were as follows: 0.8 V ≤ (t) ≤ 1.1 V; 0 ≤ mV; 100 Hz ≤ f(t) ≤ 10 kHz; −2° ≤ φ(t) ≤ +5°. Figure 7 shows the parameters’ time waveforms used in our example.
(a)
(b)
(c)
(d)
(e)
According to the intervals of the inputsignal parameters, the following macromodel parameters were set: = 10 kHz, Amax = 1.1 V, and = 5°. The value of the permissible tolerance dv was chosen as dv = V. According to (15) the maximum time step was calculated as = 101 ns. The transient analysis results are shown in Figure 8.
(a)
(b)
Δφ(t) waveform shows the phase error between the actual φ(t) and measured (t) signal phases. Noticeable, the phase error is a result of an assumption that the inputsignal parameters are constant. Consequently, the derivatives in (2) are equal to zero which, in our case, is untrue. Rapid changes in inputsignal parameters increase the phaseerror. This clearly illustrates the time waveform 4 within the timerange from 20 ms to 80 ms, when simultaneous changes in amplitude, frequency, phase, and offset voltage occur. In continuation, the time variations of the parameters were moderate and the phase error Δφ was significantly reduced. Therefore, the highest precision signal phase can only be measured under stationary conditions, when the smallest possible phase error can be expected.
5. Conclusion
The objective of this paper was to present a new macromodel of a precise signalphase meter, as required for the development of integrated circuits, implemented within linear position encoders. The macromodel is suitable for precise phase measurements whilst transient analysis is running. Namely the reduction of various unwanted effects of an encoder’s head requires very accurate signalphase measurements during the encoder development stage.
In order to simplify the concept of a macromodel, it was assumed that signal parameters are not timevariable values. As shown in Section 4, as such a signal phase was measured within an accuracy of more than ±0.02%. Simulations also showed a substantial decrease in the results’ accuracies due to the rapidtime change of the signal parameters. The signalphase measurement under such conditions is a very problematic task [13] and has not yet been fully explored. Therefore any improvement in the presented macromodel in this direction would not be meaningful. However, this macromodel concept represents an idea for realworld implementation of a precise signalphase meter.
References
 J. Rozman and A. Pleteršek, “Optical encoder head with improved linearity,” Informacije MIDEM, vol. 39, no. 1, pp. 22–27, 2009. View at: Google Scholar
 J. Rozman and A. Pleteršek, “Linear optical encoder system with sinusoidal signal distortion below 60 dB,” IEEE Transactions on Instrumentation and Measurement, vol. 59, no. 6, pp. 1544–1549, 2010. View at: Publisher Site  Google Scholar
 K. M. Ibrahim and M. A. H. AbdulKarim, “A novel digital phase meter,” IEEE Transactions on Instrumentation and Measurement, vol. 36, no. 3, pp. 711–716, 1987. View at: Google Scholar
 R. Micheletti, “Phase angle measurement between two sinusoidal signals,” IEEE Transactions on Instrumentation and Measurement, vol. 40, no. 1, pp. 40–42, 1991. View at: Publisher Site  Google Scholar
 S. M. Mahmud, “Error analysis of digital phase measurement of distorted waves,” IEEE Transactions on Instrumentation and Measurement, vol. 38, no. 1, pp. 6–9, 1989. View at: Publisher Site  Google Scholar
 S. M. Mahmud, “High precision phase measurement using adaptive sampling,” IEEE Transactions on Instrumentation and Measurement, vol. 38, no. 5, pp. 954–960, 1989. View at: Publisher Site  Google Scholar
 S. M. Mahmud, “High precision phase measurement using reduced sine and cosine tables,” IEEE Transactions on Instrumentation and Measurement, vol. 39, no. 1, pp. 56–60, 1990. View at: Publisher Site  Google Scholar
 T. Kishibe, T. Masuda, and H. Tanoue, “Phase correction circuit of encoder signal,” U.S. Patent 7 933 373, April 26, 2011. View at: Google Scholar
 B. Djokić and E. So, “Phase measurement of distorted periodic signals based on nonsynchronous digital filtering,” IEEE Transactions on Instrumentation and Measurement, vol. 50, no. 4, pp. 864–867, 2001. View at: Publisher Site  Google Scholar
 M. Sedlacek and M. Krumpholc, “Digital measurement of phase difference—a comparative study of DSP algorithms,” Metrology and Measurement Systems, vol. 12, no. 4, pp. 427–448, 2005. View at: Google Scholar
 M. Šalamon and T. Dogša, “The model of chaoticness detector,” Informacije MIDEM, vol. 39, no. 2, pp. 93–99, 2009. View at: Google Scholar
 “Intusoft Newsletter Issue 73,” February 2004, http://www.intusoft.com/nlpdf/nl73.pdf. View at: Google Scholar
 T. Dogša, M. Šalamon, B. Jarc, and M. Solar, “Optical encoder scanning head modelling and ultra small phase shift measurement,” Informacije MIDEM, vol. 39, no. 3, pp. 162–167, 2009. View at: Google Scholar
Copyright
Copyright © 2012 Matej Šalamon and Bojan Jarc. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.