Digital Signal Processing the Z-Transform

Total Page:16

File Type:pdf, Size:1020Kb

Digital Signal Processing the Z-Transform 2065-23 Advanced Training Course on FPGA Design and VHDL for Hardware Simulation and Synthesis 26 October - 20 November, 2009 Digital Signal Processing The z-transform Massimiliano Nolich DEEI Facolta' di Ingegneria Universita' degli Studi di Trieste via Valerio, 10, 34127 Trieste Italy The z-transform See Oppenheim and Schafer, Second Edition pages 94–139, or First Edition pages 149–201. 1 Introduction The z-transform of a sequence xŒn is 1 X.z/ D xŒnzn: nD1 X The z-transform can also be thought of as an operator Zfg that transforms a sequence to a function: 1 ZfxŒngD xŒnzn D X.z/: nD1 X In both cases z is a continuous complex variable. We may obtain the Fourier transform from the z-transform by making the substitution z D ej!. This corresponds to restricting jzjD 1. Also, with z D rej! , 1 1 X.rej! / D xŒn.rej! /n D .xŒnr n/ ej!n: nD1 nD1 X X That is, the z-transform is the Fourier transform of the sequence xŒnr n. For r D 1 this becomes the Fourier transform of xŒn. The Fourier transform therefore corresponds to the z-transform evaluated on the unit circle: 1 z−plane Im z D ej! ! Re Unit circle The inherent periodicity in frequency of the Fourier transform is captured naturally under this interpretation. The Fourier transform does not converge for all sequences — the infinite sum may not always be finite. Similarly, the z-transform does not converge for all sequences or for all values of z. The set of values of z for which the z-transform converges is called the region of convergence (ROC). 1 The Fourier transform of xŒn exists if the sum nD1 jxŒnj converges. However, the z-transform of xŒn is just the Fourier transform of the sequence P xŒnr n. The z-transform therefore exists (or converges) if 1 X.z/ D jxŒnr nj < 1: nD1 X This leads to the condition 1 jxŒnjjzjn < 1 nD1 X for the existence of the z-transform. The ROC therefore consists of a ring in the z-plane: 2 z−plane Im Region of convergence Re In specific cases the inner radius of this ring may include the origin, and the outer radius may extend to infinity. If the ROC includes the unit circle jzjD 1, then the Fourier transform will converge. Most useful z-transforms can be expressed in the form P.z/ X.z/ D ; Q.z/ where P.z/ and Q.z/ are polynomials in z. The values of z for which P.z/ D 0 are called the zeros of X.z/, and the values with Q.z/ D 0 are called the poles. The zeros and poles completely specify X.z/ to within a multiplicative constant. Example: right-sided exponential sequence Consider the signal xŒn D anuŒn. This has the z-transform 1 1 X.z/ D anuŒnzn D .az1/n: nD1 nD0 X X Convergence requires that 1 jaz1jn < 1; nD0 X which is only the case if jaz1j <1, or equivalently jzj > jaj. In the ROC, the 3 series converges to 1 1 z X.z/ D .az1/n D D ; jzj > jaj; 1 az1 z a nD0 X since it is just a geometric series. The z-transform has a region of convergence for any finite value of a. z−plane Im unit circle a Re 1 The Fourier transform of xŒn only exists if the ROC includes the unit circle, which requires that jaj <1. On the other hand, if jaj >1 then the ROC does not include the unit circle, and the Fourier transform does not exist. This is consistent with the fact that for these values of a the sequence anuŒn is exponentially growing, and the sum therefore does not converge. Example: left-sided exponential sequence Now consider the sequence xŒn DanuŒn 1. This sequence is left-sided because it is nonzero only for n 1. The z-transform is 1 1 X.z/ D anuŒn 1zn D anzn nD1 nD1 X1 1 X D anzn D 1 .a1z/n: nD1 nD0 X X 4 For ja1zj <1, or jzj < jaj, the series converges to 1 1 z X.z/ D 1 D D ; jzj < jaj: 1 a1z 1 az1 z a z−plane Im unit circle a Re 1 Note that the expression for the z-transform (and the pole zero plot) is exactly the same as for the right-handed exponential sequence — only the region of convergence is different. Specifying the ROC is therefore critical when dealing with the z-transform. Example: sum of two exponentials 1 n 1 n The signal xŒn D 2 uŒn C 3 uŒn is the sum of two real exponentials. The z-transform is 1 1 n 1 n X.z/ D uŒn C uŒn zn 2 3 nD1 X 1 1 n 1 1 n D uŒnzn C uŒnzn 2 3 nD1 nD1 X X 1 1 n 1 1 n D z1 C z1 : 2 3 nD0 nD0 X X From the example for the right-handed exponential sequence, the first term in this sum converges for jzj > 1=2, and the second for jzj > 1=3. The combined transform X.z/ therefore converges in the intersection of these regions, namely 5 when jzj > 1=2. In this case 1 1 2z.z 1 / X.z/ D C D 12 : 1 1 1 1 1 1 1 2 z 1 C 3 z .z 2 /.z C 3 / The pole-zero plot and region of convergence of the signal is z−plane Im unit circle 1 Re 1 1 1 3 12 2 Example: finite length sequence The signal an 0 n N 1 xŒn D 80 otherwise < has z-transform : N 1 N 1 X.z/ D anzn D .az1/n nD0 nD0 X X 1 .az1/N 1 zN aN D D : 1 az1 zN 1 z a Since there are only a finite number of nonzero terms the sum always converges when az1 is finite. There are no restrictions on a (jaj < 1), and the ROC is the entire z-plane with the exception of the origin z D 0 (where the terms in the sum are infinite). The N roots of the numerator polynomial are at j.2k=N / zk D ae ; k D 0;1;:::;N 1; 6 since these values satisfy the equation zN D aN . The zero at k D 0 cancels the pole at z D a, so there are no poles except at the origin, and the zeros are at j.2k=N / zk D ae ; k D 1;:::;N 1: 2 Properties of the region of convergence The properties of the ROC depend on the nature of the signal. Assuming that the signal has a finite amplitude and that the z-transform is a rational function: The ROC is a ring or disk in the z-plane, centered on the origin (0 rR < jzj < rL 1). The Fourier transform of xŒn converges absolutely if and only if the ROC of the z-transform includes the unit circle. The ROC cannot contain any poles. If xŒn is finite duration (ie. zero except on finite interval 1 < N1 n N2 < 1), then the ROC is the entire z-plane except perhaps at z D 0 or z D 1. If xŒn is a right-sided sequence then the ROC extends outward from the outermost finite pole to infinity. If xŒn is left-sided then the ROC extends inward from the innermost nonzero pole to z D 0. A two-sided sequence (neither left nor right-sided) has a ROC consisting of a ring in the z-plane, bounded on the interior and exterior by a pole (and not containing any poles). The ROC is a connected region. 7 3 The inverse z-transform Formally, the inverse z-transform can be performed by evaluating a Cauchy integral. However, for discrete LTI systems simpler methods are often sufficient. 3.1 Inspection method If one is familiar with (or has a table of) common z-transform pairs, the inverse can be found by inspection. For example, one can invert the z-transform 1 1 X.z/ D ; jzj > ; 1 1 1 2 z ! 2 using the z-transform pair Z 1 anuŒn ! ; for jzj > jaj: 1 az1 By inspection we recognise that 1 n xŒn D uŒn: 2 Also, if X.z/ is a sum of terms then one may be able to do a term-by-term inversion by inspection, yielding xŒn as a sum of terms. 3.2 Partial fraction expansion For any rational function we can obtain a partial fraction expansion, and identify the z-transform of each term. Assume that X.z/ is expressed as a ratio of polynomials in z1: M b zk X.z/ D kD0 k : N k PkD0 akz P 8 It is always possible to factor X.z/ as b M .1 c z1/ X.z/ D 0 kD1 k ; a N 1 0 QkD1.1 dkz / where the ck’s and dk’s are the nonzeroQ zeros and poles of X.z/. If M < N and the poles are all first order, then X.z/ can be expressed as N Ak X.z/ D 1 : 1 dkz kXD1 In this case the coefficients Ak are given by A D .1 d z1/X.z/ k k zDdk : ˇ If M N and the poles are all first order, thenˇ an expansion of the form M N N A X.z/ D B zr C k r 1 d z1 rD0 k X kXD1 can be used, and the Br ’s be obtained by long division of the numerator by the denominator.
Recommended publications
  • Detection and Estimation Theory Introduction to ECE 531 Mojtaba Soltanalian- UIC the Course
    Detection and Estimation Theory Introduction to ECE 531 Mojtaba Soltanalian- UIC The course Lectures are given Tuesdays and Thursdays, 2:00-3:15pm Office hours: Thursdays 3:45-5:00pm, SEO 1031 Instructor: Prof. Mojtaba Soltanalian office: SEO 1031 email: [email protected] web: http://msol.people.uic.edu/ The course Course webpage: http://msol.people.uic.edu/ECE531 Textbook(s): * Fundamentals of Statistical Signal Processing, Volume 1: Estimation Theory, by Steven M. Kay, Prentice Hall, 1993, and (possibly) * Fundamentals of Statistical Signal Processing, Volume 2: Detection Theory, by Steven M. Kay, Prentice Hall 1998, available in hard copy form at the UIC Bookstore. The course Style: /Graduate Course with Active Participation/ Introduction Let’s start with a radar example! Introduction> Radar Example QUIZ Introduction> Radar Example You can actually explain it in ten seconds! Introduction> Radar Example Applications in Transportation, Defense, Medical Imaging, Life Sciences, Weather Prediction, Tracking & Localization Introduction> Radar Example The strongest signals leaking off our planet are radar transmissions, not television or radio. The most powerful radars, such as the one mounted on the Arecibo telescope (used to study the ionosphere and map asteroids) could be detected with a similarly sized antenna at a distance of nearly 1,000 light-years. - Seth Shostak, SETI Introduction> Estimation Traditionally discussed in STATISTICS. Estimation in Signal Processing: Digital Computers ADC/DAC (Sampling) Signal/Information Processing Introduction> Estimation The primary focus is on obtaining optimal estimation algorithms that may be implemented on a digital computer. We will work on digital signals/datasets which are typically samples of a continuous-time waveform.
    [Show full text]
  • Z-Transformation - I
    Digital signal processing: Lecture 5 z-transformation - I Produced by Qiangfu Zhao (Since 1995), All rights reserved © DSP-Lec05/1 Review of last lecture • Fourier transform & inverse Fourier transform: – Time domain & Frequency domain representations • Understand the “true face” (latent factor) of some physical phenomenon. – Key points: • Definition of Fourier transformation. • Definition of inverse Fourier transformation. • Condition for a signal to have FT. Produced by Qiangfu Zhao (Since 1995), All rights reserved © DSP-Lec05/2 Review of last lecture • The sampling theorem – If the highest frequency component contained in an analog signal is f0, the sampling frequency (the Nyquist rate) should be 2xf0 at least. • Frequency response of an LTI system: – The Fourier transformation of the impulse response – Physical meaning: frequency components to keep (~1) and those to remove (~0). – Theoretic foundation: Convolution theorem. Produced by Qiangfu Zhao (Since 1995), All rights reserved © DSP-Lec05/3 Topics of this lecture Chapter 6 of the textbook • The z-transformation. • z-変換 • Convergence region of z- • z-変換の収束領域 transform. -変換とフーリエ変換 • Relation between z- • z transform and Fourier • z-変換と差分方程式 transform. • 伝達関数 or システム • Relation between z- 関数 transform and difference equation. • Transfer function (system function). Produced by Qiangfu Zhao (Since 1995), All rights reserved © DSP-Lec05/4 z-transformation(z-変換) • Laplace transformation is important for analyzing analog signal/systems. • z-transformation is important for analyzing discrete signal/systems. • For a given signal x(n), its z-transformation is defined by ∞ Z[x(n)] = X (z) = ∑ x(n)z−n (6.3) n=0 where z is a complex variable. Produced by Qiangfu Zhao (Since 1995), All rights reserved © DSP-Lec05/5 One-sided z-transform and two-sided z-transform • The z-transform defined above is called one- sided z-transform (片側z-変換).
    [Show full text]
  • Basics on Digital Signal Processing
    Basics on Digital Signal Processing z - transform - Digital Filters Vassilis Anastassopoulos Electronics Laboratory, Physics Department, University of Patras Outline of the Lecture 1. The z-transform 2. Properties 3. Examples 4. Digital filters 5. Design of IIR and FIR filters 6. Linear phase 2/31 z-Transform 2 N 1 N 1 j nk n X (k) x(n) e N X (z) x(n) z n0 n0 Time to frequency Time to z-domain (complex) Transformation tool is the Transformation tool is complex wave z e j e j j j2k / N e e With amplitude ρ changing(?) with time With amplitude |ejω|=1 The z-transform is more general than the DFT 3/31 z-Transform For z=ejω i.e. ρ=1 we work on the N 1 unit circle X (z) x(n) z n And the z-transform degenerates n0 into the Fourier transform. I -z z=R+jI |z|=1 R The DFT is an expression of the z-transform on the unit circle. The quantity X(z) must exist with finite value on the unit circle i.e. must posses spectrum with which we can describe a signal or a system. 4/31 z-Transform convergence We are interested in those values of z for which X(z) converges. This region should contain the unit circle. I Why is it so? |a| N 1 N 1 X (z) x(n) zn x(n) (1/zn ) n0 n0 R At z=0, X(z) diverges ROC The values of z for which X(z) diverges are called poles of X(z).
    [Show full text]
  • 3F3 – Digital Signal Processing (DSP)
    3F3 – Digital Signal Processing (DSP) Simon Godsill www-sigproc.eng.cam.ac.uk/~sjg/teaching/3F3 Course Overview • 11 Lectures • Topics: – Digital Signal Processing – DFT, FFT – Digital Filters – Filter Design – Filter Implementation – Random signals – Optimal Filtering – Signal Modelling • Books: – J.G. Proakis and D.G. Manolakis, Digital Signal Processing 3rd edition, Prentice-Hall. – Statistical digital signal processing and modeling -Monson H. Hayes –Wiley • Some material adapted from courses by Dr. Malcolm Macleod, Prof. Peter Rayner and Dr. Arnaud Doucet Digital Signal Processing - Introduction • Digital signal processing (DSP) is the generic term for techniques such as filtering or spectrum analysis applied to digitally sampled signals. • Recall from 1B Signal and Data Analysis that the procedure is as shown below: • is the sampling period • is the sampling frequency • Recall also that low-pass anti-aliasing filters must be applied before A/D and D/A conversion in order to remove distortion from frequency components higher than Hz (see later for revision of this). • Digital signals are signals which are sampled in time (“discrete time”) and quantised. • Mathematical analysis of inherently digital signals (e.g. sunspot data, tide data) was developed by Gauss (1800), Schuster (1896) and many others since. • Electronic digital signal processing (DSP) was first extensively applied in geophysics (for oil-exploration) then military applications, and is now fundamental to communications, mobile devices, broadcasting, and most applications of signal and image processing. There are many advantages in carrying out digital rather than analogue processing; among these are flexibility and repeatability. The flexibility stems from the fact that system parameters are simply numbers stored in the processor.
    [Show full text]
  • CONVOLUTION: Digital Signal Processing .R .Hamming, W
    CONVOLUTION: Digital Signal Processing AN-237 National Semiconductor CONVOLUTION: Digital Application Note 237 Signal Processing January 1980 Introduction As digital signal processing continues to emerge as a major Decreasing the pulse width while increasing the pulse discipline in the field of electrical engineering, an even height to allow the area under the pulse to remain constant, greater demand has evolved to understand the basic theo- Figure 1c, shows from eq(1) and eq(2) the bandwidth or retical concepts involved in the development of varied and spectral-frequency content of the pulse to have increased, diverse signal processing systems. The most fundamental Figure 1d. concepts employed are (not necessarily listed in the order Further altering the pulse to that of Figure 1e provides for an [ ] of importance) the sampling theorem 1 , Fourier transforms even broader bandwidth, Figure 1f. If the pulse is finally al- [ ][] 2 3, convolution, covariance, etc. tered to the limit, i.e., the pulsewidth being infinitely narrow The intent of this article will be to address the concept of and its amplitude adjusted to still maintain an area of unity convolution and to present it in an introductory manner under the pulse, it is found in 1g and 1h the unit impulse hopefully easily understood by those entering the field of produces a constant, or ``flat'' spectrum equal to 1 at all digital signal processing. frequencies. Note that if ATe1 (unit area), we get, by defini- It may be appropriate to note that this article is Part II (Part I tion, the unit impulse function in time.
    [Show full text]
  • The Scientist and Engineer's Guide to Digital Signal Processing Properties of Convolution
    CHAPTER 7 Properties of Convolution A linear system's characteristics are completely specified by the system's impulse response, as governed by the mathematics of convolution. This is the basis of many signal processing techniques. For example: Digital filters are created by designing an appropriate impulse response. Enemy aircraft are detected with radar by analyzing a measured impulse response. Echo suppression in long distance telephone calls is accomplished by creating an impulse response that counteracts the impulse response of the reverberation. The list goes on and on. This chapter expands on the properties and usage of convolution in several areas. First, several common impulse responses are discussed. Second, methods are presented for dealing with cascade and parallel combinations of linear systems. Third, the technique of correlation is introduced. Fourth, a nasty problem with convolution is examined, the computation time can be unacceptably long using conventional algorithms and computers. Common Impulse Responses Delta Function The simplest impulse response is nothing more that a delta function, as shown in Fig. 7-1a. That is, an impulse on the input produces an identical impulse on the output. This means that all signals are passed through the system without change. Convolving any signal with a delta function results in exactly the same signal. Mathematically, this is written: EQUATION 7-1 The delta function is the identity for ( ' convolution. Any signal convolved with x[n] *[n] x[n] a delta function is left unchanged. This property makes the delta function the identity for convolution. This is analogous to zero being the identity for addition (a%0 ' a), and one being the identity for multiplication (a×1 ' a).
    [Show full text]
  • Exploring Anti-Aliasing Filters in Signal Conditioners for Mixed-Signal, Multimodal Sensor Conditioning by Arun T
    Texas Instruments Incorporated Amplifiers: Op Amps Exploring anti-aliasing filters in signal conditioners for mixed-signal, multimodal sensor conditioning By Arun T. Vemuri Systems Architect, Enhanced Industrial Introduction Figure 1. Multimodal, mixed-signal sensor-signal conditioner Some sensor-signal conditioners are used to process the output Analog Digital of multiple sense elements. This Domain Domain processing is often provided by multimodal, mixed-signal condi- tioners that can handle the out- puts from several sense elements Sense Digital Amplifier 1 ADC 1 at the same time. This article Element 1 Filter 1 analyzes the operation of anti- Processed aliasing filters in such sensor- Intelligent Output Compensation signal conditioners. Sense Digital Basics of sensor-signal Amplifier 2 ADC 2 Element 2 Filter 2 conditioners Sense elements, or transducers, convert a physical quantity of interest into electrical signals. Examples include piezo resistive bridges used to measure pres- sure, piezoelectric transducers used to detect ultrasonic than one sense element is processed by the same signal waves, and electrochemical cells used to measure gas conditioner is called multimodal signal conditioning. concentrations. The electrical signals produced by sense Mixed-signal signal conditioning elements are small and exhibit nonidealities, such as tem- Another aspect of sensor-signal conditioning is the electri- perature drifts and nonlinear transfer functions. cal domain in which the signal conditioning occurs. TI’s Sensor analog front ends such as the Texas Instruments PGA309 is an example of a device where the signal condi- (TI) LMP91000 and sensor-signal conditioners such as TI’s tioning of resistive-bridge sense elements occurs in the PGA400/450 are used to amplify the small signals produced analog domain.
    [Show full text]
  • Digital Signal Processing: Theory and Practice, Hardware and Software
    AC 2009-959: DIGITAL SIGNAL PROCESSING: THEORY AND PRACTICE, HARDWARE AND SOFTWARE Wei PAN, Idaho State University Wei Pan is Assistant Professor and Director of VLSI Laboratory, Electrical Engineering Department, Idaho State University. She has several years of industrial experience including Siemens (project engineering/management.) Dr. Pan is an active member of ASEE and IEEE and serves on the membership committee of the IEEE Education Society. S. Hossein Mousavinezhad, Idaho State University S. Hossein Mousavinezhad is Professor and Chair, Electrical Engineering Department, Idaho State University. Dr. Mousavinezhad is active in ASEE and IEEE and is an ABET program evaluator. Hossein is the founding general chair of the IEEE International Conferences on Electro Information Technology. Kenyon Hart, Idaho State University Kenyon Hart is Specialist Engineer and Associate Lecturer, Electrical Engineering Department, Idaho State University, Pocatello, Idaho. Page 14.491.1 Page © American Society for Engineering Education, 2009 Digital Signal Processing, Theory/Practice, HW/SW Abstract Digital Signal Processing (DSP) is a course offered by many Electrical and Computer Engineering (ECE) programs. In our school we offer a senior-level, first-year graduate course with both lecture and laboratory sections. Our experience has shown that some students consider the subject matter to be too theoretical, relying heavily on mathematical concepts and abstraction. There are several visible applications of DSP including: cellular communication systems, digital image processing and biomedical signal processing. Authors have incorporated many examples utilizing software packages including MATLAB/MATHCAD in the course and also used classroom demonstrations to help students visualize some difficult (but important) concepts such as digital filters and their design, various signal transformations, convolution, difference equations modeling, signals/systems classifications and power spectral estimation as well as optimal filters.
    [Show full text]
  • ELEG 5173L Digital Signal Processing Ch. 3 Discrete-Time Fourier Transform
    Department of Electrical Engineering University of Arkansas ELEG 5173L Digital Signal Processing Ch. 3 Discrete-Time Fourier Transform Dr. Jingxian Wu [email protected] 2 OUTLINE • The Discrete-Time Fourier Transform (DTFT) • Properties • DTFT of Sampled Signals • Upsampling and downsampling 3 DTFT • Discrete-time Fourier Transform (DTFT) X () x(n)e jn n – (radians): digital frequency • Review: Z-transform: X (z) x(n)zn n0 j X () X (z) | j – Replace z with e . ze • Review: Fourier transform: X () x(t)e jt – (rads/sec): analog frequency 4 DTFT • Relationship between DTFT and Fourier Transform – Sample a continuous time signal x a ( t ) with a sampling period T xs (t) xa (t) (t nT ) xa (nT ) (t nT ) n n – The Fourier Transform of ys (t) jt jnT X s () xs (t)e dt xa (nT)e n – Define: T • : digital frequency (unit: radians) • : analog frequency (unit: radians/sec) – Let x(n) xa (nT) X () X s T 5 DTFT • Relationship between DTFT and Fourier Transform (Cont’d) – The DTFT can be considered as the scaled version of the Fourier transform of the sampled continuous-time signal jt jnT X s () xs (t)e dt xa (nT)e n x(n) x (nT) T a jn X () X s x(n)e T n 6 DTFT • Discrete Frequency – Unit: radians (the unit of continuous frequency is radians/sec) – X ( ) is a periodic function with period 2 j2 n jn j2n jn X ( 2 ) x(n)e x(n)e e x(n)e X () n n n – We only need to consider for • For Fourier transform, we need to consider 1 – f T 2 2T 1 – f T 2 2T 7 DTFT • Example: find the DTFT of the following signal – 1.
    [Show full text]
  • Lecture 12: Sampling, Aliasing, and the Discrete Fourier Transform Foundations of Digital Signal Processing
    Lecture 12: Sampling, Aliasing, and the Discrete Fourier Transform Foundations of Digital Signal Processing Outline • Review of Sampling • The Nyquist-Shannon Sampling Theorem • Continuous-time Reconstruction / Interpolation • Aliasing and anti-Aliasing • Deriving Transforms from the Fourier Transform • Discrete-time Fourier Transform, Fourier Series, Discrete-time Fourier Series • The Discrete Fourier Transform Foundations of Digital Signal Processing Lecture 12: Sampling, Aliasing, and the Discrete Fourier Transform 1 News Homework #5 . Due this week . Submit via canvas Coding Problem #4 . Due this week . Submit via canvas Foundations of Digital Signal Processing Lecture 12: Sampling, Aliasing, and the Discrete Fourier Transform 2 Exam 1 Grades The class did exceedingly well . Mean: 89.3 . Median: 91.5 Foundations of Digital Signal Processing Lecture 12: Sampling, Aliasing, and the Discrete Fourier Transform 3 Lecture 12: Sampling, Aliasing, and the Discrete Fourier Transform Foundations of Digital Signal Processing Outline • Review of Sampling • The Nyquist-Shannon Sampling Theorem • Continuous-time Reconstruction / Interpolation • Aliasing and anti-Aliasing • Deriving Transforms from the Fourier Transform • Discrete-time Fourier Transform, Fourier Series, Discrete-time Fourier Series • The Discrete Fourier Transform Foundations of Digital Signal Processing Lecture 12: Sampling, Aliasing, and the Discrete Fourier Transform 4 Sampling Discrete-Time Fourier Transform Foundations of Digital Signal Processing Lecture 12: Sampling, Aliasing,
    [Show full text]
  • MULTIRATE SIGNAL PROCESSING Multirate Signal Processing
    MULTIRATE SIGNAL PROCESSING Multirate Signal Processing • Definition: Signal processing which uses more than one sampling rate to perform operations • Upsampling increases the sampling rate • Downsampling reduces the sampling rate • Reference: Digital Signal Processing, DeFatta, Lucas, and Hodgkiss B. Baas, EEC 281 431 Multirate Signal Processing • Advantages of lower sample rates –May require less processing –Likely to reduce power dissipation, P = CV2 f, where f is frequently directly proportional to the sample rate –Likely to require less storage • Advantages of higher sample rates –May simplify computation –May simplify surrounding analog and RF circuitry • Remember that information up to a frequency f requires a sampling rate of at least 2f. This is the Nyquist sampling rate. –Or we can equivalently say the Nyquist sampling rate is ½ the sampling frequency, fs B. Baas, EEC 281 432 Upsampling Upsampling or Interpolation •For an upsampling by a factor of I, add I‐1 zeros between samples in the original sequence •An upsampling by a factor I is commonly written I For example, upsampling by two: 2 • Obviously the number of samples will approximately double after 2 •Note that if the sampling frequency doubles after an upsampling by two, that t the original sample sequence will occur at the same points in time t B. Baas, EEC 281 434 Original Signal Spectrum •Example signal with most energy near DC •Notice 5 spectral “bumps” between large signal “bumps” B. Baas, EEC 281 π 2π435 Upsampled Signal (Time) •One zero is inserted between the original samples for 2x upsampling B. Baas, EEC 281 436 Upsampled Signal Spectrum (Frequency) • Spectrum of 2x upsampled signal •Notice the location of the (now somewhat compressed) five “bumps” on each side of π B.
    [Show full text]
  • ESE 531: Digital Signal Processing Lecture Outline ADC Anti-Aliasing
    Lecture Outline ESE 531: Digital Signal Processing ! Data Converters " Anti-aliasing " ADC Lec 12: February 21st, 2017 " Quantization Data Converters, Noise Shaping (con’t) " Practical DAC ! Noise Shaping Penn ESE 531 Spring 2017 - Khanna Penn ESE 531 Spring 2017 - Khanna 2 Anti-Aliasing Filter with ADC ADC Penn ESE 531 Spring 2017 - Khanna 3 Penn ESE 531 Spring 2017 - Khanna 4 Oversampled ADC Oversampled ADC Penn ESE 531 Spring 2017 - Khanna 5 Penn ESE 531 Spring 2017 - Khanna 6 1 Oversampled ADC Oversampled ADC Penn ESE 531 Spring 2017 - Khanna 7 Penn ESE 531 Spring 2017 - Khanna 8 Sampling and Quantization Sampling and Quantization Penn ESE 531 Spring 2017 - Khanna 9 Penn ESE 531 Spring 2017 - Khanna 10 Effect of Quantization Error on Signal Quantization Error ! Quantization error is a deterministic function of the signal " Consequently, the effect of quantization strongly depends on the ! signal itself Model quantization error as noise ! Unless, we consider fairly trivial signals, a deterministic analysis is usually impractical " More common to look at errors from a statistical perspective " "Quantization noise” ! In that case: ! Two aspects " How much noise power (variance) does quantization add to our samples? " How is this noise distributed in frequency? Penn ESE 531 Spring 2017 - Khanna 11 Penn ESE 531 Spring 2017 - Khanna 12 2 Ideal Quantizer Ideal B-bit Quantizer ! Practical quantizers have a limited input range and a finite set of output codes ! Δ Quantization step ! E.g. a 3-bit quantizer can map onto 3 ! Quantization error
    [Show full text]