Research Article  Open Access
José R. C. Piqueira, Maurízio Q. de Oliveira, Luiz H. A. Monteiro, "Linear Approach for Synchronous State Stability in Fully Connected PLL Networks", Mathematical Problems in Engineering, vol. 2008, Article ID 364084, 13 pages, 2008. https://doi.org/10.1155/2008/364084
Linear Approach for Synchronous State Stability in Fully Connected PLL Networks
Abstract
Synchronization is an essential feature for the use of digital systems in telecommunication networks, integrated circuits, and manufacturing automation. Formerly, masterslave (MS) architectures, with precise master clock generators sending signals to phaselocked loops (PLLs) working as slave oscillators, were considered the best solution. Nowadays, the development of wireless networks with dynamical connectivity and the increase of the size and the operation frequency of integrated circuits suggest that the distribution of clock signals could be more efficient if distributed solutions with fully connected oscillators are used. Here, fully connected networks with secondorder PLLs as nodes are considered. In previous work, how the synchronous state frequency for this type of network depends on the node parameters and delays was studied and an expression for the longterm frequency was derived (Piqueira, 2006). Here, by taking the first term of the Taylor series expansion for the dynamical system description, it is shown that for a generic network with nodes, the synchronous state is locally asymptotically stable.
1. Introduction
Digital engineering technologies for communications, control, and computation require reliable clock distribution systems to guarantee the correct temporal order in the information processing by the several parts of a spatially distributed system [1–4]. Synchronization network is the general denomination of the part of the whole system responsible for this temporal order and the several possible solutions for its design are presented in [5].
Originally, masterslave architectures were used to distribute a precise clock signal generated by a master node to the other points of the systems where PLLs regenerate the phase and frequency information [1, 5].
The evolution of the telecommunication services to wireless and dynamical networks has shown the inadequacy of centralized clock distribution structures in these cases, motivating the study of fully connected architectures to generate reference signals with the phaselocked loops operating as nodes of the synchronization networks [1, 5–7].
Besides telecommunication networks, the main fields for the application of fully connected systems are time signal distribution in digital electronic circuits [2, 7–9] and wireless sensor networks [10]. Another very important application of networks is the implementation of oscillatory neuralcomputing devices, where the vector of phasedifferences amongst a group of synchronized oscillators is associated with some memory information [11, 12].
In this work, a fully connected node PLL network is studied, starting with a review of the model for a single node and the derivation of the synchronous state frequency [13]. Then, the dynamical equation for the phase differences between nodes is presented providing a model for the whole node network.
As it was shown in [13, 14], this kind of network presents a synchronous state that is reachable for any possible combination of node parameters. Consequently, the linear approximation gives practical hints about the lockin operation of the network [15]. Then, considering the first term of the Taylor series development around the synchronous state [16], it is shown that the synchronous state is locally asymptotically stable for any number of nodes.
2. Nodes in a Fully Connected Architecture and Synchronous State Frequency
In this section, the nodes of the fully connected network are analyzed, considering that they are secondorder phaselocked loops modeled as the classical analog version [17]. The results shown in [13] are summarized and the expression of synchronous state frequency is derived. This expression will be used to determine the stability conditions. PLL nodes are closed loops composed of a phase detector (PD), a lowpass filter (F), and a voltagecontrolled oscillator (VCO) [17]. The double frequency jitter in the output of the phase detector is neglected [18–20].
The nodes are geographically distributed and mutually coupled by using twoway connections. Each node exchanges clock signals with all the nodes such that . Each VCO belonging to a PLL is described by its freerunning frequency and by its phase . The whole network has its dynamics described by phase errors and frequency errors defined by and , respectively. Time delays corresponding to the propagation of signals from the node VCO output to the node PD input are considered, as shown in Figure 1, containing the model of the signal processed in node 1.
It can be noticed that in an node fully connected PLL network, the nodes need phase detectors [21] where the local VCO output signal is multiplied by the delayed signals sent by the VCO of the other nodes. As there is no integrated circuit with this architecture, the implementation of the nodes requires a combination of several PLL chips with the outputs of their phase detectors weighted and the result being the input of a single filter that feeds a single VCO. For the node 1 from Figure 1, the definitions are as follows: (i): delay from node to node ; (ii): delay from node to node ;(iii): delay from node to node ; (iv): VCO output of the node 1;(v): VCO output of the node 2 with a delay ;(vi): VCO output of the node 3 with a delay ;(vii): VCO output of the node with a delay ;(viii): output of PD 2;(ix): output of PD 3;(x): output of PD ;(xi): output of F that controls the VCO.
The instantaneous individual VCO phases and the instantaneous phase errors, and , are expressed in the forms
The node parameters related to PLL operation can be characterized by the following:(i)PD multiplying factors: , in , with ;(ii)gains of the VCOs: , in , with ;(iii) cutoff frequencies of F: , in , with . The constitutive parameters of the nodes are considered to be the same, in order to simplify the analytical reasoning. If the filters cutoff frequencies are different, but sufficient to avoid double frequency jitter [18], the results derived here are basically the same and only the acquisition times change [13]. If VCO and PD gains are considered to be different, the problem has to be numerically treated, and the expression of the synchronous state frequency derived here changes [22].
Under the assumption that all nodes have the same constitutive parameters, the output of each VCO isand signals received by the phase detector of node from node , with propagation delays , can be written aswhere is the controlled amplitude of the outputs of VCOs and PDs. Considering that each phase detector , belonging to node , is a signal multiplier [17]:replacing (2.2) and (2.3) in (2.4), and neglecting the double frequency components [18–20] as their frequency are much greater than the filter cutoff frequency, the output of each phase detector , belonging to node , is given byEach resulting signal given by (2.5) is multiplied by and added, in order to compose the filter input as below:Defining , expression (2.6) is simplified to
The filters are considered to be firstorder lowpass implying secondorder nodes. This choice is a common practice because secondorder PLLs always reach a synchronous state when submitted to phase steps and ramps [17] in spite of the doublefrequency jitter [19, 20]. If more accurate transient responses are necessary, secondorder filters are used, but complicated behaviors like bifurcation and chaos appear [23], worsening the operation.
Consequently, equations for the dynamics of the VCO phase are obtained by considering the filter transfer function [15], resulting in the expression
Replacing , given by (2.7), and the VCO control signal in (2.8), the equation for the node phase isThen, considering in (2.9) and defining and , the dynamics of each VCO phase in a fully connected network is given byEquation (2.10) is similar to the pendulum equation, containing a dissipation component , a delayed conservative term , and a forcing part [16]. Consequently, it is reasonable to suppose that the longterm solution of the system is a synchronous state, as shown in [14, 18], with the phases of all nodes oscillating with the same frequency that can be estimated. Thus, in order to estimate this frequency, the following hypotheses are considered [13]:(a),(b),(c).
Therefore,The values of are considered to be small because in the majority of the practical situations, the network is operating in the lockin mode [17]. Consequently, (2.11) can be written as a linear approximation considering . Hence, for each node ,Using (2.12) for an node network, with , and , as well as adding the resulting equations, as the sum of the terms is equal to zero because , one can writeCalculating from (2.13),Dividing (2.14) by , and replacing , the estimation of the synchronous state frequency () is obtained:Therefore, expression (2.15) is an estimation for the frequency of the synchronous state for a fully connected secondorder PLL network, depending on the individual freerunning frequencies and propagation delays. Notice that when the delays are zero, is given by the mean value of . As it was considered that all freerunning frequencies are different, there are phase shifts between the nodes in the synchronous state. In practical cases, as in communication networks, these phase differences are object of delay compensation techniques [5].
In previous work [13], numerical simulations were conducted to investigate the accuracy of expression (2.15) and to analyze how gains and delays change the behavior of the network. Here, the dynamic equations for the phase differences are derived allowing the analytic study of the local stability of the synchronous state.
3. Phase Difference Equations
In this section, the equations describing the dynamics of the phase errors, , are derived. A set of secondorder ordinary differential equations is obtained, expressing the phase differences between all the nodes and the node , taken as reference.
Starting with the individual dynamical equations for each VCO phase and expressing the differences between the nodes and node , equations for , are written as below: Using the condition and considering , then Expressing the differences with node as reference, Replacing the terms , , , and using the identity in the former expressions, the VCO phase differences become These equations give the general expression that describes the dynamical behavior of the phase differences for the network VCOs as below:The nonlinear differential equation (3.5), for smallphase deviations, can be approximately expressed by the linear term of the Taylor series expansion [16]. Then, rewriting (3.5)Expressions (3.5) and (3.6) describe the dynamics of the phase adjustments of a fully connected PLL network depending on the PLL node parameters and , the number of nodes , as well as the individual freerunning frequencies and the delays. These equations allow the research of the synchronous state stability that is conducted in the next section.
4. Synchronous State Stability
In Section 2, an expression for the synchronous state frequency for the fully connected network was derived and, in this section, the stability of the synchronous state is studied, that is, if the reachable synchronous state is robust under small perturbations. The analysis is performed considering that the solutions of the linear equation (3.6) can be topologically equivalent to the solutions of (3.5) in a small neighborhood of the synchronous state if it is a hyperbolic equilibrium point [16]. The procedure will be shown for threenode and fournode networks providing the identification of patterns in the expressions for the eigenvalues of the Jacobian matrix representing the linear equivalent system around the synchronous state.
4.1. ThreeNode Network
For a threenode network, the phase differences are and . By using (3.5), the system is described by The system composed of two secondorder equations (4.1) can be described by the state variables:Considering that , the following firstorder equations result:
Consequently, the Jacobian matrix from (4.3) isThe eigenvalues of are the roots of the characteristic polynomial which have multiplicity two, and are given by
4.2. FourNode Network
The calculation of the eigenvalues for the fournode network follows the same procedure of the former case. The phase differences are represented by , , and . By using (3.5), the system is written as The secondorder expressions (4.7) can be rewritten in terms of state variables:Considering that , , and , the following firstorder equations result:
Consequently, the Jacobian matrix from (4.9) isThe eigenvalues of are given by the roots of the characteristic polynomial: which have multiplicity , and are
4.3. NNode Network
As it was shown, for a fully connected PLL network, the Taylor series development around the synchronous state results in Jacobian matrices with a canonical Jordan form [24]. Table 1 shows the characteristic polynomials and the expressions of the eigenvalues and their multiplicities, m, corresponding to fully connected networks dynamic equations, around the synchronous state, for different number, , of nodes.

It can be observed that for any number of nodes, the eigenvalues from Table 1 are with negative real parts, implying that the synchronous state is locally asymptotically stable.
5. Conclusions
As the fully connected architecture started to be used in large scale in clock distribution systems, (2.15) can be applied by network designers as an estimation for the frequency of the synchronous state when secondorder PLLs are used to extract the timing information in the nodes.
Besides, because the synchronous state is a hyperbolic locally asymptotically stable equilibrium point, the network can recover synchronization after perturbation. The eigenvalues presented in Table 1 can be used for determining the transient behavior of such a resynchronization process.
Acknowledgment
José R. C. Piqueira is supported by FAPESP and CNPq. Luiz H. A. Monteiro is partially supported by CNPq
References
 S. Bregni, Synchronization of Digital Networks, John Wiley & Sons, Chichester, UK, 1st edition, 2002.
 H. A. Tanaka, A. Hasegawa, H. Mizuno, and T. Endo, “Synchronizability of distributed clock oscillators,” IEEE Transactions on Circuits and Systems I, vol. 49, no. 9, pp. 1271–1278, 2002. View at: Publisher Site  Google Scholar
 M.H. Hung, J. Tsai, F.T. Cheng, and H.C. Yang, “Development of an Ethernetbased equipment integration framework for factory automation,” Robotics and ComputerIntegrated Manufacturing, vol. 20, no. 5, pp. 369–383, 2004. View at: Publisher Site  Google Scholar
 P. Kellett, “Beyond the LAN: Ethernet's evolution into the public network,” Fiber and Integrated Optics, vol. 22, no. 5, pp. 289–303, 2003. View at: Publisher Site  Google Scholar
 W. C. Lindsey, F. Ghazvinian, W. C. Hagmann, and K. Dessouky, “Network synchronization,” Proceedings of the IEEE, vol. 73, no. 10, pp. 1445–1467, 1985. View at: Google Scholar
 D. J. Messerschmitt, “Synchronization in digital system design,” IEEE Journal on Selected Areas in Communications, vol. 8, no. 8, pp. 1404–1419, 1990. View at: Publisher Site  Google Scholar
 M. SaintLaurent and M. Swaminathan, “A multiPLL clock distribution architecture for gigascale integration,” in Proceedings of the IEEE Computer Society Workshop on VLSI, pp. 30–35, Orlando, Fla, USA, May 2001. View at: Publisher Site  Google Scholar
 V. Gutnik and A. P. Chandrakasan, “Active GHz clock network using distributed PLLs,” IEEE Journal of SolidState Circuits, vol. 35, no. 11, pp. 1553–1560, 2000. View at: Publisher Site  Google Scholar
 M. SaintLaurent, P. ZarkeshHa, M. Swaminathan, and J. D. Meindl, “Optimal clock distribution with an array of phaselocked loops for multiprocessor chips,” in Proceedings of the 44th IEEE Midwest Symposium on Circuits and Systems (MWSCAS '01), vol. 1, pp. 454–457, Dayton, Ohio, USA, August 2001. View at: Publisher Site  Google Scholar
 M. M. Nordman, W. E. Kozlowski, and O. Vahamaki, “A method for synchronizing low cost energy aware sensors used in industrial process monitoring,” in Proceedings of the 27th Annual Conference of the IEEE Industrial Electronics Society (IECON '01), vol. 1, pp. 100–106, Denver, Colo, USA, NovemberDecember 2001. View at: Publisher Site  Google Scholar
 F. C. Hoppensteadt and E. Izhikevich, “Pattern recognition via synchronization in phaselocked loop neural networks,” IEEE Transactions on Neural Networks, vol. 11, no. 3, pp. 734–738, 2000. View at: Publisher Site  Google Scholar
 J. R. C. Piqueira, F. M. Orsatti, and L. H. A. Monteiro, “Computing with phase locked loops: choosing gains and delays,” IEEE Transactions on Neural Networks, vol. 14, no. 1, pp. 243–247, 2003. View at: Publisher Site  Google Scholar
 J. R. C. Piqueira, M. Q. de Oliveira, and L. H. A. Monteiro, “Synchronous state in a fully connected phaselocked loop network,” Mathematical Problems in Engineering, vol. 2006, Article ID 52356, 12 pages, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 L. H. A. Monteiro, N. C. F. Canto, J. G. ChauiBerlinck, F. M. Orsatti, and J. R. C. Piqueira, “Global and partial synchronism in phaselocked loop networks,” IEEE Transactions on Neural Networks, vol. 14, no. 6, pp. 1572–1575, 2003. View at: Publisher Site  Google Scholar
 J. R. C. Piqueira and L. H. A. Monteiro, “Allpole phaselocked loops: calculating lockin range by using Evan's rootlocus,” International Journal of Control, vol. 79, no. 7, pp. 822–829, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 J. Guckenheimer and P. Holmes, Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields, vol. 42 of Applied Mathematical Sciences, Springer, New York, NY, USA, 2nd edition, 1983. View at: Zentralblatt MATH  MathSciNet
 R. E. Best, PhaseLocked Loops, McGraw Hill, New York, NY, USA, 4th edition, 1999.
 J. R. C. Piqueira, E. Y. Takada, and L. H. A. Monteiro, “Analyzing the effect of the phasejitter in the operation of second order phaselocked loops,” IEEE Transactions on Circuits and Systems II, vol. 52, no. 6, pp. 331–335, 2005. View at: Publisher Site  Google Scholar
 J. R. C. Piqueira, A. Z. Caligares, and L. H. A. Monteiro, “Doublefrequency jitter figures in masterslave PLL networks,” AEU  International Journal of Electronics and Communications, vol. 61, no. 10, pp. 678–683, 2007. View at: Publisher Site  Google Scholar
 J. R. C. Piqueira and L. H. A. Monteiro, “Considering secondharmonic terms in the operation of the phase detector for second order phaselocked loop,” IEEE Transactions on Circuits and Systems I, vol. 50, no. 6, pp. 805–809, 2003. View at: Publisher Site  Google Scholar
 R. B. Pinheiro, J. J. da Cruz, and J. R. C. Piqueira, “Robust clock generation system,” International Journal of Control, vol. 80, no. 1, pp. 35–44, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 F. M. Orsatti and J. R. C. Piqueira, “Evolutionary multiobjective optimization in mutually connected PLL synchronization networks. Modeling and control of autonomous decision support based systems,” in Proceedings of the 13rd International Workshop on Dynamics and Control, vol. 1, pp. 255–262, Shaker, Wiesensteig, Germany, May 2005. View at: Google Scholar
 L. H. A. Monteiro, D. N. F. Filho, and J. R. C. Piqueira, “Bifurcation analysis for thirdorder phaselocked loops,” IEEE Signal Processing Letters, vol. 11, no. 5, pp. 494–496, 2004. View at: Publisher Site  Google Scholar
 M. W. Hirsch, S. Smale, and R. L. Devaney, Differential Equations, Dynamical Systems, & An Introduction to Chaos, vol. 60 of Pure and Applied Mathematics, Elsevier/Academic Press, Amsterdam, The Netherlands, 2nd edition, 2004. View at: MathSciNet
Copyright
Copyright © 2008 José R. C. Piqueira et al. 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.