Phase-locked loop with multiple degrees of freedom and its design and fabrication method
09602114 · 2017-03-21
Assignee
Inventors
- Michael Pelissier (Grenoble, FR)
- Anton Korniienko (Lyons, FR)
- Mykhailo Zarudniev (Kiev, UA)
- Gèrard Scorletti (Ecully, FR)
- Olesia Mokrenko (Grenoble, FR)
- Eric Blanco (Decines-Charpieu, FR)
- Patrick Villard (Grenoble, FR)
- Gèrard Billiot (Saint-Nazaire les Eymes, FR)
Cpc classification
H03L7/099
ELECTRICITY
International classification
H03L7/06
ELECTRICITY
H04L7/00
ELECTRICITY
Abstract
A design method for a phase-locked loop comprises: a controlled-frequency oscillator; a phase comparator, to determine a phase difference between an output signal of the controlled-frequency oscillator and a reference signal; a corrector to receive as input a signal representative of the phase difference and to generate at its output a first correction signal; at least one second corrector, to receive as input a signal representative of or affected by a phase noise of the reference signal or of the output signal of the controlled-frequency oscillator and to generate at its output a second correction signal; and a circuit for generating a slaving signal for the controlled-frequency oscillator on the basis of the first and second correction signals; the method using the H-infinity method. Method for fabricating such a loop comprising a design step implementing this method. Phase-locked loop thus obtained.
Claims
1. A method of designing a phase-locked loop, the phase-looked loop being of the type comprising: a controlled-frequency oscillator; a phase comparator, configured to determine a phase difference between an output signal of the said controlled-frequency oscillator and a reference signal; a first corrector, configured to receive as input a signal representative of the said phase difference and to generate at its output a first correction signal; at least one second corrector, configured to receive as input a signal representative of or affected by a phase noise of the said reference signal or of the said output signal of the said controlled-frequency oscillator and to generate at its output a second correction signal; and a circuit for generating a slaving signal for the said controlled-frequency oscillator from the said first and second correction signals; the said method comprising a step of determining transfer functions of the said correctors allowing rejection, in one and the same frequency band, the phase noise of the said reference signal and the phase noise of the said output signal of the said controlled-frequency oscillator, and wherein said step is implemented, by means of a computer, by applying the H-infinity (H) method utilizing: at input, a first weighting function for a phase noise of the said controlled-frequency oscillator and a second weighting function for a phase noise of a reference signal, and which are determined on the basis of nominal power spectral densities of the said noise; and at output, at least one third weighting function for a phase noise of an output signal of the phase-lock loop or of an error in tracking the said reference signal, determined on the basis of a phase noise power spectral density template of the said output signal having to be complied with.
2. The method of claim 1, wherein the said step of determining transfer functions of the said correctors comprises the following sub-steps: a) determining a nominal power spectral density of the said reference signal; b) determining a nominal power spectral density of the phase noise of the said controlled-frequency oscillator; c) determining a phase noise power spectral density template of the said output signal having to be complied with; d) determining, on the basis of the said nominal power spectral densities and of the said power spectral density template, at least the said first, second and third weighting functions; e) constructing an augmented system through the said weightings; and f) applying the H-infinity (H) method to the said augmented system so as to synthesize the transfer functions of the said correctors.
3. The method of claim 2, wherein the said step d) also comprises the determination of a fourth weighting function at output for the power spectral density of the said slaving signal on the basis of a power spectral density template of the said slaving signal to be complied with.
4. The method of claim 2, also comprising a sub-step d) of simplifying the weighting functions determined during sub-step d), the weighting functions thus simplified being used during the said sub-step e), the said simplification being implemented by approximating the said weighting functions by transfer functions of lower order and of smaller modulus at least over a spectral operating span of the phase-locked loop.
5. The method of claim 2, also comprising a sub-step g) of simplifying the transfer functions synthesized during the said sub-step f), the said simplification being implemented by approximating the said transfer functions by transfer functions of lower order.
6. The method of claim 2, wherein the said step f) is implemented with an extra constraint, according to which the modulus of the transfer function going from the said output signal of the said controlled-frequency oscillator to the said reference signal exhibits, in at least one spectral span, a slope of greater than or equal to +20 db/decade and preferably of greater than or equal to 40 dB/decade.
7. The method of claim 1, wherein the said phase-locked loop comprises at least one third corrector configured to receive as input a signal representative of or affected by a phase noise generated inside the said loop, other than the said phase noise of the said output signal of the said controlled-frequency oscillator, and to generate at its output a third correction signal, the said circuit for generating a slaving signal for the said controlled-frequency oscillator being configured to generate the said slaving signal also on the basis of the said third correction signal, the said step, implemented by the said H-infinity (H) method, of determining transfer functions of the said correctors using as input also a weighting function for the said phase noise other than the said phase noise of the said output signal of the said controlled-frequency oscillator, and also determined on the basis of the power spectral density of the said noise.
8. A method for fabricating a phase-locked loop, of the type comprising: a controlled-frequency oscillator; a phase comparator for determining a phase difference between an output signal of the said controlled-frequency oscillator and a reference signal; a first corrector, configured to receive as input a signal representative of the said phase difference and to generate at its output a signal termed the first correction signal; at least one second corrector, configured to receive as input a signal representative of or affected by a phase noise of the said reference signal or of the said output signal of the said controlled-frequency oscillator and to generate at its output a signal termed the second correction signal; and a circuit for generating a signal for slaving the said controlled-frequency oscillator on the basis of the said first and second correction signals; the said method comprising: a step of design of the said phase-locked loop; and a step of physical production of the phase-locked loop thus designed; wherein the said design step is implemented by a method according to claim 1.
9. A method for fabricating a phase-locked loop, of the type comprising: a controlled-frequency oscillator; a phase comparator for determining a phase difference between an output signal of the said controlled-frequency oscillator and a reference signal; a first corrector, configured to receive as input a signal representative of the said phase difference and to generate at its output a signal termed the first correction signal; at least one second corrector, configured to receive as input a signal representative of or affected by a phase noise of the said reference signal or of the said output signal of the said controlled-frequency oscillator and to generate at its output a signal termed the second correction signal; and a circuit for generating a signal for slaving the said controlled-frequency oscillator on the basis of the said first and second correction signals; the said method comprising: a step of design of the said phase-locked loop; and a step of physical production of the phase-locked loop thus designed; wherein the said phase-locked loop comprises at least one third corrector, configured to receive as input a signal representative of or affected by a phase noise generated inside the said loop, other than the said phase noise of the said output signal of the said controlled-frequency oscillator, and to generate at its output a third correction signal, the said circuit for generating a slaving signal for the said controlled-frequency oscillator being configured to generate the said slaving signal also on the basis of the said third correction signal, the said design step being implemented by the method of claim 7.
10. A phase-locked loop comprising: a controlled-frequency oscillator; a phase comparator for determining a phase difference between an output signal of the said controlled-frequency oscillator and a reference signal; a first corrector, configured to receive as input a signal representative of the said phase difference and to generate at its output a signal termed the first correction signal; at least one second corrector, configured to receive as input a signal representative of or affected by a phase noise of the said reference signal or of the said output signal of the said controlled-frequency oscillator and to generate at its output a signal termed the second correction signal; and a circuit for generating a signal for slaving the said controlled-frequency oscillator from the said first and second correction signals; wherein the said first and second correctors exhibit non-constant transfer functions, chosen so as to allow the rejection, in one and the same frequency band, of the phase noise of the said reference signal and of the phase noise of the said output signal of the said controlled-frequency oscillator.
11. The phase-locked loop of claim 10, wherein the said second corrector is configured to receive as input: either a signal representative of the phase of the said reference signal; or an estimation of the phase noise of the said reference signal.
12. The phase-locked loop of claim 11, also comprising a circuit for determining an estimation of the phase noise of the said reference signal comprising: a delayer module, configured to generate a version of said reference signal delayed by a time (Td) relative to said reference signal; a phase-shifter module, configured to generate a version of the said reference signal exhibiting a phase shift , with
13. The phase-locked loop of claim 10, wherein the said second corrector is configured to receive as input: either a signal representative of the phase of the said output signal of the said controlled-frequency oscillator; or an estimation of the phase noise of the said output signal of the said controlled-frequency oscillator.
14. The phase-locked loop of claim 13, also comprising a circuit for determining an estimation of the phase noise of the said output signal comprising: a delayer module, configured to generate a version of said output signal delayed by a time (Td) relative to said output signal; a phase-shifter module, configured to generate a version of the said output signal exhibiting a phase shift , with
15. The phase-locked loop of claim 10, also comprising at least one third corrector, configured to receive as input: a signal representative of or affected by a phase noise generated inside the said loop, other than the said phase noise of the said output signal of the said controlled-frequency oscillator, and to generate at its output a third correction signal, the said circuit for generating a slaving signal for the said controlled-frequency oscillator being configured to generate the said slaving signal also on the basis of the said third correction signal.
16. A radiofrequency reception chain comprising: a radiofrequency preamplifier; a mixer configured to receive as input an output signal of the said radiofrequency preamplifier and a frequency conversion sinusoidal signal, and to provide as output a signal at intermediate frequency; and a circuit for generating the said frequency conversion sinusoidal signal; wherein the said circuit for generating the said frequency conversion sinusoidal signal comprises the phase-locked loop of claim 14, using the said mixer in the guise of mixer of the said circuit for determining an estimation of the phase noise of the said output signal, and an oscillator configured to generate the said reference signal.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
(1) Other characteristics, details and advantages of the invention will emerge on reading the description provided with reference to the appended drawings given by way of example and which represent, respectively:
(2)
(3)
(4)
(5)
(6)
(7)
(8)
(9)
(10)
(11)
(12)
(13)
(14)
(15)
(16)
(17)
(18)
(19)
(20)
(21)
(22)
(23)
DETAILED DESCRIPTION
(24)
(25) Under conditions where the signals are small (that is to say for .sub.ref(t)(t)), the manner of operation of the PLL can be illustrated by the diagram of
(26) The signal .sub.ref can be considered to be the useful signal to which the output of the PLL must be slaved. The corrector must therefore ensure this slaving. On the other hand, the noise signals b and b.sub.ref are parasitic signalsdue respectively to the VCO and to the oscillator producing the reference signalwhich must be filtered in such a way that their contributions to the output of the PLL is as low as possible while guaranteeing slaving. The choice of the corrector F(s) which guarantees these performance properties of the PLL will subsequently be called the optimization of the phase noise of the PLL. Other sources of phase noise exist (phase comparator noise, frequency divider noise, etc.) but their contribution is generally less significant.
(27) The slaving of the PLL can be carried out via an analogue implementation or a digital one.
(28) Given that the noise filtering requirement is naturally expressed in the form of a template in the Power Spectral Density (PSD) of the signals, it is beneficial to express the PSD of the output of the PLL with respect to the PSD of input signals considered. We therefore obtain, by considering the VCO noise, the noise of the reference and the reference signal to be uncorrelated:
(29)
where S.sub.x is the PSD of the signal x, with x being able to take the values: , .sub.ref, b, b.sub.ref, f the frequency expressed in Hz.
(30) Equation (1) features the moduli squared of the frequency responses of the transfer functions which link the reference phase signal .sub.ref and the phase noise of the reference, b.sub.ref, and of the VCO, b, with the output of the PLL. The transfer function for going from the noise of the VCO, b, to the output of the PLL , will be denoted
(31)
and will be called the sensitivity function. The transfer function for going from the phase noise of the reference b.sub.ref (or equivalently the phase signal of the reference .sub.ref) to the output of the PLL will be denoted
(32)
and will be called the complementary sensitivity function.
(33) It is easy to appreciate that:
(34)
(35) Generally the transfer function S is that of a high-pass filter so as to be able to filter the contribution of the noise b of the VCO around the carrier frequency of the PLL. Consequently, in view of relation (2), the transfer function T represents a low-pass filter, which is responsible for shaping the reference noise b.sub.ref. Moreover, it may be demonstrated that because of relation (2), and whatever corrector F(s) is chosen to ensure the stability of the PLL, it is not possible, for a given frequency, for the moduli of the frequency responses of two transfer functions S(s) and T(s) to both be low i.e. |S(j2f)|<<1 and |T(j2f)|<<1 for f[f.sub.min, f.sub.max]. If the modulus of the frequency response of the sensitivity function S(s) is low for a given frequency (|S(j2f)|<<1), it is inevitable that the modulus of the frequency response of the complementary sensitivity function T(s) is close to 1(|T(j2f)|1) so that relation (2) is complied with (see
(36)
(37) In this embodiment, the first corrector F.sub.1 acts both on the noise of the VCO and on the noise of the reference, while the second corrector F.sub.2 acts specifically on the noise of the reference.
(38)
(39) It is however difficult to produce PLLs based on these two schematic diagrams since in practice it turns out to be tricky to measure directly the phase noise signal at the input of the corrector F.sub.2(s); indeed, the signal r which appears in diagrams 4a and 5a is not, in reality, measurable in an electronic circuit. For this reason, it is generally preferable to use an extra circuit which receives as input a noisy periodic voltage signal and provides at its output an estimation of the phase noise affecting this signal (for an embodiment of such a circuit, see
(40) It is also possible to combine feed-forward control and feedback to produce a PLL with three degrees of freedom; such a PLL will be described further on with reference to
(41) It is also possible to produce phase-lock loops with more than three degrees of freedom to correct sources of noise other than the reference signal and the VCO, for example the noise introduced by a frequency divider. By way of example,
(42) The invention also proposes a systematic method for the design (synthesis of the transfer functions) of the correctors which will be described in detail hereinafter, with the aid of
(43) A general introduction to the H.sub. method and to the use of weighting functions within the framework of this methodthis constituting a significant aspect of the inventioncan be found in the following document: [Bib-1992], [ScF-2009] and [SkP-2005].
(44) In addition to guaranteeing predefined performance in terms of rejection of the phase noise, the correctors must advantageously ensure that the output signal follows the signal .sub.ref at input, and that therefore the error signal tends to zero as the time tends to infinity. In accordance with the final value theorem, assuming that the PLL is stable, the correctors must ensure:
lim.sub.t.fwdarw.+(t)=lim.sub.s.fwdarw.0s(s)=lim.sub.s.fwdarw.0sS(s).sub.ref(s)=0(3)
(45) We consider the case of a reference signal in the form of a phase ramp (the voltage signal s.sub.ref is therefore a sinusoid); its Laplace transform is therefore .sub.ref(s)=As.sup.2, where A is the slope of the ramp. We therefore obtain:
(46)
from which it is deduced that the sensitivity function must exhibit a slope of at least +40 dB/dec for low frequencies.
(47) As a variant, if only synchronization in terms of frequency is required and not in terms of frequency and phase, it is only necessary to ensure that tends to a constant value, possibly different from zero. In this case, it is possible to relax the condition on the sensitivity function slope at low frequency, which need only be greater than or equal to +20 dB/dec.
(48) Furthermore, it is desired to require that the PSD of the phase noise at output (
(49) These output constraints must be complied with while taking account of the intrinsic performance of the hardware components at the level of the phase noise of the reference oscillator, L.sub.bref(f) (
(50) The templates, represented by continuous lines in
(51) In accordance with the H.sub. method, the constraints on the moduli of the frequency responses of the transfer functions S(s) and T(s) expressed by these templates are modeled by transfer functions termed weightings. In the case of a system like that of the invention, two sorts of weightings exist: those at input W.sub.i and those at output W.sub.o. The weightings can be chosen to be whitening filters, in which case W.sub.i is the inverse of the characteristic of PSD in terms of input noise and W.sub.o is the inverse of the characteristic of PSD in terms of output noise. More generally, we consider the case of a generic transfer function H(s), wherein the modulus squared of the frequency response must comply with the constraints imposed by the noise PSD masks at input and at output L.sub.in(f) and L.sub.out(f). We then write:
|H(j2f)|.sup.2<(|W.sub.o(j2f)W.sub.i(j2f)|.sup.1).sup.2=L.sub.out(f)/L.sub.in(f) hence:
|W.sub.o(j2f)H(j2f)W.sub.i(j2f)|<1.
(52) In the case of the PLL of
(53) Simplified constraints are illustrated dashed in
(54) As a variant, it is possible to use the phase error to impose the performance constraints on the phase noise. This leads to constraining the PSD of the error with respect to that of two signals of the noise (reference and VCO) by using an equivalent constraint on L.sub./L.sub.bref (the constraint on L.sub./L.sub.b is the same as that on L.sub./L.sub.b). In this case, too, it will be appropriate to simplify the constraints.
(55)
(56) The transfer function P(s) is given, since it depends on the topology of the PLL and G(s), as well as the weightings which have been determined in such a way as to satisfy the constraints on the phase noise. The transfer function F(s) is synthesized, in accordance with the H.sub. method, in such a way as to satisfy the condition
P(s).star-solid.F(s).sub.<1(4)
where the operator .star-solid. indicates the Redheffer product and .Math..sub. indicates the H.sub. norm. The search for a vector transfer function F(s) satisfying inequality (4) is a classical problem which can be solved, with the aid of a computer, by numerical calculation algorithms well known from the prior art; see for example [ScF-2009].
(57) In accordance with the invention, the synthesis of the correctors can be subject to other constraints, for example:
(58) robust stability, which imposes |S(j2f)|.sub.dB<6 dB f;
(59) limitation of the PSD of the slaving signal, or control signal, u(t).
(60) For the choice of the weightings making it possible to satisfy these constraints it will be possible to refer to [ScF-2009].
(61)
(62) In
(63) By applying this method of synthesis of the correctors to a particular case of specifications and for the topology of
(64)
The corresponding Bode plots are represented in
(65) Although these weightings exhibit an order of not greater than 2 (they correspond to simplified constraints), they nonetheless lead to transfer functions of the correctors of significant order; curves 140 and 141 in
(66) The transfer functions of the simplified correctors are:
(67)
(68)
(69) Curves 164, 166 and 168 are also reproduced magnified in
(70)
(71)
(72) Once the correctors have been synthesized by the method described hereinabove, the design of the PLL can be finalized, and then the PLL can be produced in integrated or discrete form by conventional electronics techniques.
(73) As was mentioned above, it is difficult to measure directly the phase noise signal of the reference or output signal of the PLL which must be provided as input to the corrector F.sub.2(s). For this reason, provision is advantageously made to use an estimator circuit, adapted for receiving as input a voltage signal (reference signal, or output signal of the PLL) and for providing the corrector F.sub.2(s) with an estimation of the phase noise affecting this signal (see
(74) a first branch comprising a delayer module, introducing a delay T.sub.d of the signal V(t);
(75) a second branch, connected in parallel with the first branch and comprising a phase-shifter module, introducing a phase shift of the signal V(t), with
(76)
where .sub.ref is the frequency of the periodic signal
(77) a mixer MIX configured to multiply the delayed signal V.sub.Td arising from the first branch and the phase-shifted signal V.sub. arising from the second branch; and
(78) a so-called rejector filter F.sub.rej(s), of low-pass type receiving as input the output signal V.sub.mix of the said mixer and providing at its output a signal V.sub.EXT constituting the sought-after estimation of the phase noise affecting the input signal V.
(79) From the formal point of view, the following demonstration presents the manner of operation of the phase noise estimator.
(80) It is considered that the input signal V(t) is a sinusoid of angular frequency .sub.ref affected by phase noise b.sub.ref(t):
V(t)=A sin(.sub.reft+b.sub.ref(t))(10)
(81) We then have:
V.sub.Td(t)=A sin(.sub.ref(tT.sub.d)+b.sub.ref(tT.sub.d))(11)
V.sub.(t)=A sin(.sub.reft+b.sub.ref(t))(12)
(82) Trigonometric calculations make it possible to show that:
(83)
(84) The component of V.sub.mix at the angular frequency 2 .sub.ref is intended to be filtered by the filter F.sub.rej; furthermore, relation (9) makes it possible to simplify expression (13). We can therefore write
(85)
(86) where V.sub.mixHF(t) is the high-frequency component of V.sub.mix which is filtered.
(87) In frequency representation, and by neglecting the component V.sub.mixHF(t), we can write
(88)
(89) In order for the output signal of the circuit to provide an image of the phase noise, it would be necessary for the filter F.sub.rej to invert the transfer function H(s); this is not possible since H(s) is non-causal, but the inversion can be carried out over a limited band, or H(s) can be linearized:
H(s).fwdarw.H.sub.lin(s)=T.sub.d.Math.s(16)
(90) We can therefore take
(91)
(92) where K is a low-frequency gain and a time constant. The rejector filter is therefore a low-pass filter, this having been admitted by precedence.
(93) As a variant, H(s) can be integrated within the augmented process P(s); in this case, it is taken into account during the synthesis of the correctors. Stated otherwise, the filter F.sub.rej can be integrated within one of the correctors of the PLL, in which case the phase noise estimation circuit may comprise only the delayer module Td, the phase-shifter module and the mixer MIX.
(94) The hardware embodiment of the delay and phase-shifter modules can depend on the PLL topology considered (see infra).
(95)
(96) The filter F.sub.rej and the correctors F.sub.1 and F.sub.2 can be embodied, for example, in the form of PID (Proportional-Integral-Derivative) correctors based on operational amplifiers, as is illustrated in
(97)
(98) Finally the mixer can be embodied by logic gates of NOR or XNOR type or a Gilbert cell.
(99)
(100)
(101) The invention is not limited to the case of an analogue PLL, but may also relate to a digital PLL, and notably an All-Digital PLL (ADPLL) such as described in the articles [STA-2004] and [STA-2005], as illustrated in
(102) In this ADPLL, a digitally controlled oscillator DCO produces a signal CKV which must be slaved in phase to a reference signal of lower frequency, FREF. The ratio of the frequency of the DCO to that of the reference is not necessarily integer; it is therefore necessary to generate a re-timed reference signal CKR, used as clock signal. This signal CKR is obtained by oversampling FREF by the high-frequency signal CKV with the aid of a flip-flop ECH1.
(103) An accumulator APO counts the rising (or falling) edges of the signal CKV; the counting signal thus obtained, indicated by R.sub.v[i], is re-timed to the tempo of the clock signal CKR with the aid of a flip-flop ECH2; the re-timed signal is indicated by R.sub.v[k].
(104) The signal R.sub.R[k] is obtained by accumulating, with the aid of an accumulator APR, a binary word FCW equal to the ratio of the frequencies of CKV to FREF.
(105) The signal [k] is obtained with the aid of a Time-to-Digital Converter TDC quantizing the temporal offset between the rising (or falling) edges of the signals FREF and CKV, whose output is re-timed to the tempo of the clock signal CKR with the aid of a flip-flop ECH3.
(106) An estimation .sub.E[K] of the phase error between FREF and CKV is given by:
.sub.E[k]=R.sub.R[k]R.sub.v[k]+[k](18)
(107) It should be pointed out that, in this context, the phase of a signal is given by the number of its rising (or falling) edges, counted from an initial instant.
(108) The error signal .sub.E[k] is provided as input to a first transfer function corrector F.sub.1(z). In accordance with the invention, the ADPLL also comprises a second corrector in feed-forward control mode, F.sub.2(z), receiving as input the signal [k], and/or a third corrector in feedback mode, F.sub.3(z), receiving as input the signal R.sub.v[k]. As in an analogue embodiment, the command signals generated by the correctors F.sub.1(z), F.sub.2(z) and/or F.sub.3(z) are combined in an adder node AN to obtain a (digital) signal for driving the oscillator DCO.
REFERENCES
(109) [Bib-1992] Johen E. Bibel, Stephen Malyevac Guidelines for the selection of weighting functions for H-infinity control, AD-A251 781 Naval Surface Warfare Center, January 1992, URL: http://www.dtic.mil/dtic/tr/fulltext/u2/a251781.pdf. [ScF-2009] Scorletti G. Fromion V. Automatique frquentielle avance, polycopie de cours, Sciences de l'ingnieur, http://cel.archives-ouvertes.fr/cel-00423848. [SkP-2005] S. Skogestad and I. Postlethwaite, Multivariable Feedback Control, Analysis and Design, John Wiley and Sons Chischester, 2005. [ChM-2007] W. Chaivipas, A. Matsuzawa, Analysis and Design of Direct Reference Feed-Forward. Compensation for Fast-Settling All-Digital Phase-Locked Loop. IEICE TRANS. ELECTRON., VOL.E90-C, NO.4 April 2007 [SNI-2012] B. Indu Rani, C. K. Aravind G. Saravana Ilango, C. Nagamani, A three phase PLL with a dynamic feed forward frequency estimator for synchronization of grid connected converters under wide frequency variations. Electrical Power and Energy Systems, 41:63-70, 2012. [STA-2004] R. B. Staszewski et al., All-Digital TX Frequency Synthesizer and Discrete-time Receiver for Bluetooth Radio in 130-nm CMOS, IEEE Journal of Solid-State Circuits, Vol. 39, No. 12, December 2004 [STA-2005] R. B. Staszewski, P. T. Balsara Phase-Domain All-Digital Phase-Locked Loop, IEEE Transaction on Circuits and SystemsII: Express Briefs, Vol. 52, No. 3, March 2005.