A Beamforming Method for Blind Calibration of Time-Interleaved A/D Converters

Similar documents
Framework for functional tree simulation applied to 'golden delicious' apple trees

Optimal Control of PDEs


An Introduction to Optimal Control Applied to Disease Models

LA PRISE DE CALAIS. çoys, çoys, har - dis. çoys, dis. tons, mantz, tons, Gas. c est. à ce. C est à ce. coup, c est à ce

An Example file... log.txt

AN IDENTIFICATION ALGORITHM FOR ARMAX SYSTEMS

ETIKA V PROFESII PSYCHOLÓGA

The University of Bath School of Management is one of the oldest established management schools in Britain. It enjoys an international reputation for

4.3 Laplace Transform in Linear System Analysis

OC330C. Wiring Diagram. Recommended PKH- P35 / P50 GALH PKA- RP35 / RP50. Remarks (Drawing No.) No. Parts No. Parts Name Specifications

" #$ P UTS W U X [ZY \ Z _ `a \ dfe ih j mlk n p q sr t u s q e ps s t x q s y i_z { U U z W } y ~ y x t i e l US T { d ƒ ƒ ƒ j s q e uˆ ps i ˆ p q y

A General Procedure to Design Good Codes at a Target BER

Vectors. Teaching Learning Point. Ç, where OP. l m n

Glasgow eprints Service

Max. Input Power (W) Input Current (Arms) Dimming. Enclosure

Juan Juan Salon. EH National Bank. Sandwich Shop Nail Design. OSKA Beverly. Chase Bank. Marina Rinaldi. Orogold. Mariposa.

A Functional Quantum Programming Language

QUESTIONS ON QUARKONIUM PRODUCTION IN NUCLEAR COLLISIONS

Principal Secretary to Government Haryana, Town & Country Planning Department, Haryana, Chandigarh.

ALTER TABLE Employee ADD ( Mname VARCHAR2(20), Birthday DATE );


Pulse Shaping and ISI (Proakis: chapter 10.1, 10.3) EEE3012 Spring 2018

SOLAR MORTEC INDUSTRIES.

Pharmacological and genomic profiling identifies NF-κB targeted treatment strategies for mantle cell lymphoma

Some emission processes are intrinsically polarised e.g. synchrotron radiation.

Redoing the Foundations of Decision Theory

Manifold Regularization

New BaBar Results on Rare Leptonic B Decays

ELEN E4810: Digital Signal Processing Topic 11: Continuous Signals. 1. Sampling and Reconstruction 2. Quantization

General Neoclassical Closure Theory: Diagonalizing the Drift Kinetic Operator

F O R SOCI AL WORK RESE ARCH

ELEG-636: Statistical Signal Processing

Symbols and dingbats. A 41 Α a 61 α À K cb ➋ à esc. Á g e7 á esc. Â e e5 â. Ã L cc ➌ ã esc ~ Ä esc : ä esc : Å esc * å esc *

Towards a High Level Quantum Programming Language

DISSIPATIVE SYSTEMS. Lectures by. Jan C. Willems. K.U. Leuven

Continuing Education and Workforce Training Course Schedule Winter / Spring 2010

Connection equations with stream variables are generated in a model when using the # $ % () operator or the & ' %

B œ c " " ã B œ c 8 8. such that substituting these values for the B 3 's will make all the equations true

Planning for Reactive Behaviors in Hide and Seek

February 17, 2015 REQUEST FOR PROPOSALS. For Columbus Metropolitan Library. Issued by: Purchasing Division 96 S. Grant Ave. Columbus, OH 43215

Towards a numerical solution of the "1/2 vs. 3/2" puzzle

The Fourier series are applicable to periodic signals. They were discovered by the

Obtain from theory/experiment a value for the absorbtion cross-section. Calculate the number of atoms/ions/molecules giving rise to the line.

UNIQUE FJORDS AND THE ROYAL CAPITALS UNIQUE FJORDS & THE NORTH CAPE & UNIQUE NORTHERN CAPITALS

Cairns Hospital: Suspected Acute Coronary Syndrome Pathways. DO NOT USE if a non cardiac cause for the chest pain can be diagnosed

Plasma diagnostics and abundance determinations for planetary nebulae current status. Xiaowei Liu Department of Astronomy, Peking University

1. Allstate Group Critical Illness Claim Form: When filing Critical Illness claim, please be sure to include the following:

6. Convert 5021 centimeters to kilometers. 7. How many seconds are in a leap year? 8. Convert the speed 5.30 m/s to km/h.

1. CALL TO ORDER 2. ROLL CALL 3. PLEDGE OF ALLEGIANCE 4. ORAL COMMUNICATIONS (NON-ACTION ITEM) 5. APPROVAL OF AGENDA 6. EWA SEJPA INTEGRATION PROPOSAL

Optimization of a parallel 3d-FFT with non-blocking collective operations

Human genome contains ~30K genes Not all genes are expressed at same time

EE 521: Instrumentation and Measurements

T i t l e o f t h e w o r k : L a M a r e a Y o k o h a m a. A r t i s t : M a r i a n o P e n s o t t i ( P l a y w r i g h t, D i r e c t o r )

Queues, Stack Modules, and Abstract Data Types. CS2023 Winter 2004

Testing SUSY Dark Matter

FASTGEO - A HISTOGRAM BASED APPROACH TO LINEAR GEOMETRIC ICA. Andreas Jung, Fabian J. Theis, Carlos G. Puntonet, Elmar W. Lang

Matrices and Determinants

Examination paper for TFY4240 Electromagnetic theory

Vector analysis. 1 Scalars and vectors. Fields. Coordinate systems 1. 2 The operator The gradient, divergence, curl, and Laplacian...

$%! & (, -3 / 0 4, 5 6/ 6 +7, 6 8 9/ 5 :/ 5 A BDC EF G H I EJ KL N G H I. ] ^ _ ` _ ^ a b=c o e f p a q i h f i a j k e i l _ ^ m=c n ^

LEAD IN DRINKING WATER REPORT

Optimal data-independent point locations for RBF interpolation

Complex Analysis. PH 503 Course TM. Charudatt Kadolkar Indian Institute of Technology, Guwahati

SCOTT PLUMMER ASHTON ANTONETTI

Sensitivity of hybrid filter banks A/D converters to analog realization errors and finite word length

Damping Ring Requirements for 3 TeV CLIC

Loop parallelization using compiler analysis

Front-end. Organization of a Modern Compiler. Middle1. Middle2. Back-end. converted to control flow) Representation

A Theory of Universal AI

Using Rational Approximations For Evaluating The Reliability of Highly Reliable Systems

MARKET AREA UPDATE Report as of: 1Q 2Q 3Q 4Q

Constructive Decision Theory

PRECISE ORBIT DETERMINATION OF GPS SATELLITES FOR REAL TIME APPLICATIONS

Chapter 12 Variable Phase Interpolation

Electronic Devices and Circuits Lecture 15 - Digital Circuits: Inverter Basics - Outline Announcements. = total current; I D

Harlean. User s Guide

The Study of 3 rd Overtone Crystal Oscillator

Ed S MArket. NarROW } ] T O P [ { U S E R S G U I D E. urrrrrrrrrrrv

INTRODUCTION TO DELTA-SIGMA ADCS

First-principles calculations of insulators in a. finite electric field

CHAPTER 5 ESTABLISHING THEORETICAL / TARGET VALUES FOR DENSITY & MOISTURE CONTENT

Books. Book Collection Editor. Editor. Name Name Company. Title "SA" A tree pattern. A database instance

PART IV LIVESTOCK, POULTRY AND FISH PRODUCTION

New method for solving nonlinear sum of ratios problem based on simplicial bisection

S T A T E B U D G E T : T A X A M E N D M E N T S

Exponentially Rare Large Fluctuations: Computing their Frequency

DCSP-2: Fourier Transform

Slide Set Data Converters. Digital Enhancement Techniques


Automatic Control III (Reglerteknik III) fall Nonlinear systems, Part 3

R for beginners. 2 The few things to know before st. 8 Index 31

Emefcy Group Limited (ASX : EMC) ASX Appendix 4D Half Year Report and Interim Financial Statements for the six months ended 30 June 2016

THE ANALYSIS OF THE CONVECTIVE-CONDUCTIVE HEAT TRANSFER IN THE BUILDING CONSTRUCTIONS. Zbynek Svoboda

Tutorial Lectures on MCMC II

Unit 3. Digital encoding

On Shostak s Combination of Decision Procedures

PROGRAM DEVELOPMENTS FOR MODELING GROUNDWATER FLOW IN THREE-DIMENSIONAL HETEROGENEOUS AQUIFERS WITH MODFLOW AND MODFLOWP

Safety Mats. Presence Sensing Safety Devices MatGuard Mats. Specifications. Description. Features

DYNAMIC SAMPLED-DATA OUTPUT REGULATION: A SWITCHING ADAPTIVE APPROACH. HUANG Yaxin 1 LIU Yungang 1 ZHANG Jifeng 2,3

Transcription:

A Beamforming Method for Blind Calibration of Time-nterleaved A/D Converters Bernard C. Levy University of California, Davis Joint wor with Steve Huang September 30, 200

Outline 1 Motivation Problem Formulation Blind Calibration Method Convergence Simulations and Conclusions

Motivation 2 Applications requiring ADCs operating at high data rates Time-interleaved ADCs. Constituent ADCs have gain, offset, timing mismatches that need to be estimated and corrected in the digital domain. Correction achieved by digital filter bans operating on ADCs outputs (Johansson and Lowenborg, 2002 Prendergast et al. 2004). Requires 10 to 20 % excess samples. Timing offset estimation can be performed either with test signals or blindly. Blind methods do not lower ADC throughput and can adjust to changes online.

Outline 3 Motivation Problem Formulation Blind Calibration Method Convergence Simulations and Conclusions

Problem Formulation 4 f is a CT bandlimited signal with bandwidth, can be recovered from its samples if. nstead of using a single fast ADC, we employ at. slow ADCs operating

Problem Formulation (cont d) PSfrag replacements!!#... Analog Mux Digital Mux

for 8 + $ + % Problem Formulation (cont d) 6 Due to timing offsets, quantization errors and thermal noise, output of $-th ADC )( * =>?A@ * B C ( 0, where,.- / 21 WGN and $43 )(!6 where = sampling period of slow ADCs,!6 timing offset of $-th ADC measured wrt 1st ADC f 94: 94: <!6, analysis filter ban model for D:

Q V Q W FHS!TAU JLK M W FHO E W F O P W F O FHSLTAU JLK M V FHO E V F O P V F O Q R FHS!TAU JLK M R FHO E R F O P R F O E FHG JLK M N FHO E N F O P N F O Problem Formulation (cont d)

[ g \ \ X C j \ h i C f X Y \ \ \ a b [e\ c 9d: 94: 94: g < 6hf i C 94: < f g 6 94: g < 6f a b []\ 9Z: 9: 9: 9: ` 9: a b Problem Formulation (cont d) 8 For D, let = DTFT of exact ADC outputs, and vector of alias components of fast sampled sequence, and = vector of timing mismatches.

~ y } z z x z z x x o x z z x w qr t s / + $ + D v g f l X Y z 9 6 : h i < =>? C z z z { n 9Z: / p u 9Z: <!6 =>?@ CB 9d: / n 9d: / o 9Z: Dml 9Z: / X c 9Z: / Problem Formulation (cont d) 9 We have where and = Vandermonde matrix with

o Problem Formulation (cont d) 10 By inverting \, we can find synthesis filters 9: such that 9d: 9d: 9d: \ For small representation. s, filters admit a closed-form 1st-order Farrow Synthesis filter ban for D:

ƒ M W F O LK W FHS!TAU M V F O LK V FHSLTAU M R F O LK R FHS!TAU M N F O LK N FHS!TAU E F G Problem Formulation (cont d) 11

Outline 12 Motivation Problem Formulation Blind Calibration Method Convergence Simulations and Conclusions

DŠ y } x z x x o x z x w withl } \ X Y x \ g x f w g i C D Let z 6 z 6 z 6 z { 9: / n 9: / o 9Z: D l 9Z: / y \ 94: < f g 6h i C 9 : { 94: < f = h Blind Calibration 13 D, if is ran deficient for ˆ :, where is a reduced Vandermonde matrix. 3 = % of excess samples, the alias matrix

D 6 j 6 j j j 0 6 Œ C Œ i C Œg C i C Œ C Œ h f 9Z: / X Y 9Z: - ` j ` j ` j 6 j j 6 j j ` j < Ž f g h < Ž f g h <! Ž < Ž f g @h <! Ž g @ f Œ` B <! ˆŽ f C Blind Calibration (cont d) 14 Can find a nulling filter ban 9 : / such that for ˆ.

D 1 D a D 6 ` š 3 ( 3 ( 3 š 3 3 š 3 ( ( 3 b o - ~ b [ ` 9Z: / b n 9Z: / Blind Calibration (cont d) 1 Structure of nulling filter ban: for ˆ, =0 otherwise, where satisfies Set. Then, for small s:

% / % where D š 9 œcb ` : 6 p 0 p ~ Ÿž 9 : ideal lowpass filter of bandwidth 0 9Z: )(!6 9Z: / 9Z: /!6 <!6 =>?A@ Ÿž 9Z: Blind Calibration (cont d) 16 Consider the 1st-order Farrow approximation. Let Consider the adaptive null-steering structure

± ³ «µ ³ ± ³ «µ ³ ± ³ «µ ³ µ ³!ª ² ª! ²!ª!ª ² ² ª ² ª ª! ²!!ª ² Blind Calibration (cont d) 1

` 3 D š x x x x w w ` ¹» ¾À Á 9 / º y 3 { } ( y 3 )( ` 3 1 3 ` { } 3 ` )( 1 9 / º º )( º º ¼ 9 / º ½21 Blind Calibration (cont d) 18 Consider the objective function where = nulling filter output. We have

º Blind Calibration (cont d) 19  9 Á 9 3 / º ¾Ã / º obtained by stochastic gradient algorithm ( º º where -  step size, with initial condition Ä ~ º

Outline 20 Motivation Problem Formulation Blind Calibration Method Convergence Simulations and Conclusions

º  Convergence 21 Use ODE/stochastic averaging method. Assume WSS. small, zero-mean Write adaptive algorithm as ( ÂÆÅ º Ç b / Ç ( º where of the algorithm, ¼4È b É ½ b. Due to the stochastic gradient structure º º ½ ¾Ã º ¹ Á / / Ç» ¼ ¹ 3 so º ʺ Ê / = Lyapunov function for ODE ¹ so ODE trajectories converge to a minimum of º.

Ò 1 1 Ô Õ i h 6 Ï Ô Õ Ð Ñ Î Ó ( ( Ò º º ¹ š Ï º º º 1 0 h i 6 j 9Z: Ê = j j 1 9Z: Ê j ÐÍ 3 Ñ 3 º 3 3 ( º 1ÌË º Í 3 Î 3 3 º 3 3 Convergence (cont d) 22 For small and with where we neglect cubic terms.

º ' Ø Ù ½ ½ ( 3 Ï º Convergence (cont d) 23 For noiseless case, unique minimum of ¹for small offset and offset estimates is if -and -, so and. ¼ 3 / ¼ Ö 21 / Ö 21 must have power in band Ensures that as, ÂÆÚ / whereú= positive definite matrix.

Outline 24 Motivation Problem Formulation Blind Calibration Method Convergence Simulations and Conclusions

- - 1 1~ 1-6 ` - -~ -~ 3 ½ - ~ ÜÛ ~ ÜÛ Bandlimited WGN input 2 0.02 0.02 0.01 Simulation parameters timing estimation 0.01 0.00 0 0.00 ch2 ch3 ch4 Signal bandwidth: Â -1~, -. ¼ 3, / -. 0.01 0.01 design Ü. 0.02 0 2 4 6 8 # of Samples 10 12 14 x 10 4

- - ~ - 6 -~ - ~ ÜÛ - ~ - ÜD~ Multitone sinusoidal input 26 0.02 0.02 Simulation parameters timing estimation 0.01 0.01 0.00 0 ch2 ch3 ch4 nput frequencies:. -1~, -,,, 0.00 Â. 0.01 design Ü. 0.01 0 0. 1 1. 2 2. 3 3. 4 x 10 4

Calibrated ADC output 2 80 80 60 60 40 40 Power Spectrum 20 0 20 Power Spectrum 20 0 20 40 40 60 60 80 1 0.8 0.6 0.4 0.2 0 0.2 0.4 0.6 0.8 1 Frequency in Radian (x π) 80 1 0.8 0.6 0.4 0.2 0 0.2 0.4 0.6 0.8 1 Frequency in Radian (x π) Uncalibrated ADC, 0dB SNR Calibrated ADC, 0dB SNR

Conclusion 28 Blind calibration of time-interleaved ADCs presented, requires 10 to 20 % oversampling and intermittent excitation of certain frequency bands. Simulated for up to 16 channels, but 2 or 4 channels primary interest for today s ADC technology. Postprocessing of analog circuits with mismatched components source of interesting adaptive signal processing problems.

29 Than you!!