The present invention relates to the control of functions in a communications receiver that govern physical layer operations of equalization, synchronization (carrier frequency and phase, and timing recovery), and automatic gain control. These functions are controlled from a central signal processing and logic block to achieve robust, autonomous receiver operation.
In many digital communication systems, a source generates digital information, such as data, audio, or video, that is to be transmitted to multiple receivers. The digital information bits are divided into blocks that define a discrete alphabet of symbols. These symbols are used to modulate a radio frequency (RF) carrier's frequency, amplitude and/or phase. For example, a quadrature oscillator can be used to modulate the symbols onto the amplitude and phase of the RF carrier, and the signaling is referred to as Quadrature Amplitude Modulation (QAM). The time between adjacent symbols is referred to as the baud (or symbol) period, and inverse as baud (or symbol) rate.
The data-bearing RF carrier is amplified and transmitted through a propagation medium, for example, a cable, phone line, wireless terrestrial, satellite, underwater link, etc. As consumer appetite for information on demand continues to increase, so do requirements on data rates of present and emerging systems. Higher data rates increase the susceptibility of the data-bearing RF signal to environment-induced distortions such as fades, multipath (causing inter-symbol interference (ISI)), Doppler conditions, and additive noise processes. Communications receivers typically implement multiple stages of signal processing to individually address impairments witnessed by the RF signal en route to the receiver, and also introduced by the receiver and transmitter themselves, for example, due to implementation with finite-length filters and finite-precision arithmetic. Typical functions that govern physical layer operation include equalization, synchronization to carrier phase and frequency, synchronization to baud sampling phase and frequency, and automatic gain control.
Gain control processing adjusts the received signal to within the proper dynamic range expected by the receiver. Usually, gain control is itself done in stages, for example, at the RF level by the tuner, after RF to IF translation, and digitally within front-end and back-end signal processing. Most commercial UHF/VHF tuners include automatic gain control (AGC) circuits that tradeoff performance between adjacent and co-channel interference using average or peak power levels. AGC at the IF level usually adjusts the signal so that the full dynamic range of analog-to-digital conversion is utilized, without saturation. Sometimes, AGC utilizes pilot or reference tones or signals that are known, stored, or derived in the receiver to adjust IF AGC circuitry. In the digital domain, AGC circuitry may also use pilot or reference aids, or may operate blindly without the use of such overheard. For example, U.S. Ser. No. 60/341,931, entitled “Self-initializing decision feedback equalizer with automatic gain control,” by T. J. Endres et al., filed Dec. 17, 2001, which is incorporated herein by this reference, describes an automatic gain control (AGC) circuit that is nested with a decision feedback equalizer (DFE) structure, which operates on novel error signals that are blind (do not use training, reference, or pilot aids). Similarly, U.S. Ser. No. 10/246,084, entitled, “Adaptive expanded information capacity for communications systems,” by C. Long et al., filed Sep. 8, 2002, which is incorporated herein by reference, describes a similar feedback AGC circuit that is nested with an equalizer and uses decision-directed (DD) error adjustment.
Timing recovery refers to the process of synchronization to correct baud sampling frequency and phase. Sometimes, the oscillator clock that adjusts the analog-to-digital converter (ADC) at the input to the receiver is adjusted in frequency and phase. Alternatively, the sampled data in the receiver can be interpolated to achieve proper baud sampling. This second approach can reduce complexity by eliminating costly oscillator circuitry needed to control the ADC in the first approach. Timing recovery methods sometimes try to detect zero-crossings of the received signal in the time domain. As such, some methods apply filtering techniques to upper and lower band edges of the data spectrum. For example, see U.S. Pat. No. 5,872,815 by C. Strolle et al, entitled “Apparatus for generating timing signals for a digital receiver” which is incorporated by this reference. Alternatively, another approach uses a decision-directed (DD) technique, usually requiring feedback from a decision device, or slicer, or quantizer. The DD methods described by Mueller and Muller in “Timing recovery in digital synchronous data receivers,” IEEE Transactions on Communications, vol. COM-24, no. 5, May 1976, which is incorporated by this reference, are widely applied. Other blind algorithms, usually applied solely to equalization, have lately been successfully applied to timing recovery, too. For example, Guglielmi et al. in “Joint Clock Recovery and Baseband Combining for the Diversity Radio Channel,” IEEE Transactions on Communications, vol. 44. p. 114-117, January 1996, which is incorporated by this reference, jointly applies a Constant Modulus Algorithm (CMA), originally proposed by D. N. Godard in “Self-recovering equalization in two-dimensional data communication systems,” IEEE Transactions on Communications, vol. 28, no. 11, pp. 1867-1875, October 1980, which is also incorporated by this reference, for equalization, to joint optimization of equalization and timing recovery. This joint optimization is further analyzed by Chung et al. in “Timing Recovery Based on Dispersion Minimization,” Proceedings of the 2001 Conference on Information Sciences and Systems, March 2001 which is incorporated by reference.
Equalization in a digital communications receiver is analogous to an equalizer on a stereo system: the equalizer filters the distorted, received waveform (from a tape head or RF antenna) and tries to mitigate distortions and restore the signal properties of the original source. In a digital communications receiver, filter mismatches in transmitter and receiver, finite-precision implementation, and propagation channel effects induce inter-symbol interference (ISI), in which the receiver processes a signal that contains multiple, delayed and weighted copies of the transmitted signal. For example, reflections of the RF signal from a large building will induce ISI (or multipath), and reflections from an airplane will induce multipath distortion with a time-varying Doppler component. Since the exact channel characteristics are not known apriori at the receiver, the equalizer is usually implemented with adaptive methods. Adjustment of filter coefficients can be done with trained equalization methods, relying on the embedding of a pre-determined training sequence in the transmitted data. Usually, equalizer coefficient convergence relies on multiple transmissions of the training sequence, and the channel characteristics are also time varying, requiring periodic re-training. The Least Mean Squares (LMS) algorithm, which was proposed by Widrow, McCool, and Ball, in The Proceedings of the IEEE, vol. 63, no. 4, pp. 719-720, April 1975, which is incorporated by reference, minimizes a Mean Squared Error (MSE) cost function and is a stochastic gradient descent update rule utilizing the training sequence. Alternatively, blind methods do not rely on a reference signal, or derive a reference signal from the data itself, and are therefore desirable, since in the absence of a training signal, revenue-generating user data can instead be transmitted. A common blind equalization method replaces the reference signal in the LMS algorithm with the receiver's best guess at the data, and is referred to as Decision Directed LMS (DD-LMS), as proposed in a paper entitled “Techniques for adaptive equalization of digital communication systems,” by R. W. Lucky, in the Bell Systems Technical Journal, vol. 45, no. 2, pp. 255-286, February 1966 which is incorporated by reference. Unfortunately, DD-LMS needs a reasonably low percentage of incorrect decisions to prevent algorithm divergence, and is therefore impractical from a cold-start initialization. Godard's CMA is an attractive alternative that provides robust signal acquisition in harsh environments. These algorithms are sometimes applied to linear, finite-impulse response (FIR) filters, or decision feedback equalizer (DFE) structures, in which an FIR filter is embedded in a feedback loop (so that the overall impulse response is infinite) that processes hard decision samples from a decision device, slicer, or quantizer. U.S. Ser. No. 60/341,931 by T. J. Endres et al. entitled “Self-initializing decision feedback equalizer with automatic gain control,” which is incorporated by reference, uses novel, adaptive combining techniques of soft and hard decision samples in a DFE structure with CMA, LMS, and DD-LMS for robust, blind acquisition and re-acquisition.
Synchronization to carrier frequency and phase is usually done in stages, including, for example, RF to IF translation, IF to passband (or near baseband), and translation to precise baseband. Since the receiver in a coherent communications system must ultimately make hard decisions for symbol estimates, the data-bearing RF signal must ultimately be brought down to precise baseband. Synchronous or quasi-synchronous detectors can be used for intermediate translations. Translation to precise baseband can be done in a DD fashion, by nesting a phase-locked loop with the slicer and equalizer. For example, the instantaneous phase offset between slicer input and output is detected, filtered, and used to drive an oscillator in “Digital Communication” by E. A. Lee and D. G. Messerschmitt, Kluwer Academic Publishers, 1994, which is incorporated by reference. The PLL can be closed at various points within the equalizer. For example, Strolle et al describe alternative DFE architectures for reception of digital television signals that include a DD carrier loop, deriving its input from slicer input and output, but derotating the received signal to precise baseband at various points within the DFE structure. (See: Feasibility of reliable 8-VSB reception,” C. H. Strolle, S. N. Hulyalkar, T. J. Endres, Proceedings of the NAB Broadcast Engineering Conference, Las Vegas, Nev., pp. 483-488, Apr. 8-13, 2000, which is incorporated by reference.)
The present invention relates to the joint, adaptive, control of physical layer functions that govern equalization, baseband synchronization to precise carrier phase and frequency, baud sampling phase and frequency, and automatic gain control in the digital domain.
In accordance with certain aspects and embodiments of the present invention, physical layer functions governing equalization, translation to precise baseband, synchronization of sampling phase and frequency, and automatic gain control in the digital domain are jointly controlled and adapted with a central signal processing and logic block.
Other aspects, features, and advantages of various aspects and embodiments of the present invention will become more fully apparent from the following detailed description, the appended claims, and the accompanying drawings in which:
The receiver 150 receives the RF signal from propagation medium 147 via receiver antenna 149. Receiver antenna 149 is coupled to tuner 160. Tuner 160 is set to receive the RF signal in the desired frequency range, while rejecting signals in nearby or adjacent frequency ranges. Tuner 160 may provide automatic gain control at the RF frequency and also downconvert the received signal to an intermediate frequency (IF) before passing the signal to the analog IF Processing block, 162. Analog IF processing prepares the signal for analog-to-digital conversion, applying filtering and gain control, before passing the signal to the Front End Processing block 165. Front End Processing block 165 samples the signal with an analog-to-digital converter and may contain further automatic gain control and filtering, digital downconversion in frequency, and quadrature demodulation to split the signal into in-phase (I) and quadrature-phase (Q) samples. Front End Processing block 165 is coupled to Timing Recovery module 170 that determines a correct sampling phase. Timing Recovery module 170 may adjust the sampling phase by interpolating the data samples, or adjusting the phase and sampling frequency of the analog-to-digital converter in Front End Processing block 165. Timing Recovery module 170 is coupled to Equalizer 175, which is used to mitigate the distortions, such as inter-symbol interference and noise, that are introduced by the propagation medium 147, transmitter 100, Tuner 160, IF Processing block 162, receiver Front End Processing block 165, and receiver Timing Recovery module 170. Equalizer 175 is coupled to Carrier Recovery module 180, which detects residual offset in frequency and phase. The detected carrier offset in Carrier Recovery module 180 may be supplied back to the Equalizer 175 for translation of equalized samples to precise baseband, or used to adjust the downconversion process in Front End Processing block 165, or both. Also coupled to Equalizer 175 is Feedback AGC module 177, which works in conjunction with Equalizer 175 to detect and correct rapid amplitude variations in the received signal. The output of Equalizer 175 is coupled to Error Correction module 185, which detects and corrects bit errors in the recovered bit stream. The Error Correction module 185 is coupled to Decoder 190, which decodes the bit stream in accordance with the standard or protocol used in the Encoder 120 of Transmitter 100. The decoded bits from Decoder 190 represent the recovered information source, consisting of data, audio, or video, and are supplied to a user interface 195. Certain aspects and embodiments of the present invention are manifested in the joint, adaptive, blind control of Equalizer 175, Carrier Recovery module 180, Timing Recovery module 170, and Feedback AGC module 177.
Feedback AGC
A primary goal of the Feedback AGC is to provide quick, robust gain compensation of the signal that is input into the Slicer. This gain compensation may be required as the result of dynamic multipath conditions or flat fades. Though the equalizer itself provides gain compensation of the signal by adjusting the scale of the adaptive filters' coefficients, such compensation can be slow, since a plurality of equalizer coefficients (sometimes hundreds) must move to scale the data appropriately. Therefore, the Feedback AGC provides a single, strictly positive gain value, and is designed to provide gain compensation more quickly than the equalizer is able to do. Two different Feedback AGC architectures are disclosed.
Feedback AGC Preferred Embodiment
The preferred embodiment of the Feedback AGC in the present invention utilizes techniques described in USPTO application No. 60/341,931, entitled “Self-initializing decision feedback equalizer with automatic gain control,” by T. J. Endres et al., filed Dec. 17, 2001. These techniques rely on minimization of a specified cost function over the choice of gain value, and are implemented in an adaptive way with a stochastic gradient descent update rule.
The automatic gain control signal α(n) is a real, strictly positive scalar, that is calculated at each baud instance by
where ρα is a leakage factor used in practice to mitigate divergence due to finite-precision effects or quantization noise, and is chosen less than or equal to unity, but close to unity, J is the cost function to be minimized by choice of α(n), and μα is a real-valued, positive stepsize, chosen less than unity, and governs algorithm convergence rate, tracking capabilities, and stochastic jitter.
Different cost functions can be used, and these different cost functions will define different modes of operation that will be controlled by the Signal Processing, Logic and Control Block 260 in
The first operational mode of the Feedback AGC uses a Mean Squared Error (MSE)-like cost function, expressed as
J=E{(|α(n−1)·{tilde over (w)}(n−1)·e−θ(n)|q−|ŵ(n−1)|q)2}
where q is a positive integer and is set to one for the preferred embodiment. This cost function penalizes the squared difference in magnitudes between the slicer input and output, and is analogous to a mean squared error (MSE) cost function. The partial derivative calculation, denoted as {tilde over (ξ)}(n−1)≡∂J/∂α(n−1), assuming correct decisions and neglecting the expectation, for q=1, results in
{tilde over (ξ)}(n−1)=[α(n−1)·|{tilde over (w)}(n−1)|−|{circumflex over (w)}(n−1)|]·|{tilde over (w)}(n−1)|
where we have absorbed a factor of two into the stepsize μα, and used the fact that sign(α(n−1)) is always one since the automatic gain control signal is strictly positive by definition.
Since the automatic gain control signal, α(n), is applied to multiplier 225 in
An alternative embodiment of the present invention uses arbitrary positive integer q in the MSE-like cost function. See USPTO application No. 60/341,931 for a description of the update term for this case.
The second operational mode of the Feedback AGC uses a Constant Modulus (CM)-like cost function, expressed as
J=E{(|α(n−1)·{tilde over (w)}(n−1)·e−jθ(n)|q−γ)2}
This cost function has the advantage that it does not rely on hard decisions from the output of the Slicer. Letting q=2, taking the partial derivative and neglecting the expectation (as in the previous case for the MSE-like cost function) results in an error term for this cost function expressed as
{tilde over (ξ)}(n−1)=α(n−1)·|{tilde over (w)}(n−1)|2·(|α(n−1)·{tilde over (w)}(n−1)|2−γ)
The constant γ is calculated as γ=E{|s|4}/E{|s|2}, i.e., analogously to the Godard radius used in adaptation of the equalizer coefficients.
An alternative CM-like cost function uses q=1. The error term from the partial derivative is found as
{tilde over (ξ)}(n−1)=|{tilde over (w)}(n−1)|·(|α(n−1)·{tilde over (w)}(n−1)|−γ)
In this case with q=1, the Godard radius is calculated as γ=E{|s|2}/E{|s|}.
The third operational mode of the Feedback AGC uses a linear combination of MSE-like and CM-like cost functions. In this case, the error term derived from a MSE-like cost function is weighted by 1−λ(n), while the error term derived from a CM-like cost function is weighted by λ(n), and the results are summed to produce an error term used to update the Feedback AGC gain in the stochastic gradient descent rule. This combining is expressed as
{tilde over (ξ)}(n−1)=(1−λ(n))·{tilde over (ξ)}MSE(n−1)+λ(n)·{tilde over (ξ)}CM(n−1)
The combining weight λ(n) is time-varying and can be chosen among a variety of methods. For example, the combining weight can be set based on the signal-to-noise ratio, error rate, cluster variance, or number of symbols processed. In the preferred embodiment of the present invention, the combining weight is calculated in an adaptive method at each baud instance by the techniques discussed in USPTO application No. 60/341,931, by T. J. Endres et al., entitled “Self-initializing decision feedback equalizer with automatic gain control”. Alternatively, the combining weight can be determined from a look-up-table, which is addressed by the array index that is described in the future sections of the present document.
The fourth operational mode of the Feedback AGC halts adaptation of the Feedback AGC gain value. In this case, the error term is set to zero, {tilde over (ξ)}(n−1)=0, and leakage is set to one, πα=1. Hence, the new gain value is equal to the old gain value.
The instantaneous error term {tilde over (ξ)}(n−1), calculated in accordance with one of the four operational modes previously described, can be averaged before being applied to the stochastic gradient descent update rule to calculate the automatic gain control signal α(n), to induce memory in and reduce the variance of the error signal. For example, the error term used in the stochastic gradient update can be calculated as
ξ(n−1)=(1−ρage)·ξ(n−2)+ρage·{tilde over (ζ)}(n−1)
with automatic gain control signal updated by
α(n)=ραα(n−1)−μαξ(n−1)
where ρage is chosen greater than or equal to zero, but less than or equal to one, and is a leakage factor. Selecting ρage=1 represents no leakage and induces no memory in the error term. In this case, the error term relies purely on the unity-delayed samples and ξ(n−1)={tilde over (ξ)}(n−1).
To reduce undesired interaction between the Equalizer and Feedback AGC, a penalty term can be added to the cost function. This extra term penalizes gain values that are different from a target value, Γ. For example, the cost function including the penalty term is expressed as
J+β·(α(n−1)−Γ)2
where β is a small, non-negative weighting factor for the additional penalty term. The automatic gain control value is then updated according to
α(n)=ραα(n−1)−μαξ(n−1)−μα·β·(α(n−1)−Γ).
The product μα·β can be collapsed into a single value, to ease implementation.
The instantaneous error term {tilde over (ξ)}(n−1) is applied to a leaky integrator, comprised of multiplier 315, multiplier 325, and delay element 330, to induce averaging. The stepsize μα is applied in multiplier 335 to the averaged error term ξ(n−1). Adder 340 sums the update error term, penalty term, and previous AGC gain value to produce an updated AGC gain value, α(n). Leakage to avoid finite precision effects can be applied in multiplier 345 with leakage factor ρα to the previous AGC gain value, α(n−1), supplied by delay element 360. The previous AGC gain value, α(n−1), is also applied to summer 350 where a difference is formed with target threshold Γ. This result is multiplied by the weighting factor −μα·β in multiplier 355. The updated AGC gain value from adder 340, α(n), is also fed back to the Error Term Generator 310 to form the next instantaneous error term.
Feedback AGC Alternative Embodiment
An alternative embodiment of the Feedback AGC uses a bandpass control loop to adjust the AGC gain value. The bandpass characteristic is designed so that the rise (or attack) time of the Feedback AGC is very fast (about a few symbol periods) compared to the decay time of the Feedback AGC (hundreds or thousands of symbol periods). The gain value of the Feedback AGC decays to a prescribed threshold, nominally unity, as time goes to infinity. Hence, in a dynamic environment the overall gain of the received signal is quickly adjusted by the Feedback AGC and seamlessly transferred to the Equalizer.
The control loop derives the gain value from an error penalizing the difference of the ratio of the magnitudes of soft to hard decisions from a target value, for example, unity. The complex magnitude (square root of sum of real and imaginary components squared) of the equalizer output sample, w(n), is calculated in the CABS block 405. Look-up-table 410 stores inverse magnitudes of constellation points, and is addressed using the hard decision samples, ŵ(n), or indices associated with these hard decision samples, coming from the Slicer. Multiplier 415 forms the ratio of magnitudes of soft to hard decision samples, and this is termed the magnitude error and denoted as
Adder 420 takes the difference of the magnitude error from a target threshold, F, nominally unity, and forms the loop error, denoted as
Esym(n)=(Γ−εs(n))
The loop error, Esym(n), is first filtered by a low pass section that has transfer function
Adder 425 sums the loop error Esym(n) with a delayed and scaled version of the low pass output. Delay element 435 delays the low pass signal fL. Control signal reset is used to clear the contents of delay element 435 and output a zero-value from mux 440 instead of the output of delay 435. The output of mux 440 is multiplied with low pass scale αL in multiplier 430. Low pass scale αL is used to set the breakpoint of the low pass filter section. The output of multiplier 430 is summed with the input to the low pass section in adder 425. The low pass signal fL is scaled by low pass gain, βL, in multiplier 445, producing the output of the low pass section, fLβ, also used as the input to the high pass section.
The high pass section has transfer function
Signal fLβ is summed with a delayed and scaled version of signal fhint in adder 450, producing signal fhint. Signal fhint is delayed in delay element 460. Control signal reset is used to clear the contents of delay element 460 and output a zero-value from mux 465 instead of the output of delay element 460. The output of mux 465 is multiplied with high pass scale αH in multiplier 455. High pass scale αH is used to set the breakpoint of the high pass filter section. The output of mux 465 is also subtracted from signal fhint in adder 470, producing high pass signal fh. High pass signal fh is scaled by high pass gain, βH, in multiplier 475, producing the high pass output signal, fHβ.
The high pass output fHβ is added back to the target threshold, Γ, in adder 480, producing the calculated Feedback AGC gain value. Upon reset, however, mux 485 selects a unity-valued gain signal instead of the calculated AGC gain value. Using control signal FbagcFixedScaleEn, mux 490 sets the AGC gain value to programmable gain value FbagcFixedScale instead of the calculated gain value.
For a baud rate of approximately 613 KHz, nominal gain and scale values can be set to:
The Signal Processing, Logic and Control Block 260 in
When the Feedback AGC is operated in a decision directed (DD) mode, the magnitude error is found as
and the filtering is as described above for
When the Feedback AGC is operated in a Freeze mode, the system clocking of the circuit in
When the Feedback AGC is in Fixed Scale mode, all internal calculations are performed. However, mux 490 in
When the Feedback AGC is in Reset mode, the accumulators 435 and 460 are set to zero, mux's 440 and 465 output a zero-valued sample, and mux 485 outputs a unity-valued sample for the AGC gain value.
Carrier Loop
The goal of the final stage of carrier recovery is to translate the received signal to precise baseband in the synchronous demodulator. A precise baseband signal is needed to form hard decision samples, or symbol estimates, in the Slicer. A Carrier Recovery loop is used to track out any residual carrier frequency or phase offset, and phase noise introduced in the demodulation process. In
The Carrier Loop 250 is described in more detail in
Decision-directed phase and frequency recovery is discussed in the textbook by Lee and Messerchmitt, entitled Digital Communication, Appendix 17-B, Kluwer Academic Publishers, Norwell, Mass., Second Edition, 1994, which is incorporated by this reference.
In
Coarse Tune Hard Decision Operating Mode 1, Linear Combination of Soft and Hard:
In this preferred embodiment of the Coarse Tune Hard Decision Block 505, the alternative hard decision sample, ŵA, is formed by taking linear combinations of soft decision samples from the equalizer and hard decision samples from the Slicer,
ŵA=(1−λ)·w+λ·ŵ
where the real-valued combing weight λ is greater than or equal to zero but less than or equal to unity, and is selected adaptively. Methods to determine the combining weight are described in U.S. Ser. No. 60/341,931, entitled “Self-initializing decision feedback equalizer with automatic gain control,” by T. J. Endres et al., filed Dec. 17, 2001, which is incorporated by reference, and will also be described in the sequel.
Coarse Tune Hard Decision Operating Mode 2, Quadrant Mapping:
An alternative embodiment of the present invention uses Quadrant Mapping of the soft decision sample. The basic goal of generating an alternative hard decision sample using Quadrant Mapping is to estimate the quadrant in I/Q space in which the soft decision sample belongs, and scale the sample depending on the modulation format. This simple operation is expressed by selecting the alternative hard decision sample, ŵA, according to
ŵA=Δ·(sign(wI)+j·sign(wQ))
where Δ is a real-valued, positive scalar that can be selected according to the modulation format.
Coarse Tune Hard Decision Operating Mode 3, Reduced Constellation:
Certain alternative embodiments of the present invention quantize the soft decision sample not to the nearest element of the source alphabet (as the Slicer does), but to the nearest element of a subset of the source alphabet. This is referred to as a reduced constellation since only part of the source constellation is used to generate the alternative hard decision sample.
Coarse Tune Hard Decision Operating Mode 4, Reduced Precision:
Suppose that the baseband soft decision sample, w, is represented by M bits, while the Slicer's hard decision sample needs only N bits to be exactly represented, with M>N. An alternative hard decision sample is found by rounding and saturating the end points of the soft decision sample, to a bit-width P, where M>P>N. For example, suppose the soft decision sample is represented by M=16 bits, but only N=7 bits are needed to precisely represent the Slicer hard decision sample on the constellation grid. Then an alternative hard decision sample can be found by formatting the soft decision down to P=9 or P=10 bits.
Coarse Tune Hard Decision Operating Mode 5, Arbitrary Nonlinearity:
The Coarse Tune Hard Decision Block 505 can also generate alternative hard decision samples according to an arbitrary nonlinear curve, without memory. For example, a hyperbolic tangent curve can be approximated.
These operating modes are selected using input signal “control 3” to Coarse Tune Hard Decision Block 505. Input signal “control 3” is supplied from the Signal Processing, Logic, and Control Block (260 in
The alternative hard decision sample that is output from The Coarse Tune Hard Decision Block 505 in
The Phase Detector 515 derives an instantaneous estimate of the phase error between soft and hard decisions, and is described in Lee and Messerchmitt, Digital Communication, Appendix 17-B, Kluwer Academic Publishers, Norwell, Mass., Second Edition, 1994, which is incorporated by reference. Using a small angle approximation, the phase error estimate is calculated according to
where w=wI+j·wQ is the baseband soft decision sample, and ŵ=ŵI+j·ŵQ is the hard decision sample supplied from mux 510. In practice, the denominator is sometimes omitted.
For the preferred embodiment of the Coarse Tune Hard Decision Block 505, using Operating Mode 1, the hard decision sample used in the phase error calculation is the weighted sum of soft decision samples from the Equalizer and hard decision samples from the Slicer, using adaptive combining weight λ. When λ=0, signal quality is very good, and the instantaneous phase estimate is calculated with a hard decision sample that is equal to the Slicer hard decision sample, as intended in Lee and Messerchmitt, Digital Communication, Appendix 17-B, Kluwer Academic Publishers, Norwell, Mass., Second Edition, 1994. However, when λ=1 and signal quality is poor, the instantaneous phase estimate, φ, collapses to zero, which corresponds to no update to the phase and frequency produced by the PLL. In an acquisition state of the demodulator, when signal quality is poor, and hard decision samples are more unreliable, it is desirable to halt adjustment of frequency and phase that could slow or event prevent convergence of the loop. Hence, adaptively combining soft decisions from the Equalizer and hard decisions from the Slicer, the PLL is automatically transitioned among acquisition and tracking states, and halted (or slowed) when signal quality worsens and its inputs become unreliable.
The output of Phase Detector 515 is supplied to mux 520. Mux 520 selects from among Phase Detector 515 output or a zero-value. Mux 520 is controlled by signal “control 2”, that is supplied from the Signal Processing, Logic, and Control Block (260 in
The output of Mux 520 is supplied to Loop Filter 599. Multiplier 530 applies proportional control by multiplication with τ1, a programmable real-valued scalar, while multiplier 525 applies integral control by multiplication with τ2, a programmable real-valued scalar. The output of multiplier 525 is then integrated in adder 535 and delay element 540. The output of adder 535 is the integral control, and is summed in adder 545 with the proportional control from multiplier 530. Adder 545 supplies the Loop Filter 599 output to the NCO 598.
Numerically-controlled oscillator (NCO) 598 integrates the signal from adder 545 (Loop Filter 599 output) in adder 550 and delay element 555. The output of adder 550, representing instantaneous phase of the desired sinusoid, is supplied to a Sin/Cos Generator 560. Sin/Cos Generator 560 uses prior art techniques, such as LUT storage of sinusoid, in whole or part, or direct calculation of approximation by Taylor (or other) series, for example, to provide sine and cosine terms. These sine and cosine terms are used to de-rotate the Equalizer soft decision sample by e−j·θ(n), and also to re-rotate the Slicer hard decision sample by e+j·θ(n).
Equalization
One preferred embodiment of the Equalizer portion of the present invention is extracted from
In
Forward filter 220 is a finite impulse response (FIR) filter, computing its output according to the convolution sum
x(n)=fo(n)·r(n)+f1(n)·r(n−1)+f2(n)·r(n−2)+ . . . +fL
where r(n) is the sample sequence input to forward filter 220, x(n) is the output sample sequence of forward filter 220, fi are the forward filter coefficients (or parameters,) and Lf is the number of forward filter coefficients. Note that the forward filter coefficients are also shown with time index n to indicate that the forward filter 220 is adaptive.
The feedback filter 230 is not multi-channel, and is a FIR filter that calculates its output according to the convolution sum
y(n)=g0(n)·v(n)+g1(n)·v(n−1)+g2(n)·v(n−2)+ . . . +gL
where v(n) is the sample sequence input to feedback filter 230, y(n) is the output sample sequence of feedback filter 230, gi are the feedback filter coefficients (or parameters,) and Lg is the number of feedback filter coefficients. Note that the feedback filter coefficients are also shown with time index n to indicate that the feedback filter 230 is adaptive. Though the feedback filter 230 is a FIR filter, it is embedded in a feedback loop, so that the equalizer has an overall impulse response that is infinite.
Adaptation of the forward filter 220 coefficients and feedback filter 230 coefficients uses a stochastic gradient descent update rule:
fi(n+1)=fi(n)−μf·χ*(n)·e(n)
gi(n+1)=gi(n)−μg·φ*(n)·e(n)
where (•)* represents complex conjugation, and μf and μg are small, positive stepsizes governing algorithm convergence rate, tracking capabilities and stochastic jitter. Using simplified updates, the data used in the adaptation equations are set to χ(n)=r(n) and φ(n)=v(n). The error term, e(n), depends on the operating mode of the Equalizer and is controlled by the Signal Processing, Logic and Control Block 260 (in
Adder 215 combines the outputs of forward filter 220 and feedback filter 230, x(n) and y(n), respectively, to form sample sequence {tilde over (w)}(n). Sample sequence {tilde over (w)}(n) is referred to as passband soft decision samples. Multiplier 225 applies de-rotation from the Carrier Loop 250 by e−j·θ(n) and gain scaling from the Feedback AGC 240 by α(n). The output of multiplier 225 is denoted w(n) and referred to as baseband soft decision samples. These basband soft decision samples are input to Slicer 270. Slicer 270 is a nearest-element decision device that outputs a hard decision, ŵ(n), corresponding to the source alphabet member with closest Euclidean distance to its input sample. The hard decisions, ŵ(n), from slicer 270 are re-rotated to the passband by application of Carrier Loop 250 signal e+j·θ(n), and re-scaled by Feedback AGC 240 signal α−1(n) in multiplier 235.
The output of multiplier 235 is weighted by 1−λ(n) in multiplier 245, and the signal {tilde over (w)}(n) is weighted by λ(n) in multiplier 265. Signals 1−λ(n) and λ(n) are real-valued, adaptive combining weights provided by Signal Processing, Logic, and Control Block 260 (in
The operational mode of the Equalizer is determined by error term, e(n), whose selection is governed by the Signal Processing, Logic, and Control Block 260 (in
The Least Mean Squares (LMS) algorithm minimizes a Mean Squared Error (MSE) cost function that penalizes the squared difference between source sequence s(n) (in
elms(n)=(w(n)−s(n))·e+j·θ(n)
so that e(n) is set to elms(n). The baseband error term is re-rotated back to passband by multiplying the baseband error term with the sine and cosine terms from Carrier Loop 250. The training sequence is usually periodically resent, wasting valuable bandwidth.
In the absence of a training sequence, the output of Slicer 270 is used instead. Replacing s(n) with ŵ(n) in the LMS error term is called Decision Directed (DD) LMS and uses the error term
elms(n)=(w(n)−ŵ(n))·e+j·θ(n)
so that e(n) is set to elms(n).
When this operational mode is selected, training sequence data is used in the LMS calculation; else, hard decision samples are used in the LMS calculation. It is understood that modifications to these LMS error terms, using MSE techniques, to suit the specific application, are obvious to one skilled in the art.
The Constant Modulus Algorithm (CMA) minimizes the Constant Modulus (CM) cost function and calculates an error term without the use of a training sequence. Hence, it is classified as a blind algorithm. The CM criterion penalizes the dispersion of the magnitude squared of the equalizer output sample from a predetermined scalar. The preferred embodiment of the present invention calculates the CMA error term according to
ecma(n)={tilde over (w)}(n)·(|{tilde over (w)}(n)|2−γ)
where γ is a predetermined, real-valued scalar referred to as the Godard radius, or dispersion constant, and is usually determined from the source alphabet as γ=E{|s(n)|4}/E{|s(n)|2}. Notice that the passband soft decision sample {tilde over (w)}(n) is used in the CMA error term calculation, so that there is no multiplication required with sine and cosine terms. This sample is prior to Feedback AGC gain scaling in multiplier 225 (
An alternative embodiment of the present invention calculates the CMA error term according to
ecma(n)=w(n)·(|w(n)|2−γ)·e+j·θ(n)
In this case, e(n) is set to ecma(n).
These CMA error terms described are of order 2, as described by Godard. It is understood that modifications to these CMA error terms, using constant modulus techniques, to suit the specific application, are obvious to one skilled in the art. Furthermore, related Bussgang techniques, as summarized in chapter 2 of “Blind Deconvolution,” Prentice Hall, written by S. Bellini, edited by S. Haykin, 1994, which work is incorporated by this reference, are generalizations of constant modulus techniques, and can be substituted for the described CMA error terms by one skilled in the art.
Using the adaptive combing weights 1−λ(n) and λ(n) provided by Signal Processing, Logic, and Control Block 260 (in
e(n)=λ(n)·elms(n)+(1−λ(n))·ecma(n).
In this operational mode, the equalizer coefficients remain unchanged, so that the error term is
e(n)=0.
The error terms described above are for the Equalizer in
Timing Recovery
The purpose of the Timing Recovery Loop (170 in
To adjust the sampling epoch, most systems either adjust the sampling phase and frequency of the analog-to-digital converter (ADC), or interpolate the sampled data somewhere in the signal processing chain.
In
Phase Detector 720 receives the sampled data r(n), and computes a phase error responsive to signals from the Signal Processing, Logic and Control Block 260 in accordance with certain aspects and embodiments of the present invention. Loop Filter 730 uses proportional and integral control to filter high-frequency signal components and integrate the phase error for sampling frequency adjustment. The Reference Generator 740 uses the loop error signal from Loop Filter 730 to provide control signals for sampling frequency and phase adjustment in Sample Generator 710.
In accordance with certain aspects and embodiments of the present invention, the operating mode of the Timing Recovery Loop 170 is determined from the signals provided by the Signal Processing, Logic, and Control Block 260 by selecting the criterion by which the phase error signal is calculated in the Phase Detector 720.
Band edge timing recovery methods are prior art techniques basically based on the maximization of Signal to Noise Ratio (SNR) of sampler output. Band edge timing recovery schemes such as in U.S. Pat. No. 5,872,815 by C. Strolle et al, entitled “Apparatus for generating timing signals for a digital receiver,” which is incorporated by this reference, are implemented in pass band by exploiting the band edges of the data spectrum in order to obtain the timing phase maximizing the sampler output energy.
Phase Detector 720 computes the base band equivalent band edge timing phase error signal by differentiating the following base band output energy cost function
where eTL(n) is an error term (in this case derived from r(n)) and
is approximated in a variety or prior art methods, for example, using FIR filters in Lee and Messerchmitt, Digital Communication, Appendix 17-B, Kluwer Academic Publishers, Norwell, Mass., Second Edition, 1994 which is incorporated by reference.
With the cost function defined as
the band-edge phase error is
eTL(n)=r(n).
This operating mode is selected by the Signal Processing, Logic, and Control Block 260 through the four control bits DTL0-3. Bit DTL0 instructs Phase Detector 720 to calculate a phase error signal with its own internal data, rather than using data provided from Equalizer 175. The remaining control bits, DTL1-3, instruct Phase Detector 720 to calculate the specific error term selected by Signal Processing, Logic and Control Block 260.
Band edge timing recovery techniques work well for blind acquisition, provided there is sufficient signal power in the band edges. However, a single echo multipath ray can attenuate the signal power in the band edges, and prevent phase lock.
An alternative blind approach to timing recovery minimizes a CM cost function, JCM(r,τ), of the sampled data r(n) with respect to timing phase τ, as proposed by Guglielmi et al. in “Joint Clock Recovery and Baseband Combining for the Diversity Radio Channel,” IEEE Transactions on Communications, vol. 44. pp 114-117, January 1996, which is incorporated by reference, and further described by Chung et al. in “Timing Recovery Based on Dispersion Minimization,” Proceedings of the 2001 Conference on Information Sciences and Systems, March 2001, which is incorporated by this reference.
The phase error is found by differentiating the cost function with respect to timing phase, τ. Using calculus, this differentiation is usually decomposed into
where eTL(n) is an error term (in this case derived from r(n)) and
is approximated in a variety or prior art methods, for example, using FIR filters in Lee and Messerchmitt, Digital Communication, Appendix 17-B, Kluwer Academic Publishers, Norwell, Mass., Second Edition, 1994 which is incorporated by reference. (This same decomposition is used for each of the phase error signals that we will describe and are derived by minimizing a cost function.) The CMA error term is calculated as
eTL(n)=r(n)·(|r(n)|2−γ)
The CMA error term calculated by Phase Detector 720 is selected by Signal Processing, Logic, and Control Block 260 using control bits DTL0-3. Bit DTL0 instructs Phase Detector 720 to supply a phase error signal that is calculated with an error term using its own internal data, rather than using the error term provided by the equalizer. In
A Mean Squared Error (MSE) style cost function JMSE(r,τ) is minimized over the sampled data r(n) with respect to timing phase τ by Chung in “Blind Parameter Estimation for Data Acquisition in Digital Communications Systems”. Ph.D. Dissertation, Cornell University, Ithaca, N.Y., August 2002 which is incorporated by reference.
The phase error is found by differentiating the cost function with respect to timing phase, τ, and is of the form
eTL(n)=|{circumflex over (r)}(n)−r(n)|
where {circumflex over (r)}(n) is a locally generated or stored hard decision estimate of the current received sample, for example, from a training sequence, a local (to Timing Recovery Loop 170) slicer, a reduced constellation mapping, or other method.
The LMS error term calculated by Phase Detector 720 is selected by Signal Processing, Logic, and Control Block 260 using control bits DTL0-3. Bit DTL0 instructs Phase Detector 720 to supply a phase error signal that is calculated with an error term using its own internal data, rather than using the error term provided by the equalizer. In
Using the adaptive combining weight, λ(n), provided by the Signal Processing, Logic, and Control Block 260, a combination error term, comprised of CMA and LMS error terms calculated in Operating Modes 2 and 3, is calculated by Phase Detector 720 according to
eTL(n)=(1−λ(n))·eCMA(n)+λ(n)·eLMS(n)
In
In this operating mode, the phase error is set to zero. Control bit DTL0 selects a phase error in mux 835 that is supplied by multiplier 825, while control bits DTL1-3 select a zero-valued error term in mux 830.
In this operating mode, the error term calculated in Equalizer 175 is used in the phase error calculation. In
Joint, Adaptive Control of Equalization, Synchronization, and AGC
A block diagram of the Signal Processing, Logic, and Control Block 260 in accordance with certain aspects and embodiments of the present invention is shown in
To calculate candidate error terms, baseband and passband Equalizer 175 (
The combining weight is chosen at each baud instance by comparing the distance of the baseband soft decision, w(n), to its nearest element in the source constellation, and normalizing by the size of the decision region. This idea is illustrated in
The left-hand-side of
To add memory or induce averaging to the instantaneous combining weight {tilde over (λ)}(n), a leaky integrator is used, and the value of combining weight λ(n) is calculated as
λ(n)=(1−ρλ)·λ(n−1)+ρλ·{tilde over (λ)}(n)
where ρλ is the leakage term and is chosen less than or equal to one and greater than or equal to zero. Also, the combining weight λ(n) may be compared to two thresholds, TU and TL. If λ(n)>TU, then λ(n) is set to a register value, for example, unity; if λ(n)<TL, then λ(n) is set to another register value, for example, zero.
A circuit used to calculate the combining weight λ(n) in accordance with certain aspects of the present invention is shown in
Multiplier 1130, adder 1135, delay 1140, and multiplier 1145 comprise a leaky integrator, or averaging circuit. Multiplier 1130 scales the clipped value from block 1125 by the leakage amount ρλ. This product is summed with the output of multiplier 1145 in adder 1135. The result is delayed in delay element 1140, scaled by the leakage amount 1−ρλ in multiplier 1145 and then supplied to adder 1135.
The output of adder 1135 is first compared to upper threshold TU in comparator 1155. If this first comparison is satisfied, the combining weight is set to programmable upper value TU-value in assignment block 1165. If this first comparison is not satisfied, the output of adder 1135 is next compared to lower threshold TL. If this second comparison is satisfied, the combining weight is set to programmable lower value TL-value in assignment block 1170. If this second comparison is not satisfied, the output of adder 1135 remains unchanged in assignment block 1175.
An alternative embodiment of the present invention calculates a combining weight λ(n) by indexing a programmable array of stored combining weight values, indexed by array index provided by Array Index Update block 975.
Referring again to
A second control signal used in the Signal Processing, Logic, and Control Block 260 is generated based on the position of the hard decision sample in the source constellation, conditioned upon a prescribed template. For example, Slicer 270 supplies the symbol number (0 . . . M for M-ary QAM) prescribed to the hard decision sample in the constellation bit-to-symbol mapping (usually a type of grey code) to Region LUT 920. Region LUT 920 stores the desired template that is used to discern the second control signal. For example,
A third control signal used in the Signal Processing, Logic, and Control Block 260 is a training sequence indicator and is denoted by TR. This control signal is valid when the current sample is part of a training sequence sent by the transmitter, instead of bearing user data. It is apparent to one skilled in the art how to derive a training sequence indicator, for example using correlation techniques.
Control signals in_box and TR are logically compared in AND gate 950, and the result is supplied to an accumulator 970 that increments when the result is true. The accumulator 970 operates until a clear signal is provided by Array Index Update block 975, at which time accumulator 970 supplies its accumulated count, designated as JT, to Array Index Update block 975, and resets its accumulator value to zero.
Control signals in_box and in_region are logically compared in AND gate 955, and the result is supplied to AND gate 960, which logically compares this result with the complement of control signal TR. The result of AND gate 960 is supplied to accumulator 965 that increments when the result is true. The accumulator 965 operates until a clear signal is provided by Array Index Update block 975, at which time accumulator 965 supplies its accumulated count, designated as JD, to Array Index Update block 975, and resets its accumulator value to zero.
Array Index Update block 975 calculates array index ranging from 0 . . . 15 in the preferred embodiment of the present invention, according to the flow diagram in
The circuit in
Decision block 1320 checks to see that the current value of array index does not exceed the minimum value, 0. If >0 is satisfied, then decision block 1325 is entered; else, output block 1340 is entered. Decision block 1325 compares JD from accumulator 970 to the value of decrement array AD that is indexed by the current value of array index , denoted by AD(). If JD<AD() then array index is decremented in decrement block 1335 and output block 1340 is entered; else, output block 1340 is entered directly.
Output block 1340 sends “clear” signals to accumulators 965 and 970, resetting the count to zero in each, and loads another value of array index , which has either been incremented, decremented, or is unchanged.
Referring again to the Signal Processing, Logic, and Control Block 260 in
For the alternative embodiment of Feedback AGC 240, Array Block 980 stores a length-P array of freeze bits. The freeze bit selected by array index is provided to logical OR gate 990.
Moreover, Array Block 980 stores a length-P array of godard radius values, γ, one of which is used for CMA error term calculation. The elements stored in the array of godard radius values can be pre-calculated in accordance with the templates that are stored in Region LUT 920 corresponding to the current value of array index . For example, since arbitrary templates of constellation points can be designed for inclusion in Region LUT 920, the CMA error term should be calculated based on the effective source constellation it sees (a subset of the full QAM source constellation), which could be the values stored in the present template of Region LUT 920, or their complement. This choice depends on the design of the Adaptation Select LUT 915. The godard radius selected by array index is provided to Error Term and Combining Weight Calculator 905, and used in the CMA error term calculation.
Array Block 980 can also store arrays for other demodulator parameters, like stepsizes for adaptive filters 220 and 230. Each array is indexed by array index , and a new parameter is provided by Array Block 980 for demodulator operation.
The three control signals, in_box, in_region, and TR are used also as address bits to Adaptation Select LUT 915. Adaptation Select LUT 915 stores a programmable LUT that contains various control bits for demodulator operation, whose selected values are contingent upon the current values of control signals in_box, in_region, and TR. For example, Adaptation Select LUT 915 outputs two control bits that govern mux 910, which selects the error term used in update of Equalizer 175 coefficients, choosing from among LMS, CMA, combined LMS and CMA (using combining weight λ(n)), and zero, at each baud instance. An example adaptation strategy may be to use an LMS error term, derived from training data, when TR is true; CMA when in_region is true and in_box is not; a combination error term when both in_region and in_box are true, and a zero-valued error term otherwise. The Adaptation Select LUT 915 can be programmed to accommodate a variety of adaptation strategies, not limited to this example. Furthermore, the Adaptation Select LUT 915 can be reloaded over the course of demodulator operation, for example, contingent on the value of array index , an estimate of demodulated SNR, or after a prescribed amount of time, so that the adaptation strategy is itself changed over the course of demodulator operation.
The Adaptation Select LUT 915 contains control bits for other demodulator operations, too. For example, referring again to
Adaptation Select LUT 915 also contains control bits for Feedback AGC 177 operation. For the preferred embodiment of Feedback AGC 177, four operating modes have been identified, and can be selected on a symbol-by-symbol basis using the control bits for from Adaptation Select LUT 915. The association of operating mode to control signals in_box, in_region, and TR depends upon the templates stored in Region LUT 920. For example, if Region LUT 920 stores an annulus (as described in the left subplot of
If the alternative embodiment of Feedback AGC 177, illustrated in
Adaptation Select LUT 915 also contains control bits for Timing Recovery Loop 170, that determine the operating mode from among the six operating modes previously discussed in
The methods described in
One skilled in the art would understand that the equations described herein may include scaling, change of sign, or similar constant modifications that are not shown for simplicity. One skilled in the art would realize that such modifications can be readily determined or derived for the particular implementation. Thus, the described equations may be subject to such modifications, and are not limited to the exact forms presented herein.
Certain aspects and embodiments of the present invention have been described using Quadrature Amplitude Modulation (QAM) signals with complex signal processing, unless specifically noted. However, one skilled in the art would realize that the techniques described herein may be applied to a receiver processing Phase-Shift Keyed (PSK), Pulse Amplitude Modulation (PAM), or other signals.
As would be apparent to one skilled in the art, the various functions of equalization, signal combining, automatic gain control, carrier recovery, and timing recovery may be implemented with circuit elements or may also be implemented in the digital domain as processing steps in a software program. Such software may be employed in, for example, a digital signal processor, micro-controller, or general-purpose computer.
The present invention can be embodied in the form of methods and apparatuses for practicing those methods. The present invention can also be embodied in the form of program code embodied in tangible media, such as floppy diskettes, CD-ROMs, hard drives, or any other machine-readable storage medium, wherein, when the program code is loaded into and executed by a machine, such as a computer, the machine becomes an apparatus for practicing the invention. The present invention can also be embodied in the form of program code, for example, whether stored in a storage medium, loaded into and/or executed by a machine, or transmitted over some transmission medium, such as over electrical wiring or cabling, through fiber optics, or via electromagnetic radiation, wherein, when the program code is loaded into and executed by a machine, such as a computer, the machine becomes an apparatus for practicing the invention. When implemented on a general-purpose processor, the program code segments combine with the processor to provide a unique device that operates analogously to specific logic circuits.
It will be further understood that various changes in the details, materials, and arrangements of the parts which have been described and illustrated in order to explain the nature of this invention may be made by those skilled in the art without departing from the principle and scope of the invention as expressed in the following claim.
This application is a continuation application of commonly assigned U.S. patent application Ser. No. 10/782,316, filed on Feb. 19, 2004 now U.S. Pat. No. 7,580,482, and entitled “JOINT, ADAPTIVE CONTROL OF EQUALIZATION, SYNCHRONIZATION, AND GAIN IN A DIGITAL COMMUNICATIONS RECEIVER,” which claims priority to U.S. Provisional Patent Application Ser. No. 60/448,643 filed on Feb. 19, 2003, both of which are hereby incorporated herein in their entirety by reference.
Number | Name | Date | Kind |
---|---|---|---|
2138330 | Thompson | Nov 1938 | A |
2977551 | Gibson et al. | Mar 1961 | A |
3386033 | Copeland et al. | May 1968 | A |
3703685 | Simopoulos et al. | Nov 1972 | A |
3714659 | Firman | Jan 1973 | A |
3774218 | Fowler | Nov 1973 | A |
3838444 | Loughlin et al. | Sep 1974 | A |
3873771 | Kleinerman et al. | Mar 1975 | A |
3924060 | Bedford | Dec 1975 | A |
3974449 | Falconer | Aug 1976 | A |
3984624 | Waggener | Oct 1976 | A |
3999005 | Dickinson | Dec 1976 | A |
4081497 | Tokumitsu et al. | Mar 1978 | A |
4115778 | Snow | Sep 1978 | A |
4128880 | Cray, Jr. | Dec 1978 | A |
4302626 | Streeter | Nov 1981 | A |
4310920 | Hayes | Jan 1982 | A |
4316215 | Yasumoto et al. | Feb 1982 | A |
4322842 | Martinez | Mar 1982 | A |
4337479 | Tomimoto et al. | Jun 1982 | A |
4379947 | Warner | Apr 1983 | A |
4459595 | Kramer et al. | Jul 1984 | A |
4459681 | Ohtsuka | Jul 1984 | A |
4469437 | Yuasa et al. | Sep 1984 | A |
4476484 | Haskell | Oct 1984 | A |
4513415 | Martinez | Apr 1985 | A |
4523225 | Masuda et al. | Jun 1985 | A |
4535352 | Haskell | Aug 1985 | A |
4551011 | Yuasa et al. | Nov 1985 | A |
4575225 | Yuasa et al. | Mar 1986 | A |
4576470 | Yuasa et al. | Mar 1986 | A |
4587495 | Osawa et al. | May 1986 | A |
4589011 | Rzeszewski | May 1986 | A |
4602260 | Lindenmeier et al. | Jul 1986 | A |
4607230 | Kaku et al. | Aug 1986 | A |
4704025 | Yuasa et al. | Nov 1987 | A |
4750036 | Martinez | Jun 1988 | A |
4754233 | Pickett | Jun 1988 | A |
4821097 | Robbins | Apr 1989 | A |
4870489 | Ducret | Sep 1989 | A |
4879606 | Walter et al. | Nov 1989 | A |
4882614 | Kageyama et al. | Nov 1989 | A |
4882725 | Noda et al. | Nov 1989 | A |
4907218 | Inoue et al. | Mar 1990 | A |
4926244 | Isnardi | May 1990 | A |
4928177 | Martinez | May 1990 | A |
4944032 | Kageyama et al. | Jul 1990 | A |
4958230 | Jonnalagadda et al. | Sep 1990 | A |
4985769 | Yasumoto et al. | Jan 1991 | A |
5006926 | Tsinberg | Apr 1991 | A |
5019830 | Harada | May 1991 | A |
5036386 | Yasumoto et al. | Jul 1991 | A |
5038402 | Robbins | Aug 1991 | A |
5087975 | Citta et al. | Feb 1992 | A |
5093718 | Hoarty et al. | Mar 1992 | A |
5103295 | Uwabata et al. | Apr 1992 | A |
5103310 | Gibson et al. | Apr 1992 | A |
5126998 | Stem, Jr. | Jun 1992 | A |
5132988 | Fisher et al. | Jul 1992 | A |
5134464 | Basile et al. | Jul 1992 | A |
5142353 | Kageyama et al. | Aug 1992 | A |
5151783 | Faroudja | Sep 1992 | A |
5172126 | Naito | Dec 1992 | A |
5177604 | Martinez | Jan 1993 | A |
5214501 | Cavallerano et al. | May 1993 | A |
5220420 | Hoarty et al. | Jun 1993 | A |
5235619 | Beyers, II et al. | Aug 1993 | A |
5243423 | DeJean et al. | Sep 1993 | A |
5247351 | Cho | Sep 1993 | A |
5270817 | Hayashi et al. | Dec 1993 | A |
5276507 | Uwabata et al. | Jan 1994 | A |
5287180 | White | Feb 1994 | A |
5291289 | Hulyalkar et al. | Mar 1994 | A |
5293633 | Robbins | Mar 1994 | A |
5321514 | Martinez | Jun 1994 | A |
5327460 | Batruni | Jul 1994 | A |
5357284 | Todd | Oct 1994 | A |
5379324 | Mueller et al. | Jan 1995 | A |
5386239 | Wang et al. | Jan 1995 | A |
5412720 | Hoarty | May 1995 | A |
5430661 | Fisher et al. | Jul 1995 | A |
5442403 | Yasumoto et al. | Aug 1995 | A |
5448299 | Yang et al. | Sep 1995 | A |
5461426 | Limberg et al. | Oct 1995 | A |
5534933 | Yang | Jul 1996 | A |
5550578 | Hoarty et al. | Aug 1996 | A |
5550579 | Martinez | Aug 1996 | A |
5557316 | Hoarty et al. | Sep 1996 | A |
5559559 | Jungo et al. | Sep 1996 | A |
5563664 | Yang et al. | Oct 1996 | A |
5579055 | Hamilton et al. | Nov 1996 | A |
5585975 | Bliss | Dec 1996 | A |
5586121 | Moura et al. | Dec 1996 | A |
5587743 | Montgomery et al. | Dec 1996 | A |
5589872 | Martinez | Dec 1996 | A |
5596361 | Martinez | Jan 1997 | A |
5600573 | Hendricks et al. | Feb 1997 | A |
5617148 | Montgomery | Apr 1997 | A |
5666168 | Montgomery et al. | Sep 1997 | A |
5673293 | Scarpa et al. | Sep 1997 | A |
5694424 | Ariyavisitakul | Dec 1997 | A |
5714965 | Taguchi | Feb 1998 | A |
5734853 | Hendricks et al. | Mar 1998 | A |
RE35774 | Moura et al. | Apr 1998 | E |
5742902 | Shore | Apr 1998 | A |
5809086 | Ariyavisitakul | Sep 1998 | A |
5812184 | Martinez | Sep 1998 | A |
5818441 | Throckmorton et al. | Oct 1998 | A |
5818845 | Moura et al. | Oct 1998 | A |
5828655 | Moura et al. | Oct 1998 | A |
5859852 | Moura et al. | Jan 1999 | A |
5881302 | Omata | Mar 1999 | A |
5909253 | Jonnalagadda et al. | Jun 1999 | A |
5946047 | Levan | Aug 1999 | A |
5946048 | Levan | Aug 1999 | A |
5946322 | Moura et al. | Aug 1999 | A |
5946351 | Ariyavisitakul et al. | Aug 1999 | A |
5956346 | Levan | Sep 1999 | A |
5959660 | Levan | Sep 1999 | A |
5959997 | Moura et al. | Sep 1999 | A |
5961603 | Kunkel et al. | Oct 1999 | A |
5994891 | Hubbell | Nov 1999 | A |
5995168 | Yagi | Nov 1999 | A |
5995563 | Ben-Efraim et al. | Nov 1999 | A |
6005850 | Moura et al. | Dec 1999 | A |
6012161 | Ariyavisitakul et al. | Jan 2000 | A |
6016316 | Moura et al. | Jan 2000 | A |
6018526 | Liu et al. | Jan 2000 | A |
6018764 | Field et al. | Jan 2000 | A |
6034678 | Hoarty et al. | Mar 2000 | A |
6047159 | Powell et al. | Apr 2000 | A |
6066993 | Yamamoto et al. | May 2000 | A |
6073030 | Nair et al. | Jun 2000 | A |
6100883 | Hoarty | Aug 2000 | A |
6104727 | Moura et al. | Aug 2000 | A |
6192384 | Dally et al. | Feb 2001 | B1 |
6201536 | Hendricks et al. | Mar 2001 | B1 |
6226323 | Tan et al. | May 2001 | B1 |
6240133 | Sommer et al. | May 2001 | B1 |
6249180 | Maalej et al. | Jun 2001 | B1 |
6268774 | Soumyanath | Jul 2001 | B1 |
6305020 | Hoarty et al. | Oct 2001 | B1 |
6310548 | Stephens, Jr. et al. | Oct 2001 | B1 |
6314134 | Werner et al. | Nov 2001 | B1 |
6341360 | Abdelilah et al. | Jan 2002 | B1 |
6356586 | Krishnamoorthy et al. | Mar 2002 | B1 |
6366613 | Sommer et al. | Apr 2002 | B2 |
6370160 | Knutson et al. | Apr 2002 | B1 |
6370571 | Medin, Jr. | Apr 2002 | B1 |
6377529 | Lee et al. | Apr 2002 | B1 |
6385237 | Tsui et al. | May 2002 | B1 |
6426973 | Madhow et al. | Jul 2002 | B1 |
6433835 | Hartson et al. | Aug 2002 | B1 |
6501310 | Takami | Dec 2002 | B2 |
6510553 | Hazra | Jan 2003 | B1 |
6563862 | Knutson et al. | May 2003 | B1 |
6574293 | Khayrallah et al. | Jun 2003 | B1 |
6665790 | Glossner, III et al. | Dec 2003 | B1 |
6680971 | Tazebay et al. | Jan 2004 | B1 |
6754277 | Heinzelman et al. | Jun 2004 | B1 |
6754293 | Henriksson et al. | Jun 2004 | B1 |
6754715 | Cannon et al. | Jun 2004 | B1 |
6791995 | Azenkot et al. | Sep 2004 | B1 |
6816548 | Shiue et al. | Nov 2004 | B1 |
6829298 | Abe et al. | Dec 2004 | B1 |
6834039 | Kelly et al. | Dec 2004 | B1 |
6904098 | Isaksen et al. | Jun 2005 | B1 |
6980602 | Kleinerman et al. | Dec 2005 | B1 |
6985549 | Biracree et al. | Jan 2006 | B1 |
7027500 | Casas et al. | Apr 2006 | B1 |
7027528 | Liu et al. | Apr 2006 | B2 |
7031405 | Touzni et al. | Apr 2006 | B1 |
7158567 | Wang et al. | Jan 2007 | B2 |
7580482 | Endres et al. | Aug 2009 | B2 |
20010022813 | Tan et al. | Sep 2001 | A1 |
20020056140 | Oishi et al. | May 2002 | A1 |
20020059597 | Kikinis et al. | May 2002 | A1 |
20020104083 | Hendricks et al. | Aug 2002 | A1 |
20020186764 | Amin et al. | Dec 2002 | A1 |
20030035468 | Corbaton et al. | Feb 2003 | A1 |
20030095590 | Fuller et al. | May 2003 | A1 |
20030156603 | Rakib et al. | Aug 2003 | A1 |
20030202568 | Choi et al. | Oct 2003 | A1 |
20040042566 | Eidson et al. | Mar 2004 | A1 |
20060181797 | Sugawara et al. | Aug 2006 | A1 |
Number | Date | Country |
---|---|---|
0308241 | Mar 1989 | EP |
0577351 | Jan 1994 | EP |
1156634 | Nov 2001 | EP |
WO-9222983 | Dec 1992 | WO |
WO-9423470 | Oct 1994 | WO |
WO-9501676 | Jan 1995 | WO |
WO-9955087 | Oct 1999 | WO |
WO-0072446 | Nov 2000 | WO |
WO-0207425 | Jan 2002 | WO |
WO-02084965 | Oct 2002 | WO |
Number | Date | Country | |
---|---|---|---|
20100061488 A1 | Mar 2010 | US |
Number | Date | Country | |
---|---|---|---|
60448643 | Feb 2003 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 10782316 | Feb 2004 | US |
Child | 12546605 | US |