This invention relates to the field of precision timing, and in particular to a method of compensating for non-linear effects during phase interpolation in digital-to-time conversion, and a clock synthesizer for implementing the method.
A variety of electronic circuits require the synthesis of extremely accurate clock signals at a precise frequency. One way of achieving this result is to lock a digital phase-locked loop employing a software digital controlled oscillator (SDCO) to generate phase and frequency values from a reference signal, which after suitable conversion, drive a hardware digital or numerical controlled oscillator (DCO/NCO), referred to hereafter as a DCO. This comprises an accumulator that increments by a certain amount on each system pulse. When the accumulator overflows, it generates a carry output, which represents the output pulse. The remainder left in the accumulator when a carry is output represents the phase. The DCO produces clock pulses at the desired overall frequency, but since these pulses are coincident with the system pulses when the carry occurs, they are not generally correctly phased and form a so-called gapped clock signal because when there is no carry no output occurs on some of the system clock pulses.
A digital-to-time converter (DTC) makes use of the phase information provided by the remainder to interpolate the output of the hardware DCO and produce evenly spaced clock pulses in their correct phase locations. Such an arrangement is described in U.S. Pat. No. 8,692,599, the contents of which are herein incorporated by reference.
A problem arises in the prior art due to the fact that the DTC, which is implemented in analog form, is not perfectly linear. The non-linear space between the interpolation points in each system clock pulse period creates output clock jitter that is significantly larger than the resolution of the DTC interpolation as determined by the number M of interpolation points k. Pre-calibration is generally used to eliminate this non-linear distortion, but this has to be constantly adjusted to compensate for temperature variations.
Embodiments of this invention employ an accurate compensation method that can quickly and adaptively adjust the non-linear distortion in the DTC interpolation and thus create a low jitter output clock. In accordance with embodiments of the invention, two paths are employed, each comprising a DCO and a DTC. A deliberate initial delay is introduced at the input of the second path, which is compensated for by an analog delay introduced into the output of the first path. Varying the initial delay moves the interpolation point for the second path but has no effect on the first path. By averaging measurements for multiple matching pairs of initial and compensating delays, the invention produces discrete average integral non-linear distortion values for each interpolation point k in the first path, and an average integral non-linear distortion value for all M interpolation points over a system clock period in the second path. This average does not vary much with the initial phase value input to the DCO, and therefore can be treated as a constant, which can be allowed for in the final circuitry.
According to the present invention there is provided a method of compensating for integral nonlinear interpolation (INL) distortion in a clock synthesizer driven by a system clock running at a frequency fsys, comprising:
Nominally the analog delay I*dt and nominal digital delay I*D are set to the same values to that they offset each other, but in practice the analog delay cannot be precisely set and will have a small offset δ from the nominal delay giving an actual delay of I*dt+δ. The actual digital delay I*D+Δ matching this actual analog delay has to be found. The matching operation can be done in advance for all pairs of analog/digital delays I*dt/I*D, or in one embodiment it can be done concurrently with the determination of the INL compensation values.
It will be appreciated that, as generally the case in an adaptive algorithm, both I and k can be changed randomly so that statistically all pairs of k and I are covered. The sequencing of steps (iii) and (iv) is not important. For example, the system could maintain a fixed value of I and until all values of k are covered, or continually change both I and k so as to eventually cover all combinations of I and k.
In one embodiment, at each iteration, the contents x(k) of the second compensation module are added to a value y(k) stored in a memory according to the equation:
y(k)=y(k)+x(k−Ib)
where Ib is an index addition bias, and the final summation is divided by N to derive the compensation values for at least one of said first and second paths
The index addition bias Ib may be set to zero to obtain the compensation values for the second DTC and to the same value as the initial phase of the second DCO to obtain the compensation values for the first DTC.
According to a second aspect of the invention there is provided a clock synthesizer with nonlinear digital-to-time conversion compensation driven by system clock pulses, comprising:
According to yet another aspect of the invention there is provided a clock synthesizer with nonlinear digital-to-time conversion compensation driven by system clock pulses, comprising:
This invention will now be described in more detail, by way of example only, with reference to the accompanying drawings, in which:
The phase locked loop (PLL) shown in
The DCO 20 is shown in more detail in
The overall output frequency matches the desired frequency, but the phase is inaccurate if when overflow occurs there is a remainder value in the accumulator register. This represents the phase (also essentially just a number) information that is used by digital-to-time converter (DTC) 28 (
Since it is critical to understanding the invention, a simplified example of the operation of the accumulator may be appropriate. If we assume for the sake of argument that the accumulator is a 4-bit counter (in practice it is larger than this) and contains a number, say 7, this represents the phase value because it identifies the point in the accumulator cycle where the carry would occur if the system did not have to wait for the next system clock pulse to generate an output. If a frequency value of say 9 is added to the accumulator on a system clock pulse, the accumulator will reset to zero and output a carry pulse representing the output of the DTC 28. The output pulse will be correctly positioned coincident with the system clock pulse. If the frequency value 9 is added on the next system clock pulse, the phase value in the accumulator will be 9, and no output will be generated. If the frequency value 9 is added again on the next system clock pulse, an output will be generated but the value in the accumulator will roll over to 2 since the accumulator is counting modulo 16. It can thus be seen that the output pulses are uneven (there are gaps between them as shown in
However, as noted above, due to the nonlinearity of the analog circuitry in the DTC, the DTC circuit 28 may not produce a linearly spaced clock output, and the actual phase of an output pulse at an interpolated location k will be
Phase/Freq*1/fsys+INL(k)
where INL(k) is the integral Non-Linear (INL) distortion at interpolated location k. This is a function of phase/Freq, namely the correct position of the output pulse between adjacent clock pulses. INL(k) adds jitter noise to the noise resulting from the limitations imposed by the resolution of DTC interpolation.
However, if INL(k) can be estimated, it can be compensated for as shown in
In a typical PLL application, INL(k) may be constantly changing with temperature, voltage and the DCO frequency. In accordance with embodiments of the invention the values of INL(k) can be adaptively estimated on the basis of principles explained with reference to
In the block diagram shown in
Since both DCO0200 and DCO1201 are driven by the same DCO frequency Freq and the two DCOs 200 201 start with different initial phases with phase difference of I*D, at any given time, the phase of DCO0200 is phasek and the phase of DCO1201 is phasek+I*D.
If M is the number of steps of DTC interpolation in a single system clock period 1/fsys, the DTC resolution is 1/(M*fsys). Therefore, INL has M values (INL(k), k=0, 1, . . . , M−1), one value for each interpolation step, k.
It should be noted that INL(k+M)=INL(k−M)=INL(k) because INL only represents the non-linear distortion for the fractional portion of the phase (i.e., remainder of the DCO phase), so the value of INL at the same interpolation point k between successive system pulses will always be the same.
Now, if we let k=└phasek*M/Freq┘ be defined as the integer of (phasek*M/Freq) (this is true because the DTC only interpolates to discrete points), where phasek<Freq and k<M, and d=D*M/Freq, the phase difference between the output pulses of the two paths, namely the output of the delay estimation circuit will be:
eI(k)=[(phasek+I*D)/Freq*1/fsys+INL1(k+└I*d┘)]−[phasek/Freq*1/fsys+INL0(k)]
where INL1(k+└I*d┘) and INL0(k) represent INL values for DTC1 and DTC0 respectively.
Since the initial phase I*D/Freq*1/fsys is known, the above equation be simplified by offsetting the initial phase with a delay circuit in the output of the DTC in that path without the initial delay to give:
eI(k)=INL1(k+└I*d┘)−INL0(k)
Over a series of system clock pulses, as the value Freq is successively added to the accumulator in the DCOs 200, 201, the DCO will pass through all phases k, giving a set of averaged values eI(k) where k=0, 1, . . . , M, with a given initial phase difference I*D.
If the same procedure is repeated for different initial phase settings I, the average values with respect to I are
With a large value N,
will be close to a constant and will vary very little with respect to k. If N=M, then because the condition INL(k+M)=INL(k) holds, the summation INL(k)+INL(k+1)+ . . . +INL(k+M) is independent of value k and is the same as INL(0)+INL(1)+ . . . +INL(M−1). If N<M, the summation will vary slightly but will be close to a constant provided N is large enough.
The average of eI(k) will give an INL value INL0(k) with a reverse sign and a small constant bias. The constant bias does not affect the jitter performance. This average value can be applied to the compensation module 30 to cancel the INL for DTC0.
In a similar manner, if instead of summing with respect to eI(k), we sum with respect to eI(k−└I*d┘), we have
By similar reasoning the average value of
remains substantially constant. As a result the INL values for DTC1 can be obtained to permit DTC1 to be compensated.
It is noted that in the above reasoning the initial delay is subtracted from the equation. This is achieved in theory, as shown in
When the two clocks are aligned in
eI(k)=[phasek/Freq*1/fsys+INL1(k+└I*d┘)−I*dt]−[(phasek−I*D)/Freq*1/fsys+INL0(k)]
If the initial phase I*D/Freq*1/fsys in DCO0 cancelled the analog delay I*dt, we have the same result as obtained in connection with
eI(k)=INL1(k+└I*d┘)−INL0(k)
The problem with the arrangement shown in
A suitable circuit for matching the digital delay introduced into one path with the analog delay introduced into the other path is shown in
The Q, Q−1 outputs of DFF 40 are input to counter 44, whose output is applied to two threshold comparators 46, 48 coupled respectively to up/down inputs of up/down module 31. This adjusts the delay I*D in the delay memory 42 by small increments of Δ.
With a given actual analog delay I*dt+δ, both DCO0 and DCO1 are driven by the same DCO frequency with the same initial phase values. DCO0 will move its phase backwards or forwards based on the output of DFF 40 so as to align the two output clocks.
The output of DFF 40 will be either positive or negative depending on whether the D input is ahead or behind the driven clock. The output of the DFF 40 will be sent to the counter 44, which will add or subtract 1 (+/−1) based on DFF output. The counter 44 will be read on a controlled edge of the system clock. If the value in the counter is larger than a given threshold, it means that DCO0 is behind DCO1 and DCO0 will move its phase forward. If the counter is smaller than threshold, it means that DCO0 is ahead of DCO1 and DCO0 will move its phase backward.
If the counter value is between two thresholds, no change is made to the phase of DCO0. Gradually, the two clock phases will be aligned at the inputs of DFF 40. After alignment is completed, the adjusted DCO0 phase I*D+Δ will reflect the analog delay I*dt+δ and be stored in delay memory 42. The same procedure is be repeated for another setting of nominal analog delay I*dt to build a library of actual digital delays I*D+Δ corresponding to nominal analog delays I*dt (actual analog delays I*dt+δ).
When all the nominal analog delays I*dt have been premeasured, the compensation procedure can begin as shown in
The counter output average module 44 receives the Q and Q−1 outputs of the D-type flop 40 and increments up or down depending on whether the first path leads or lags the second path. When the counter 44 exceeds a preset threshold, determined by comparator 46, it adjusts the compensation memory of DTC0 down, and when it falls below a preset threshold determined by comparator 48, it adjusts the compensation memory of DTC0 up so that clock outputs from both paths are aligned to each other.
After convergence, i.e. when the outputs of the paths 290 and 291 are aligned, the contents of the compensation memory x(k), which are in the form of a vector representing the INL values for each k, are dumped into a circular accumulator memory provided by the average block 50 with the index addition bias Ib in accordance with the equation:
y(k)=y(k)+x(k−Ib),
where y(k) is the accumulated contents of the memory 50 and Ib. If DTC1 is to be compensated, the index addition bias is I*d, so the equation becomes:
y(k)=y(k)+x(k−└I*d┘)
If DTC0 is to be compensated, the index addition bias is zero, and the equation becomes:
y(k)=y(k)+x(k)
After finishing the alignment for the first delay setting, the second analog delay is selected and the previous procedure is repeated, and so on for all N values of I.
The circular memory provided by average block 50 contains the sum of all the compensation memory values for different delays. By dividing by N into the final summation and setting the first element to zero and subtracting all elements by the first element (assuming INL(0)=0), we have the average compensation value, which represents INL for DTC0 if the addition index bias is set to zero and INL for DTC1 if the addition index bias is set to the same as the initial phase of DCO0.
In the embodiment shown in
During the compensating process, the compensated synthesizer (the corresponding DCO and its DTC) runs normally without interruption. When the compensation process is completed, its DTC compensation memory is updated and next synthesizer is selected with the multiplexer 38.
The embodiment shown in
In the previous description, the dual DTC nonlinear compensation was performed in two separate steps: delay estimation and DTC compensation with variable delay setting. These two steps can be preformed concurrently to provide a combined estimate of analog delay and DTC nonlinear compensation as shown in
As in the case of the embodiment of
In
The process will continue until convergence when the output pulses are aligned. During the convergence procedure, the controller 52 rotates the analog delay 0 to (N−1)*dt regularly after the INL values for each value of k have been obtained. The controller 52 also selects the corresponding value I*D in the delay memory 42, which is adaptively updated by the scaling factor Scale2 output from multiplier 58. Likewise, the compensation memory 300 will be adaptively updated by the scaling factor Scale1 output by the multiplier 56. Both I and k can be changed randomly as is common in an adaptive algorithm. I can be held fixed until all values of k have been covered, or the system can continually change I and k to cover all combinations of I and k. If I is changed at random, statistically all combinations of I and k will eventually be covered.
As with any adaptive algorithm, the final error will approach a stable level. In this case, the DFF 40 output will toggle between +/−1. The final value in the compensation memory 300 will be the average INL for the corresponding DTC, and the delay memory 42 will contain corresponding actual digital I delay value D+Δ for the selected nominal analog delay I*dt.
Non-limiting aspects of the invention include a novel architecture for DTC nonlinear compensation, DTC nonlinear compensation using dual PLL architecture, real time DTC nonlinear compensation to compensate for temperature variations, dual PLL compensation with added variable analog delay, multi-DPLL DTC real time compensation, joint/adaptive DTC compensation and analog delay measurement, DTC compensation for INL cancellation, and PI calibration for INL cancellation.
It should be appreciated by those skilled in the art that any block diagrams herein represent conceptual views of illustrative circuitry embodying the principles of the invention. For example, a processor may be provided through the use of dedicated hardware as well as hardware capable of executing software in association with appropriate software. When provided by a processor, the functions may be provided by a single dedicated processor, by a single shared processor, or by a plurality of individual processors, some of which may be shared. Moreover, explicit use of the term “processor” should not be construed to refer exclusively to hardware capable of executing software, and may implicitly include, without limitation, digital signal processor (DSP) hardware, network processor, application specific integrated circuit (ASIC), field programmable gate array (FPGA), read only memory (ROM) for storing software, random access memory (RAM), and non volatile storage. Other hardware, conventional and/or custom, may also be included. The functional blocks or modules illustrated herein may in practice be implemented in hardware or software running on a suitable processor, and the terms circuit or circuitry include collections of functional blocks implemented in software.
This application claims the benefit under 35 USC 119(e) of U.S. provisional application No. 62/447,583, filed, Jan. 18, 2017, the contents of which are herein incorporated by reference.
Number | Name | Date | Kind |
---|---|---|---|
8692599 | Gong et al. | Apr 2014 | B2 |
Entry |
---|
CircuitVision Analysis of the Crystal Oscillator, DCO Circuitry and Output Buffers on the Silicon Labs Si5345 Clock Generator (TechInsights 2014). |
Number | Date | Country | |
---|---|---|---|
20180205370 A1 | Jul 2018 | US |
Number | Date | Country | |
---|---|---|---|
62447583 | Jan 2017 | US |