Digital phase locked loop with integer channel mitigation

Abstract
An embodiment of the present invention provides a phase locked loop that operates on clock signals derived from an RF clock signal generated by the phase locked loop. A frequency reference input provides a reference clock. A controllable oscillator generates the RF clock signal with a plurality of phases. A switch is coupled to receive the RF clock, and is operative to select one of the plurality of phases. A phase detection circuit is coupled to the switch and is operable to receive a selected phase and to provide digital phase error samples indicative of a time difference between the reference clock and the selected phase.
Description
FIELD OF THE INVENTION

This invention generally relates to the field of control systems and data communications. In particular, it relates to cellular telephony and communication devices such as Bluetooth, WLAN, etc. using all-digital radio frequency (RF) circuitry.


BACKGROUND OF THE INVENTION

With each successive cellular phone handset generation, users demand more features in a smaller form factor. Some recent examples include cell phones with integrated Bluetooth, GPS, digital camera, and MP3 functionality. Process shrinks help deliver a cost and size advantage for digital designs with relative ease. However, for analog/RF designs, the immaturity of advanced processes comes with design challenges that may outweigh the intended advantage. In a typical handset, 30 to 40% of handset board space is occupied by analog/RF functionality which cannot be re-designed or migrated to the newer process/technology nodes easily, inhibiting vendor ability to cost effectively add features and reduce footprint.


Digital radio has recently allowed the replacement of space consuming analog RF circuitry with much more compact digital circuitry, thereby facilitating the ability to port designs rapidly to more advanced lithographies. Texas Instruments (TI) has proven this concept with its Digital RF Processor (DRP™) architecture, which it has successfully implemented in production versions of its Bluetooth BRF6xxx transceivers, GSM/GPRS LoCosto TCS23xx transceivers among other chips. DRP implementation is consistent with the on-going trend toward RF-CMOS in the cellular area, making it attractive in terms of power consumption, cost, and the integration of multiple radios.


Oscillators are a key component in the design of radio frequency (RF) communication systems. The estimation and calibration of the modulation gain of an RF oscillator is currently an area of active research. Accurate knowledge of this gain significantly reduces the complexity and increases the performance of the phase-locked loop (PLL) as well as the transmit frequency modulation path. It is particularly beneficial in systems implemented in deep submicron CMOS and based on orthogonal frequency/phase and amplitude (i.e. polar) topology. Estimation of RF oscillator frequency-modulation gain is especially important in low-cost dominantly digital high-volume transceivers. In such systems, the phase locked loop sets the loop bandwidth while the transmitter sets the transfer function of the direct frequency modulation path wherein the acceptable gain estimation error ranges from less than 1% for LTE/WCDMA to several percents for EDGE, GSM and Bluetooth, for example.


An all-digital frequency synthesizer architecture built around a digitally controlled oscillator (DCO) that is tuned in response to a digital tuning word (OTW) is described in U.S. Pat. No. 7,046,098 entitled “All-digital frequency synthesis with capacitive re-introduction of dithered tuning information” and is incorporated by reference in its entirety herein. A gain characteristic (KDCO) of the digitally controlled oscillator can be determined by observing a digital control word before and after a known change Δ.fmax) in the oscillating frequency. This has been described in U.S. patent application Ser. No. 11/460,221 entitled “Hybrid Stochastic Gradient Based Digital Controlled Oscillator Gain KDCO Estimation”. A portion (TUNE_TF) of the tuning word can be dithered and the resultant dithered portion can then be applied to a control input of switchable devices within the digitally controlled oscillator.


SUMMARY OF THE INVENTION

An embodiment of the present invention provides a phase locked loop that operates on clock signals derived from an RF clock signal generated by the phase locked loop. A frequency reference input provides a reference clock. A controllable oscillator generates the RF clock signal with a plurality of phases. A switch is coupled to receive the RF clock, and is operative to select one of the plurality of phases. A phase detection circuit is coupled to the switch and is operable to receive a selected phase and to provide digital phase error samples indicative of a time difference between the reference clock and the selected phase.





BRIEF DESCRIPTION OF THE DRAWINGS

Particular embodiments in accordance with the invention will now be described, by way of example only, and with reference to the accompanying drawings:



FIG. 1 is a block diagram of a single-chip radio with an all-digital local oscillator and transmitter and a discrete-time receiver;



FIG. 2 is a more detailed block diagram illustrating the interpolative all-digital phase locked loop (iADPLLiADPLL) based polar transmitter of FIG. 1 constructed in accordance with the present invention;



FIG. 3 is a block diagram of the digitally controlled oscillator (DCO) of FIG. 2 illustrating in more detail tracking bits with dynamic element matching (DEM) of the integer part and sigma-delta (SD) dithering of the fractional part;



FIG. 4 is a normalized abstraction layer of the DCO of FIG. 3;



FIG. 5 illustrates a general block diagram of a phase detection mechanism used in the iADPLL of FIG. 2;



FIG. 6 is a block diagram of the time-to-digital converter (TDC) of FIG. 2;



FIGS. 7A and 7B illustrate normalization and edge skipping operation of the TDC;



FIGS. 8A and 8B illustrate a conceptual view of data resampling using a sample rate converter;



FIG. 9A is a magnitude response plot and FIG. 9B is a phase response plot for the zero-order hold (ZOH) and first-order hold (FOH) resamplers of the iADPLL of FIG. 2;



FIG. 10A is a simple block diagram and FIG. 10B is a more detailed schematic of the phase error resampler of FIG. 2;



FIG. 10C illustrates Mu generation between CKVD32 and FREF in the resampler of FIG. 10A;



FIG. 11 is a block diagram of the loop filter for the iADPLL of FIG. 2;



FIG. 12A is a z-domain block diagram for the iADPLL with all iADPLL operations shown at fVD;



FIGS. 12B and 12C are plots of z-domain phase response with open-loop amplitude and phase transfer functions of the iADPLL of FIG. 12A with default loop settings, α=2−7, ρ=2−16, λ=[2−3 2−5 2−5 2−5] and integral gain operating on resampled PHE signal;



FIGS. 12D and 12E are plots of z-domain phase response with open-loop amplitude and phase transfer functions of the iADPLL of FIG. 12A with default loop settings, α=2−7, ρ=2−16, λ=[2−3 2−5 2−5 2−5], and integral gain operating on filtered PHE signal;



FIG. 12F is a block diagram of the iADPLL loop filter with alpha gear-shifting;



FIG. 12G is a plot showing the closed loop iADPLL transfer function magnitude response using both a raw PHE signal and a filtered PHEF signal for integral control;



FIG. 12H is a plot showing the closed loop iADPLL transfer function phase response using both a raw PHE signal and a filtered PHEF signal for integral control;



FIG. 13 is a block diagram of a hybrid stochastic gradient (HSG) KDCO normalizing factor calibration/compensation algorithm;



FIG. 14 is an illustration of coarse dithering applied to the slicer of the reference frequency generator;



FIG. 15 is a block diagram of circuit for providing both coarse and fine dithering to the slicer input;



FIG. 16 is a more detailed block diagram of the coarse dither circuit of FIG. 15;



FIG. 17 is a more detailed block diagram of the fine dither circuit of FIG. 15;



FIG. 18 is a more detailed block diagram of the sigma-delta generator of the fine dither circuit of FIG. 17;



FIG. 19 is a model of open loop FREF dithering of the iADPLL of FIG. 2;



FIG. 20 is a plot of phase error vs. dither resolution (fraction of inverter delay) for different inverter delay values;



FIG. 21 is a plot illustrating a comparison of phase error spectrums for FREF dithering in a low noise environment;



FIG. 22 is a plot illustrating effects of fine dithering on integer and half-integer channels with neighboring channels for high band (DCS, PCS);



FIG. 23A is a conceptual block diagram illustrating quadrature CKV rotation in the iADPLL of FIG. 2;



FIG. 23B is a plot illustrating quadrature CKV rotation when a frequency control work (FCW) has a value of integer five;



FIG. 24 is an alternative embodiment of an iADPLL;



FIG. 25 is another alternative embodiment of an iADPLL;



FIG. 26 is another alternative embodiment of an iADPLL; and



FIG. 27 is a block diagram of a digital system with an embodiment of an iADPLL within a digital radio transceiver.





DETAILED DESCRIPTION OF EMBODIMENTS OF THE INVENTION

The notation is used throughout this document is listed in Table 1.









TABLE 1







Glossary








Term
Definition





ADPLL
All Digital Phase Locked Loop


ASIC
Application Specific Integrated Circuit


CDMA
Code Division Multiple Access


CKR
Retimed Reference Clock


CKV
Variable Oscillator Clock, (~2 GHz) equals channel



frequency in HB and twice the channel frequency in LB


CMOS
Complementary Metal Oxide Semiconductor


DCO
Digitally Controlled Oscillator


DECT
Digital Enhanced Cordless Telecommunications


DPLL
Digital Phase Locked Loop


DRP
Digital RF Processor or Digital Radio Processor


EDGE
Enhanced Data rates for GSM Evolution


EDR
Enhanced Data Rate


FCW
Frequency Command Word is the ratio between variable



clock frequency and the reference frequency,










FCW
=


Variable





clock





frequency


Reference





Frequency











FDD
Frequency Division Duplex


FPGA
Field Programmable Gate Array


FREF
Frequency Reference


FSK
Frequency Shift keying


GFSK
Gaussian Frequency Shift Keying


GPRS
General Packet Radio Service


GSM
Global System for Mobile communications


HB
High Band


HDL
Hardware Description Language


IEEE
Institute of Electrical and Electronic Engineers


IIR
Infinite Impulse Response


KDCOINV
DCO Gain Normalization multiplier is the



reciprocal of the DCO gain step normalized



with the input sampling frequency,










KDCOINV
=


Datapath





Sampling





Frequency


Estimated





DCO





Gain











LB
Low Band


LDO
Low Drop-Out regulator


LFSR
linear feedback shift register


MOSCAP
Metal Oxide Semiconductor Capacitor


OTW
Oscillator Tuning Word is the normalization of the DCO



modulation frequency ratio by the DCO gain,










OTW
=


Modulation





Frequency


Estimated





DCO





Gain











PERINV
Period Normalization multiplier is the ratio between the



nominal DCO clock frequency and an



average inverter delay of the TDC,










PERINV
=


Variable





clock





time


-


period


TDC





Inverter





Delay











PHE
Phase Error


PLL
Phase Locked Loop


PVT
Process, Voltage, Temperature


RF
Radio Frequency


RTL
Register Transfer Logic


TDC
Time to Digital Converter


UMTS
Universal Mobile Telecommunications System


VCO
Voltage Controlled Oscillator


WCDMA
Wideband Code Division Multiple Access


WiMAX
Worldwide Interoperability for Microwave Access


WLAN
Wireless Local Area Network









An Interpolative All-Digital Phase Locked Loop (iADPLL) is an enhanced version of the All-digital PLL (ADPLL) frequency synthesizer, which is described in U.S. Pat. No. 7,145,399 entitled “Type-II all digital phase locked loop (PLL)” as well as U.S. Pat. No. 7,183,860 entitled “Gain Calibration of a Digital Controlled Oscillator” each of which are incorporated by reference herein in their entirety. Also, U.S. Patent Publication No. 2006/0038710A1, published Feb. 23, 2006, Staszewski et al., entitled “Hybrid Polar/Cartesian Digital Modulator” and U.S. Pat. No. 6,809,598, to Staszewski et al., entitled “Hybrid of Predictive And Closed-Loop Phase-Domain Digital PLL Architecture,” each of which are incorporated herein by reference in their entirety.


This advanced all-digital frequency synthesizer possesses two-point digital, i.e., reference and direct, frequency modulation capability. iADPLL, like its predecessor, uses digital signal processing and efficient circuit design techniques. A 4 or 8 GHz digitally-controlled oscillator (DCO) is at the heart of iADPLL, where fine frequency resolution is achieved through high-speed sigma-delta (ΣΔ) dithering. This setup allows for iADPLL loop control circuitry to be implemented in a fully digital manner using Infinite Impulse Response (IIR) filters and a digital Proportional-Integral (PI) controller. For the feedback of DCO phase in the time-domain, very fine time resolution (˜20 ps) is achieved by employing a Time-to-Digital Converter (TDC) comprising a chain of digital inverters.


The iADPLL architecture is distinct in that the loop filter and bulk of the logic operates on an RF derived variable clock domain. Due to the common RF based clock source for the major digital components, convenient digital handoff is achieved in the signal paths. The architecture uses two resamplers (one linear and the other a simple zero-order hold) for the exchange of data from and to the reference clock domain to the variable RF clock domain. In the variable clock domain, availability of several edge-aligned clock domains also facilitate physical register transfer logic (RTL) hardware sharing by time slicing. iADPLL can be embodied by more than one arrangement of the resampler and the other digital blocks. While several embodiments will be described herein, these descriptions are not intended to be limiting.


iADPLL architecture enhancements included in the various embodiments that will be described in more detail below aim at improving the spectral and the phase-domain performance of the original ADPLL architecture. In particular, the architectural changes use the DRP silicon area more efficiently as well as tackle the following critical aspects of the ADPLL design:


1. A reduction in the magnitude of the reference spurs is achieved by executing the loop filter in the iADPLL on a channel frequency derived clock. This significantly reduces the circuitry directly clocked by the reference clock or its retimed versions. Further, all the circuits interfacing to the DCO are being clocked at frequencies which are both higher and asynchronous to the reference frequency. Note that this arrangement necessitates the use of resampling circuits to deal with the multi-rate operation of the iADPLL. Mathematical analysis ensures the stability of iADPLL with these resamplers.


2. By using a direct point injection rate of channel frequency divided by 16, any sampling replicas are pushed to frequencies greater than 100 MHz from the carrier, where they are greatly attenuated by the DCO phase noise and loop skirt. Essentially these spurs are below the noise floor and do not appear in the simulations.


3. The use of a higher modulation injection rate also results in decreased spectral contribution of an iADPLL based transmitter in the GSM receive bands as described above. Further use of a higher direct-point injection frequency (>100 MHz) theoretically reduces the possibility of multiple spurs in the RX band.


4. A KDCO inverse calibration algorithm based on the Hybrid Stochastic Gradient (HSG) has been incorporated into iADPLL. This new algorithm ensures better than 1% accuracy of KDCO inverse estimate after initial compensation before each GSM burst under PVT. This aspect is described in more detail in U.S. patent application Ser. No. 11/460,221 entitled “Hybrid Stochastic Gradient Based Digital Controlled Oscillator Gain KDCO Estimation”, filed Jul. 7, 2006 and incorporated by reference herein in its entirety.


5. Improved dithering mechanisms are used to counter the TDC quantization dead-band impediments in the iADPLL phase feedback path, especially in case of integer-N channels, i.e., when channel frequency is a multiple of reference frequency, as well as half-integer-N channels. Dithering has been improved to provide uniform phase-domain characteristics for all possible GSM channel frequencies.


6. A mechanism for quadrature phase rotation of the ADPLL feedback DCO variable frequency (CKV) at the input of the TDC. This forces the ADPLL to be operated in an offset feedback mode, where the frequency command word (FCW) for integer channel frequencies will not be integer any more (due to the introduced offset).


To aid in understanding the principles of the present invention, a description is provided in the context of a digital RF processor (DRP) transmitter and receiver that may be adapted to comply with a particular wireless communications standard such as GSM, Bluetooth, WCDMA, etc. It is appreciated, however, that the invention is not limited to use with any particular communication standard and may be used in control, optical, wired and wireless applications. Further, the invention is not limited to use with a specific modulation scheme but is applicable to any modulation scheme including both digital and analog modulation.


Note that throughout this document, the term communications device is defined as any apparatus or mechanism adapted to transmit, or transmit and receive data through a medium. The communications device may be adapted to communicate over any suitable medium such as RF, wireless, infrared, optical, wired, microwave, etc. In the case of wireless communications, the communications device may comprise an RF transmitter, RF receiver, RF transceiver or any combination thereof. The notation DRP is intended to denote either a Digital RF Processor or Digital Radio Processor. References to a Digital RF Processor infer a reference to a Digital Radio Processor and vice versa.


A block diagram illustrating a single chip radio incorporating an interpolative all-digital local oscillator based polar transmitter and digitally-intensive receiver is shown in FIG. 1. For illustration purposes only, the transmitter, as shown, is adapted for the GSM/EDGE/WCDMA cellular standards. It is appreciated, however, that one skilled in the communication arts can adapt the transmitter illustrated herein to other modulations and communication standards as well without departing from the spirit and scope of the present invention.


The radio circuit, generally referenced 130, comprises a radio integrated circuit 136 coupled to a crystal 152, antenna front end module 176 connected to antenna 180 and battery management circuit 132. The radio chip 136 comprises a script processor 146, digital baseband (DBB) processor 144, memory 142 (e.g., static RAM), transmit (TX) block 148, receiver (RX) block 150, digitally controlled crystal oscillator (DCXO) 154, slicer 156, power management unit 138, RF built-in self test (BIST) 140. Battery 134 and battery management circuit 132 are connected to radio chip 136 for providing power. The TX block comprises high speed and low speed digital logic block 158 including ΣΔ modulators 160, 162, digitally controlled oscillator (DCO) 164, digitally controlled power amplifier (DPA) 174 or pre power amplifier (PPA), time-to-digital converter (TDC) circuit 170 and TDC quantization noise shaping block 166. The iADPLL and transmitter generate various radio frequency signals. The RX block comprises a low noise trans conductance amplifier 182, current sampler 184, discrete time processing block 186, analog to digital converter (ADC) 188 and digital logic block 190.


The radio also comprises TDC quantization noise shaping block 166 operative to reduce the quantization noise contribution of the TDC. It is noted that the TDC quantization noise shaping mechanism is especially applicable in an ADPLL circuit.


The interpolated all-digital phase-locked loop (iADPLL) based transmitter employs a polar architecture with all digital phase/frequency and amplitude modulation paths. The receiver employs a discrete-time architecture in which the RF signal is directly sampled and processed using analog and digital signal processing techniques.


A key component is the digitally controlled oscillator (DCO) 164, which avoids any analog tuning controls. A digitally-controlled crystal oscillator (DCXO) generates a high-quality base station-synchronized frequency reference such that the transmitted carrier frequencies and the received symbol rates are accurate to within 0.1 ppm. Fine frequency resolution is achieved through high-speed ΣΔ dithering of its varactors. Digital logic built around the DCO realizes an interpolated all-digital PLL (iADPLL) that is used as a local oscillator for both the transmitter and receiver. The polar transmitter architecture utilizes the wideband direct frequency modulation capability of the iADPLL and a digitally controlled power amplifier (DPA) 174 for the amplitude modulation. The DPA operates in near-class-E mode and uses an array of nMOS transistor switches to regulate the RF amplitude and acts as a digital-to-RF amplitude converter (DRAC). It is followed by a matching network and an external antenna front-end module 176, which comprises a power amplifier (PA), a transmit/receive switch for the common antenna 180 and RX surface acoustic wave (SAW) filters. Fine amplitude resolution is achieved through high-speed ΣΔ dithering of the DPA NMOS transistors.


The receiver 150 employs a discrete-time architecture in which the RF signal is directly sampled at the Nyquist rate of the RF carrier and processed using analog and digital signal processing techniques. The transceiver is integrated with a script processor 146, dedicated digital base band processor 144 (i.e. ARM family processor and DSP) and SRAM memory 142. The script processor handles various TX and RX calibration, compensation, sequencing and lower-rate data path tasks and encapsulates the transceiver complexity in order to present a much simpler software programming model.


The frequency reference (FREF) is generated on-chip by a 38.4 MHz (but could be 26.0 MHz or other in another embodiment) digitally controlled crystal oscillator (DCXO) 154 coupled to slicer 156. An integrated power management (PM) system is connected to an external battery management circuit 132 that conditions and stabilizes the supply voltage. The PM comprises a switched mode power supply (SMPS) as well as multiple low drop out (LDO) regulators that provide internal supply voltages and also isolate supply noise between circuits, especially protecting the DCO. The SMPS is used for efficient conversion of the battery voltage to a level that can be used by on-chip LDOs. The RF built-in self-test (RFBIST) 140 performs autonomous phase noise and modulation distortion testing, various loopback configurations for bit-error rate measurements and implements various DPA calibration and BIST procedures. The transceiver is integrated with the digital baseband, SRAM memory in a complete system-on-chip (SoC) solution. Almost all the clocks on this SoC are derived from and are synchronous to the RF oscillator clock. This helps to reduce susceptibility to the noise generated through clocking of the massive digital logic.


The transmitter comprises a polar architecture in which the amplitude and phase/frequency modulations are implemented in separate paths. Transmitted symbols generated in the digital baseband (DBB) processor are first pulse-shape filtered in the Cartesian coordinate system. The filtered in-phase (I) and quadrature (Q) samples are then converted through a COordinate Rotation DIgital Computer (CORDIC) algorithm into amplitude and phase samples of the polar coordinate system. The phase is then differentiated to obtain frequency deviation. The polar signals are subsequently conditioned through signal processing to sufficiently increase the sampling rate in order to reduce the quantization noise density and lessen the effects of the modulating spectrum replicas.



FIG. 2 is a more detailed block diagram of an iADPLL 200 used in the transceiver of FIG. 1 and constructed in accordance with the present invention. For illustration purposes only, the transmitter of the present embodiment is adapted for the GSM/EDGE cellular standard. It is appreciated, however, that one skilled in the communication arts can adapt the transmitter illustrated herein to other modulations and communication standards as well without departing from the spirit and scope of the present invention. For example, the transmitter illustrated in FIG. 1 can be extended for performing an arbitrary quadrature modulation scheme.


A description of the iADPLL, generally referenced 200, including the frequency/phase modulation path is provided herein below. The core of the iADPLL is a digitally controlled oscillator (DCO) 228 adapted to generate the RF oscillator clock CKV. The oscillator core (not shown) operates at twice the 1.6-2.0 GHz high frequency band or four times the 0.8-1.0 GHz low frequency band. The output of the DCO is then divided for precise generation of RX quadrature signals, and for use as the transmitter's carrier frequency. For GSM/EDGE transceivers, a single DCO is shared between transmitter and receiver and is used for both the high frequency bands (HB) and the low frequency bands (LB). However, for modern 3G (WCDMA) or other duplex transmission systems, separate local oscillators might be needed to supply TX and RX carrier frequencies.


A digitally-controlled oscillator (DCO) lies at the heart of the interpolated all-digital PLL (iADPLL) frequency synthesizer. It deliberately avoids any analog tuning voltage controls and is realized as an ASIC cell with truly digital inputs and outputs. The DCO comprises tunable switchable varactor elements, cross-coupled pairs of NMOS transistors, and a biasing circuit. The DCO varactors may be realized as n-poly/n-well MOS capacitor (MOSCAP) devices that operate in the flat regions of their C-V curves. Current advanced CMOS process lithography allows creation of extremely small-size but well-controlled varactors. The switchable capacitance of the finest differential TB varactor is in tens of attofarads. This resolution, however, is still too coarse for wireless applications and requires high-speed ΣΔ dithering to enhance the time-averaged frequency resolution, which is described in the following sections. The output of the DCO is input to the RF high band power amplifier 234. It is also input to the RF low band power amplifier 232 after divide by two in divider 230.


In case of transmit modulation, the symbols, (for example Gsm EDGE, wcdma, etc), in the form of in-phase and quadrature data streams are received from the digital baseband (DBB) circuit, not shown in this figure. The GSM symbols are passed through a pulse-shaping filter (PSF) within processor 212 that converts it to phase modulation. This phase modulation is interpolated in transmit data (DTX) processing circuit 250 and then passed on to the iADPLL after differentiation at the CKVD16 clock rate using differentiator 252. CKV is the iADPLL RF output digital variable clock in case of high-bands (HB>1 GHz) or twice the RF output clock in case of low-band (LB, <1 GHz).


For the case of EDGE, WCDMA, etc. the complex modulation I/Q data streams are fed to a COordinate Rotation Digital Computer (CORDIC) within processor 212, which converts it from Cartesian to polar representation. The amplitude modulation signal is passed through sigma-delta amplitude (SAM) signal processing blocks 214 before they are passed onto the on-chip digital pre-power amplifier (DPA) 234, while the phase modulation output of the cordic is passed onto iADPLL after the necessary interpolation and signal processing, which performs the phase modulation of the DCO.


Under no modulation conditions, iADPLL digitally controls the DCO to produce a stable variable clock (CKV) in the targeted RF frequency band. In the feedback path, CKV is used for phase detection and reference retiming. The time to digital phase conversion in the feedback happens using a TDC inverter chain 242.


The channel and data frequency control words are in the frequency command word (FCW) format, which is defined as the fractional frequency division ratio N, with a fine frequency resolution limited only by the FCW word-length. For example, with 24 fractional FCW bits, the frequency granularity using a 38.4 MHz reference frequency is 38.4 MHz/224 ≡2.29 Hz. In this embodiment, the direct point frequency injection is at the CKVD16 (which is 1×HB/2×LB channel frequency divided by 16, i.e., CKVD16=fV/16) rate, so the possible DCO frequency resolution is in the range of 6˜7.5 Hz (computed as fV/16/224).


The frequency reference (FREF) clock contains the only reference timing information for the RF frequency synthesizer to which phase and frequency of the RF output are to be synchronized. The RF output frequency (fV) is related to the reference frequency fR according to the following formula.

fV=N·fR   (1)


where, N=fV/fR≡FCW.


Synchronous Phase-Domain Operation


The iADPLL operates in a digitally-synchronous fixed-point phase domain. The variable phase RV[i] is determined by counting the number of rising clock transitions of the DCO oscillator clock CKV using variable phase accumulator 236. The variable phase RV[i] is sampled via sampler 238 to yield sampled FREF variable phase RV[k], where k is the index of the FREF edge activity. The sampled FREF variable phase RV[k] is fixed-point concatenated with the normalized time-to-digital converter (TDC) 242 output ε[k]. The TDC measures and quantizes the time differences between the frequency reference FREF and the DCO clock edges. The sampled differentiated (via block 240) variable phase is subtracted from the frequency command word (FCW) by a synchronous arithmetic phase detector 218. The reference phase RR[k] (equivalent to PHR in FIG. 5) is conceptually obtained by accumulating FCW with every cycle of the retimed frequency reference (FREF) clock input.


The frequency error fE[k] samples are accumulated via the frequency error accumulator 220 to create the phase error φE[k] samples. The digital phase error φE[k] is filtered by a digital loop filter 222 and then normalized by the DCO gain normalization circuit 270 in order to correct the DCO phase/frequency in a negative feedback manner. The loop behavior due to its digital nature is independent of process, voltage and temperature variations. The FREF retiming quantization error ε[k] is determined by the time-to-digital converter (TDC) 242 and the DCO period normalization multiplier 244. The TDC is built as a simple array of cascaded inverter delay elements and flip-flops, which produces time conversion resolution finer than 25 ps in the design process.


It must be recognized that the two clock domains, FREF and DCO, are not entirely synchronous and it is difficult to physically compare the two digital phase values without having to face meta-stability problems. During the frequency acquisition, their edge relationship is not known and during the phase lock the edges will exhibit rotation if the fractional FCW is non-zero. Consequently, the digital-word phase comparison is performed in the same clock domain. The synchronous operation is achieved by over-sampling the FREF clock using a higher-rate DCO derived clock (typically CKVD8) in reference retiming circuit 246. The resulting retimed CKR clock is thus stripped of the FREF timing information and is used throughout the system. This ensures that the massive digital logic is clocked after the quiet interval of the phase error detection by the TDC.


The main advantage of representing the phase information in fixed-point digital numbers is that, after the conversion, it cannot be further corrupted by noise. Consequently, the phase detector could be simply realized as an arithmetic subtracter that performs an exact digital operation. Thus, having a single conversion place, where the continuously-valued clock transition edge delay is quantized within the TDC, the susceptibility to noise and quantization errors is minimized and well controlled. It should be emphasized that it is very advantageous to operate in the phase domain for several reasons. First, the phase detector used is not a conventional correlative multiplier generating reference spurs. DRP architecture uses an arithmetic subtracter 218, which does not introduce any spurs into the loop. Second, the dynamic range of the phase error could be made arbitrarily large simply by the increasing word-length of the phase/frequency accumulators. Conventional three-state phase/frequency detectors are typically limited to only ±2π of the compare rate. Third, the phase domain operation is more amenable to digital implementations, contrary to the conventional approach.


High-Speed Direct Frequency Modulation Capability


As shown in FIG. 2, the oscillating frequency deviation Δf is dynamically controlled by directly modulating the DCO frequency in a feed-forward manner. The iADPLL loop compensates by effectively removing the loop dynamics from the modulating transmit path (using the reference modulation injection). The remainder of the loop, including all error sources, operates under the normal closed-loop regime. This method is similar to the conventional two-point direct modulation scheme but because of the digital nature, it is exact and does not require any analog component matching, except for the DCO gain KDCO=Δf/ΔOTW calibration, which is achieved in using a robust hybrid stochastic-gradient algorithm implemented in digital domain, where OTW is the oscillator tuning word and is analogous to the voltage tuning of a VCO.


The fixed-point frequency modulating data FCW is over sampled in resampler 254 by the iADPLL DCO injection frequency fV/16 and normalized in multiplier 262 to the value of iADPLL DCO injection frequency fV/16. Using the direct injection of the normalized FCW directly at the DCO impacts the oscillating frequency. The PLL loop will try to correct this perceived frequency perturbation integrated over the update period of 1/fR, which is then interpolated to the iADPLL operational frequency of fV/32 in resampling interpolater 256. This corrective action is compensated by the other (compensating) reference feed that is integrated by the reference phase accumulator. If the estimated DCO gain is accurate, i.e., {circumflex over (K)}DCO≅KDCO, then the loop response to the modulation is flat from dc to fV/64 (or half of iADPLL operational frequency fV/32). The immediate and direct DCO frequency control, made possible by accurate prediction of the DCO transfer function, is combined with the phase compensation of the PLL loop response. The two factors constitute the hybrid of predictive/closed PLL loop modulation method.


Advantages of Using Higher Rate Direct-Point Modulation Injection


One of the key advantages of using a direct point injection rate (say of channel frequency divided by 16) is that the phase modulation can be presented to the DCO with a finer resolution. For example, the phase modulation in GSM has a BW of 200.00 kHz, while for a polar TX, in EDGE mode the phase modulation BW is approx. 2.0 MHz (LB) and 1.0 MHz (HB). The CKVD16 rate corresponds to an injection frequency range of 103-124 MHz, which is at least three times higher than an FREF of 38.4 MHz, and 4 times higher than an FREF of 26 MHz. This implies that the phase modulation data update using a CKVD16 rate will be 3-4 times finer than the FREF rate used in the previous generations of ADPLL.


Furthermore, the data injection into the DCO comprises an integer and fractional parts, described in more detail below. The injection rate creates an effective zero order hold (ZOH) at resampler 254. The ZOH operation does not provide a large attention to the sampling replicas, which is only 13 dB lower for 2nd harmonic and approx. 17 dB for 3rd harmonic. As CKVD16 frequency is much higher than FREF, these replicas are correspondingly at 3-4 times higher frequency for CKVD16 (>100 MHz) vs. FREF (26-38.4 MHz). The DCO phase noise beyond the flicker corner of 1-2 MHZ has a 20 dB/decade slope, which implies that the residual sampling replicas after ZOH 254 sync filtering will receive an additional attenuation of 12 dB using CKVD16 injection rate as compared to FREF. In short, use of CKVD16 for direct point phase modulation injection results in pushing any sampling replicas to frequencies greater than 100 MHz from the carrier, where they are greatly attenuated by the DCO phase noise and the spectral skirt of the loop filter. Essentially these signal processing spurs are below the noise floor and can not be seen in simulations or measurements.


Another important benefit of using CKVDx, where x=16 or 8 for direct point injection is that the quality of phase modulation injection becomes independent of the FREF frequency. The same iADPLL when used with different FREF's, say 26 , 38.4 or 52 MHz will exhibit the same direct point injection fidelity. However, note that there are other noise scaling terms that will be impacted by the FREF frequency change, but the iADPLL loop filters, modulation injection rates etc. will maintain their resolution across multiple possible reference frequencies.


It has been observed in previous versions of the ADPLL that the current spikes caused by clocking of bulk of the logic can be a source of spurious emissions. This is especially true for highly integrated transceivers targeted using DRP technology. For iADPLL, a significant part of the loop filter and DCO interface logic executes on the LO derived clock domain. Since most of these frequencies are chosen to be higher than FREF, any such spurious products will have a larger intra-spur distance than FREF. For example, using CKVD32, the spurs (if present) will be 52-62 MHz apart as compared to FREF frequencies. In retrospect, the current spikes due to the modulation injection rate into DCO have the highest impact, as the rush current to the boundary level-shifters might be supplied by the same LDO supply regulator, which powers DCO. The most critical among these spurs are the ones that appear in the corresponding GSM/EDGE RX band during transmission. The widest GSM RX band is 65 MHz, and using CKVD16 at the interface at most one spur may appear in the RX band due to these parasitic supply regulation issues. Therefore, the use of a higher direct-point injection frequency (>100 MHz) theoretically reduces the possibility of multiple spurs in the RX band.



FIG. 3 is a block diagram of the digitally controlled oscillator (DCO) of FIG. 2 illustrating tracking bits with dynamic element matching (DEM) of the integer part and sigma-delta (SD) dithering of the fractional part. The digitally-controlled quantized capacitance of the LC tank is split into four major varactor banks that are sequentially activated during frequency locking. Large 1.8-2.5 MHz steps are performed during a process-voltage-temperature (PVT) calibration modes using a MIM-capacitor PVT bank (PB). Smaller 250-450 kHz steps of the acquisition bank (AB) are used during a channel select. The finest tens of kHz steps of the tracking bank (TB) are used during the actual transmit and receive. The unit weighted tracking bank is further partitioned into 32 or more integer and 3 fractional varactors. The fractional varactors undergo high-speed ΣΔ dithering via a 1st and 2nd order digital ΣΔ dithering circuit 312.



FIG. 3 illustrates an oscillator tuning mechanism 227 that improves the DCO 228 frequency resolution beyond the basic 100 kHz frequency step of the TB varactors. The 15 fixed-bit oscillator tuning word (OTW) 302 is split into seven integer 304 and eight fractional bits 306. In other embodiments, based on the configuration of the DCO varactor banks, the integer bit can be scaled to be five or six, for example. The LSB of the integer part (TB varactor) corresponds to the minimum frequency step of the DCO. Within the DCO, the integer part is thermometer encoded to control the same-size DCO varactors. This guarantees monotonicity and helps to achieve adequate linearity. The switch matrix 308, together with the row and column select logic, operates as a binary-to-unit-weighted encoder in response to the integer part of the TB tuning word. To minimize the impact on DCO linearity due to TB varactor mismatches, dynamic element matching (DEM) mechanisms that perform rotation of the TB row varactors have been built into the design.


The fractional part, on the other hand, employs a time-averaged dithering that produces a high-rate integer stream whose time-averaged value equals the lower-rate fractional TB input. The spurs due to sigma-delta (ΣΔ) idle tones are randomized using an linear feedback shift register (LFSR) 310 dithering of the fractional word. Use of the second order ΣΔ 312 and the high speed modulator clock speed (as high as CKVD1, i.e., 1.8 GHz), makes the in-band oscillator phase noise degradation almost immeasurable. With eight fractional bits and an integer bit size of 100 kHz, the effective open-loop DCO resolution will be 100 kHz/28=390.625 Hz, which is sufficient for GSM applications.


The fractional path of the DCO tracking bits is entirely separated from the lower-rate integer part. It even has a dedicated DCO input just to avoid “contamination” of the rest of the tracking bits with frequent transitions. The ΣΔ modulator 312 is responsive to only the fractional part of the tracking tuning word. Under certain operational conditions, especially when the input of the ΣΔ modulator 312 is a small but constant fraction, the ΣΔ output may exhibit spurs due to the cyclic pattern being generated by the modulator. These spurious products are called idle tones, which can be avoided by enabling an LFSR 310 output which gets added to raw input to create a dithered input to the ΣΔ modulator. The ΣΔ output becomes free of these idle tones by using the described mechanism.


The actual merging of both parts is performed inside the oscillator through time-averaged capacitance summation at the LC tank. Thus the critical high-speed arithmetic operations are performed in the “analog domain” through the additions of capacitance inside the DCO.



FIG. 4 is a normalized abstraction layer of the DCO of FIG. 3. At a higher level of abstraction, the DCO oscillator 228 together with the equivalent DCO gain normalization multiplication, results in a gain of:









f
V

/
16



K
^


D





C





O



=




f
V

/
16


f
R


·


f
R



K
^


D





C





O








The above equation logically comprises the normalized DCO (nDCO) 400, as illustrated in FIG. 4 for the data modulating path. The DCO gain normalization circuitry 270 (referring to FIG. 2) conveniently decouples the phase and frequency information throughout the system from the process, voltage and temperature variations that normally affect the KDCO. The frequency information is normalized to the value of the CKVD16 (1×HB/2×LB CKV frequency divided by 16, fV/16) injection frequency from the DTX 250/252. The DCO interface is updated at the CKVD16 rate (which is also the direct point frequency modulation injection rate) via multiplier 258 at the adder 224. The decomposition of the








f
V

/
16



K
^


D





C





O







DCO normalization factor, shown above, accounts for the interpolation already done on the digital frequency modulation data.


For clarity, note that the iADPLL phase accumulation 220 happens at the CKR (or FREF) rate, the remainder of the loop filters operates at the CKVD32 rate, see FIG. 2. The phase-domain resampling 256 from the CKR to CKVD32 rate serves to translate the data from one clock domain to the other using embedded interpolative filtering. This does not change the phase accumulation rate (which is CKR), which requires the correct DCO normalization to be








f
R



K
^


D





C





O



.





However, in order to derive this reference signal normalization from the








f
V

/
16



K
^


D





C





O







factor, which is precisely estimated by adaptation, a scaled version of








f
V

/
16



K
^


D





C





O







is used for the scaling 226 of the iADPLL loop correction applied to the DCO. For a scaling 260 by 2, the missing factor of







f
R



f
V

/
32






is considered as a part of the PI controller tuning parameters (i.e., proportional gain, α and integral gain, ρ). This avoids the use of an additional multiplier in the architecture.


The digital input to the Normalized DCO (nDCO) is a fixed-point normalized tuning word (NTW), whose integer part LSB bit corresponds to CKVD16. The quantity KDCO should be contrasted with the process-temperature-voltage-independent oscillator gain KnDCO which is defined as the frequency deviation (in Hz units) of the DCO in response to the 1 LSB change of the integer part of the NTW input. If the DCO gain estimate is exact KnDCO=fV/(16×LSB), where fV/16 is the direct point modulation injection rate. If there is a KDCO estimation error, then










K
nDCO

=




f
V


16
×
LSB


·


K
DCO



K
^

DCO



=



f
V


16
×
LSB


·
r






(
2
)







The dimensionless ratio r=KDCO/{circumflex over (K)}DCO is a measure of the accuracy of the DCO gain estimate.



FIG. 5 illustrates a general block diagram of the phase detection mechanism that can be mathematically captured as:

ΦE[k]=RR[k]−RV[k]+ε[k]  (3)


The operation consists of the phase detector 218 itself (see also FIG. 2, which shows a mathematically equivalent version), which operates on the three phase sources: reference phase RR[k] from adder 216, variable phase RV[k] from sampler 238, and the fractional error correction ε[k] from DCO period normalizer 244. The actual variable phase RV[l] from accumulator 236 is clocked by the CKV clock of index I and it must be resampled 238 by the CKR clock of index k. After the PHV resampling, all the three phase sources are synchronous to the CKR clock which guarantees the resulting phase error φE[k] to be also synchronous. An extra output bit from the fractional phase error correction comprising TDC and PF is due to metastability avoidance and is explained in the following sections.


The measurement of variable phase (i.e., the phase of the DCO output, CKV) is carried out in two steps. The integer part (i.e., the integer number of the CKV clock cycles) is determined using a non-resetable CKV edge counter called Variable Phase Accumulator 236. The remainder is the fractional part of the phase, which is the estimation of the sub-CKV clock period estimation of phase between FREF and the nearest CKV edge. This step is carried out by the time-to-digital converter (TDC) 242 described in the next section.


Integer Variable Phase Accumulation


The integer part of the variable phase comprises a count of the complete CKV clock cycles. The variable phase accumulator 236 implements the DCO clock count incrementing with the rollover effect as described in the following equation.









R
V



(

i
·
Tv

)





R
V



[
i
]



=




I
=
0

i


1





The deep submicron CMOS process is fast enough to perform an 8-bit binary incrementer at 2 GHz clock in one cycle using a simple carry-ripple structure. Critical timing of this operation would comprise a chain of seven half-adders and an inverter. However, for a commercial application it was necessary to add an extra timing margin in order to guarantee robust operation with acceptable yield for all the process and environmental conditions, as well as anticipated clock distribution skew statistics. This extra margin was obtained by increasing the maximum operational speed through topological means. The carry-ripple binary incrementer was transformed into two separate smaller incrementers, not shown. The first high-speed incrementer operates on the two lower-order bits and triggers the higher order increment whenever its count reaches “11”. The second incrementer operates on the same CKV clock, but the 6-bit increment operation is allowed now to take four clock cycles. The long critical path of the 8-bit carry-ripple incrementer has thus been split into smaller parts allowing for the necessary timing margin.


Time-To-Digital Converter (TDC)



FIG. 6 is a block diagram of the time-to-digital converter (TDC) 242 of FIG. 2. The TDC estimates the fractional part of the variable phase between the reference frequency edge and the next nearest edge of the DCO CKV clock. The TDC operates by passing the DCO clock through a chain of inverters 602. The delayed clock vector is then sampled by the FREF clock using an array of registers 604 whose outputs form a pseudo-thermometer code. The decoded binary TDC output is normalized 244 by the DCO clock period TV before feeding it to the loop. The combination of the arithmetic phase detector and the TDC is considered to be a replacement of the conventional phase/frequency detector. The number of TDC taps, L=56, has been determined as the count of inverters needed to cover the full DCO period under the strong process corner (min tinv=15 ps) plus some margin. Other embodiments may possess a different TDC topology and have fewer or more taps, depending on process parameters.


TDC resolution is a single inverter delay, Δtinv, which in this deep-submicron CMOS process is considered the most stable logic—level regenerative delay and in the 65 nm process node is approximately 20-25 ps. This results in a quality phase detection mechanism, as evidenced by the close-in and rms phase noise measurement results of the DRP architecture. While other TDC architectures can achieve the TDC resolution that is better than one inverter delay, they are quite complex and analog intensive. They simply appear not needed for GSM applications in this deep-submicron CMOS process when Δtinv=25 ps can be easily achieved


The phase quantization resolution of the variable phase accumulator 236, as described in the preceding section, is limited to ±½ of the DCO CKV clock cycle, TV. For wireless applications, a finer phase resolution is required. This is achieved using time-to-digital converter (TDC). The TDC measures the fractional delay difference ε between the reference clock and the next rising edge of the DCO clock. Using the TDC, the integer clock-domain quantization error is corrected by means of the fractional error correction term computed by the TDC, where the time between the rising and falling edges of CKV and FREF follow the relationship below











T
V

2

=

{





Δ






t
r


-

Δ






t
f







Δ






t
r




Δ






t
f









Δ






t
f


-

Δ






t
r





otherwise








(
4
)







Therefore, the total variable phase in the feedback path is RV[k]−ε[k]. As the iADPLL phase comparator is implemented in the frequency domain, which is mathematically equivalent to the pure phase domain operation, the variable phase is differentiated 240 to estimate the variable frequency error, i.e.,

feV[k]=(RV[k]−ε[k])−(RV[k−1]−ε[k−1])  (5)



FIG. 7A and 7B illustrate normalization and edge skipping operation of the TDC. The dimensionless TDC normalizing factor output from inversion circuit 702 is represented as follows:










N
TDC

=


2

W
F






T
_

V

/
Δ







t
res







(
6
)







NTDC is a fixed-point representation of the inverse of the expected DCO period TV in units of inverter delays Δtres=tinv. It is obtained through long-term averaging 704 of

TV=2|Δtr−Δtf|

(in inverter units) followed by inversion 702.


The averaging time constant could be as slow as the expected drift of the inverter delay, possibly due to temperature and supply voltage variations. The instantaneous value of the clock period TV=2|Δtr−Δtf| is an integer but averaging results in addition of significant fractional bits to the integer estimate with longer operations.











T
_

V

=


1

N
avg







k
=
1


N
avg





T
V



[
k
]








(
7
)







Note that each doubling of the accumulation length Navg would add one bit to the TV resolution.


It was experimentally confirmed that accumulating 128 clock cycles would produce accuracy within 1 ps of the inverter delay. The length of the operation is chosen to be a power of 2 since the division by the number of samples Navg could now be replaced with a simple shift-right operation. Of course, other embodiments using a different technology or process point may use a different number of accumulations.


Resampling Operation in the iADPLL



FIGS. 8A and 8B illustrate a conceptual view of data resampling using a sample rate converter 802. This section describes the resampling operations done in the iADPLL on the reference path and the phase error feeding into the loop filters. The sample rate conversion (SRC) or resampling is the mechanism with which data can be handed-off between digital systems (or circuits) which can have independent asynchronous clocks, such as example system-1804 and system-2806. The choice of the resampling mechanism determines the fidelity of the signal as it is passed from one clock domain to the other. The resampling operation can be as simple as a having a set of boundary registers on both the clock-domains, known as zero-order hold, ZOH; may contain linear polynomial interpolation, known as first-order hold; or more advanced polynomial filtering schemes such as Gardner, Lagrange, or Cubic Spline interpolations to name a few. For a resampler, the maximum difference in sample times indicated by 808 and the rate of change in the sampled signal determine how faithfully the output signal dOUT[k] tracks the input signal dIN[n].


In case of iADPLL, the two clock domains around the resampling domains are either the reference clock (FREF) or the retimed reference clock (CKR) and the CKV-derived clock. Although these two clocks are in general asynchronous to each other, the CKV being the output of iADPLL has a known relationship with FREF, which is tracked by the variable phase in the feedback path of the iADPLL. This information is used to efficiently resample the data between the above mentioned clock domains.



FIG. 9A is a magnitude response plot and FIG. 9B is a phase response plot for the zero-order hold (ZOH) 254 and first-order hold (FOH) 256 resamplers of the iADPLL of FIG. 2. The resampling operation in a closed loop control system such as the iADPLL, designed with a focus on cellular RF systems is challenging as sophisticated higher-order interpolation schemes can become very expensive to implement. On the other hand, a control system utilizing a resampler essentially becomes a multi-rate system, whose mathematical analysis is quite complicated. Therefore, for the iADPLL design, different resampling schemes were carefully analyzed and it was determined that due to very heavy filtering in the iADPLL loop filter, linear resampling (or FOH) is adequate for the iADPLL phase domain resampling. FIG. 9A shows the frequency response of both the ZOH 254 and the FOH 256 resamplers for a data injection rate of FREF (38.4 MHz in this case). The FOH has been implemented as a symmetric delayed interpolater with its co-efficient η=½. It can be seen that the FOH provides much more in-band filtering as compared to ZOH. The 3-dB cut-off frequency for ZOH is approx 17 MHz, while for FOH, it is approx 12 MHz. This reduces the possible impact of any aliasing that might be present in the iADPLL phase error due to the many noise sources that can potentially contaminate the DCO spectrum, such as the processor clock, current spikes in the power management system, etc. Moreover the minimum rejection above FREF is −13.26 dB for ZOH and −26.5 dB for FOH. Please, note that the IIR filters which are a part of the iADPLL loop filter are typically set in the range of 1-2 MHz, therefore for a spectral offset of 10 MHz, the phase error noise gets suppressed by 18.3 dB for ZOH and 19.5 dB for FOH, while the in-band aliased noise will be only 12 dB down for ZOH>24 dB lower for FOH.


Resampling of the Reference Feed Modulation


In this two point modulation scheme, the reference modulation input is the phase modulation compensation input, i.e., the reference modulation input cancels the phase modulation in the DCO output being feedback to the phase detector. Recall that the DCO was modulated using the direct point input. The reference modulation FCW is received from the DTX 250 at the LO-derived clock rate and needs to be resampled to the FREF domain at which the frequency/phase detector of the iADPLL operates, see FIG. 2. After careful analysis and to preserve area, it was deemed adequate that ZOH resampling is adequate on the FCW signal.


In various embodiments, two different ZOH implementations may be used. In option I, the resampling is done in the phase domain, the impact of differentiation (which is a high-pass filter to convert phase to frequency) is evident in the output spectrum. This can be problematic in case the phase modulation has tonal content, which may get amplified. Therefore option II, in which the ZOH resampling 254 takes place after differentiation 252 was adopted for the present embodiment. In this case the resampler output spectrum is flat. Note that use of ZOH at the reference modulation input introduces some aliasing. Note that the iADPLL variable phase accumulation also integrates the DCO variable phase over a period of reference frequency, which also produces aliasing in the feedback to the phase detector. The two aliased signals cancel each other substantially up to Fref/2, beyond which the residual aliasing is below the sensitivity floor for the iADPLL and gets attenuated by the phase error signal processing of the iADPLL.


Resampling of the iADPLL Phase Error



FIG. 10A is a simple block diagram and FIG. 10B is a more detailed schematic of phase error resampler 256 of FIG. 2. The frequency error detector and the phase error accumulator operate on fixed rate retimed reference clock i.e. CKR. The iADPLL loop filter 222 operates on a LO-derived clock domain, ie CKV. Therefore, the phase error needs to be resampled from this fixed rate clock domain into an RF-derived variable rate clock domain (CKV) before it is sent to the digital loop filter. The iADPLL output carries phase modulation in a polar TX. Phase resampler 256 helps to accomplish this functionality by using PLL variable phase information. Variable phase signals 1014 from the PLL carries the phase relationship information between the reference clock FREF and variable clock CKV. In Mu generator 1002, this information is used to calculate Mu values required for interpolating the phase error to a higher frequency variable rate clock, as illustrated in FIG. 10C.


The phase resampler performs 1st order linear interpolation on phase error signal 1006 in interpolater 1004. The phase error resampler is implemented using a pair of FIFO's 1016, 1018 to correctly interpolate phase error from fixed rate reference clock (CKR) 1010 to RF derived variable rate clock (CKVDx) 1012. FIFO depth is set appropriately to support multiple variable rate clock domains for LOOP_FILTER blocks of CKV/32, CKV/64, CKV/128 and also to support any reference clock frequency between 8 MHz and 38.4 MHz. Using FIFO's also eliminates the effects of CKR jitter on iADPLL performance.


MEM_DLO_RATE is set according to the desired division factor of the variable rate clock frequency for the LOOP_FILTER blocks. MEM_DLO_RFCW is a normalization which needs to be calculated and set for the resampler to correctly calculate Mu values, based on the measured variable phase. This is a fractional value which represents the ratio between the divide factor of loop filter clock rate from CKV and FCW. For example, for a division factor of 32, this value will be 32/FCW.


This linear interpolation resampler can be configured to operate as a zero order hold resampler by setting MEM_RESAMPLER_EN=0.


When the spectrum of phase error carries a tone, resampling using a ZOH resampler may result in spectral re-growth. This spectral re-growth will be somewhat mitigated by the IIR filters in the higher-order iADPLL loop filter. However, as mentioned above, the cumulative attenuation using ZOH is approximately 13 dB less than FOH. On the contrary, using FOH for phase error resampling there is no spectral re-growth for the same input signal. Moreover the phase error noise also gets more attenuation in this mode. Note that after resampling of the PHE, the sampling frequency (or rate) for the signal gets translated to an RF-derived clock domain. Since PHE signal was originally accumulated at the FREF (or CKR) rate, its magnitude needs to be appropriately adjusted for the change of sampling rate. An explicit scaling to this effect is avoided in the resampler to save area. This scaling is incorporated later into the proportional and integral gain coefficients of the PI controller in the iADPLL.



FIG. 11 is a block diagram of loop filter 222 for the iADPLL of FIG. 2. The Loop Filter's main input is the PHASE_ERR signal 1008 output from the resampler block 256. Its main function is to first filter the phase error signal and then feed the filtered phase error signal to the proportional-integral controller (equivalent of type-II operation in an analog PLL) of the DCO. Note that the applied proportional and integral gains include DCO gain normalization for each of the DCO capacitor banks. The loop filter outputs TUNE_PLL, which is the tuning word contribution due to the PLL. This is added (not shown on FIG. 11, but shown on FIG. 2 by adder 224) with TUNE_TX (the feed-forward modulation injection) in the DCO interface (DCOIF) block to produce the composite OTW word (as shown in FIG. 2).


The proportional part 1192 of the iADPLL loop has a gain alpha, which is used in both type-I (used for the DCO PVT and Acquisition bank tuning) and type-II (used for DCO tracking bank tuning) operational modes of iADPLL. iADPLL can have different programmable alpha coefficients that can be used based on the operational mode.


The four IIR filters 1196 exist to filter the noisy PHASE_ERR signal. All four IIR stages are similar with following first-order transfer function:

y[n]=2−λ*x[n]+(1−2−λ)*y[n−1]  (8)


where lambda is a programmable value and can take any value from 0 to 7.


The integral loop 1194 is intended to only work in tracking mode once the proportional loop is already quite close to the desired channel frequency. This is followed to minimize the overall settling time of iADPLL, otherwise, there is no hardware limitations forcing this to only be used during tracking. By adding an integrator to the loop, the loop will work to zero out the phase error. To prevent this long settling, the integral loop samples the current IIR filtered PHASE_ERR and uses the result to subtract 1198 from subsequent phase error samples. This results in a new offset error signal that can be forced to zero deviation by the integral control. The offset error value is then accumulated by the digital integrator and multiplied by the integral gain, rho, before it is added to the proportional correction.



FIG. 12A shows a z-domain, frequency mode model of the iADPLL, i.e., the input and the output of the iADPLL are in terms of frequency units. The model includes phase detector module 1218, loop filter module 1222, normalized DCO module 1270 and feedback TDC/PVL module 1242. The loop filter 1222 for iADPLL comprises four independently controlled IIR stages and a proportional-integral (PI) controller. Eq. (9) is a linearized z-domain model that includes the four cascaded single-stage IIR filters, each with an attenuation factor λi, where i=0 . . . 3 as well as the proportional and integral parts of the loop. The loop filter operates at the fVD≅fV/32 rate, where fV is the DCO output frequency.











H
loop



(
z
)


=



α
_

·




t
=
0

3




λ
i


1
+


(

1
-

λ
i


)



z

-
1







+


ρ
_

·


z

-
1



1
-

z

-
1










(
9
)








FIGS. 12B and 12C are open-loop amplitude and phase transfer functions of the iADPLL of FIG. 12A with default loop settings, α=2−7, ρ=2−16, λ=[2−3 2−5 2−5 2−5] and integral gain operating on resampled PHE signal.


Alternately, iADPLL also supports a mode in which the integral part of the PI controller also operates on filtered phase error computed by operating a set of cascaded IIR filters on the resampled PHE signal. The main advantage of this mode is that as the integral part operates on an error signal, which has undergone low pass filtering after resampling (identical to the proportional part) resulting in better group delay equalization of the two paths. Note that in this mode the phase margin of the iADPLL under similar conditions will be sacrificed. In this mode, the loop filter can be expressed as follows:











H
loop



(
z
)


=


(


α
_

+


ρ
_

·


z

-
1



1
-

z

-
1






)

·




t
=
0

3




λ
i


1
+


(

1
-

λ
i


)



z

-
1











(
10
)








FIGS. 12D and 12E are plots of z-domain phase response with open-loop amplitude and phase transfer functions of the iADPLL of FIG. 12A with default loop settings, α=2−7, ρ=2−16 , λ=[2−3 2−5 −5 2−5] and integral gain operating on filtered PHE signal.








ϕ
R
n


2

π


,



ϕ
V
n


2

π







and







ϕ
TDC
n


2

π








are the reference, variable (DCO) and the TDC error source contributions, respectively. For simplicity, the feedforward path comprises of the DCO represented only as the frequency scaling factor KDCO. For simplicity, the phase accumulation of the DCO has been modeled as a discrete integrator running at fVD rate in the feedback path. Note that the DCO integration can also be represented using the Bilinear/Tustin transformation. This is equivalent to running the DCO at the fV rate and then scaling the phase accumulation at the output to the fVD rate.











H
DCO



(
z
)


=



z

-
1



1
-

z

-
1




·

1

f
VD







(
11
)







The simplified feedforward transfer function of the loop is











H
FF



(
z
)


=



z

-
1



1
-

z

-
1




·

(



α
_

·




t
=
0

3




λ
i


1
+


(

1
-

λ
i


)



z

-
1







+


ρ
_

·


z

-
1



1
-

z

-
1






)

·


f
VD


f
R


·
r





(
12
)







where r is the dimensionless ratio representing the DCO gain estimation error.


The transfer function of the feedback components is











H
FB



(
z
)


=





z

-
1



1
-

z

-
1




·

1

f
VD


·
1

-

z

-
1



=


z

-
1



f
VD







(
13
)







The closed loop transfer function for the reference is lowpass with the gain multiplier N=FCW, i.e.,











H

cl
,
REF




(
z
)


=

N
·



H
FF



(
z
)



1
+



H
FF



(
z
)


·


H
FB



(
z
)










(
14
)







The closed loop transfer function for the TDC is lowpass. Neglecting the accumulation and the differentiation operations in the feedback path, we can write the closed-loop transfer function for TDC as











H

cl
,
TDC




(
z
)


=



H
FF



(
z
)



1
+


H
FF



(
z
)








(
15
)







The closed loop transfer function for the direct-point injection into the DCO is highpass in nature and is given by











H

cl
,
DCO




(
z
)


=

1

1
+



H
FF



(
z
)


·


H
FB



(
z
)









(
16
)








FIG. 12F is a block diagram of the iADPLL loop filter 1222 with alpha gear-shifting. The amplitude and phase transfer function differences using the iADPLL loop filters in the configurations as captured in Eq. (9) and (10), also see FIGS. 12B-12E, impacts the ADPLL closed loop performance as well. FIG. 12G is a plot showing the closed loop iADPLL transfer function magnitude response when using a raw PHE signal 1280 for integral control, and when using a filtered PHEF signal 1282 for integral control. FIG. 12H is a plot showing the closed loop iADPLL transfer function phase response using both a raw PHE signal and a filtered PHEF signal for integral control;


Essentially the use of filtered phase error for the integral control part of the iADPLL loop filter eliminates the amplitude and phase transfer function abruption (introduced due to resampling in raw PHE signal). However, note that physically the kink 1284 shown in the reference iADPLL transfer function is 80 dB below the PLL unity gain and as such does not produce any stability concerns for iADPLL loop bandwidth extending up to a few hundreds of kilohertz. Note, that under typical operating conditions, iADPLL loop bandwidth will be in the range of 20-60 kHz only.


Note that all the transfer functions are a function of the DCO frequency fV or its derivatives, such as CKVDx


The primary advantage of operating the iADPLL on the fVD frequency is the improved spectrum and better rejection of the reference spurs. However, this results in the iADPLL loop response to vary as a function of DCO frequency for fixed iADPLL tuning parameters. Although this iADPLL BW variation is quite small, it can be taken care of by optimal tuning of the iADPLL loop parameters as a function of the DCO output frequency. Note that for the entire GSM/EDGE tuning range this variation is less than 4% of PLL BW and less than 2% of the iADPLL phase margin.


Some operating points computed for GSM modulation using the model presented in equations (9)˜(16) are shown in Table 2.









TABLE 2







iADPLL Loop Parameters - Channel Frequency = 1800 MHz
























0 dB loop
Gain
Phase
50 kHz
400 kHz
3 MHz








BW
Margin
Margin
Loop
Loop
Loop


lambda 1
lambda 2
lambda 3
lambda 4
alpha
rho
[kHz]
[dB]
[deg]
Gain [dB]
Gain [dB]
Gain [dB]





















4
6
6
6
8
18
15.57
17.64
45.00
−10.22
−66.53
−98.31


2
4
4
6
7
16
30.89
20.28
45.23
−2.04
−36.02
−87.14


5
6
2
0
7
16
30.87
21.09
45.05
−2.06
−36.92
−93.12


3
5
5
5
7
16
31.15
17.38
44.81
−1.70
−38.27
−86.31


1
5
6
6
6
19
44.16
9.96
45.48
2.53
−41.59
−106.55


1
3
4
6
6
15
54.24
16.68
45.30
2.60
−28.40
−85.65


3
3
5
0
6
14
61.66
20.02
45.10
3.36
−21.80
−86.01


2
3
3
6
5
16
85.54
13.29
45.48
1.35
−21.81
−78.03


2
3
4
4
5
13
110.90
12.56
44.72
2.27
−13.40
−72.09


2
2
3
3
5
12
124.06
16.51
44.76
2.24
−9.71
−57.03










Proportional and Integral Gain Gear Shifting


Both the proportional and integral loops have the capability to change their bandwidths via gear shifting. When a gear shift happens, the gear shift value increments (or decrements) the corresponding gain. A gear shift mechanism is described in US Patent Publication N. 2007-0085622 entitled “CONTINUOUS REVERSIBLE GEAR SHIFTING MECHANISM” and which is incorporated by reference herein in its entirety.


For the dynamic proportional gain gear-shifting, the feedback register needs to sample a new value for the feedback path. Gear shifting can occur during any of the iADPLL operational modes; however, it is only sensible during tracking. Although the iADPLL loop bandwidth can be adjusted independently using either proportional (coarse change) or integral (finer change) gains, the possibility to simultaneously gear-shift both alpha and rho allows the iADPLL loop bandwidth to change while maintaining the phase and stability margins for the loop. Example operating points are shown in Table 3.









TABLE 3







Example iADPLL gear-shift combinations

























Gain
Phase
50 kHz
400 kHz
3 MHz








0 dB loop
Margin
Margin
Loop Gain
Loop Gain
Loop Gain


lambda 1
lambda 2
lambda 3
lambda 4
alpha
rho
BW [kHz]
[dB]
[deg]
[dB]
[dB]
[dB]





















2
3
4
4
5
13
110.90
12.56
44.72
2.27
−13.40
−72.09


2
3
4
4
6
14
61.58
18.58
46.52
3.18
−20.25
−78.11


2
3
4
6
5
19
82.72
11.43
44.94
1.18
−23.42
−82.08


0
3
4
6
6
15
54.17
17.02
45.58
2.55
−28.34
−88.52


2
2
4
6
6
15
54.14
17.28
45.85
2.51
−28.10
−86.26


2
2
4
6
7
16
30.57
23.30
46.72
−2.49
−34.25
−92.28


0
5
6
6
7
17
27.49
15.07
47.88
−2.50
−47.94
−92.90


5
5
6
6
8
18
15.49
18.23
46.06
−10.10
−63.64
−98.31










FIG. 13 is a block diagram of a hybrid stochastic gradient (HSG) KDCO normalizing factor calibration/compensation algorithm. A new KDCO adaptation algorithm is used in iADPLL 200 using the hybrid stochastic gradient (HSG) algorithm. The new algorithm is faster and more efficient as compared to the algorithm in ADPLL, with slightly higher hardware cost/complexity. The algorithm is primarily used for initial calibration and compensation using a pilot square wave training sequence. However, the algorithm because of its nature can also be tuned to work with arbitrary modulation waveform.


The new algorithm computes a stochastic sampled gradient between the filtered frequency modulation input (FCW) from FCW filter 1302 and the filtered zero-mean raw phase error (PHE) from PHE filter 1304 to compute the correction to the KDCO scaling value used in both reference and direct modulation paths of the iADPLL. The algorithm has been efficiently implemented using first order IIR filters on the incoming FCW and PHE inputs. The filters not only reduce the noise on these signals especially the raw phase error, but are also used to time align the filtered FCW and PHE signals. The correction gradient is computed to be the product of these IIR filters, which is then scaled appropriately using a shift operation 1306 to compute the correction to be applied to the KDCO normalizing factor(s). An optional momentum circuitry 1308 (or hysteresis filter) can be employed to further smoothen the computed value of the KDCO normalizing factor. This step is only required if the algorithm is used for background calibration with modulation.


Integer Channel Mitigation by Reference Clock Dithering


The purpose of reference dithering is to combat two independent issues which cause degradation of phase error of transmitted symbols on integer and to a lesser extent half-integer channels:


Fine FREF Dithering: Performance degradation due to the quantization error of the TDC. In short, the TDC quantization error causes the frequency/phase error between the desired frequency and the frequency produced by the DCO to periodically accumulate before it gets corrected by the loop, thus causing tonal content at the TX output.


Coarse FREF dithering is applied to desensitize iADPLL from the coupling of RF and other interference signals onto FREF.


Both fine and course reference dithering are controlled independently. Therefore they can be running concurrently. However, the system will be impacted by the overall dither amount. Dithering is described in more detail in U.S. patent application Ser. No. 11/853,182 filed Oct. 1, 2007 entitled “Adaptive spectral noise shaping to improve time to digital converter quantization resolution using dithering” and is incorporated herein by reference in its entirety.



FIG. 14 is an illustration of coarse dithering applied to inputs of the slicer of the reference frequency generator, such as DCXO 154 of FIG. 1. RF signal coupling on FREF is a parasitic phenomenon that potentially corrupts the reference signal. The dirty reference-clock can be further magnified by TDC quantization non-linearity in ADPLL loop dynamics. Coarse dithering is applied 1404, 1405 on the reference signal at slicer inputs as a DC offset to the sinusoidal reference signal. The DC offset moves the zero crossing of the reference signal. The slicer 1406 converts the sinusoidal reference wave into rectangular wave with its edges corresponding to the new zero-crossings. Since the coarse dither signal is much stronger than the RF interference signal the reference signal gets phase modulated by the dither signal with some dithering around it caused by the RF interference. If the coarse dither signal is chosen to have high frequency compared to ADPLL loop filter, the modulation caused by coarse dithering as well as RF coupling will be filtered out by the loop filter.


The TDC (Time to Digital Converter) converts the time delay between FREF clock and CKV clock to a digital word. TDC as a quantizer can only resolve delays specified by its resolution, in this case the inverter delay. When the CKV clock edges with respect to FREF clock edge is idle the quantizer is inadequate to represent the delay offset between the two. For integer channel operation non-modulated CKV clock aligns with FREF clock with some offset. Due to low frequency noise in the system the delay may drift within an inverter delay while TDC is unable to detect such changes. The delay range of the inverter delay is like a dead zone that ADPLL is unable to correct for. The undetected error can worsen the overall ADPLL performance. The effect of TDC quantization error is related to the frequency contents of this noise because it is filtered by the ADPLL's 4th order IIR filter. High frequency content of the TDC quantization noise is much more desirable over low frequency content of the quantization noise due to frequency selective nature of ADPLL. While CKV clock edge is trapped in TDC dead zone, the state of the ADPLL can nurture certain oscillation within this dead-zone. These oscillations are called limit-cycles that show up as undesirable, “mysterious” idle tones at RF output. The purpose of dithering is to rescue CKV clock from dead zone by causing random perturbation in FREF clock edges. Dithering is an intentional noise injected in the iADPLL loop, and thus it will introduce additional phase noise at iADPLL output. However this noise can be high pass shaped so that the overall iADPLL noise due to dithering is minimum. Higher order high pass shaping of the dither noise is not necessarily a good idea. The amount of high pass shaping required is dependent on the iADPLL loop filters. For example, the iADPLL has a 4th order IIR filter, therefore a reasonable choice of high pass shaping is 3rd order. The fine dithering signals are designed to be either noise shaped by ΣΔ or short (high frequency) periodic pattern.


Generation of Dithering Sequences



FIG. 15 is a block diagram of circuit for providing both course and fine dithering to the slicer input. This section describes the various types of reference clock dithering schemes implemented in the iADPLL. Control registers 1502 are connected to a system control processor, such as processor 146 of FIG. 1, and can be set with various parameters as needed to control dithering. Coarse dither circuit 1504 generates a square wave dither signal, as discussed above. Fine dither circuit 1506 generates various types of fine dithering signals, as will be described in more detail below.



FIG. 16 is a more detailed block diagram of the coarse dither circuit 1504 of FIG. 15. The purpose of coarse dithering is to tackle the specific problem of RF coupling onto digitally-controlled crystal oscillator (DCXO) and its buffer. The strong coarse dither signal applied on RF coupled reference signal desensitizes the FREF signal from RF coupling, as discussed above. Further, the frequency of coarse dither is high compared to ADPLL loop band-width. Therefore, the interference applied is filtered by the ADPLL loop filters. The coarse dither signal is a simple square wave signal with controllable frequency. It is implemented with a 1st order ΣΔ structure with a 5-bit accumulation stage 1604. The input DC value 1602 to the ΣΔ controls the frequency of the square wave. In conventional ΣΔ the output is the carry over or carry-out bit, however, in this case the dither bit is the MSB 1606 of the residue. The MSB of the residue has 50% duty cycle for power of 2 inputs. For example, if the input is 2(5−p) then the frequency of the dither is going to be FREF/2p with 50% duty cycle. For 5-bit implementation of ΣΔ the range of achievable frequency is FREF/2 to FREF/24. Small duty cycle variations can be obtained with non-power of two inputs. The carry over bit is kept available for future use.



FIG. 17 is a more detailed block diagram of the fine dither circuit 1506 of FIG. 15. Fine dithering is primarily applied to integer channels to improve rms phase-error due to TDC quantization noise. For integer channels the TDC quantization can cause wide range of rms phase error depending on the initial state of ADPLL. Dithering will ensure that TDC is kept busy enough so that CKV clock edges are not trapped in TDC dead zones. The dithering mechanism can be characterized into two categories: 1) Short periodic pattern, 2) Random sequence with noise shaping.


Theoretically, noise injected at TDC gets low-pass filtered by the ADPLL loop transfer function. It is desirable that the dither sequence has less low frequency contents. Short periodic sequence are good in that regard since there is no energy near DC, in fact, the only frequency contents are at the fundamental frequency and its harmonics and the periodicity being short puts these harmonic frequencies farther away. However if the dither amount is large enough due to high dither resolution step sizes the harmonic frequencies may show up at the iADPLL output spectrum. For this case a random noise shaped sequence may be more desirable. Based on test results, when the system has low noise random noise shaped sequence will perform better compared to a noisy system. Therefore both kinds of dithers are accommodated to combat different system scenarios. However for a reasonable noisy iADPLL system the short periodic pattern will perform better.


Using appropriate fine dither step (or resolution), i.e., dither step≧Tinv/N, where N is the number of dither elements, it is observed that the iADPLL rms phase error performance of integer channels is as good as non-integer channels plus the minor degradation caused by the dithering itself. In short, it is observed that using a short dither sequence, the peak-to-peak dither amplitude (dither resolution times dither sequence value) needs to be greater than 2 inverter delays, to improve the rms phase error of integer channels.


Fine dithering has a wide range of operating modes to accommodate different types of dithering. In FIG. 17, sequence generator 1702 uses a programmable LFSR to generate fixed alternating patterns. Sigma-delta generator 1704 is a 3rd order 5-bit sigma-delta that provides a signed integer output. The outputs of both generators are multiplexed and then scaled by multiplier 1706 that performs a 2GAIN shift. Weighting is performed and the output is then synchronized according to the reference fine dither mode. Table 4 lists various example fine dither modes that can be programmed via a set of control registers 1502.









TABLE 4





Fine dither modes


















0000
Use MEM_REFDTHRF_VAL



0001
Constant toggle



0010
3-bit LFSR



0101
1st order ΣΔ without LFSR-16



0110
2nd order ΣΔ without LFSR-16



0111
3rd order ΣΔ without LFSR-16



1101
1st order ΣΔ with LFSR-16



1110
2nd order ΣΔ with LFSR-16



1111
3rd order ΣΔ with LFSR-16











FIG. 18 is a more detailed block diagram of sigma-delta generator 1704 of the fine dither circuit of FIG. 17. Cascaded accumulators 1704a-c allow selection of first, second or third order ΣΔ, respectively.



FIG. 19 is a model of open loop FREF dithering of the iADPLL. The dithering mechanism is analyzed using iADPLL simulation in Matlab. Due to the nonlinear effect of TDC the overall performance of the dither mechanism is too complicated to analyze theoretically. It is important to evaluate the closed loop performance of the dither mechanisms with simulation. The effect of dithering in the closed loop system is shown in the model 1900 of iADPLL. Dither sequence generator 1506 generates a dither value that is added 1904 to the phase of the FREF clock supplied by the reference system fed to the variable phase generation mechanism. TDC 1942 then generates the FREF quantization error using this dithered signal and provides it to the loop which then operates as described above.


It has been shown through simulation that short periodic pattern fine dithering is an effective dithering mechanism for breaking TDC quantization limit cycles. In one embodiment, a triangle waveform is used for dithering. This waveform is inspired from the fact that non integer channel produces quantization noise that will have similar characteristics. Thus, with the help of dithering non-integer channel like quantization noise is induced on TDC. For ease of implementation triangular dithering is replaced by 3-bit LFSR dithering. In general short periodic patterns are chosen such that they have no low frequency contents to directly affect the rms phase error; however, their high frequency nature makes the TDC busy enough that it is not trapped in a dead zone.



FIG. 20 is a plot of RMS phase error vs. dither delay resolution (fraction of inverter delay) vs. dither range.(multiples of inverter delay). In Matlab, rms phase error is computed for varying dither delay resolutions and inverter delays. For good results, the dither resolution has to be large enough to span at least 1 inverter delay. However, as shown in FIG. 20, a larger amount of dithering beyond a single inverter delay does not impact rms phase error significantly. This may be helpful, but the increased amplitude of dither implies increased harmonic levels for the short dither pattern at the output RF spectrum. Unnecessarily large dither amplitude may cause spurious violations and co-existence issues with other surrounding radios. Furthermore, if the dither repetition period is large, the harmonic content at lower frequencies might not be adequately suppressed by the iADPLL loop filters. Therefore it is desirable to make the dominant frequency of the short pattern to be a higher frequency so that iADPLL loop filter attenuates it adequately below DCO noise contribution. In this simulation DCO phase noise is chosen to be −158 dBc/Hz.


3-bit LFSR dithering performance analysis is similar to that of triangular dithering.


A short periodic pattern can also be generated from a ΣΔ with proper use of input. When the input to the ΣΔ represents a fraction with small denominator the output of the ΣΔ will produce short periodic patterns. For example, all power of 2 inputs will produce periodic patterns. Simulation shows that, as long as the dither amount covers more than an inverter delay, the rms phase error is improved.


For random ΣΔ dithering, short periodic pattern has dominant tones multiples of its fundamental frequency. The assumption is that the loop filter will attenuate them below DCO noise level. However for low noise system or large dither amounts these tones can show up at the ADPLL RF output. This problem is mitigated by using a dither sequence that is random but high-pass noise shaped. The noise shaping is accomplished with a ΣΔ with random inputs. The random binary sequence is generated with a 16-bit LFSR. The order of the noise shaping is dictated by the order of low-pass filtering in iADPLL. This embodiment of the iADPLL has 4th order IIR filter, therefore, the highest choice of noise shaping is an order of 3.



FIG. 21 is a plot illustrating a comparison of phase error spectrums for FREF dithering in a low noise environment. The TDC quantization noise (or other nonlinear interference) has much worse effect on the overall ADPLL performance when ADPLL noise sources are small. Therefore it is important to investigate the effect of FREF dithering in a very low noise system.


Both ΣΔ generated or 3-bit LFSR generated short periodic patterns have similar performance in terms of rms phase error. However, the spectral characteristics of them are consistently different. This can be seen from the spectral plots of phase error at iADPLL output for both these dither sequences. The 2104 plots show the phase error at RF output due to ΣΔ dithering and the 2102 plots show the same for LFSR dithering. Many spectral plots are generated with different initial conditions as well as different seed for random noise sources. It is evident that the spectral contents of phase error for the two dither sequences are different.


For low-noise system, this difference in spectral characteristics is more evident. The same spectrums are plotted with all the noise sources turned off. ΣΔ generated short sequence has worse limit cycle effect compared to that of LFSR short sequence. The tones (limit cycles) generated near the natural frequency of the ADPLL has nothing to do with the frequency of the dithering, rather displays characteristics of the ADPLL system. It is to be noted that, although LFSR dithering is better in terms of limit cycles, ΣΔ dithering produces much lower tones at its fundamental frequency.



FIG. 22 is a plot illustrating effects of fine dithering on integer and half-integer channels with neighboring channels shown for high band (DCS, PCS). The plot shows simulations of integer and half-integer channels with the fine dithering solution turned on and off. The dither sequence shown is 3-bit LFSR with gain=1, which shows good performance in reducing the RMS phase error in peaking integer or half-integer channels. The DCXO slicer is modeled as a hard limiter with variable delay which is dependent on the number of fine dither bits that are turned on. The TDC delay is assumed to be 25 ps.


Integer Channel Mitigation by CKV Phase Rotation



FIG. 23A is a conceptual block diagram illustrating quadrature CKV rotation in the iADPLL of FIG. 2. At integer-N channels (i.e., channels at which the ratio between LO frequency and FREF is an integer), at near integer-N channels, and to a lesser extend at half-integer channels (i.e., FCW ratios which are a multiple of 0.5 but not exact integers), phase noise of ADPLL may degrade since TDC is in a dead-band condition. In this case, the timing relationship between the two clocks that the TDC receives is substantially constant and various TDC codes are not exercised. The well known large-signal approximation used in the field of the ADC design to model the quantization noise as a uniformly distributed noise is not longer valid there. This can give rise to idle tones that can get modulated. Another issue of operating at integer-N channels is a possibility of the LO clock coupling back into the TDC circuitry and affecting the FREF edges. The method described below can be used in addition to dithering of the reference clock edges as explained in the section above, which tries to relieve the TDC out of dead-band by injection of noise. This embodiment of iADPLL also provides another method to improve RMS phase error of the transmitter under similar conditions by performing of quadrature phase rotation of CKV.


This method effectively rotates the phase of the variable clock to PLL feedback path (phase/frequency detector) by 90 degrees every FREF edge to knock the TDC out of its dead-band. CKV phase can be either advanced or retarded by using the digital LO sequencing control signals. It should be noted that this operation takes place only for the feedback CKV clock to the iADPLL, the variable clock to the rest of the system including the amplitude path and the RF output remains to be the regular single phase CKV.


The simple trick that has been exploited in this scheme is that

FCW=fV/fR


is perceived by the user of the ADPLL, but

FCW′=fV′/fR=fV/fR±¼


is visible only to the extended TDC, where the sign depends on the direction of CKV rotation.


CKV phase rotation operates as follows. “Extended TDC” is defined as the core TDC with the decoding circuitry, the normalizing gain multiplier, and the CKV edge counters. It produces the fixed-point variable phase (RV[k]−ε[k]) having integer and fractional parts. After every FREF cycle (alternatively after a number of CKV clock cycles), the CKV clock supplied to the extended TDC gets delayed (alternatively advanced) by a quarter of the CKV clock cycle by quad switch 2302. Phase select controller 2304 controls the phase selection sequencing of quad switch 2302 through the select control signal SEL. Since the TDC operates on the rising FREF edge events, it virtually disregards any CKV activity not immediately preceding the FREF events. To properly maintain the appearance of the frequency multiplication ratio FCW to the external user, the internal value of FCW′ needs to be appropriately adjusted lower.


The phase rotation does not have to be in general increment/decrement or rotation, although rotation is a convenient scheme to implement. Other embodiments may use other schemes such as a random phase selection, for example. Likewise, the phase relationship does not have to be four phases with a quadrature or the 90 degree relationship. Other embodiments may use phases other than quadrature, such as three phases with 120° phases, for example.


In the illustrated example in FIG. 23A, the FCW value is an integer of 5. Every 5 CKV clock cycles or every FREF clock cycle, the CKV′ that is connected to the TDC gets delayed. The TDC perceives it as if the DCO period is (5 +¼)/5 of the actual DCO period. In general, the perceived DCO period is

TV′=TV(FCW+¼)/FCW=TV(1+1/(4*FCW))


The quad switch change could be done safely away from the rising FREF edges. The switching perturbations on the CKV′ are acceptable if the integer counting part of the extended TDC is powered down. The counting of CKV edges for the purpose of calculating the integer part of the variable phase is not needed once the iADPLL output frequency is locked. The CKVDx (i.e., CKV divided down by x, where x=8 for CKVD8) is based on CKV clock, which is not perturbed.


For proper operation with CKV phase rotation, the following sequence must be followed:

  • a) Disable use of integer part of variable phase. This also needs to be taken care of in calculation of script processor computations, such as FREQ_ERR (differentiated value equivalent of the phase error PHE).
  • b) Freeze the loop (PHE calculation) for at least 2 FREF cycles, and test for proper TDC operation
  • c) Enable CKV phase rotation synchronously to FREF/CKR but safely away from FREF rising edge; and at about the same time, as indicated at 2304
  • d) Update FCW to FCW′


Equation for new FCW (visible to ADPLL)

FCW′=FCW±¼  (17)


For example, fV=1742 MHz, fR=26 MHz, then Under normal operational conditions, FCW=67; but, FCW′=66.75 for +90 deg rotation, and FCW′=67.25 for −90 deg rotation


Variable Phase Integer Freeze Capability


LO phase rotation can cause ADPLL to get out of lock unless it supports a mechanism to handle the large known phase shifts caused by such rotation. This capability to handle large phase shifts in LO can be realized by variable phase integer freeze/disable capability.


Interpolative ADPLL operates by sampling reference clock phase using LO clock signal. This generates variable phase information which is represented digitally with a wide integer range and fine fractional resolution. Proposed solution to the above described issue is to create an ADPLL mode, in which iADPLL can handle large but known phase shifts.

    • a) This is done by first locking the ADPLL to the desired RF frequency and then engage in a mode of operation in which higher-magnitude phase information can be ignored and only lower-magnitude phase information is used to keep the ADPLL frequency locked.
    • b) Once the ADPLL is locked, then under normal operation, the higher integer magnitude phase is very predictable.
    • c) For a short duration, this higher integer phase can also be generated (or frozen to a known value) internally if required for the correct operation of the interpolative phase resampler of the iADPLL.


Practically, this implies that the integer part of the iADPLL variable phase information is ignored at every reference clock edge. This would not cause the iADPLL to react to larger magnitude phase shifts in LO. This is an enabler to the LO phase rotation capability which helps with TDC dead-band issues while dynamically keeping the ADPLL locked.



FIG. 24 is an alternative embodiment of an iADPLL 2400. In this embodiment, input to the iADPLL is differentiated phase (or in frequency domain) from digital transmit filter (DTX) 2450. The complete feedforward path of the iADPLL operates in the variable RF clock-domain, which also renders the reference path scaling to be more precise. In this embodiment, divider 2460 performs a divide by two and multiplier 862 performs a multiply by two. Zero order hold resampler 2454 is clocked by both CKVD32 and CKVD16. The linear resampler 2456 has been moved to the variable-phase (feedback) path.



FIG. 25 is another alternative embodiment of an iADPLL 2500. Input to iADPLL 2500 is in phase-domain from digital transmit filter (DTX) 2550. The reference point is computed using a resampler 2570 and a differentiator 2572 in the “reference” clock domain. The phase error (PHE) is computed at CKR rate and then linearly resampled by resampler 2556 to CKVD32.



FIG. 26 is another alternative embodiment of an iADPLL 2600. Input to iADPLL 2600 is in phase-domain from digital transmit filter (DTX) 2650. The compensating path operates on variable RF clock domain with resampler 2670 being clocked by CKVD32. The frequency domain operational regions of other iADPLL options have been eliminated in this option.


The implementation of the phase error block in this mode is practically more complex. Resampler 2656 buffers the phase error as with option 2400, but the error accumulator and phase error differentiator blocks have been removed and replaced with reference frequency accumulator 2680.



FIG. 27 is a block diagram of a digital system with an embodiment of an iADPLL within a digital radio transceiver, as described above. Digital system 1100 a representative cell phone that is used by a mobile user. Digital baseband (DBB) unit 1102 is a digital processing processor system that includes embedded memory and security features. In this embodiment, DBB 1102 is an open media access platform (OMAP™) available from Texas Instruments designed for multimedia applications. Some of the processors in the OMAP family contain a dual-core architecture consisting of both a general-purpose host ARM™ (advanced RISC (reduced instruction set processor) machine) processor and one or more DSP (digital signal processor). The digital signal processor featured is commonly one or another variant of the Texas Instruments TMS320 series of DSPs. The ARM architecture is a 32-bit RISC processor architecture that is widely used in a number of embedded designs.


Analog baseband (ABB) unit 1104 performs processing on audio data received from stereo audio codec (coder/decoder) 1109. Audio codec 1109 receives an audio stream from FM Radio tuner 1108 and sends an audio stream to stereo headset 1116 and/or stereo speakers 1118. In other embodiments, there may be other sources of an audio stream, such a compact disc (CD) player, a solid state memory module, etc. ABB 1104 receives a voice data stream from handset microphone 1113a and sends a voice data stream to handset mono speaker 1113b. ABB 1104 also receives a voice data stream from microphone 1113a and sends a voice data stream to mono headset 1114b. Previously, ABB and DBB were separate ICs but here are integrated into one IC. In most embodiments, ABB does not embed a programmable processor core, but performs processing based on configuration of audio paths, filters, gains, etc being setup by software running on the DBB. In an alternate embodiment, ABB processing is performed on the same OMAP processor that performs DBB processing. In another embodiment, a separate DSP or other type of processor performs ABB processing.


RF transceiver 1106 is a digital radio processor and includes a receiver for receiving a stream of coded data frames from a cellular base station via antenna 1107 and a transmitter for transmitting a stream of coded data frames to the cellular base station via antenna 1107. At the heart of transceiver 1106 lies a digitally controlled oscillator (DCO), which deliberately avoids any analog tuning controls. Fine frequency resolution is achieved through high-speed dithering of its varactors. Digital logic built around the DCO realizes an interpolative all-digital PLL (iADPLL) that is used as a local oscillator for both the transmitter and receiver and operates as described above. The polar transmitter architecture utilizes the wideband direct frequency modulation capability of the iADPLL and a digitally controlled power amplifier (DPA) for the power ramp and amplitude modulation. In this embodiment, a single transceiver supports both GSM and WCDMA operation but other embodiments may use multiple transceivers for different transmission standards. Other embodiments may have transceivers for a later developed transmission standard with appropriate configuration. RF transceiver 1106 is connected to DBB 1102 which provides processing of the frames of encoded data being received and transmitted by cell phone 1100.


The basic WCDMA DSP radio consists of control and data channels, rake energy correlations, path selection, rake decoding, and radio feedback. Interference estimation and path selection is performed by instructions stored in memory 1112 and executed by DBB 1102 in response to signals received by transceiver 1106. Programmable features of the iADPLL within transceiver 1106 are controlled by instructions executed by DBB 1102.


DBB unit 1102 may send or receive data to various devices connected to USB (universal serial bus) port 1126. DBB 1102 is connected to SIM (subscriber identity module) card 1110 and stores and retrieves information used for making calls via the cellular system. DBB 1102 is also connected to memory 1112 that augments the onboard memory and is used for various processing needs. DBB 1102 is connected to Bluetooth baseband unit 1130 for wireless connection to a microphone 1132a and headset 1132b for sending and receiving voice data.


DBB 1102 is also connected to display 1120 and sends information to it for interaction with a user of cell phone 1100 during a call process. Display 1120 may also display pictures received from the cellular network, from a local camera 1126, or from other sources such as USB 1126.


DBB 1102 may also send a video stream to display 1120 that is received from various sources such as the cellular network via RF transceiver 1106 or camera 1126. DBB 1102 may also send a video stream to an external video display unit via encoder 1122 over composite output terminal 1124. Encoder 1122 provides encoding according to PAL/SECAM/NTSC video standards.


As used herein, the terms “applied,” “connected,” and “connection” mean electrically connected, including where additional elements may be in the electrical connection path. “Associated” means a controlling relationship, such as a memory resource that is controlled by an associated port. The terms assert, assertion, de-assert, de-assertion, negate and negation are used to avoid confusion when dealing with a mixture of active high and active low signals. Assert and assertion are used to indicate that a signal is rendered active, or logically true. De-assert, de-assertion, negate, and negation are used to indicate that a signal is rendered inactive, or logically false.


While the invention has been described with reference to illustrative embodiments, this description is not intended to be construed in a limiting sense. Various other embodiments of the invention will be apparent to persons skilled in the art upon reference to this description. This invention applies to all scheduled communication systems which perform power control and channel sounding across multiple resource blocks. This invention applies in uplink and downlink. The embodiments of this invention apply for all modulation strategies, which include but are not limited to, OFDMA, CDMA, DFT-spread FDMA, SC-OFDMA, and others. Embodiments of this invention can be applied in most if not all emerging wireless standards, including EUTRA.


While a mobile user equipment device has been described, embodiments of the invention are not limited to mobile devices. Desk top equipment and other stationary equipment being served by a cellular network may also embody an iADPLL as described herein.


Although the invention finds particular application to Digital Signal Processors (DSPs), implemented, for example, in an Application Specific Integrated Circuit (ASIC), it also finds application to other forms of processors. An ASIC may contain one or more megacells which each include custom designed functional circuits combined with pre-designed functional circuits provided by a design library.


While the invention has been described with reference to illustrative embodiments, this description is not intended to be construed in a limiting sense. Various other embodiments of the invention will be apparent to persons skilled in the art upon reference to this description.


It is therefore contemplated that the appended claims will cover any such modifications of the embodiments as fall within the true scope and spirit of the invention.

Claims
  • 1. A digital system comprising a phase locked loop (PLL), wherein said PLL comprises: a frequency reference input for receiving a reference clock;a controllable oscillator for generating a radio frequency (RF) clock with a plurality of phases;a switch coupled to receive said RF clock, said switch operative to select one of said plurality of phases;a phase detection circuit coupled to said switch operable to receive a selected phase, said phase detection circuit providing digital phase error samples indicative of a time difference between said reference clock and said selected phase; anda phase select controller for changing input selection of said switch, wherein said changing circuitry executes every cycle of said reference clock and wherein said phase select controller is activated only on a subset of RF clock frequencies.
  • 2. The PLL of claim 1 further comprising a loop filter coupled between said phase detection circuit and said controllable oscillator, said loop filter operable to provide a filtered frequency tuning adjustment to said controllable oscillator.
  • 3. The PLL of claim 1 further comprising an oscillator gain normalization circuit coupled to said controllable oscillator for providing tuning adjustment to said controllable oscillator.
  • 4. The PLL of claim 1, wherein said phase detection circuit comprises a time-to-digital converter (TDC), said TDC providing digital samples indicative of a quantized value of said time difference between said reference clock and said selected phase.
  • 5. The PLL of claim 1, wherein said phase select controller selects a phase that is immediately advanced or immediately delayed from the previously selected phase.
  • 6. The PLL of claim 1, wherein said plurality of phases comprise quadrature phases.
  • 7. The PLL of claim 1, wherein said digital phase error samples are indicative of a phase difference between said reference clock and said RF clock.
  • 8. A digital system comprising a phase locked loop (PLL), wherein said PLL comprises: a frequency reference input for receiving a reference clock;a controllable oscillator for generating a radio frequency (RF) clock with a plurality of phases;a switch coupled to receive said RF clock, said switch operative to select one of said plurality of phases; anda phase detection circuit coupled to said switch operable to receive a selected phase, said phase detection circuit providing digital phase error samples indicative of a time difference between said reference clock and said selected phase, wherein said digital phase error samples are a difference between measured and expected timestamps of said RF clock, wherein said expected timestamps are calculated by accumulation of a frequency command word (FCW).
  • 9. A digital system comprising a phase locked loop (PLL), wherein said PLL comprises: a frequency reference input for receiving a reference clock;a controllable oscillator for generating a radio frequency (RF) clock with a plurality of phases;a switch coupled to receive said RF clock, said switch operative to select one of said plurality of phases;a phase detection circuit coupled to said switch operable to receive a selected phase, said phase detection circuit providing digital phase error samples indicative of a time difference between said reference clock and said selected phase; anda frequency command word (FCW) input for receiving an FCW, and circuitry for adjusting the received FCW.
  • 10. The PLL of claim 1, wherein a frequency of said RF clock is related to a frequency of said reference clock through a frequency multiplication ratio.
  • 11. A digital system comprising a phase locked loop (PLL), wherein said PLL comprises: a frequency reference input for receiving a reference clock;a controllable oscillator for generating a radio frequency (RF) clock with a plurality of phases, wherein a frequency of said RF clock is related to a frequency of said reference clock through a frequency multiplication ratio;a switch coupled to receive said RF clock, said switch operative to select one of said plurality of phases;a phase detection circuit coupled to said switch operable to receive a selected phase, said phase detection circuit providing digital phase error samples indicative of a time difference between said reference clock and said selected phase;a phase select controller for changing input selection of said switch, wherein said changing circuitry executes every cycle of said reference clock; andcircuitry for adjusting said frequency multiplication ratio.
  • 12. A digital transmitter comprising a phase locked loop (PLL), wherein said PLL comprises: a frequency reference input for receiving a reference clock;a controllable oscillator for generating a radio frequency (RF) clock with a plurality of phases;a switch coupled to receive said RF clock, said switch operative to select one of said plurality of phases;a phase detection circuit coupled to said switch operable to receive a selected phase, said phase detection circuit providing digital phase error samples indicative of a time difference between said reference clock and said selected phase; anda data modulating input for providing data deviation commands coupled to said PLL, the data deviation commands operable to adjust operation of said PLL.
  • 13. The digital system of claim 1 being a cellular telephone, wherein said PLL forms a local oscillator with a frequency modulation capability, said PLL further comprising a data modulating input for providing data deviation commands operable to adjust the frequency of said controllable oscillator.
  • 14. A cellular telephone comprising: a reference frequency generator for forming a reference clock; anda transmitter connected to receive said reference clock, comprising: a phase locked loop (PLL), wherein said PLL comprises:a controllable oscillator for generating a radio frequency (RF) clock with a plurality of phases;a switch coupled to receive said RF clock, said switch operative to select one of said plurality of phases;a phase detection circuit coupled to said switch operable to receive a selected phase and to compare said reference clock and said selected phase, said phase detection circuit providing digital phase error samples; anda data modulating input for providing data deviation commands operable to adjust the frequency of said controllable oscillator.
  • 15. A digital system comprising a phase locked loop (PLL), wherein said PLL comprises: a frequency reference input for receiving a reference clock;means for generating a radio frequency (RF) clock with a plurality of phases;means for selecting one of said plurality of phases;means for providing digital phase error samples by receiving a selected phase and comparing said reference clock and said selected phase; andmeans for changing input selection of said switch, wherein said means for changing executes every cycle of said reference clock and wherein said means for changing is activated only on a subset of RF clock frequencies.
  • 16. The PLL of claim 15, wherein said means for changing selects a phase that is immediately advanced or immediately delayed from the previously selected phase.
  • 17. A method for operating a phase locked loop (PLL), comprising: receiving a reference clock;generating a radio frequency (RF) clock with a plurality of phases;selecting one of said plurality of phases;providing digital phase error samples by receiving a selected phase and comparing said reference clock and said selected phase; andchanging to select a different one of said plurality of phases, wherein said changing occurs every cycle of said reference clock and wherein changing is activated only on a subset of RF clock frequencies.
  • 18. A cellular telephone comprising: a reference frequency generator for forming a reference clock; anda transmitter connected to receive said reference clock, comprising: a phase locked loop (PLL), wherein said PLL comprises:a controllable oscillator for generating a radio frequency (RF) clock with a plurality of phases;a switch coupled to receive said RF clock, said switch operative to select one of said plurality of phases; anda phase detection circuit coupled to said switch operable to receive a selected phase and to compare said reference clock and said selected phase, said phase detection circuit providing digital phase error samples; anda phase select controller for changing input selection of said switch, wherein said changing circuitry executes every cycle of said reference clock and wherein said phase select controller is activated only on a subset of RF clock frequencies.
  • 19. The cellular telephone of claim 18, wherein said phase select controller selects a phase that is immediately advanced or immediately delayed from the previously selected phase.
  • 20. A method for operating a phase locked loop (PLL), comprising: receiving a reference clock;generating a radio frequency (RF) clock with a plurality of phases;selecting one of said plurality of phases; andproviding digital phase error samples by receiving a selected phase and comparing said reference clock and said selected phase, wherein said digital phase error samples are a difference between measured and expected timestamps of said RF clock, wherein said expected timestamps are calculated by accumulation of a frequency command word (FCW).
  • 21. A method for operating a phase locked loop (PLL), comprising: receiving a reference clock;generating a radio frequency (RF) clock with a plurality of phases;selecting one of said plurality of phases;providing digital phase error samples by receiving a selected phase and comparing said reference clock and said selected phase; andreceiving a frequency command word (FCW): andadjusting the received FCW.
  • 22. A method for operating a phase locked loop (PLL), comprising: receiving a reference clock;generating a radio frequency (RF) clock with a plurality of phases;selecting one of said plurality of phases;providing digital phase error samples by receiving a selected phase and comparing said reference clock and said selected phase; andadjusting operation of said PLL in response to receiving a data deviation command.
  • 23. The method of claim 17, wherein changing to select comprises selecting a phase that is immediately advanced or immediately delayed from the previously selected phase.
CLAIM OF PRIORITY UNDER 35 U.S.C. 119(e)

The present application claims priority to and incorporates by reference U.S. Provisional Application No. 60/945,818 filed Jun. 22, 2007, entitled “Interpolative All-Digital Phase Locked Loop”.

US Referenced Citations (24)
Number Name Date Kind
4568888 Kimura et al. Feb 1986 A
5943369 Knutson et al. Aug 1999 A
6809598 Staszewski et al. Oct 2004 B1
6833875 Yang et al. Dec 2004 B1
7046098 Staszewski May 2006 B2
7061409 Jantti et al. Jun 2006 B1
7127022 Dieguez Oct 2006 B1
7145399 Staszewski et al. Dec 2006 B2
7183860 Staszewski et al. Feb 2007 B2
7205924 Vemulapalli et al. Apr 2007 B2
7496168 Leonowich et al. Feb 2009 B2
7570182 Sheba et al. Aug 2009 B2
20020080901 Ham, III Jun 2002 A1
20050243163 Ozasa et al. Nov 2005 A1
20060038710 Staszewski et al. Feb 2006 A1
20060078079 Lu Apr 2006 A1
20060165204 Shumarayev et al. Jul 2006 A1
20070075785 Kossel et al. Apr 2007 A1
20070085622 Wallberg et al. Apr 2007 A1
20070085779 Smith et al. Apr 2007 A1
20070164835 Co Jul 2007 A1
20080056426 Si et al. Mar 2008 A1
20080192877 Eliezer et al. Aug 2008 A1
20080317188 Staszewski et al. Dec 2008 A1
Related Publications (1)
Number Date Country
20080317188 A1 Dec 2008 US
Provisional Applications (1)
Number Date Country
60945818 Jun 2007 US