FAST RECOVERY ALGORITHMS FOR TIME ENCODED BANDLIMITED SIGNALS. Ernő K. Simonyi
|
|
- Jeremy Henderson
- 6 years ago
- Views:
Transcription
1 FAST RECOVERY ALGORITHS FOR TIE ENCODED BANDLIITED SIGNALS Aurel A Lazar Dept of Electrical Engineering Columbia University, New York, NY 27, USA aurel@eecolumbiaedu Ernő K Simonyi National Council of Hungary for Information and Communication Technology, H-73 Budapest, Hungary simonyi@nhithu László T Tóth Dept of Telecom and edia Informatics, Budapest Univ of Technology and Economics, H-7 Budapest, Hungary tothl@tmitbmehu ABSTRACT Time encoding is a real-time asynchronous mechanism of mapping amplitude information into a time sequence We investigate fast algorithms for the recovery of time encoded bandlimited signals and construct an algorithm that has provably low computational complexity We also devise a fast algorithm that is parameter-insensitive INTRODUCTION Time encoding is a real-time asynchronous mechanism of mapping amplitude information into a time sequence In [4] time encoding and irregular sampling were shown to be largely equivalent modalities of information representation and recovery For communications applications, however, irregular sampling requires the transmission of both amplitude and sampling time information Time encoding requires only the transmission of the time sequence Thus, capacity requirements for time encoded versus irregular sampled bandlimited signals are lower by a factor of two [5] For the case of irregular sampling, fast algorithms for signal recovery have been extensively studied in the literature (eg, [2], [3], [6]) While these algorithms are an excellent starting point to search for fast algorithms for recovery of time encoded bandlimited signals, they can not be directly applied to the time encoding case The reasons are purely technical: the method developed for the irregular sampling case calls for reducing the solution of a linear systems of equations to a Toeplitz matrix inversion Our method for obtaining a fast algorithm for the time encoding case is based on the observation that the indefinite integral of the same signal can be directly recovered from its amplitude values sampled at instances provided by the time sequence In the process we find that the complexity of recovering the indefinite integral of an arbitrary bandlimited signal from irregular samples is essentially the same as the complexity of recovering the same signal from its time encoded sequence In addition, we find that the condition numbers of the pseudo-inverse matrices that arise in both formulations can be chosen to fall in the same range This paper is organized as follows A brief overview of the classical recovery algorithm for irregular sampling followed by a review of an efficient recovery for irregular sampling is presented in section 2 In section 3 we derive a fast recovery algorithm for time encoded bandlimited signals In section 4, a parameter-insensitive reconstruction algorithm is devised Simulation results are presented in section 5 2 FAST RECOVERY ALGORITH FOR IRREGULAR SAPLING 2 The Classical Recovery Algorithm for Irregular Sampling Given a set of irregular sampling times (t k ), k Z, we shall represent the Ω-bandlimited signal x = x(t), t R, as x(t) = g T c = g(t t k ), () where c = [ ] is a vector of weights, T denotes the transpose and g = [g(t t k )] with g(t) = sinωt πt is the vector of t k -shifted impulse response of an ideal lowpass filter with bandwidth Ω Also, R and Z above denote the set of real numbers and integers, respectively Denoting by: [q] l = x( ) and [G] l,k = g( t k ), (3) for all k, k Z, and l, l Z, it is easy to see that (2) q = Gc (4) Proceedings of ICASSP 25, arch 8-23, 24, Philadelphia c IEEE 25
2 If the average density of the s k s is at or above the Nyquist rate then, c can be evaluated as c = G + q, (5) where G + is the pseudo-inverse (oore-penrose) of G The practical solution, however, is challenging because G is typically ill-conditioned, and, q, c and G are infinite dimensional In summary, the bandlimited signal x can be approximately represented by equation (3) with the weighting coefficients evaluated following () The low complexity of the algorithm () is due to the fact that T is a Hermitian Toeplitz matrix 3 FAST RECOVERY ALGORITH FOR TIE ENCODING 22 A Fast Recovery Algorithm for Irregular Sampling For irregular sampling, let us consider with T = π/ω and α = (2 + )T g(t) = α e jn Ω t = α sin ( 2+ 2 Ωt) sin ( ) (6) Ωt 2 instead of the impulse response defined in (2) When tends to infinity this function converges to the original impulse response At the same time, the latter choice of g(t) is both Ω-bandlimited and periodic with a period 2 T In analogy to (), we define the function f = f(t), t R, by f(t) = α e jn Ω (t t k) (7) Assuming an appropriate norm, f is an approximation of x as provided that f( ) = x( ) for all l, l Z [3] Evaluating the above equality at t = we obtain f( ) = α e jn Ω ( t k ), (8) for all l, l Z By denoting [q] l = x( ) and [S] m,l = e jm Ω, equation (8) above after multiplication with S becomes Sq = αss H Sc, (9) or d = αt + Sq, () 3 Classical Recovery Algorithm for Time Encoding In the time encoding case an Ω-bandlimited signal x = x(t), t R, is represented as a discrete strictly increasing time sequence (t k ), k Z The time sequence (t k ), k Z, is generated using a Time Encoding achine (TE) [4] An example of a TE is depicted in Figure [4] It consists of an ideal integrator and a noninverting Schmitt trigger in a feedback arrangement The output z(t) takes the values b or b at transition times denoted by t k It can be shown [4] that this circuit is described by for all l, l Z, by the recursive equation x(u)du = ( ) k [2κδ b(t k+ t k ] (4) Integrator Z dt κ where d = αsc and T = αss H () [q] l = x(u)du and [G] l,k = g(u s k )du, For the record (6) [d] n = α e jn Ω t k and [T] m,n = α where s k = (t k+ + t k )/2, it is easy to see that e j(n m) Ω q = Gc and c = G + q (7) l Z (2) Note that T is both Toeplitz and Hermitian since [T] m,n = 32 Reformulation of the Classical Recovery Algorithm [T] m n and T H = T (the superscript H indicates conjugate transposition) Finally with [d] n = d n, x(u)du, t If x is a bandlimited function, then so is t R, and therefore: t f(t) = d n e jn Ω t, (3) x(u)du = g T c = g(t t k ), (8) 2 x(t) + y(t) Noninverting Schmitt trigger δ b z b Fig An Example of a Time Encoding achine δ y z(t) The recovery of the signal x is achieved via the representation x(t) = g T c = g(t s k ) (5) with an appropriate set of weights, k Z Denoting by:
3 where g = [g(t t k )], g(t) is given by (2) and c = [ ] is an appropriate set of coefficients Since we have x(u)du = [g(+ t k ) g( t k )], (9) q = PGc, (2) where [P] l,k = δ l+,k δ l,k (using Kronecker s notation) and [G] l,k = [g( t k )] (same as in (3)) and thus with [P ] i,k = c = G + P q (2) { if i k if i > k (22) Remark Note that by multiplying both sides of equation (2) with +, we obtain: (+ )[P q] l = (+ ) sin Ω( t k ) π( t k ) }{{} [ ] l,k (23) In the time encoding case, [q] l = g(u s k )du = (+ ) sin Ω(ξ (24) l s k ) π(ξ l s k ) for some ξ l [t tl, + ] Therefore, the elements of the matrix identified by the lower brace in equation (23) are approximately equal to the elements of the G matrix for time encoding This points to the close relationship between the recovery algorithm of a time encoded bandlimited signal x and the recovery of an irregularly sampled integrated signal t x(u)du For both recovery methods the same time sequence is used 33 Fast Recovery Algorithm Now, (23) can be transformed into a Toeplitz system in the same way as carried out in Sec 22 Replacing the sinc term by the approximation introduced in (6) gives: (+ )[P q] l = α(+ ) e jn Ω ( t k ) Denoting by D = diag(+ ), l Z, we have in matrix form DP q = αds H Sc ultiplying both sides by S from the left gives: Equivalently with we have SDP q = αsds H Sc T = αsds H, d = αsc, (25) d = αt + SDP q (26) The matrix T is both Toeplitz and Hermitian since [T] n,m = α (t k+ t k )e j(m n) Ω t k Finally, the original signal x is approximated the function f given by f(t) = jω nd n e jn Ω t, (27) with the vector d given by (26) 4 PARAETER-INSENSITIVE RECONSTRUCTION The original reconstruction of a time encoded bandlimited function depends on the TE parameters κ and δ In [4] the Compensation Principle was used to address the insensitivity of the recovery algorithm with respect to the parameters of the TE A simple argument shows, however, that the method employed in [4] does not directly apply here Therefore, a different technique is needed Parameter insensitivity can be achieved if the termp q in equation (25) does not depend on κ and δ Carrying out P q by using (22) gives: [P q] l = κδ [ + ( ) i] + b L ( ) k (t k+ t k ) k=i where L denotes the size of G (ideally L ) Therefore, as seen, P q will be independent of κδ for odd values of i If - in the last column of every second line of P is changed to zero then each integral is carried out over even number of subintervals This can be achieved by adding the vector q = }{{} a [ ] q }{{} b T 3
4 to P q and defining p = P q + ab T q, (28) where for illustration purposes and in finite dimensions q q 2 p = q 3 q 4 q 5 q 6 Note that p can be calculated without using parameters δ and κ Since p = P q + ab T PP q and P q = αs H Sc, we have: p = αs H Sc + αab T PS H Sc ultiplying both sides above by SD gives: SDp = αsds H Sc + αsdab T PS H Sc Using the matrix T and vector d defined in (25) and (26), respectively, and denoting by we have: and therefore u = αsda, v T = b T PS H, (29) (T + uv T )d = SDp, d = (T + uv T ) + SDp (3) Finally, employing the result of [] (page 5, Corollary 33), the pseudo-inversion in (3) can be determined based on the pseudo-inverse of T as (T + uv T ) + = T + T+ uv T T + + v T T + u provided that +v T T + u Thus, to solve (3) we again have to calculate the pseudo- inverse of a Hermitian Toeplitz matrix and the parameter-insensitivity is also guaranteed The reconstructed signal is again given by (27) 5 SIULATION RESULTS With Ω = π 8 khz and the corresponding Nyquist period T = π/ω = 25 µs the bandlimited signal x(t) was generated by its Shannon-representation x(t) = 65 n= 35 x(nt)sinc(ω(t nt)), x(t) t (ms) Fig 2 Overall input signal x(t) and the simulation range (dashed box) for time encoding where the samples x( 35T),,x(65T) were randomly selected in the range ( 3, 3) and sinc(t) = sin t/t if t and sinc() = This is shown in Fig 2 in the range t [ 5T, 8T] The dashed box in the figure shows the simulation range t [ µs, 2538 µs] for the TE With parameters δ = 7 6, κ = /2, and c = 3 the TE simulation produced 35 trigger times t =, t,, t 35 Fig 3 shows the error signal defined as the difference between the right-hand-side and the left-hand-side of (27) in an even further reduced range t [2525 µs, µs] to decrease the boundary effects [4] e(t) 9 e RS = 543 db t (µs) Fig 3 Error signal using the original formulation in the reduced range to reduce the boundary effects Finally, Fig 4 shows the results for the RS error and the condition number of T (based on using 2 ) in terms of different values of by using the new reconstruction technique It can be seen that as the condition number increases the accuracy of the reconstruction improves Fig 5 shows the simulation results with the original (squares) and the parameter-insensitive reconstruction technique (stars) The error with the original method is the same as that shown in Fig 4 4
5 db cond(g) e RS original cond(t) e RS () Fig 4 Simulation results for the condition numbers and RS error e RS (db) Fig 5 Simulation results for the RS error 7 REFERENCES [] Campbell, SL and eyer, CD Jr, Generalized Inverses of Linear Transformations, Dover Publications, 979 [2] Feichtinger, HG, Gröchenig, K and Strohmer, T, Efficient Numerical ethods in Non-uniform Sampling Theory, Numerische athematik, Vol 69, pp , 995 [3] Gröchenig, K, Irregular Sampling, Toeplitz atrices, and the Approximation of Entire Functions of Exponential Type, athematics of Computation, Vol 68, #226, pp , April 999 [4] Lazar, AA and Tóth, LT, Perfect Recovery and Sensitivity Analysis of Time Encoded Bandlimited Signals, IEEE Transactions on Circuits and Systems-I: Regular Papers, Vol 5, No, pp , October 24, [5] Lazar, AA, Roska, T, Simonyi, EK and Toth, LT, A Time Decoding Realization with a CNN, Proceedings of Neurel 24, pp 97-2, September 23-25, 24, Belgrade [6] Strohmer, T, Numerical Analysis of the Non-uniform Sampling Problem, J Comp Appl ath, Vol 2, #-2, pp , 2 Using the same parameters as before we found that increasing the accuracy improves The nonzero eigenvalues of T are approaching to those of G in (23) and practically zero eigenvalues are introduced if is increased In particular, having 27 trigger times the error for = (2 + = 23), = 3 (2 + = 27), = 5 (2 + = 3), the corresponding error turned out to be -86 db, db, The accuracy cannot be improved below -543 db corresponding to the case when the original g(t) is used 6 CONCLUSIONS We have shown that the indefinite integral of an arbitrary bandlimited signal can be directly recovered from the time encoded sequence associated with the same bandlimited signal This simple observation enabled us to devise a fast algorithm for signal recovery We have also demonstrated that our approach can be extended to parameter insensitive signal recovery Taken together, these results shed further light on the close relationship between irregular sampling and time encoding 5
TIME ENCODING AND PERFECT RECOVERY OF BANDLIMITED SIGNALS
TIME ENCODING AND PERFECT RECOVERY OF BANDLIMITED SIGNALS Aurel A. Lazar Department of Electrical Engineering Columbia University, New York, NY 27 aurel@ee.columbia.edu László T. Tóth Department of Telecommunications
More informationChapter 5 Frequency Domain Analysis of Systems
Chapter 5 Frequency Domain Analysis of Systems CT, LTI Systems Consider the following CT LTI system: xt () ht () yt () Assumption: the impulse response h(t) is absolutely integrable, i.e., ht ( ) dt< (this
More informationChapter 5 Frequency Domain Analysis of Systems
Chapter 5 Frequency Domain Analysis of Systems CT, LTI Systems Consider the following CT LTI system: xt () ht () yt () Assumption: the impulse response h(t) is absolutely integrable, i.e., ht ( ) dt< (this
More informationHomework 4. May An LTI system has an input, x(t) and output y(t) related through the equation y(t) = t e (t t ) x(t 2)dt
Homework 4 May 2017 1. An LTI system has an input, x(t) and output y(t) related through the equation y(t) = t e (t t ) x(t 2)dt Determine the impulse response of the system. Rewriting as y(t) = t e (t
More informationDistributed Noise Shaping of Signal Quantization
1 / 37 Distributed Noise Shaping of Signal Quantization Kung-Ching Lin Norbert Wiener Center Department of Mathematics University of Maryland, College Park September 18, 2017 2 / 37 Overview 1 Introduction
More informationELEN E4810: Digital Signal Processing Topic 11: Continuous Signals. 1. Sampling and Reconstruction 2. Quantization
ELEN E4810: Digital Signal Processing Topic 11: Continuous Signals 1. Sampling and Reconstruction 2. Quantization 1 1. Sampling & Reconstruction DSP must interact with an analog world: A to D D to A x(t)
More informationAmplitude Adaptive ASDM without Envelope Encoding
Amplitude Adaptive ASDM without Envelope Encodin Kaspars Ozols and Rolands Shavelis Institute of Electronics and Computer Science 14 Dzerbenes Str., LV-16, Ria, Latvia Email: kaspars.ozols@edi.lv, shavelis@edi.lv
More informationE2.5 Signals & Linear Systems. Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & 2)
E.5 Signals & Linear Systems Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & ) 1. Sketch each of the following continuous-time signals, specify if the signal is periodic/non-periodic,
More informationChapter 3. LQ, LQG and Control System Design. Dutch Institute of Systems and Control
Chapter 3 LQ, LQG and Control System H 2 Design Overview LQ optimization state feedback LQG optimization output feedback H 2 optimization non-stochastic version of LQG Application to feedback system design
More informationModule 3 : Sampling and Reconstruction Lecture 22 : Sampling and Reconstruction of Band-Limited Signals
Module 3 : Sampling and Reconstruction Lecture 22 : Sampling and Reconstruction of Band-Limited Signals Objectives Scope of this lecture: If a Continuous Time (C.T.) signal is to be uniquely represented
More information18th European Signal Processing Conference (EUSIPCO-2010) Aalborg, Denmark, August 23-27, University of Pittsburgh Pittsburgh, PA, USA
8th European Signal Processing Conference (EUSIPCO-) Aalborg, Denmar, August 3-7, ADAPTIVE LEVEL-CROSSING SAMPLING AND RECONSTRUCTION Seda Senay, Luis F. Chaparro, Mingui Sun and Robert J. Sclabassi, Department
More informationDigital Signal Processing
COMP ENG 4TL4: Digital Signal Processing Notes for Lecture #21 Friday, October 24, 2003 Types of causal FIR (generalized) linear-phase filters: Type I: Symmetric impulse response: with order M an even
More informationEE 16B Final, December 13, Name: SID #:
EE 16B Final, December 13, 2016 Name: SID #: Important Instructions: Show your work. An answer without explanation is not acceptable and does not guarantee any credit. Only the front pages will be scanned
More informationProblem Sheet 1 Examples of Random Processes
RANDOM'PROCESSES'AND'TIME'SERIES'ANALYSIS.'PART'II:'RANDOM'PROCESSES' '''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''Problem'Sheets' Problem Sheet 1 Examples of Random Processes 1. Give
More informationGrades will be determined by the correctness of your answers (explanations are not required).
6.00 (Fall 2011) Final Examination December 19, 2011 Name: Kerberos Username: Please circle your section number: Section Time 2 11 am 1 pm 4 2 pm Grades will be determined by the correctness of your answers
More informationQ. 1 Q. 25 carry one mark each.
Q. Q. 5 carry one mark each. Q. Consider a system of linear equations: x y 3z =, x 3y 4z =, and x 4y 6 z = k. The value of k for which the system has infinitely many solutions is. Q. A function 3 = is
More informationReview of Linear Time-Invariant Network Analysis
D1 APPENDIX D Review of Linear Time-Invariant Network Analysis Consider a network with input x(t) and output y(t) as shown in Figure D-1. If an input x 1 (t) produces an output y 1 (t), and an input x
More informationVarious signal sampling and reconstruction methods
Various signal sampling and reconstruction methods Rolands Shavelis, Modris Greitans 14 Dzerbenes str., Riga LV-1006, Latvia Contents Classical uniform sampling and reconstruction Advanced sampling and
More informationDiscrete-Time Signals and Systems. Efficient Computation of the DFT: FFT Algorithms. Analog-to-Digital Conversion. Sampling Process.
iscrete-time Signals and Systems Efficient Computation of the FT: FFT Algorithms r. eepa Kundur University of Toronto Reference: Sections 6.1, 6., 6.4, 6.5 of John G. Proakis and imitris G. Manolakis,
More informationLECTURE 12 Sections Introduction to the Fourier series of periodic signals
Signals and Systems I Wednesday, February 11, 29 LECURE 12 Sections 3.1-3.3 Introduction to the Fourier series of periodic signals Chapter 3: Fourier Series of periodic signals 3. Introduction 3.1 Historical
More informationIB Paper 6: Signal and Data Analysis
IB Paper 6: Signal and Data Analysis Handout 5: Sampling Theory S Godsill Signal Processing and Communications Group, Engineering Department, Cambridge, UK Lent 2015 1 / 85 Sampling and Aliasing All of
More informationLinear Algebra and ODEs review
Linear Algebra and ODEs review Ania A Baetica September 9, 015 1 Linear Algebra 11 Eigenvalues and eigenvectors Consider the square matrix A R n n (v, λ are an (eigenvector, eigenvalue pair of matrix A
More informationSignals and Systems I Have Known and Loved. (Lecture notes for CSE 3451) Andrew W. Eckford
Signals and Systems I Have Known and Loved (Lecture notes for CSE 3451) Andrew W. Eckford Department of Electrical Engineering and Computer Science York University, Toronto, Ontario, Canada Version: December
More informationA Family of Nyquist Filters Based on Generalized Raised-Cosine Spectra
Proc. Biennial Symp. Commun. (Kingston, Ont.), pp. 3-35, June 99 A Family of Nyquist Filters Based on Generalized Raised-Cosine Spectra Nader Sheikholeslami Peter Kabal Department of Electrical Engineering
More informationChap 4. Sampling of Continuous-Time Signals
Digital Signal Processing Chap 4. Sampling of Continuous-Time Signals Chang-Su Kim Digital Processing of Continuous-Time Signals Digital processing of a CT signal involves three basic steps 1. Conversion
More informationLecture Schedule: Week Date Lecture Title
http://elec34.org Sampling and CONVOLUTION 24 School of Information Technology and Electrical Engineering at The University of Queensland Lecture Schedule: Week Date Lecture Title 2-Mar Introduction 3-Mar
More informationA Second Course in Elementary Differential Equations
A Second Course in Elementary Differential Equations Marcel B Finan Arkansas Tech University c All Rights Reserved August 3, 23 Contents 28 Calculus of Matrix-Valued Functions of a Real Variable 4 29 nth
More informationContinuous Time Signal Analysis: the Fourier Transform. Lathi Chapter 4
Continuous Time Signal Analysis: the Fourier Transform Lathi Chapter 4 Topics Aperiodic signal representation by the Fourier integral (CTFT) Continuous-time Fourier transform Transforms of some useful
More informationGrades will be determined by the correctness of your answers (explanations are not required).
6.00 (Fall 20) Final Examination December 9, 20 Name: Kerberos Username: Please circle your section number: Section Time 2 am pm 4 2 pm Grades will be determined by the correctness of your answers (explanations
More informationRepresentation of Signals and Systems. Lecturer: David Shiung
Representation of Signals and Systems Lecturer: David Shiung 1 Abstract (1/2) Fourier analysis Properties of the Fourier transform Dirac delta function Fourier transform of periodic signals Fourier-transform
More informationPROBLEMS In each of Problems 1 through 12:
6.5 Impulse Functions 33 which is the formal solution of the given problem. It is also possible to write y in the form 0, t < 5, y = 5 e (t 5/ sin 5 (t 5, t 5. ( The graph of Eq. ( is shown in Figure 6.5.3.
More informationChapter 3 : Linear Differential Eqn. Chapter 3 : Linear Differential Eqn.
1.0 Introduction Linear differential equations is all about to find the total solution y(t), where : y(t) = homogeneous solution [ y h (t) ] + particular solution y p (t) General form of differential equation
More informationDigital Baseband Systems. Reference: Digital Communications John G. Proakis
Digital Baseband Systems Reference: Digital Communications John G. Proais Baseband Pulse Transmission Baseband digital signals - signals whose spectrum extend down to or near zero frequency. Model of the
More informationVII. Bandwidth Limited Time Series
VII. Bandwidth Limited Time Series To summarize the discussion up to this point: (1) In the general case of the aperiodic time series, which is infinite in time and frequency, both the time series and
More informationEFFICIENT CONSTRUCTION OF POLYPHASE SEQUENCES WITH OPTIMAL PEAK SIDELOBE LEVEL GROWTH. Arindam Bose and Mojtaba Soltanalian
EFFICIENT CONSTRUCTION OF POLYPHASE SEQUENCES WITH OPTIAL PEAK SIDELOBE LEVEL GROWTH Arindam Bose and ojtaba Soltanalian Department of Electrical and Computer Engineering University of Illinois at Chicago
More informationDigital Signal Processing
COMP ENG 4TL4: Digital Signal Processing Notes for Lecture #20 Wednesday, October 22, 2003 6.4 The Phase Response and Distortionless Transmission In most filter applications, the magnitude response H(e
More informationLine Codes and Pulse Shaping Review. Intersymbol interference (ISI) Pulse shaping to reduce ISI Embracing ISI
Line Codes and Pulse Shaping Review Line codes Pulse width and polarity Power spectral density Intersymbol interference (ISI) Pulse shaping to reduce ISI Embracing ISI Line Code Examples (review) on-off
More informationLaplace Transform Analysis of Signals and Systems
Laplace Transform Analysis of Signals and Systems Transfer Functions Transfer functions of CT systems can be found from analysis of Differential Equations Block Diagrams Circuit Diagrams 5/10/04 M. J.
More informationExample: Bipolar NRZ (non-return-to-zero) signaling
Baseand Data Transmission Data are sent without using a carrier signal Example: Bipolar NRZ (non-return-to-zero signaling is represented y is represented y T A -A T : it duration is represented y BT. Passand
More informationECE 301 Fall 2010 Division 2 Homework 10 Solutions. { 1, if 2n t < 2n + 1, for any integer n, x(t) = 0, if 2n 1 t < 2n, for any integer n.
ECE 3 Fall Division Homework Solutions Problem. Reconstruction of a continuous-time signal from its samples. Consider the following periodic signal, depicted below: {, if n t < n +, for any integer n,
More informationDifference between Reconstruction from Uniform and Non-Uniform Samples using Sinc Interpolation
IJRRES INERNAIONAL JOURNAL OF RESEARCH REVIEW IN ENGINEERING SCIENCE & ECHNOLOGY Difference between Reconstruction from Uniform and Non-Uniform Samples using Sinc Interpolation *Apurva H. Ghate, **Dr.
More information3.2 Complex Sinusoids and Frequency Response of LTI Systems
3. Introduction. A signal can be represented as a weighted superposition of complex sinusoids. x(t) or x[n]. LTI system: LTI System Output = A weighted superposition of the system response to each complex
More informationEE301 Signals and Systems In-Class Exam Exam 3 Thursday, Apr. 19, Cover Sheet
EE301 Signals and Systems In-Class Exam Exam 3 Thursday, Apr. 19, 2012 Cover Sheet Test Duration: 75 minutes. Coverage: Chaps. 5,7 Open Book but Closed Notes. One 8.5 in. x 11 in. crib sheet Calculators
More informationMultirate signal processing
Multirate signal processing Discrete-time systems with different sampling rates at various parts of the system are called multirate systems. The need for such systems arises in many applications, including
More informationReview: Continuous Fourier Transform
Review: Continuous Fourier Transform Review: convolution x t h t = x τ h(t τ)dτ Convolution in time domain Derivation Convolution Property Interchange the order of integrals Let Convolution Property By
More informationSignal and systems. Linear Systems. Luigi Palopoli. Signal and systems p. 1/5
Signal and systems p. 1/5 Signal and systems Linear Systems Luigi Palopoli palopoli@dit.unitn.it Wrap-Up Signal and systems p. 2/5 Signal and systems p. 3/5 Fourier Series We have see that is a signal
More informationQuestion Paper Code : AEC11T02
Hall Ticket No Question Paper Code : AEC11T02 VARDHAMAN COLLEGE OF ENGINEERING (AUTONOMOUS) Affiliated to JNTUH, Hyderabad Four Year B. Tech III Semester Tutorial Question Bank 2013-14 (Regulations: VCE-R11)
More information2.161 Signal Processing: Continuous and Discrete Fall 2008
MIT OpenCourseWare http://ocw.mit.edu 2.6 Signal Processing: Continuous and Discrete Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. MASSACHUSETTS
More informationMassachusetts Institute of Technology
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.011: Introduction to Communication, Control and Signal Processing QUIZ 1, March 16, 2010 ANSWER BOOKLET
More informationCompressibility of Infinite Sequences and its Interplay with Compressed Sensing Recovery
Compressibility of Infinite Sequences and its Interplay with Compressed Sensing Recovery Jorge F. Silva and Eduardo Pavez Department of Electrical Engineering Information and Decision Systems Group Universidad
More information5742 IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 55, NO. 12, DECEMBER /$ IEEE
5742 IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 55, NO. 12, DECEMBER 2009 Uncertainty Relations for Shift-Invariant Analog Signals Yonina C. Eldar, Senior Member, IEEE Abstract The past several years
More informationESE 531: Digital Signal Processing
ESE 531: Digital Signal Processing Lec 8: February 12th, 2019 Sampling and Reconstruction Lecture Outline! Review " Ideal sampling " Frequency response of sampled signal " Reconstruction " Anti-aliasing
More information2. SPECTRAL ANALYSIS APPLIED TO STOCHASTIC PROCESSES
2. SPECTRAL ANALYSIS APPLIED TO STOCHASTIC PROCESSES 2.0 THEOREM OF WIENER- KHINTCHINE An important technique in the study of deterministic signals consists in using harmonic functions to gain the spectral
More informationLecture 6 January 21, 2016
MATH 6/CME 37: Applied Fourier Analysis and Winter 06 Elements of Modern Signal Processing Lecture 6 January, 06 Prof. Emmanuel Candes Scribe: Carlos A. Sing-Long, Edited by E. Bates Outline Agenda: Fourier
More informationNetwork Parameters of a Two-Port Filter
Contents Network Parameters of a Two-Port Filter S Parameters in the s-domain Relation between ε and ε R 2 2 Derivation of S 22 3 2 Network Parameters in the jω-domain 3 2 S Parameters in the jω-domain
More informationNonlinear and Non-Ideal Sampling: Theory and Methods
Nonlinear and Non-Ideal Sampling: 1 Theory and Methods Tsvi G. Dvorkind, Yonina C. Eldar, Senior Member, IEEE, and Ewa Matusiak Abstract We study a sampling setup where a continuous-time signal is mapped
More informationQuantization and Compensation in Sampled Interleaved Multi-Channel Systems
Quantization and Compensation in Sampled Interleaved Multi-Channel Systems The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation
More informationStabilization of persistently excited linear systems
Stabilization of persistently excited linear systems Yacine Chitour Laboratoire des signaux et systèmes & Université Paris-Sud, Orsay Exposé LJLL Paris, 28/9/2012 Stabilization & intermittent control Consider
More informationHomework: 4.50 & 4.51 of the attachment Tutorial Problems: 7.41, 7.44, 7.47, Signals & Systems Sampling P1
Homework: 4.50 & 4.51 of the attachment Tutorial Problems: 7.41, 7.44, 7.47, 7.49 Signals & Systems Sampling P1 Undersampling & Aliasing Undersampling: insufficient sampling frequency ω s < 2ω M Perfect
More informationAPPM 2360: Midterm exam 3 April 19, 2017
APPM 36: Midterm exam 3 April 19, 17 On the front of your Bluebook write: (1) your name, () your instructor s name, (3) your lecture section number and (4) a grading table. Text books, class notes, cell
More informationEncoding of Multivariate Stimuli with MIMO Neural Circuits
20 IEEE International Symposium on Information Theory Proceedings Encoding of Multivariate Stimuli with MIMO Neural Circuits Aurel A. Lazar Department of Electrical Engineering Columbia University, New
More informationUniversity of Connecticut Lecture Notes for ME5507 Fall 2014 Engineering Analysis I Part III: Fourier Analysis
University of Connecticut Lecture Notes for ME557 Fall 24 Engineering Analysis I Part III: Fourier Analysis Xu Chen Assistant Professor United Technologies Engineering Build, Rm. 382 Department of Mechanical
More informationSobolev Duals of Random Frames
Sobolev Duals of Random Frames C. Sinan Güntürk, Mark Lammers, 2 Alex Powell, 3 Rayan Saab, 4 Özgür Yılmaz 4 Courant Institute of Mathematical Sciences, New York University, NY, USA. 2 University of North
More informationSignals and Systems. Lecture 11 Wednesday 22 nd November 2017 DR TANIA STATHAKI
Signals and Systems Lecture 11 Wednesday 22 nd November 2017 DR TANIA STATHAKI READER (ASSOCIATE PROFFESOR) IN SIGNAL PROCESSING IMPERIAL COLLEGE LONDON Effect on poles and zeros on frequency response
More information16.362: Signals and Systems: 1.0
16.362: Signals and Systems: 1.0 Prof. K. Chandra ECE, UMASS Lowell September 1, 2016 1 Background The pre-requisites for this course are Calculus II and Differential Equations. A basic understanding of
More informationChapter 4 The Fourier Series and Fourier Transform
Chapter 4 The Fourier Series and Fourier Transform Fourier Series Representation of Periodic Signals Let x(t) be a CT periodic signal with period T, i.e., xt ( + T) = xt ( ), t R Example: the rectangular
More informationLecture 28 Continuous-Time Fourier Transform 2
Lecture 28 Continuous-Time Fourier Transform 2 Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/6/14 1 Limit of the Fourier Series Rewrite (11.9) and (11.10) as As, the fundamental
More informationETSF15 Analog/Digital. Stefan Höst
ETSF15 Analog/Digital Stefan Höst Physical layer Analog vs digital Sampling, quantisation, reconstruction Modulation Represent digital data in a continuous world Disturbances Noise and distortion Synchronization
More informationChapter 12 Variable Phase Interpolation
Chapter 12 Variable Phase Interpolation Contents Slide 1 Reason for Variable Phase Interpolation Slide 2 Another Need for Interpolation Slide 3 Ideal Impulse Sampling Slide 4 The Sampling Theorem Slide
More informationFROM ANALOGUE TO DIGITAL
SIGNALS AND SYSTEMS: PAPER 3C1 HANDOUT 7. Dr David Corrigan 1. Electronic and Electrical Engineering Dept. corrigad@tcd.ie www.mee.tcd.ie/ corrigad FROM ANALOGUE TO DIGITAL To digitize signals it is necessary
More informationG52IVG, School of Computer Science, University of Nottingham
Image Transforms Fourier Transform Basic idea 1 Image Transforms Fourier transform theory Let f(x) be a continuous function of a real variable x. The Fourier transform of f(x) is F ( u) f ( x)exp[ j2πux]
More informationContents. 0.1 Notation... 3
Contents 0.1 Notation........................................ 3 1 A Short Course on Frame Theory 4 1.1 Examples of Signal Expansions............................ 4 1.2 Signal Expansions in Finite-Dimensional
More informationMultirate Digital Signal Processing
Multirate Digital Signal Processing Basic Sampling Rate Alteration Devices Up-sampler - Used to increase the sampling rate by an integer factor Down-sampler - Used to decrease the sampling rate by an integer
More information22 APPENDIX 1: MATH FACTS
22 APPENDIX : MATH FACTS 22. Vectors 22.. Definition A vector has a dual definition: It is a segment of a a line with direction, or it consists of its projection on a reference system xyz, usually orthogonal
More informationShift-Variance and Nonstationarity of Generalized Sampling-Reconstruction Processes
Shift-Variance and Nonstationarity of Generalized Sampling-Reconstruction Processes Runyi Yu Eastern Mediterranean University Gazimagusa, North Cyprus Web: faraday.ee.emu.edu.tr/yu Emails: runyi.yu@emu.edu.tr
More informationModule 4. Related web links and videos. 1. FT and ZT
Module 4 Laplace transforms, ROC, rational systems, Z transform, properties of LT and ZT, rational functions, system properties from ROC, inverse transforms Related web links and videos Sl no Web link
More informationMassachusetts Institute of Technology
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.11: Introduction to Communication, Control and Signal Processing QUIZ 1, March 16, 21 QUESTION BOOKLET
More informationFilter Analysis and Design
Filter Analysis and Design Butterworth Filters Butterworth filters have a transfer function whose squared magnitude has the form H a ( jω ) 2 = 1 ( ) 2n. 1+ ω / ω c * M. J. Roberts - All Rights Reserved
More information06/12/ rws/jMc- modif SuFY10 (MPF) - Textbook Section IX 1
IV. Continuous-Time Signals & LTI Systems [p. 3] Analog signal definition [p. 4] Periodic signal [p. 5] One-sided signal [p. 6] Finite length signal [p. 7] Impulse function [p. 9] Sampling property [p.11]
More informationNoise Reduction in Oversampled Filter Banks Using Predictive Quantization
IEEE TRANSACTIONS ON INFORMATION THEORY, VOL 47, NO 1, JANUARY 2001 155 Noise Reduction in Oversampled Filter Banks Using Predictive Quantization Helmut Bölcskei, Member, IEEE, Franz Hlawatsch, Member,
More informationECE 3084 QUIZ 2 SCHOOL OF ELECTRICAL AND COMPUTER ENGINEERING GEORGIA INSTITUTE OF TECHNOLOGY APRIL 2, Name:
ECE 3084 QUIZ 2 SCHOOL OF ELECTRICAL AND COMPUTER ENGINEERING GEORGIA INSTITUTE OF TECHNOLOGY APRIL 2, 205 Name:. The quiz is closed book, except for one 2-sided sheet of handwritten notes. 2. Turn off
More informationSpin Glass Approach to Restricted Isometry Constant
Spin Glass Approach to Restricted Isometry Constant Ayaka Sakata 1,Yoshiyuki Kabashima 2 1 Institute of Statistical Mathematics 2 Tokyo Institute of Technology 1/29 Outline Background: Compressed sensing
More informationLecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design.
ISS0031 Modeling and Identification Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design. Aleksei Tepljakov, Ph.D. September 30, 2015 Linear Dynamic Systems Definition
More informationMEDE2500 Tutorial Nov-7
(updated 2016-Nov-4,7:40pm) MEDE2500 (2016-2017) Tutorial 3 MEDE2500 Tutorial 3 2016-Nov-7 Content 1. The Dirac Delta Function, singularity functions, even and odd functions 2. The sampling process and
More informationLecture 3 - Design of Digital Filters
Lecture 3 - Design of Digital Filters 3.1 Simple filters In the previous lecture we considered the polynomial fit as a case example of designing a smoothing filter. The approximation to an ideal LPF can
More informationThe integrating factor method (Sect. 1.1)
The integrating factor method (Sect. 1.1) Overview of differential equations. Linear Ordinary Differential Equations. The integrating factor method. Constant coefficients. The Initial Value Problem. Overview
More informationQuadrature-Mirror Filter Bank
Quadrature-Mirror Filter Bank In many applications, a discrete-time signal x[n] is split into a number of subband signals { v k [ n]} by means of an analysis filter bank The subband signals are then processed
More informationECE 301. Division 2, Fall 2006 Instructor: Mimi Boutin Midterm Examination 3
ECE 30 Division 2, Fall 2006 Instructor: Mimi Boutin Midterm Examination 3 Instructions:. Wait for the BEGIN signal before opening this booklet. In the meantime, read the instructions below and fill out
More informationCOMPUTATION OF FOURIER TRANSFORMS FOR NOISY B FOR NOISY BANDLIMITED SIGNALS
COMPUTATION OF FOURIER TRANSFORMS FOR NOISY BANDLIMITED SIGNALS October 22, 2011 I. Introduction Definition of Fourier transform: F [f ](ω) := ˆf (ω) := + f (t)e iωt dt, ω R (1) I. Introduction Definition
More informationAdvanced Analog Building Blocks. Prof. Dr. Peter Fischer, Dr. Wei Shen, Dr. Albert Comerma, Dr. Johannes Schemmel, etc
Advanced Analog Building Blocks Prof. Dr. Peter Fischer, Dr. Wei Shen, Dr. Albert Comerma, Dr. Johannes Schemmel, etc 1 Topics 1. S domain and Laplace Transform Zeros and Poles 2. Basic and Advanced current
More informationLecture 15: Time and Frequency Joint Perspective
WAVELETS AND MULTIRATE DIGITAL SIGNAL PROCESSING Lecture 15: Time and Frequency Joint Perspective Prof.V.M.Gadre, EE, IIT Bombay Introduction In lecture 14, we studied steps required to design conjugate
More informationarxiv: v1 [cs.it] 26 Sep 2018
SAPLING THEORY FOR GRAPH SIGNALS ON PRODUCT GRAPHS Rohan A. Varma, Carnegie ellon University rohanv@andrew.cmu.edu Jelena Kovačević, NYU Tandon School of Engineering jelenak@nyu.edu arxiv:809.009v [cs.it]
More informationOn the Frequency-Domain Properties of Savitzky-Golay Filters
On the Frequency-Domain Properties of Savitzky-Golay Filters Ronald W Schafer HP Laboratories HPL-2-9 Keyword(s): Savitzky-Golay filter, least-squares polynomial approximation, smoothing Abstract: This
More informationFourier transform representation of CT aperiodic signals Section 4.1
Fourier transform representation of CT aperiodic signals Section 4. A large class of aperiodic CT signals can be represented by the CT Fourier transform (CTFT). The (CT) Fourier transform (or spectrum)
More informationLaplace Transforms and use in Automatic Control
Laplace Transforms and use in Automatic Control P.S. Gandhi Mechanical Engineering IIT Bombay Acknowledgements: P.Santosh Krishna, SYSCON Recap Fourier series Fourier transform: aperiodic Convolution integral
More informationWavelet Bi-frames with Uniform Symmetry for Curve Multiresolution Processing
Wavelet Bi-frames with Uniform Symmetry for Curve Multiresolution Processing Qingtang Jiang Abstract This paper is about the construction of univariate wavelet bi-frames with each framelet being symmetric.
More informationBridge between continuous time and discrete time signals
6 Sampling Bridge between continuous time and discrete time signals Sampling theorem complete representation of a continuous time signal by its samples Samplingandreconstruction implementcontinuous timesystems
More informationChapter XI Novanion rings
Chapter XI Novanion rings 11.1 Introduction. In this chapter we continue to provide general structures for theories in physics. J. F. Adams proved in 1960 that the only possible division algebras are at
More informationSignals & Systems. Chapter 7: Sampling. Adapted from: Lecture notes from MIT, Binghamton University, and Purdue. Dr. Hamid R.
Signals & Systems Chapter 7: Sampling Adapted from: Lecture notes from MIT, Binghamton University, and Purdue Dr. Hamid R. Rabiee Fall 2013 Outline 1. The Concept and Representation of Periodic Sampling
More informationLOPE3202: Communication Systems 10/18/2017 2
By Lecturer Ahmed Wael Academic Year 2017-2018 LOPE3202: Communication Systems 10/18/2017 We need tools to build any communication system. Mathematics is our premium tool to do work with signals and systems.
More information