System for highly repeatable clock parameter recovery from data modulated signals

Abstract: A system for recovering parameters of an underlying clock signal from data modulated signals (12) includes a time digitizer (22), or time stamping device, for sampling and locating edges of the data modulated signals, thereby selectively providing time stamp readings identifying incoming edges of the data signals. The time stamp readings are used, along with a user-input nominal clock frequency (36), to determine (30) an array of actual clock cycle stamps (46), each representing an integral number of actual clock cycles occurring between adjacent time stamps (48). The array is then used in conjunction with the time stamps to characterize parameters (46,48) of the underlying clock data. Once the clock parameters are recovered, a wide variety of information, for example, frequency and phase information, including jitter, can be determined (20) with a very modest rate of sampling compared to the high data bandwidth.

Description: This invention relates to a method and system for recovering parameters of a clock, for example to a system for highly repeatable clock parameter extraction from data modulated signals. In digital transmission and communications, a wide variety of modulation schemes are used to encode data. For example, in a random NRZ data signal, data are represented by a constant signal level for the full duration of a bit interval, i.e., one clock cycle. By contrast, RZ (return-to-zero) data signals return to zero with each bit interval. In each case, all data transitions occur at clock transitions. However, not all clock transitions generate data transitions in NRZ data. To decode the digital data under NRZ, RZ, or other modulation schemes, it is necessary, in general, to know the correct instances to sample the data. Preferably, the clock used to encode the data is available for decoding. However, in many cases, the encoding clock is not available, and a new clock must be synthesized at the receiving end. US-A-5260978 discloses a method of recovering the timing clock of a network source input signal. At the end of each time stamp period a time stamp is transmitted from the source, and at a destination node there is determined the number of network clock cycles in the corresponding time stamp period to recover the timing clock. The present invention seeks to provide improved recovery of clock parameters. According to an aspect of the present invention, there is provided a method of recovering parameters of an underlying clock as specified in claim 1. According to another aspect of the present invention, there is provided a system for recovering data clock parameters as specified in claim 7. It is possible with preferred embodiments of the invention to satisfy a need of system operators and system component manufacturers, who are not directly interested in the exact content of the data messages being sent and received, to know or monitor the frequency and phase behavior, both average and instantaneous, of the underlying clock of the data stream. Such measurements are particularly desirable in plexiosynchronous digital systems in which deliberate phase shifts are occasionally introduced to compensate for small clock frequency differences at different points of the system. Preferably, such measurements should be made in "live traffic", that is, made without interrupting the system operation. Furthermore, it is possible to make such measurements at relatively low speed, commensurate with the clock jitter bandwidth, in the presence of much higher data modulation bandwidth. For these measurements, it is not necessary literally to recover a physical clock. Also, knowledge of the underlying clock frequency and phase can be used to monitor and control system performance and to predict bit error rate (BER), especially when the rate is too low for convenient measurement. Systems may be specified to tolerate a certain amount of jitter at various frequency components (jitter tolerance), and, therefore, the actual jitter at various system points should be measured (jitter generation). Jitter is considered as phase deviation above 10 Hz bandwidth. To guard against magnification of jitter when cascading components, the transfer of jitter from input to output should be characterized (jitter transfer). Ideally, these measurements are made during data modulation to be more realistic and, if possible, during actual system operation. When incompatibility does arise in system operation, a repeatable standard method of measurement can be provided to settle disputes between vendors and buyers, system operators, and service providers, as to whether each piece of equipment is within specifications, and, indeed, to determine if the specifications are adequate to ensure correct operation. Although widely used, the conventional methods using phase lock loop (PLL) and surface acoustic wave (SAW) oscillators suffer from lack of accepted standards in PLL design regarding loop components, and loop dynamic filter and SAW detail performance parameters. In practice, the physical clocks synthesized likely differ from one design to another, even for identical data streams, thereby producing non-unique and non-repeatable measurements. This lack of repeatability is true even if the more serious PLL errors such as loss-of-lock or false locking have been overcome. The preferred embodiment can be provide a system for highly repeatable measurement based on an accurate representation of the underlying clock, independent of any particular circuit design detail or parametric performance of various components. In arrangements according to the invention, there is provided a system for highly repeatable clock parameter recovery from data modulated signals. The system utilizes readings of a conventional time stamping device for selectively sampling and locating edges of the data modulated signals in the system. The time stamp readings are then used, along with a user-input nominal clock frequency, to calculate an array of actual clock cycle stamps, each representing an integral number of actual clock cycles occurring between adjacent time stamps. The array is used, in conjunction with the time stamps, to characterize the underlying data clock. Once the clock parameters are recovered, frequency and phase information, including jitter, can be determined. Since the outcome of the cycle estimation process comprises whole numbers (integers), barring gross errors analogous to false locking and loss-of-lock in a phase lock loop (PLL), the outcome is the same, independent of any particular implementation detail. Any non-uniqueness comes from the time stamping process, but these are typically very small, being on the order of 50 picoseconds or less. Once a clock is effectively recovered, jitter and other phase-frequency information from the underlying clock can be obtained. The process can be accomplished with very modest sampling rates compared to the high data bandwidth. An embodiment of the present invention is described below, by way of example only, with reference to the accompanying drawings, in which: FIG. 1 is a general block diagram of an embodiment of a system for measurements based on highly repeatable clock parameter recovery; FIG. 2 is a block diagram of a data clock parameter recovery section of the system of FIG. 1; FIG. 3 is a graphical representation of the time stamp and cycle stamp parameters used to characterize the data clock; FIG. 4 is a block diagram of an alternate dual stage embodiment of the clock parameter recovery section of FIG. 2; FIG. 5 is a block diagram of an embodiment of the processing section of the system of FIG. 1; and FIG. 6 is a hardware representation of a real time measurement system based on FIG. 2. Referring to FIG. 1, the highly repeatable clock parameter recovery system shown includes a clock parameter recovery section 10 for extracting parameters of an underlying clock from data modulated signals 12 input into the system. The recovered clock parameters are fed to a processing section 20 for determining a wide variety of measurements through digital filtering techniques. These measurements include, but are not limited to, phase deviation (jitter), frequency deviation, and average frequency, phase-spectral response, wander, TVAR, jitter transfer, etc. The measurements may then be displayed or used, either manually or automatically, to adjust the overall system application through system display/controls 21. Referring to FIG. 2, the clock parameter recovery section 10 includes a conventional time interval analyzer 22, also known as a time digitizer, or time stamping device. A suitable time interval analyzer, commercially available from Hewlett-Packard Co., Palo Alto, California, under the product names HP1725A, HP53310, and HP5372, sequentially samples and provides a high precision time stamp of selected data edges of the incoming data signal. The data signals are received by the device at input 24, and the time stamp information is delivered at an output 26. By time stamping, it is meant that the analyzer samples the data signals, locates at least selected edges of the data signal, and identifies when each such edge occurs. The time stamping device includes a built-in counter (not shown) referenced to a highly accurate reference time base oscillator 28 which is more precise than the underlying data clock. The oscillator may be, for example, a quartz oscillator. The time stamping device can be controlled to selectively sample the rising or falling edges of the data signal, or both. The sampling rate of the device can also be controlled based on a user-defined pace setting or cycle time. Pace is one method of controlling the sampling rate of the device by having it count a specific number of edges between samples. The sampling rate of other time stamping devices can also be regulated on a real time basis, e.g., by a user-defined cycle time, during which the device is enabled to sample data after a set time interval. Referring to both FIGS. 2 and 3, an example of the operation of the time stamping device is shown on the graph of the incoming data modulated signal 12, in this case, an NRZ signal. The time stamping device identifies every other rising edge of the data signal and records the occurrence of such edges in time stamps t i for i=1 to N, where N = the total number of time stamps taken. The time stamps t i are (shown graphically in FIG. 3) made on every other rising edge of the data signal according to the pace setting (P=2) of the device. The time stamps have values corresponding to the precise reference time based oscillator of the stamping device. (As shown in FIG. 3, the time stamps have values t 1 =0ns, t 2 =223.45ns, t 3 =558.75ns, and t 4 =804.80ns.) The set of time stamps t i to t N are then collected and recorded, as are a series of data stamps e i . In clock parameter extraction, no use is made of the data stamps e i , which, in this case, are regular with an increment of P per stamp. (As shown in FIG. 3, the data stamps have values of e 1 =1, e 2 =3, e 3 =5, e 4 =7, etc.) Instead, the gathered time stamps t i are used to determine an array of actual clock cycle stamps z i , each actual clock cycle stamp representing an integral number of actual clock cycles occurring between adjacent time stamps. The actual clock cycle stamps may be calculated and compiled in an array by any suitable hardware or software-driven computer. To calculate the array of actual clock cycle stamps, the actual clock cycle stamp change ▵z i (phase change in unit interval or u.i.) is first computed on the basis of the time stamps t i , t i+1 , and a nominal clock frequency f nom provided through input 36, e.g., a keyboard. The actual clock cycle stamp change ▵z i is derived by inserting the nearest integral number of nominal clock-cycles between time stamps according to the following formula: Δz i = INT[f nom · (t i+1 -t i )+0.5] The actual clock cycle stamp change ▵z i is equal to the phase change in unit intervals (of the counter) between one time stamp t i and an adjacent stamp t i+1 . The integer (INT) function truncates any non-integer value (e.g., 4.2 becomes 4; 4.8 becomes 4). By adding 0.5 to the equation, the INT function rounds the operand value to the nearest integer, thereby converting the truncating into a round-off function. (As shown in FIG. 3, the z i values are all integers.) This process of calculating the actual clock cycle stamp changes is performed for all of the ▵z i values from i=1 to N, synthesizing an overall z i array from the t i array of time stamps. The z i values can be determined by sequentially adding the particular ▵z i values to an arbitrary reference value z o ; in this example, z o is equal to 1. (As shown in FIG. 3, the actual clock cycle stamps have values of z 1 =1, z 2 =11, z 3 =26, and z 4 =37). The two arrays of parameters t i and z i sufficiently characterize the underlying clock by providing clock cycle data t i extracted directly from the data signal and comparing this data to whole number (integer) actual clock cycle stamp Z i reference values. The system described is highly repeatable and operable over a wide and well-defined range. The system is highly repeatable because the outcome of the underlying clock cycle estimation process comprises only whole numbers, and barring gross errors analogous to false locking and loss of loop in phase lock loops, the outcome is the same independent of any particular implementation detail. Any non-uniqueness arises from limitations in timing resolution of the time stamps, but this is typically under 50 picoseconds. To avoid clock parameter recovery failures analogous to lose-lock in a phase lock loop, and for the system to operate properly, the nominal frequency must be sufficiently close to the actual frequency (f act ) of the underlying clock sought to be recovered by the system. More particularly, a necessary condition for successful operation is that the following error formula be satisfied: |▵f·▵t| < ½ where ▵f is the frequency error |f nom -f act | between the nominal and actual frequencies, and Δt is the time interval (t i+1 -t i ) being fitted by a actual clock cycle stamp change ▵z i . When the error formula is violated, this means that ▵z i may be rounded off to the wrong integer. As is evident from the formula, short time intervals can tolerate larger frequency errors. Conversely, large modulation on the clock, and/or an insufficient sampling rate, may lead to clock parameter recovery failure. In order to meet the requirements of the error formula, one can choose a smaller measurement time (smaller pace setting) and larger sampling size (for a higher sampling rate). The pace setting of the time stamping device may also be adjusted to change ▵t. Therefore, one should set the pace value of the time stamping device based on the characteristics of the underlying data clock to be recovered. A pace value of one is best for accurate results, but suffers from the constraint of available memory to store the array of measured time stamps and calculated actual clock cycle stamps. Therefore, such a low pace value should be used when the system involves a fluctuating or unstable underlying data clock, or when the nominal frequency input cannot be adequately characterized to meet the error formula. If the data clock is relatively stable, and the nominal frequency can be accurately predicted to be sufficiently close to the actual frequency, a low pace value is not required. In this case, one should choose the pace value simply to meet the requirements of Nyquist's theorem, that is, the sampling rate of the time stamping device should be at least twice the actual jitter bandwidth of the underlying clock. This can be significantly less than the clock frequency, thousands and even millions of times less. Once the array of actual clock cycle stamps z i is calculated, the combination of z i and t i can be used in a computer or other hardware means to provide a variety of information concerning the data signal. Instantaneous frequency F i , for example, is obtained by dividing the phase change ▵z i by the time difference ▵t i between adjacent stamps according to the following formula: F i = ΔZ i Δt i In essence, the instantaneous frequency F i is a corrected frequency derived from a predictor frequency F nom . Rewriting F i as F corrected , and F nom as F predictor , the following equations are derived: Δt = t i+1 -t i , F corrected = INT[F predictor ·Δt+0.5] Δt i ΔZ i = Δt i ·F corrected The number of actual clock cycles is obtained by rounding off to the nearest integer the number of nominal clock cycles occurring within Δt i . It should be noted that ▵z i , the actual clock cycle stamp changes, are integers, each change representing the nearest whole number of actual clock cycles occurring within a time interval Δt i . It should also be noted that other frequencies besides an arbitrary nominal frequency value f nom can be used as a predictor frequency, for example, an adjacent instantaneous frequency F i-1 may effectively be used. A dual stage data clock parameter recovery section as illustrated in FIG. 4 embodies a system for predicting frequencies for longer time intervals which may be prone to cause violation of the error formula and instability of the system. First, the data modulated signal is time stamped by the time stamping device 22. The time stamp intervals Δt between adjacent time stamps are then computed and sorted according to length of the interval. The relatively shorter intervals ▵t i(short) are suitable for actual clock cycle stamp changes ▵z i based on a nominal clock frequency obtained from input 36. A set of instantaneous frequencies F i is then determined based on the shorter values of ▵t i and corresponding values of ▵z i . These, or interpolated or extrapolated versions of these, are then, in turn, used as nominal clock predictor frequencies to predict relatively moderate and, finally, relatively longer time intervals ▵t i(long) . The actual clock cycle stamp changes ▵z i are then sequentially accumulated to obtain a complete array of z i values for the entire data signal. This bootstrap method, therefore, allows successful clock recovery under larger modulation than can be obtained using a user-input nominal clock frequency f nom alone as a predictor. Referring to FIG. 5, the phase and frequency deviation measurement section (one exemplary embodiment of the processing section 20 of the system), receives the arrays of time stamp (t i ) and actual clock cycle stamp (z i ) data and uses the arrays to calculate phase and frequency information of the data modulated signals 12. By using the effectively recovered clock parameters, this information concerning the underlying clock can be obtained with a very modest rate of sampling compared to the high data bandwidth, and no literal physical clock need be recovered. Average frequency (f avg ) of the data clock is determined by first synthesizing an overall array of actual clock cycle stamps z i by sequentially summing all of the ▵z i values as follows: The initial value of z i is arbitrary and of no significance. Once the array of actual clock cycle stamps is determined, a curve fitter 46 is used to determine average frequency through linear regression, for example. Other curve fitting schemes may be used. The curve fitter simply fits the best straight line in the min-square sense to the z i and t i data over the data range of N points. The slope of the line yields the least square average frequency. Phase deviation, or jitter,  is determined on the basis of the average frequency (f avg ) according to the formula:  i = z i -[ o +f avg ·t i ] where  i = phase deviation, z i = the actual clock cycle stamp,  o = the initial phase deviation, and t i is the corresponding time stamp reading. The calculated phase deviation  i is then passed through a digital filter 48 in order to produce a more precise reading. Frequency deviation is determined by merely differentiating the phase deviation measurements over time using a differentiator 50. A wide variety of other measurements may also be made in the digital processor 48 as a function of the recovered underlying clock parameters using digital filtering techniques. For example, a high pass filter may be used to separate jitter and wander measurements. (Jitter is phase deviation above 10 Hz bandwidth; wander is phase deviation below 10 Hz bandwidth.) Other measurements may also be made including, for example, phase-spectral density, Allan variance, TVAR, and TDEV. Fourier transforms may be used to obtain other useful data. Some or all of the recovered data can be displayed or used to control the overall system through controls 21 (FIG. 1). Although described thus far in terms of characterizing parameters of an underlying clock for determining phase and frequency of the data signal, the system may also be used during the course of decoding the received data signals using the measured array of time stamps t i and the calculated array of actual clock cycle stamps z i . To decode, however, every edge of the data signal must be sampled to satisfy the Nyquist rate by adjusting the pace value of the time stamping device. Once every edge is sampled, and characterized as rising or falling, the data signal can be reconstructed, and the data read by a decoder 48. In addition to decoding, when every edge of the data modulated signal is sampled, clock-to-clock jitter and cycle-to-cycle jitter measurements can be made. It should be noted that the clock recovery system can operate from time stamps stored in memory or in real time. Using the clock parameter recovery system in situations where the time stamps are stored in memory, the system can be used to scan the signal and determine an accurate predictor frequency based on a small ▵t. Then, the data can be read backward or forward, as the case may be, to sample the rest of the data. This is a significant advantage because such computations yield more accurate results (less failure) than simply a forward-looking method. The system of FIG. 4 is suitable for calculations based on time stamps stored in memory 41. This is the preferred embodiment. In real time clock parameter extraction, only a limited number of consecutive time stamps are available at any given time in a dynamic fashion. FIG. 2 illustrates an embodiment whereby only two adjacent time stamps t i (current) and t i-1 (previous) from time stamping device 22 are available for subtraction to generate Δt i =t i -t i-1 and for fitting actual clock count difference Δz i (34) via the use of blocks 30 and 36. After summing Δz i to generate z i (46) by accumulator 40, both t i (48) and z i (46) are outputted for digital signal processing. In the next sample, the index i is incremented by one. The current time stamp then becomes the previous time stamp and the newly acquired time stamp then becomes the current time stamp. Then the process repeats. The system still can be useful to calculate underlying clock parameters and phase-frequency information within a relatively modest sampling rate compared to the high data bandwidth. The system of FIG. 2 is suitable for calculations based on real time data signal input. Referring to FIG. 6, in the hardware embodiment shown a data modulated signal, such as an NRZ signal, is paced by a pacer, e.g., a programmable hardware frequency divider 104 by a pace factor P. The paced signal is used to trigger a time digitizer 108 which produces an n-bit digital output T i , which is a digital representation of the physical time t i . The least significant bit of T i represents the time-resolution of the time digitizer t nor , typically well below 1 ns. The mathematical relationship between t i and T i is simply t i =T i ·t nor . Register 112 is a buffer memory which records the result of the previous measurement T i-1 . Both current result T i and previous result T i-1 are fed into accumulator 120 which subtracts the two consecutive results to form an m-bit difference ΔT i =T i -T i-1 . The size of m is chosen to accommodate the largest time difference allowed. The m-bit difference is fed into a multiplier 126 which multiplies the difference ΔT i by a number equal to the product of the nominal frequency F nom and t nor , and generates a q-bit result. The effect of this multiplication is the same as multiplying the physical time Δt i by the frequency F nom , since ΔT i · (F nom ·t nor ) =(ΔT i ·t nor ) · F nom =Δt i ·F nom . The least significant bit of this result should be 1/2, and further refinement has no benefit. The round-off process is represented by adding 1/2 to the q-bit result by adder 130 and then truncating the least significant bit to give a q-1 bit result of Δz i =z i -z i-1 . This Δz i difference is accumulated by an accumulator 136 to give the sequence z i . A buffer 142 sequences the z results for the proper timing for accumulation. Just prior to being triggered, the buffer feeds the result z i-1 to accumulator 136, but shortly after trigger, the result z i will appear at the output of the buffer. Both T i and z i are outputs of the clock parameter recovery process. The time stamps T i , strictly speaking, should be multiplied by t nor before being outputted as physical time stamps t i . But in practice, T i is used in further processing for various applications rather than t i . Normalization is done when all processes are completed.

Claims: A method of recovering parameters of an underlying clock from data signals (12) comprising the steps of: selectively sampling the data signals (22) to produce a plurality of time stamps; inputting a nominal clock cycle frequency (36); and determining a nearest integral number of actual underlying clock cycles (30) between adjacent time stamps as a function of the nominal clock cycle frequency and time stamps and characterizing parameters of the underlying clock by determining (40) an array of the time stamps (48) and an array of actual clock cycle stamps (46), each actual clock cycle stamp representing an integral number of actual clock cycles occurring between adjacent time stamps. A method according to claim 1, comprising the step of determining the average frequency (20) of the underlying clock by fitting a curve (46) through the array of characterized parameters. A method according to claim 1 or 2, comprising the steps of: determining instantaneous frequency (20) of the underlying clock as a function of the characterized clock parameters; and inputting the instantaneous frequency (20) as the nominal clock cycle frequency for determining additional elements of the array of actual clock cycle stamps. A method according to claim 3, comprising the step of computing instantaneous frequencies (20,41) determined from short time intervals and extrapolating the computed instantaneous frequencies before the inputting step of determining actual clock cycle stamps of long time intervals. A method according to any preceding claim, comprising the step of controlling system performance (21) as a function of recovered clock parameters. A system for recovering data clock parameters from a data modulated signal comprising: time stamping means (108) for producing time stamps of selected portions of an incoming data modulated signal; first accumulator means (120) for accumulating the time stamps and determining time intervals between selected time stamps; counter means (126,130) for determining a number of digital counts occurring within selected time intervals; and second accumulator means (136) for accumulating digital count data and time stamps and having an output (142) for producing actual clock cycle stamps for characterizing the data clock. A system according to claim 6, wherein the time stamping means comprises pace setting means (104) for controlling the sampling rate of the incoming data modulated signal. A system according to claim 6 or 7, wherein counter means comprises a multiplier (126) and a nominal frequency input. A system according to claim 6, 7 or 8, wherein the counter means comprises rounding means (130) for determining the nearest whole number of underlying clock cycles occurring within each time interval. A system according to any of claims 6 to 9, comprising buffer means (112) for buffering incoming data coupled to the time stamp means. A system according to any of claims 6 to 10, comprising decoding means (48) for decoding the underlying clock data as a function of the accumulated digital count data and time stamps. A system according to any one of claims 6 to 11, comprising processing means (48,50) for characterizing phase-frequency information of the underlying clock as a function of the accumulated digital count data and time stamps.

Inventors: CHU DAVID CHAU-KWONG, CHU, DAVID CHAU-KWONG
Applicants: HEWLETT PACKARD CO, HEWLETT-PACKARD COMPANY
Application year: 1995
Main classification: H04L7
IPC codes: H04L7, H04L1, G01R23