DIGITAL SIGNAL PROCESSING UNIT III INFINITE IMPULSE RESPONSE DIGITAL FILTERS. 3.6 Design of Digital Filter using Digital to Digital
|
|
- Melvyn Summers
- 5 years ago
- Views:
Transcription
1 DIGITAL SIGNAL PROCESSING UNIT III INFINITE IMPULSE RESPONSE DIGITAL FILTERS Contents: 3.1 Introduction IIR Filters 3.2 Transformation Function Derivation 3.3 Review of Analog IIR Filters Butterworth Filter Chebyshev Filter 3.4 Design of IIR Digital Filters Impulse Invariant Technique Bilinear Transformation 3.5 Pre Warping 3.6 Design of Digital Filter using Digital to Digital Transformation 3.7 Realization of IIR Filters 3.1 Introduction: IIR FILTERS Infinite impulse response (IIR) is a property of signal processing systems. Systems with this property are known as IIR systems or, when dealing with electronic filter systems, as IIR filters. IIR systems have an impulse response function that is non-zero over an infinite length of time. This is in contrast to finite impulse response filters (FIR), which have fixed-duration impulse responses. The simplest analog IIR filter is an RC filter made up of a single resistor (R) feeding into 1
2 a node shared with a single capacitor (C). This filter has an exponential impulse response characterized by an RC time constant. IIR filters may be implemented as either analog or digital filters. In digital IIR filters, the output feedback is immediately apparent in the equations defining the output. Design of digital IIR filters is heavily dependent on that of their analog counterparts because there are plenty of resources, works and straightforward design methods concerning analog feedback filter design while there are hardly any for digital IIR filters. As a result, usually, when a digital IIR filter is going to be implemented, an analog filter (e.g. Chebyshev filter, Butterworth filter, Elliptic filter) is first designed and then is converted to a digital filter by applying discretization techniques such as Bilinear transform or Impulse invariance. Example for IIR filters: Chebyshev filter, Butterworth filter, and the Bessel filter. An infinite impulse response filter (IIR) is one whose impulse response goes on forever. Recursive filters are closely associated with IIR filters but are not exclusively IIR. The transfer function of an IIR filter is shown below. For it to be IIR, at least one of the a's in the denominator must be non zero. M, the number of zeros should be less than N, the number of poles. To obtain the frequency response of the filter, substitute the following relation between z and omega. T is the sampling time. This results in a Frequency response with the general form. 2
3 To obtain the time sequence of the output or the impulse response use the following equivalent difference equation. The Linear Predictive model is a special case of an IIR filter and shown in figure The general IIR filter is given by: If p = 0 then the system represents a finite impulse response (FIR) filter. If p is not zero, then the system is an infinite impulse response (IIR) filter. Review Questions: 1. What is an Infinite Impulse Response filter? 3
4 2. Give examples for IIR filter. 3. Define Transfer function of IIR filter. 3.2 TRANSFER FUNCTION DERIVATION Digitals filters are often described and implemented in terms of the difference equation that defines how the output signal is related to the input signal: where: is the feedforward filter order are the feedforward filter coefficients is the feedback filter order are the feedback filter coefficients is the input signal is the output signal. A more condensed form of the difference equation is: which, when rearranged, becomes: 4
5 To find the transfer function of the filter, we first take the Z-transform of each side of the above equation, where we use the time-shift property to obtain: We define the transfer function to be: Considering that in most IIR filter designs coefficient is 1, the IIR filter transfer function takes the more traditional form: Stability The transfer function allows us to judge whether or not a system is bounded-input, bounded-output (BIBO) stable. To be specific, the BIBO stability criteria requires the ROC of the system include the unit circle. For example, for a causal system, all poles of the transfer function have to have an absolute value smaller than one. In other words, all poles must be located within a unit circle in the z-plane. The poles are defined as the values of z which make the denominator of H(z) equal to 0: 3.3 Review of Analog IIR Filters 3.3.1BUTTERWORTH FILTER: The Butterworth filter is one type of electronic filter design. It is designed to have a frequency response which is as flat as mathematically possible in the passband. Another name for it is maximally flat magnitude filter. 5
6 The transfer function Plot of the gain of Butterworth low-pass filters of orders 1 through 5. Note that the slope is 20n db/decade where n is the filter order. Like all filters, the typical prototype is the low-pass filter, which can be modified into a high-pass filter, or placed in series with others to form band-pass and bandstop filters, and higher order versions of these. The gain G(ω) of an n-order Butterworth low pass filter is given in terms of the transfer function H(s) as: where n = order of filter ω c = cutoff frequency (approximately the -3dB frequency) G 0 is the DC gain (gain at zero frequency) It can be seen that as n approaches infinity, the gain becomes a rectangle function and frequencies below ω c will be passed with gain G 0, while frequencies above ω c will be suppressed. For smaller values of n, the cutoff will be less sharp. We wish to determine the transfer function H(s) where s = σ + jω. Since H(s)H(-s) evaluated at s = jω is simply equal to H(jω) 2, it follows that: The poles of this expression occur on a circle of radius ω c at equally spaced points. The transfer function itself will be specified by just the poles in the negative real half-plane of s. The k-th pole is specified by: 6
7 The transfer function may be written in terms of these poles as: The denominator is a Butterworth polynomial in s. Normalized Butterworth polynomials The Butterworth polynomials may be written in complex form as above, but are usually written with real coefficients by multiplying pole pairs which are complex conjugates, such as s 1 and s n. The polynomials are normalized by setting ω c = 1. The normalized Butterworth polynomials then have the general form: 7
8 Digital implementation Digital implementations of Butterworth filters often use bilinear transform or matched z-transform to discretize an analog filter. For higher orders, they are sensitive to quantization errors. Comparison with other linear filters Here is an image showing the gain of a discrete-time Butterworth filter next to other common filter types. All of these filters are fifth-order. The Butterworth filter rolls off more slowly around the cutoff frequency than the others, but shows no ripples CHEBYSHEV FILTER Chebyshev filters are analog or digital filters having a steeper roll-off and more passband ripple (type I) or stopband ripple (type II) than Butterworth filters. Chebyshev filters have the property that they minimize the error between the idealized filter characteristic and the actual over the range of the filter, but with ripples in the passband. This type of filter is named in honor of Pafnuty Chebyshev 8
9 because their mathematical characteristics are derived from Chebyshev polynomials. Because of the passband ripple inherent in Chebyshev filters, filters which have a smoother response in the passband but a more irregular response in the stopband are preferred for some applications. Type I Chebyshev filters The frequency response of a fourth-order type I Chebyshev low-pass filter with These are the most common Chebyshev filters. The gain (or amplitude) response as a function of angular frequency ω of the nth order low pass filter is where ε is the ripple factor, ω 0 is the cutoff frequency and T n () is a Chebyshev polynomial of the nth order. The passband exhibits equiripple behavior, with the ripple determined by the ripple factor ε. In the passband, the Chebyshev polynomial alternates between 0 and 1 so the filter gain will alternate between maxima at G = 1 and minima at. At the cutoff frequency ω 0 the gain again has the value but continues to drop into the stop band as the frequency increases. The ripple is often given in db: Ripple in db = so that a ripple amplitude of 3 db results from Type II Chebyshev filters The frequency response of a fifth-order type II Chebyshev low-pass filter with 9
10 Also known as inverse Chebyshev, this type is less common because it does not roll off as fast as type I, and requires more components. It has no ripple in the passband, but does have equiripple in the stopband. The gain is: In the stop band, the Chebyshev polynomial will oscillate between 0 and 1 so that the gain will oscillate between zero and and the smallest frequency at which this maximum is attained will be the cutoff frequency ω 0. The parameter ε is thus related to the stopband attenuation γ in decibels by: Review Questions: 1. What is an analog Butterworth filter? 2. Define the transfer function of Butterworth filter in terms of its gain. 3. Mention the types of Chebyshev filter. 4. Compare the gain of Butterworth filter with other types. 3.4 DESIGN OF IIR DIGITAL FILTERS: Methods: i) Impulse Invariant Technique ii) Bilinear Transformation 10
11 3.4.1 i) IMPULSE INVARIANT TECHNIQUE Output from a digital filter is made up from previous inputs and previous outputs, using the operation of convolution: Two convolutions are involved: one with the previous inputs, and one with the previous outputs. In each case the convolving function is called the filter coefficients. If such a filter is subjected to an impulse (a signal consisting of one value followed by zeroes) then its output need not necessarily become zero after the impulse has run through the summation. So the impulse response of such a filter can be infinite in duration. Such a filter is called an Infinite Impulse Response filter or IIR filter. Note that the impulse response need not necessarily be infinite: if it were, the filter would be unstable. In fact for most practical filters, the impulse response will die away to a negligibly small level. One might argue that mathematically the response can go on for ever, getting smaller and smaller: but in a digital world once a level gets below one bit it might as well be zero. The Infinite Impulse Response refers to the ability of the filter to have an infinite impulse response and does not imply that it necessarily will have one: it serves as a warning that this type of filter is prone to feedback and instability. 11
12 The filter can be drawn as a block diagram: The filter diagram can show what hardware elements will be required when implementing the filter: The left hand side of the diagram shows the direct path, involving previous inputs: the right hand side shows the feedback path, operating upon previous outputs. If the phase response is not specified, one prefers to use IIR digital filter. In case of an IIR filter design, the most common practice is to convert the digital filter specifications to analog low pass prototype filter specifications, to determine the analog low pass transfer function meeting these specifications, and then to transform it into desired digital filter transfer function. This methods is used for the following reasons: 1. Analog filter approximation techniques are highly advanced. 2. They usually yield closed form solutions. 3. Extensive tables are available for analog-design. 4. Many applications require the digital solutions of analog filters. 12
13 The transformations generally have two properties (1) the imaginary axis of the s-plane maps into unit circle of the z-plane and (2) a stable continuous time filter is transformed to a stable discrete time filter. Filter design by impulse invariance Here is a recipe for designing an IIR digital filter: decide upon the desired frequency response design an appropriate analogue filter calculate the impulse response of this analogue filter sample the analogue filter's impulse response use the result as the filter coefficients This process is called the method of impulse invariance. The method of impulse invariance seems simple: but it is complicated by all the problems inherent in dealing with sampled data systems. In particular the method is subject to problems of aliasing and frequency resolution. In the impulse variance design procedure the impulse response of the impulse response of the discrete time system is proportional to equally spaced samples of the continues time filter, i.e., where T d represents a sampling interval, since the specifications of the filter are given in discrete time domain, it turns out that T d has no role to play in design of 13
14 the filter. From the sampling theorem we know that the frequency response of the discrete time filter is given by Since any practical continuous time filter is not strictly bandlimited there issome aliasing. However, if the continuous time filter approaches zero at high frequencies, the aliasing may be negligible. Then the frequency response of the discrete time filter is We first convert digital filter specifications to continuous time filter specifications. Neglecting aliasing, we get specification by applying the relation where is transferred to the designed filter H(z), we again use equation (9.2) and the parameter T d cancels out. Let us assume that the poles of the continuous time filter are simple, then The corresponding impulse response is Then The system function for this is 14
15 We see that a pole at in the s-plane is transformed to a pole at T d in the z-plane. If the continuous time filter is stable, that is, then the magnitude of will be less than 1, so the pole will be inside unit circle. Thus the causal discrete time filter is stable. The mapping of zeros is not so straight forward II) BILINEAR TRANSFORMATION The method of filter design by impulse invariance suffers from aliasing. The aliasing will be a problem if the analogue filter prototype's frequency response has significant components at or beyond the Nyquist frequency. The problem of aliasing arises because we wrap an infinitely long, straight frequency axis around a circle. So the frequency axis wraps around and around, and any components above the Nyquist frequency get wrapped back on top of other components. The bilinear transform is a method of squashing the infinite, straight analogue frequency axis so that it becomes finite. To avoid squashing the filter's desired frequency response too much, the bilinear transform squashes the far ends of the frequency axis the most - leaving the middle portion relatively unsquashed: 15
16 The infinite, straight analogue frequency axis is squashed so that it becomes finite - in fact just long enough to wrap around the unit circle once only. This is also sometimes called frequency warping If is the continues time transfer function the discrete time transfer function is detained by replacing s with Rearranging terms in equation (9.3) we obtain. Substituting, we get If, it is then magnitude of the real part in denominator is more than that of the numerator and so. Similarly if, than for all. Thus poles in the left half of the s-plane will get mapped to the poles inside the unit circle in z-plane. If then So,, writing we get rearranging we get 16
17 or The compression of frequency axis represented by (9.5) is nonlinear. This is illustrated in figure
18 Because of the nonlinear compression of the frequency axis, there is considerable phase distortion in the bilinear transformation. Review Questions : 1. What is meant by impulse invariant technique? 2. List out the steps to design a digital IIR filter using impulse invariant technique. 3. What is meant by Bilinear Transformation? 4. List out the steps to design a digital IIR filter using Bilinear Transformation. 3.5 PRE WARPING Frequency warping does change the shape of the desired filter frequency response. In particular, it changes the shape of the transition bands. This is a pity, since we went to a lot of trouble designing an analogue filter prototype that gave us the desired frequency response and transition band shapes. One way around this is to warp the analogue filter design before transforming it to the sampled data z plane Argand diagram: this warping being designed so that it will be exactly undone by the frequency warping later on. This is called prewarping. This technique avoids the problem of aliasing by mapping one revaluation of the unit circle in the z-plane. axis in the s-plane to Some frequently used analog filters In the previous two examples we have used Butterworth filter. The Butterworth filter of order n is described by the magnitude square frequency response of It has the following properties 18
19 is monotonically decreasing function of As n gets larger, approaches an ideal low pass filter are zero at is called maximally flat at origin, since all order derivative exist and they The poles of a Butterworth filter lie on circle of radius in s-plane. There are two types of Chebyshev filters, one containing ripples in the passband (type I) and the other containing a ripple in the stopband (type II). A Type I low pass normalizer Chebyshev filter has the magnitude squared frequency response. where is n th order Chebyshev polynomial. We have the relationship with Chebyshev filters have the following properties 1. The magnitude squared frequency response oscillates between 1 and within the passband, the so called equiripple and has a value of at, the normalized cut off frequency. 2. The magnitude response is monotonic outside the passband including transitionand stopband. 3. The poles of the Chebysher filter lie on an ellipse in s-plane. An elliptic filter has ripples both in passband and in stopband. The square magnitude frequency response is given by 19
20 where is Chebyshev rational function of O determined from specified ripple characteristics. An n th order Chebyshev filter has sharper cutoff than a Butterworth filter, that is, has a narrower transition bandwidth. Elliptic filter provides the smallest transition width. 3.6 DESIGN OF DIGITAL FILTER USING DIGITAL TO DIGITAL TRANSFORMATION There exists a set of transformation that takes a low pass digital filter and turn into highpass, bandpass, bandstop or another lowpass digital filter. These transformations are given in table 9.1. The transformations all take the form of replacing the in by some function of. Type From To Transformation Design Formula Low pass cutoff Low pass cutoff LPF HPF LPF BPF LPF BSF 20
21 Starting with a set of digital specifications and using the inverse of the design equation given in table 9.1, a set of lowpass digital requirements can be established. A LPF digital prototype filter is then selected to satisfy these requirements and the proper digital to digital transformation is applied to give the desired. Review Questions: 1. What is aliasing? 2. What do you mean by prewarping? 3. How a digital LPF can be transformed into HPF,BPF or BSF? 3.7 REALIZATION OF IIR FILTERS So far, we have studied the IIR filters by analysis of transfer functions or impulse responses. In this section we want to face the problem of implementing these filters as computational structures that can be directly coded using sound processing languages or real-time sound processing environments. i) Direct Form 1 ii) Direct Form 2 or Transposed form iii) Parallel Form iv) Cascade Form i) Direct form 1 So called because it can be drawn direct from the filter equation. Consider a second-order filter with two poles and two zeros, which is rep- Resented by the transfer function { δ(n) =1 ; n=0 }with N = M = 2. This can be realized by the signal flowgraph of fig., where the nodes having converging edges are considered as points of addition, and the nodes having diverging edges are considered as branching points. Such a realization is called Direct Form I. fig.1 21
22 Signal flowgraphs can be manipulated in several ways, thus leading to alternative realizations having different numerical properties and, possibly, more computationally efficient. For instance, if we want to implement a filter as a cascade of second-order cells such as that of fig.1 we can share, between two contiguous cells, the unit delays that are on the output stage of the first cell, with the unit delays that are on the input stage of the second cell, thus saving a number of memory accesses. This particular diagram is called the direct form 1 because the diagram can be drawn directly from the filter equation. The filter diagram can show what hardware elements will be required when implementing the filter: The left hand side of the diagram shows the direct path, involving previous inputs: the right hand side shows the feedback path, operating upon previous outputs. ii) Direct form 2 so called because it can be derived by changing the diagram of direct form 1 22
23 A first transformation comes from the observation that the structure of fig.1 is formed by the cascade of two blocks, each being linear and time invariant. Therefore, the two blocks can be commuted without altering the input-output behavior. Moreover, from the block exchange we get a flowgraph with two sideto-side stages of pure delays, and these stages can be combined in one only. The realization of these transformations is shown in fig. and it is called Direct Form II. fig.2 The filter diagram for direct form 1 can be drawn direct from the filter equation: The block diagram is in two halves: and since the results from each half are simply added together it does not matter in which order they are calculated. So the order of the halves can be swapped: 23
24 Now, note that the result after each delay is the same for both branches. So the delays down the centre can be combined: This is called direct form 2. Its advantage is that it needs less delay elements. And since delay elements require hardware (for example, processor registers) the direct form 2 requires less hardware and so is more efficient than direct form I. direct form 2 is also called canonic, which simply means 'having the minimum number of delay elements'. The transposition theorem says that if we take a filter diagram and reverse all the elements - swapping the order of execution for every element, and reversing the direction of all the flow arrows - then the result is the same: if everything is turned back to front, it all works just the same This means that the direct form 1 diagram can be obtained by transposition of the direct form 2 diagram: 24
25 For this reason, direct form 1 is often called transposed direct form 2. Transposed so called because it is obtained by transposition of direct form 2 - but really, this is just direct form 1 Another transformation that can be done on a signal flow graph without altering its input-output behavior is the transposition The transposition of a signal flow graph is done with the following operations: Inversion of the direction of all the edges Transformation of the nodes of addition into branching nodes, and vice versa The transposition of a realization in Direct Form II leads to the Transposed Form II, which is shown in fig.2 Similarly, the Transposed Form I is obtained by transposition of the Direct Form I. fig.3 25
26 Canonic so called because it has the minimum number of delay elements - but really, it is just direct form 2 iii) Parallel Form: Because IIR filters are very sensitive to quantisation errors, they are usually implemented as second order sections. The parallel form is simple: The outputs from each second order section are simply added together. If scaling is required, this is done separately for each section. It is possible to scale each section appropriately, and by a different scale factor, to minimise quantisation error. In this case another extra multiplier is required for each section, to scale the individual section outputs back to the same common scale factor before adding them. The order in which parallel sections are calculated does not matter, since the outputs are simply added together at the end. iv) Cascade form: In the cascade form, the output of one section forms the input to the next: 26
27 Mathematically, it does not matter in which order the sections are placed - the result will be the same. This assumes that there are no errors. In practice, the propagation of errors is crucial to the success of an IIR filter so the order of the sections in the cascade form is vital. In the cascade form, the output of one section forms the input to the next: In practice, the propagation of errors is crucial to the success of an IIR filter so the order of the sections in the cascade, and the selection of which filter coefficients to group in each section, is vital: Sections with high gain are undesirable because they increase the need for scaling and so increase quantisation errors It is desirable to arrange sections to avoid excessive scaling Review Questions: 1. Mention the types of realizing a digital IIR filter. 2. How a digital filter can be realized using cascade form? 27
EE482: Digital Signal Processing Applications
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu EE482: Digital Signal Processing Applications Spring 2014 TTh 14:30-15:45 CBC C222 Lecture 05 IIR Design 14/03/04 http://www.ee.unlv.edu/~b1morris/ee482/
More informationECE 410 DIGITAL SIGNAL PROCESSING D. Munson University of Illinois Chapter 12
. ECE 40 DIGITAL SIGNAL PROCESSING D. Munson University of Illinois Chapter IIR Filter Design ) Based on Analog Prototype a) Impulse invariant design b) Bilinear transformation ( ) ~ widely used ) Computer-Aided
More informationDigital Signal Processing
COMP ENG 4TL4: Digital Signal Processing Notes for Lecture #24 Tuesday, November 4, 2003 6.8 IIR Filter Design Properties of IIR Filters: IIR filters may be unstable Causal IIR filters with rational system
More informationMITOCW watch?v=jtj3v Rx7E
MITOCW watch?v=jtj3v Rx7E The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources for free. To
More informationDesign of IIR filters
Design of IIR filters Standard methods of design of digital infinite impulse response (IIR) filters usually consist of three steps, namely: 1 design of a continuous-time (CT) prototype low-pass filter;
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 informationDIGITAL SIGNAL PROCESSING. Chapter 6 IIR Filter Design
DIGITAL SIGNAL PROCESSING Chapter 6 IIR Filter Design OER Digital Signal Processing by Dr. Norizam Sulaiman work is under licensed Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International
More informationPS403 - Digital Signal processing
PS403 - Digital Signal processing 6. DSP - Recursive (IIR) Digital Filters Key Text: Digital Signal Processing with Computer Applications (2 nd Ed.) Paul A Lynn and Wolfgang Fuerst, (Publisher: John Wiley
More informationDigital Signal Processing Lecture 8 - Filter Design - IIR
Digital Signal Processing - Filter Design - IIR Electrical Engineering and Computer Science University of Tennessee, Knoxville October 20, 2015 Overview 1 2 3 4 5 6 Roadmap Discrete-time signals and systems
More informationUNIT - III PART A. 2. Mention any two techniques for digitizing the transfer function of an analog filter?
UNIT - III PART A. Mention the important features of the IIR filters? i) The physically realizable IIR filters does not have linear phase. ii) The IIR filter specification includes the desired characteristics
More informationDigital Control & Digital Filters. Lectures 21 & 22
Digital Controls & Digital Filters Lectures 2 & 22, Professor Department of Electrical and Computer Engineering Colorado State University Spring 205 Review of Analog Filters-Cont. Types of Analog Filters:
More informationDigital Signal Processing IIR Filter Design via Bilinear Transform
Digital Signal Processing IIR Filter Design via Bilinear Transform D. Richard Brown III D. Richard Brown III 1 / 12 Basic Procedure We assume here that we ve already decided to use an IIR filter. The basic
More informationChapter 7: IIR Filter Design Techniques
IUST-EE Chapter 7: IIR Filter Design Techniques Contents Performance Specifications Pole-Zero Placement Method Impulse Invariant Method Bilinear Transformation Classical Analog Filters DSP-Shokouhi Advantages
More informationFilters and Tuned Amplifiers
Filters and Tuned Amplifiers Essential building block in many systems, particularly in communication and instrumentation systems Typically implemented in one of three technologies: passive LC filters,
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 informationLike bilateral Laplace transforms, ROC must be used to determine a unique inverse z-transform.
Inversion of the z-transform Focus on rational z-transform of z 1. Apply partial fraction expansion. Like bilateral Laplace transforms, ROC must be used to determine a unique inverse z-transform. Let X(z)
More information1 1.27z z 2. 1 z H 2
E481 Digital Signal Processing Exam Date: Thursday -1-1 16:15 18:45 Final Exam - Solutions Dan Ellis 1. (a) In this direct-form II second-order-section filter, the first stage has
More informationStability Condition in Terms of the Pole Locations
Stability Condition in Terms of the Pole Locations A causal LTI digital filter is BIBO stable if and only if its impulse response h[n] is absolutely summable, i.e., 1 = S h [ n] < n= We now develop a stability
More informationCast of Characters. Some Symbols, Functions, and Variables Used in the Book
Page 1 of 6 Cast of Characters Some s, Functions, and Variables Used in the Book Digital Signal Processing and the Microcontroller by Dale Grover and John R. Deller ISBN 0-13-081348-6 Prentice Hall, 1998
More informationDHANALAKSHMI COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRICAL AND ELECTRONICS ENGINEERING EC2314- DIGITAL SIGNAL PROCESSING UNIT I INTRODUCTION PART A
DHANALAKSHMI COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRICAL AND ELECTRONICS ENGINEERING EC2314- DIGITAL SIGNAL PROCESSING UNIT I INTRODUCTION PART A Classification of systems : Continuous and Discrete
More informationResponses of Digital Filters Chapter Intended Learning Outcomes:
Responses of Digital Filters Chapter Intended Learning Outcomes: (i) Understanding the relationships between impulse response, frequency response, difference equation and transfer function in characterizing
More informationLecture 19 IIR Filters
Lecture 19 IIR Filters Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/5/10 1 General IIR Difference Equation IIR system: infinite-impulse response system The most general class
More informationIT DIGITAL SIGNAL PROCESSING (2013 regulation) UNIT-1 SIGNALS AND SYSTEMS PART-A
DEPARTMENT OF ELECTRONICS AND COMMUNICATION ENGINEERING IT6502 - DIGITAL SIGNAL PROCESSING (2013 regulation) UNIT-1 SIGNALS AND SYSTEMS PART-A 1. What is a continuous and discrete time signal? Continuous
More informationEE 521: Instrumentation and Measurements
Aly El-Osery Electrical Engineering Department, New Mexico Tech Socorro, New Mexico, USA November 1, 2009 1 / 27 1 The z-transform 2 Linear Time-Invariant System 3 Filter Design IIR Filters FIR Filters
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 informationINFINITE-IMPULSE RESPONSE DIGITAL FILTERS Classical analog filters and their conversion to digital filters 4. THE BUTTERWORTH ANALOG FILTER
INFINITE-IMPULSE RESPONSE DIGITAL FILTERS Classical analog filters and their conversion to digital filters. INTRODUCTION 2. IIR FILTER DESIGN 3. ANALOG FILTERS 4. THE BUTTERWORTH ANALOG FILTER 5. THE CHEBYSHEV-I
More informationTime Series Analysis: 4. Linear filters. P. F. Góra
Time Series Analysis: 4. Linear filters P. F. Góra http://th-www.if.uj.edu.pl/zfs/gora/ 2012 Linear filters in the Fourier domain Filtering: Multiplying the transform by a transfer function. g n DFT G
More informationLecture 7 Discrete Systems
Lecture 7 Discrete Systems EE 52: Instrumentation and Measurements Lecture Notes Update on November, 29 Aly El-Osery, Electrical Engineering Dept., New Mexico Tech 7. Contents The z-transform 2 Linear
More informationTime Series Analysis: 4. Digital Linear Filters. P. F. Góra
Time Series Analysis: 4. Digital Linear Filters P. F. Góra http://th-www.if.uj.edu.pl/zfs/gora/ 2018 Linear filters Filtering in Fourier domain is very easy: multiply the DFT of the input by a transfer
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 informationLecture 16: Filter Design: Impulse Invariance and Bilinear Transform
EE58 Digital Signal Processing University of Washington Autumn 2 Dept. of Electrical Engineering Lecture 6: Filter Design: Impulse Invariance and Bilinear Transform Nov 26, 2 Prof: J. Bilmes
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 informationImplementation of Discrete-Time Systems
EEE443 Digital Signal Processing Implementation of Discrete-Time Systems Dr. Shahrel A. Suandi PPKEE, Engineering Campus, USM Introduction A linear-time invariant system (LTI) is described by linear constant
More informationChapter 7: Filter Design 7.1 Practical Filter Terminology
hapter 7: Filter Design 7. Practical Filter Terminology Analog and digital filters and their designs constitute one of the major emphasis areas in signal processing and communication systems. This is due
More informationLecture 9 Infinite Impulse Response Filters
Lecture 9 Infinite Impulse Response Filters Outline 9 Infinite Impulse Response Filters 9 First-Order Low-Pass Filter 93 IIR Filter Design 5 93 CT Butterworth filter design 5 93 Bilinear transform 7 9
More informationVALLIAMMAI ENGINEERING COLLEGE. SRM Nagar, Kattankulathur DEPARTMENT OF INFORMATION TECHNOLOGY. Academic Year
VALLIAMMAI ENGINEERING COLLEGE SRM Nagar, Kattankulathur- 603 203 DEPARTMENT OF INFORMATION TECHNOLOGY Academic Year 2016-2017 QUESTION BANK-ODD SEMESTER NAME OF THE SUBJECT SUBJECT CODE SEMESTER YEAR
More informationSpeaker: Arthur Williams Chief Scientist Telebyte Inc. Thursday November 20 th 2008 INTRODUCTION TO ACTIVE AND PASSIVE ANALOG
INTRODUCTION TO ACTIVE AND PASSIVE ANALOG FILTER DESIGN INCLUDING SOME INTERESTING AND UNIQUE CONFIGURATIONS Speaker: Arthur Williams Chief Scientist Telebyte Inc. Thursday November 20 th 2008 TOPICS Introduction
More informationUNIVERSITI SAINS MALAYSIA. EEE 512/4 Advanced Digital Signal and Image Processing
-1- [EEE 512/4] UNIVERSITI SAINS MALAYSIA First Semester Examination 2013/2014 Academic Session December 2013 / January 2014 EEE 512/4 Advanced Digital Signal and Image Processing Duration : 3 hours Please
More informationINF3440/INF4440. Design of digital filters
Last week lecture Today s lecture: Chapter 8.1-8.3, 8.4.2, 8.5.3 INF3440/INF4440. Design of digital filters October 2004 Last week lecture Today s lecture: Chapter 8.1-8.3, 8.4.2, 8.5.3 Last lectures:
More informationFilter structures ELEC-E5410
Filter structures ELEC-E5410 Contents FIR filter basics Ideal impulse responses Polyphase decomposition Fractional delay by polyphase structure Nyquist filters Half-band filters Gibbs phenomenon Discrete-time
More informationEECE 301 Signals & Systems Prof. Mark Fowler
EECE 3 Signals & Systems Prof. ark Fowler Note Set #28 D-T Systems: DT Filters Ideal & Practical /4 Ideal D-T Filters Just as in the CT case we can specify filters. We looked at the ideal filter for the
More informationPoles and Zeros in z-plane
M58 Mixed Signal Processors page of 6 Poles and Zeros in z-plane z-plane Response of discrete-time system (i.e. digital filter at a particular frequency ω is determined by the distance between its poles
More informationAnalysis of Finite Wordlength Effects
Analysis of Finite Wordlength Effects Ideally, the system parameters along with the signal variables have infinite precision taing any value between and In practice, they can tae only discrete values within
More informationDigital Signal Processing Lecture 9 - Design of Digital Filters - FIR
Digital Signal Processing - Design of Digital Filters - FIR Electrical Engineering and Computer Science University of Tennessee, Knoxville November 3, 2015 Overview 1 2 3 4 Roadmap Introduction Discrete-time
More informationVel Tech High Tech Dr.Ranagarajan Dr.Sakunthala Engineering College Department of ECE
Subject Code: EC6502 Course Code:C302 Course Name: PRINCIPLES OF DIGITAL SIGNAL PROCESSING L-3 : T-1 : P-0 : Credits 4 COURSE OBJECTIVES: 1. To learn discrete Fourier transform and its properties 2. To
More informationV. IIR Digital Filters
Digital Signal Processing 5 March 5, V. IIR Digital Filters (Deleted in 7 Syllabus). (dded in 7 Syllabus). 7 Syllabus: nalog filter approximations Butterworth and Chebyshev, Design of IIR digital filters
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 informationLet H(z) = P(z)/Q(z) be the system function of a rational form. Let us represent both P(z) and Q(z) as polynomials of z (not z -1 )
Review: Poles and Zeros of Fractional Form Let H() = P()/Q() be the system function of a rational form. Let us represent both P() and Q() as polynomials of (not - ) Then Poles: the roots of Q()=0 Zeros:
More informationDISCRETE-TIME SIGNAL PROCESSING
THIRD EDITION DISCRETE-TIME SIGNAL PROCESSING ALAN V. OPPENHEIM MASSACHUSETTS INSTITUTE OF TECHNOLOGY RONALD W. SCHÄFER HEWLETT-PACKARD LABORATORIES Upper Saddle River Boston Columbus San Francisco New
More informationQuestion Bank. UNIT 1 Part-A
FATIMA MICHAEL COLLEGE OF ENGINEERING & TECHNOLOGY Senkottai Village, Madurai Sivagangai Main Road, Madurai -625 020 An ISO 9001:2008 Certified Institution Question Bank DEPARTMENT OF ELECTRONICS AND COMMUNICATION
More informationThe Approximation Problem
EE 508 Lecture 3 The Approximation Problem Classical Approximating Functions - Thompson and Bessel Approximations Review from Last Time Elliptic Filters Can be thought of as an extension of the CC approach
More information(Refer Slide Time: 02:11 to 04:19)
Digital Signal Processing Prof. S. C. Dutta Roy Department of Electrical Engineering Indian Institute of Technology, Delhi Lecture - 24 Analog Chebyshev LPF Design This is the 24 th lecture on DSP and
More informationNotes on L (Optimal) Filters
Notes on L (Optimal) Filters by C. Bond, 20 Background In 959 A. Papoulis published a paper which completed the description of a new class of filters with optimal properties, originally reported in 958.
More information-Digital Signal Processing- FIR Filter Design. Lecture May-16
-Digital Signal Processing- FIR Filter Design Lecture-17 24-May-16 FIR Filter Design! FIR filters can also be designed from a frequency response specification.! The equivalent sampled impulse response
More informationBasic Design Approaches
(Classic) IIR filter design: Basic Design Approaches. Convert the digital filter specifications into an analog prototype lowpass filter specifications. Determine the analog lowpass filter transfer function
More informationAll-Pole Recursive Digital Filters Design Based on Ultraspherical Polynomials
RADIOENGINEERING, VOL. 3, NO. 3, SEPTEMBER 4 949 All-Pole Recursive Digital Filters Design Based on Ultraspherical Polynomials Nikola STOJANOVIĆ, Negovan STAMENKOVIĆ, Vidosav STOJANOVIĆ University of Niš,
More informationButterworth Filter Properties
OpenStax-CNX module: m693 Butterworth Filter Properties C. Sidney Burrus This work is produced by OpenStax-CNX and licensed under the Creative Commons Attribution License 3. This section develops the properties
More informationOversampling Converters
Oversampling Converters David Johns and Ken Martin (johns@eecg.toronto.edu) (martin@eecg.toronto.edu) slide 1 of 56 Motivation Popular approach for medium-to-low speed A/D and D/A applications requiring
More informationAnalog and Digital Filter Design
Analog and Digital Filter Design by Jens Hee http://jenshee.dk October 208 Change log 28. september 208. Document started.. october 208. Figures added. 6. october 208. Bilinear transform chapter extended.
More informationElectronic Circuits EE359A
Electronic Circuits EE359A Bruce McNair B26 bmcnair@stevens.edu 21-216-5549 Lecture 22 578 Second order LCR resonator-poles V o I 1 1 = = Y 1 1 + sc + sl R s = C 2 s 1 s + + CR LC s = C 2 sω 2 s + + ω
More informationComputer-Aided Design of Digital Filters. Digital Filters. Digital Filters. Digital Filters. Design of Equiripple Linear-Phase FIR Filters
Computer-Aided Design of Digital Filters The FIR filter design techniques discussed so far can be easily implemented on a computer In addition, there are a number of FIR filter design algorithms that rely
More informationAn Iir-Filter Example: A Butterworth Filter
An Iir-Filter Example: A Butterworth Filter Josef Goette Bern University of Applied Sciences, Biel Institute of Human Centered Engineering - microlab JosefGoette@bfhch February 7, 2017 Contents 1 Introduction
More informationFilter Design Problem
Filter Design Problem Design of frequency-selective filters usually starts with a specification of their frequency response function. Practical filters have passband and stopband ripples, while exhibiting
More informationELEG 305: Digital Signal Processing
ELEG 305: Digital Signal Processing Lecture : Design of Digital IIR Filters (Part I) Kenneth E. Barner Department of Electrical and Computer Engineering University of Delaware Fall 008 K. E. Barner (Univ.
More informationCITY UNIVERSITY LONDON. MSc in Information Engineering DIGITAL SIGNAL PROCESSING EPM746
No: CITY UNIVERSITY LONDON MSc in Information Engineering DIGITAL SIGNAL PROCESSING EPM746 Date: 19 May 2004 Time: 09:00-11:00 Attempt Three out of FIVE questions, at least One question from PART B PART
More informationR13 SET - 1
R13 SET - 1 III B. Tech II Semester Regular Examinations, April - 2016 DIGITAL SIGNAL PROCESSING (Electronics and Communication Engineering) Time: 3 hours Maximum Marks: 70 Note: 1. Question Paper consists
More informationELEG 5173L Digital Signal Processing Ch. 5 Digital Filters
Department of Electrical Engineering University of Aransas ELEG 573L Digital Signal Processing Ch. 5 Digital Filters Dr. Jingxian Wu wuj@uar.edu OUTLINE 2 FIR and IIR Filters Filter Structures Analog Filters
More informationUNIVERSITY OF OSLO. Faculty of mathematics and natural sciences. Forslag til fasit, versjon-01: Problem 1 Signals and systems.
UNIVERSITY OF OSLO Faculty of mathematics and natural sciences Examination in INF3470/4470 Digital signal processing Day of examination: December 1th, 016 Examination hours: 14:30 18.30 This problem set
More informationTransform analysis of LTI systems Oppenheim and Schafer, Second edition pp For LTI systems we can write
Transform analysis of LTI systems Oppenheim and Schafer, Second edition pp. 4 9. For LTI systems we can write yœn D xœn hœn D X kd xœkhœn Alternatively, this relationship can be expressed in the z-transform
More informationAPPLIED SIGNAL PROCESSING
APPLIED SIGNAL PROCESSING DIGITAL FILTERS Digital filters are discrete-time linear systems { x[n] } G { y[n] } Impulse response: y[n] = h[0]x[n] + h[1]x[n 1] + 2 DIGITAL FILTER TYPES FIR (Finite Impulse
More informationw n = c k v n k (1.226) w n = c k v n k + d k w n k (1.227) Clearly non-recursive filters are a special case of recursive filters where M=0.
Random Data 79 1.13 Digital Filters There are two fundamental types of digital filters Non-recursive N w n = c k v n k (1.226) k= N and recursive N M w n = c k v n k + d k w n k (1.227) k= N k=1 Clearly
More information2.161 Signal Processing: Continuous and Discrete Fall 2008
MT 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 informationTheory and Problems of Signals and Systems
SCHAUM'S OUTLINES OF Theory and Problems of Signals and Systems HWEI P. HSU is Professor of Electrical Engineering at Fairleigh Dickinson University. He received his B.S. from National Taiwan University
More informationSPEECH ANALYSIS AND SYNTHESIS
16 Chapter 2 SPEECH ANALYSIS AND SYNTHESIS 2.1 INTRODUCTION: Speech signal analysis is used to characterize the spectral information of an input speech signal. Speech signal analysis [52-53] techniques
More information# FIR. [ ] = b k. # [ ]x[ n " k] [ ] = h k. x[ n] = Ae j" e j# ˆ n Complex exponential input. [ ]Ae j" e j ˆ. ˆ )Ae j# e j ˆ. y n. y n.
[ ] = h k M [ ] = b k x[ n " k] FIR k= M [ ]x[ n " k] convolution k= x[ n] = Ae j" e j ˆ n Complex exponential input [ ] = h k M % k= [ ]Ae j" e j ˆ % M = ' h[ k]e " j ˆ & k= k = H (" ˆ )Ae j e j ˆ ( )
More informationCMPT 889: Lecture 5 Filters
CMPT 889: Lecture 5 Filters Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University October 7, 2009 1 Digital Filters Any medium through which a signal passes may be regarded
More informationDFT & Fast Fourier Transform PART-A. 7. Calculate the number of multiplications needed in the calculation of DFT and FFT with 64 point sequence.
SHRI ANGALAMMAN COLLEGE OF ENGINEERING & TECHNOLOGY (An ISO 9001:2008 Certified Institution) SIRUGANOOR,TRICHY-621105. DEPARTMENT OF ELECTRONICS AND COMMUNICATION ENGINEERING UNIT I DFT & Fast Fourier
More informationDigital Filters. Linearity and Time Invariance. Linear Time-Invariant (LTI) Filters: CMPT 889: Lecture 5 Filters
Digital Filters CMPT 889: Lecture 5 Filters Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University October 7, 29 Any medium through which a signal passes may be regarded as
More information(Refer Slide Time: 01:28 03:51 min)
Digital Signal Processing Prof. S. C. Dutta Roy Department of Electrical Engineering Indian Institute of Technology, Delhi Lecture 40 FIR Design by Windowing This is the 40 th lecture and our topic for
More informationLINEAR-PHASE FIR FILTERS DESIGN
LINEAR-PHASE FIR FILTERS DESIGN Prof. Siripong Potisuk inimum-phase Filters A digital filter is a minimum-phase filter if and only if all of its zeros lie inside or on the unit circle; otherwise, it is
More informationDSP Configurations. responded with: thus the system function for this filter would be
DSP Configurations In this lecture we discuss the different physical (or software) configurations that can be used to actually realize or implement DSP functions. Recall that the general form of a DSP
More informationOptimum Ordering and Pole-Zero Pairing of the Cascade Form IIR. Digital Filter
Optimum Ordering and Pole-Zero Pairing of the Cascade Form IIR Digital Filter There are many possible cascade realiations of a higher order IIR transfer function obtained by different pole-ero pairings
More informationHomework Assignment 11
Homework Assignment Question State and then explain in 2 3 sentences, the advantage of switched capacitor filters compared to continuous-time active filters. (3 points) Continuous time filters use resistors
More informationDSP-CIS. Chapter-4: FIR & IIR Filter Design. Marc Moonen
DSP-CIS Chapter-4: FIR & IIR Filter Design Marc Moonen Dept. E.E./ESAT-STADIUS, KU Leuven marc.moonen@esat.kuleuven.be www.esat.kuleuven.be/stadius/ PART-II : Filter Design/Realization Step-1 : Define
More informationThe Approximation Problem
EE 508 Lecture The Approximation Problem Classical Approximating Functions - Elliptic Approximations - Thompson and Bessel Approximations Review from Last Time Chebyshev Approximations T Type II Chebyshev
More informationSecond and Higher-Order Delta-Sigma Modulators
Second and Higher-Order Delta-Sigma Modulators MEAD March 28 Richard Schreier Richard.Schreier@analog.com ANALOG DEVICES Overview MOD2: The 2 nd -Order Modulator MOD2 from MOD NTF (predicted & actual)
More informationREAL TIME DIGITAL SIGNAL PROCESSING
REAL TIME DIGITAL SIGNAL PROCESSING www.electron.frba.utn.edu.ar/dplab Digital Filters FIR and IIR. Design parameters. Implementation types. Constraints. Filters: General classification Filters: General
More informationECE503: Digital Signal Processing Lecture 5
ECE53: Digital Signal Processing Lecture 5 D. Richard Brown III WPI 3-February-22 WPI D. Richard Brown III 3-February-22 / 32 Lecture 5 Topics. Magnitude and phase characterization of transfer functions
More informationLECTURE NOTES DIGITAL SIGNAL PROCESSING III B.TECH II SEMESTER (JNTUK R 13)
LECTURE NOTES ON DIGITAL SIGNAL PROCESSING III B.TECH II SEMESTER (JNTUK R 13) FACULTY : B.V.S.RENUKA DEVI (Asst.Prof) / Dr. K. SRINIVASA RAO (Assoc. Prof) DEPARTMENT OF ELECTRONICS AND COMMUNICATIONS
More informationLecture 5 - Assembly Programming(II), Intro to Digital Filters
GoBack Lecture 5 - Assembly Programming(II), Intro to Digital Filters James Barnes (James.Barnes@colostate.edu) Spring 2009 Colorado State University Dept of Electrical and Computer Engineering ECE423
More informationECGR4124 Digital Signal Processing Exam 2 Spring 2017
ECGR4124 Digital Signal Processing Exam 2 Spring 2017 Name: LAST 4 NUMBERS of Student Number: Do NOT begin until told to do so Make sure that you have all pages before starting NO TEXTBOOK, NO CALCULATOR,
More informationMatched Second Order Digital Filters
Matched Second Order Digital Filters Martin Vicanek 14. February 016 1 Introduction Second order sections are universal building blocks for digital filters. They are characterized by five coefficients,
More information2.161 Signal Processing: Continuous and Discrete
MIT OpenCourseWare http://ocw.mit.edu 2.161 Signal Processing: Continuous and Discrete Fall 28 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. MASSACHUSETTS
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 Series Representation of
Fourier Series Representation of Periodic Signals Rui Wang, Assistant professor Dept. of Information and Communication Tongji University it Email: ruiwang@tongji.edu.cn Outline The response of LIT system
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 informationResearch Article Design of One-Dimensional Linear Phase Digital IIR Filters Using Orthogonal Polynomials
International Scholarly Research Network ISRN Signal Processing Volume, Article ID 8776, 7 pages doi:.54//8776 Research Article Design of One-Dimensional Linear Phase Digital IIR Filters Using Orthogonal
More informationDesign of Stable IIR filters with prescribed flatness and approximately linear phase
Design of Stable IIR filters with prescribed flatness and approximately linear phase YASUNORI SUGITA Nagaoka University of Technology Dept. of Electrical Engineering Nagaoka city, Niigata-pref., JAPAN
More informationLAB 6: FIR Filter Design Summer 2011
University of Illinois at Urbana-Champaign Department of Electrical and Computer Engineering ECE 311: Digital Signal Processing Lab Chandra Radhakrishnan Peter Kairouz LAB 6: FIR Filter Design Summer 011
More informationRoll No. :... Invigilator s Signature :.. CS/B.Tech (EE-N)/SEM-6/EC-611/ DIGITAL SIGNAL PROCESSING. Time Allotted : 3 Hours Full Marks : 70
Name : Roll No. :.... Invigilator s Signature :.. CS/B.Tech (EE-N)/SEM-6/EC-611/2011 2011 DIGITAL SIGNAL PROCESSING Time Allotted : 3 Hours Full Marks : 70 The figures in the margin indicate full marks.
More information