Multirate signal processing
|
|
- Arron Sherman Collins
- 6 years ago
- Views:
Transcription
1 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 digital audio, communication, image processing, and HDTV. Conceptually, the process of sampling rate conversion consists of two steps: 1) reconstruction of a continuous-time signal, and 2) resampling at a different rate. In practice, however, the conversion is implemented by means of discrete-time signal processing. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 1/38
2 Downsampling by an integer factor A sampling rate compressor or simply compressor is a systems that reduced the sampling rate of an input discrete-time signal by an integer factor of D. The spectrum of the sampled discrete-time signal x D [m] is given by X D (e jωt D ) = 1 X c (j(ω lω D )) = T D l Z = 1 DT l Z ( X c (j Ω l 2π ) ), DT where Ω D = 2π/T D = Ω s /D is the reduced sampling rate. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 2/38
3 Downsampling by an integer factor (cont.) Expressing l = m + kd, with m = 0, 1,..., D 1, we obtain X D (e jωt D ) = 1 D D 1 m=0 [ 1 T ( X c (j Ω m 2π DT k 2π )) T k Z }{{} X(e j(ω 2πm/DT )T ) We note that the above observation is true since 1 ( X c (j Ω k 2π )) = X ( e jωt ). T T k Z Consequently, we obtain X D (e jωt D ) = 1 D D 1 m=0 ( ) X e j(ω mω D)T, which relates directly the spectra of x[n] and x D [n]. ]. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 3/38
4 Downsampling by an integer factor (cont.) If we define the normalized frequency variable as ω = ΩT D, and note that ΩT = ΩT D /D = ω/d, we obtain X D (e jω ) = 1 D D 1 m=0 X (e j(ω 2πm)/D) Hence, the process of computing X D (e jω ) directly from X(e jω ) consists of: 1 Stretching X(e jω ) by a factor of D to obtain X(e jω/d ); 2 Placing D copies of X(e jω/d ) at the frequencies ω = 2πm, where m = 0, 1,..., D 1; 3 Adding the D stretched and shifted replicas and then dividing the result by D to obtain X D(e jω ). Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 4/38
5 Example: D = 3 In the case of downsampling by a factor of D = 3, we have X D (e jω ) = 1 [ ] X(e jω/3 ) + X(e j(ω/3+2π/3) ) + X(e j(ω/3+4π/3) ). 3 Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 5/38
6 Sampling rate conversion To establish a closer analogy between CT and DT sampling, we define the discrete-time sampling rate as f s = 1/D or ω s = 2π/D, which simply means keep one out of every D samples. Then, to avoid aliasing, we should have the following conditions If X(e jω ) = 0, ω H ω π then ω s = 2π D 2ω H. If the above condition is not satisfied, we will have aliasing distortion. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 6/38
7 Sampling rate conversion (cont.) To avoid aliasing distortion we precede the downsampling by a lowpass filter with cutoff frequency ω c = π/d. The combined system is called a decimator and the corresponding process is known as decimation. If we use an FIR lowpass filter of order M, the output of the decimator is given by x D [m] = v[md] = M h[k]x[md k]. k=0 How is it different from convolution? Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 7/38
8 Example: D = 2 Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 8/38
9 Back reconstruction In the absence of aliasing, x c (t) can be recovered from x D [m] according to x c (t) = m Z x D [m] Evaluating x c (t) at t = nt results in x[n] = m Z sin [π(t mdt )/DT ]. π(t mdt )/DT sin [π(n md)/d] x D [m] π(n md)/d, which provides perfect reconstruction of x[n] directly from x D [m]. Clearly, the ideal discrete-time interpolator is not practically realizable, which necessitates the use of approximations. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 9/38
10 Upsampling by an integer factor Increasing the sampling rate of x[n] by an integer factor I requires the insertion of (I 1) samples between consecutive samples of x[n]. Thus, upsampling is an information preserving operation. The interpolated signal x I [n] can be obtained using the formula x I [n] = m Z x[m] sin [π(n mi)/i], π(n mi)/i which is a special case of a more general interpolation formula x I [n] = m Z x[m]g r [n mi], where g r [n] is an interpolating filter (aka interpolator ). Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 10/38
11 Upsampling by an integer factor (cont.) In the domain of DTFT, x I [n] = m Z x[m]g r[n mi] becomes X I (e jω ) = G r (e jω )X(e jωi ), where X(e jωi ) is the DTFT of sequence x u [n] obtained by inserting (I 1) zeros between consecutive samples of x[n], that is, by { x[n/i], n is a multiple of I, x u [n] = 0, otherwise, with X u (e jω ) = X(e jωi ). The above system is called upsampler or sampling rate expander. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 11/38
12 Upsampling by an integer factor (cont.) Thus, the interpolation can be alternatively implemented as X I (e jω ) = G r (e jω )X u (e jω ), which corresponds to an LTI system described by x I [n] = x u [m]g r [n m]. m Z Schematically, we have Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 12/38
13 Example: I = 3 In practice, the ideal interpolators are replaced by their realizable approximations. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 13/38
14 Example: Linear interpolation Due to its computational efficiency, linear interpolation has long been a default choice of many multirate systems. Let m and n = mi + k (with k = 0, 1,..., I 1) be a slow and a fast time index, respectively. Then, linear interpolation is defined as x I [(m 1)I + k] = x[m 1] + (x[m] x[m 1]) (k/i), for k = 0, 1,..., I 1. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 14/38
15 Example: Linear interpolation (cont.) It is easy to verify that, in this case, x I [n] is obtained by convolving x u [n] with the triangular interpolator { 1 n /I, I < n < I, g lin [n] = 0, otherwise. The frequency response of the linear interpolation filter is given by [ ] 2 sin(ωi/2). G lin (e jω ) = 1 I ω/2 Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 15/38
16 Example: Linear interpolation (cont.) Despite the limited ability of linear interpolation to cope with aliasing, its performance is often deemed satisfactory. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 16/38
17 Sampling rate change by a rational factor Interpolation and decimation can be combined to change the sampling rate by a rational factor I/D. Let s consider the cascade connection of an interpolator and a decimator. The effective sampling rate of the output signal is T 0 = T D/I. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 17/38
18 Sampling rate change by a rational factor (cont.) Since the interpolation and decimation filters are in cascade and operate at the same sampling rate, they can be replaced by a single lowpass filter H(e jω ) = { I, 0 ω ω c 0, ω c ω π { π where ω c = min I, π }. D When I > D (sampling rate increase), H(e jω ) acts as a postfilter to remove the replicas introduced by the upsampling operation. When I < D (sampling rate decrease), H(e jω ) acts as an antialiasing filter for the downsampling operation. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 18/38
19 Half-band FIR filters Decimation and interpolation by a factor of two require lowpass filters with ω c = π/2, that is, ideal half-band filters. The impulse response of this filter is given by h[n] = ω c π sin(ω c n) ω c n = { 1/2, n = 0, 0, n = ±2, ±4,..., where, for simplicity, we assumed a zero-phase response. In general, any filter with h[n] = h[ n] that satisfies h[n] = 0, for n = ±2, ±4,..., and h[0] = 1/2 is called a half-band filter. Such filters can be shown to satisfy or, on the unit circle, H(z) + H( z) = 1, H(e jω ) + H(e j(ω π) ) = 1. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 19/38
20 Half-band FIR filters (cont.) H(e jω ) can be shown to have odd symmetric about ω = ±π/2. As a result, we have δ p = δ s and ω p = π ω s. Moreover, M/2 must be an odd number. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 20/38
21 Equiripple half-band filter design 1 Given the specifications ω s, A p, and A s of the half-band filter, compute δ = min{δ s, δ p }, ω p = π ω s. 2 Design a single band Type II FIR filter G(z) of order M/2 (odd) with ω p = 2ω p, ω s = π, and δ = 2δ using minimax optimization (i.e., equiripple design). Note that since G is Type II, G(e jπ ) = 0. 3 Scale the impulse response g[n] by 1/2, upsample the result by a factor of 2, and set the middle coefficient to 1/2. 4 The result is an impulse response h[n] with system function H(z) given by H(z) = 1 [ ] z M/2 + G(z 2 ), 2 which is a half-band filter with passband cutoff frequency ω p. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 21/38
22 Filter banks A filter bank is a collection of filters with a common input or a common output. The analysis filter bank splits x[n] into K sub-band signals v k [n], using the analysis filters H k (z). The synthesis filter bank consists of K synthesis filters G k (z) that combine K signals s k [n] into a signal y[n]. If s k [n] = v k [n], for all k, we would like to have y[n] = x[n]. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 22/38
23 Two-channel filter banks Maximally decimated two-channel filter banks use K = 2. In this case, y[n] can be shown to be given by Y (z) = 1 [T (z)x(z) + A(z)X( z)], 2 where T (z) = H 0 (z)g 0 (z) + H 1 (z)g 1 (z), A(z) = H 0 ( z)g 0 (z) + H 1 ( z)g 1 (z). Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 23/38
24 Two-channel filter banks (cont.) Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 24/38
25 Two-channel filter banks (cont.) The system function A(z) (between X( z) and Y (z)) causes aliasing. The system function T (z) (between X(z) and Y (z)) may cause magnitude and phase distortion. That the alias-free filter bank is obtained by enforcing A(z) = 0, in which case we obtain an LTI system with a transfer function T (z)/2. In fact, the filters G 0 (z) (lowpass) and G 1 (z) (highpass) are designed to eliminate the aliasing term. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 25/38
26 Perfect reconstruction Perfect reconstruction is possible if and only if A(z) = H 0 ( z)g 0 (z) + H 1 ( z)g 1 (z) = 0, T (z) = H 0 (z)g 0 (z) + H 1 (z)g 1 (z) = 2z n D, where n D is an arbitrary integer delay. The condition can also be written as [ ] [ ] [ H0 (z) H 1 (z) G0 (z) 2z n D = H 0 ( z) H 1 ( z) G 1 (z) 0 ]. Given H 0 (z) & H 1 (z), the above system can be solved for G 0 (z) & G 1 (z) as long as m (z) = H 0 (z)h 1 ( z) H 0 ( z)h 1 (z) 0. The resulting synthesis filters guarantee cancellation of aliasing error for any choice of analysis filters. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 26/38
27 Perfect reconstruction (cont.) In particular, defining G 0 (z) = 2z n D H 1( z) m (z) and G 1 (z) = 2z n D H 0( z) m (z), the condition for perfect reconstruction becomes R(z) + R( z) = 2 where R(z) is a product filter defined as R(z) = z n D H 0 (z)g 0 (z). The product filter R(z) plays a crucial role in analyzing and designing filter banks. It can be shown that R(z) has the form of R(z) = 1 + z 1 R 1 (z 2 ), and, therefore, it must be a half-band filter. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 27/38
28 Perfect reconstruction (cont.) Thus, to design a two-channel perfect reconstruction filter bank, it is necessary and sufficient to 1 obtain an R(z) = 1 + z 1 R 1(z 2 ), 2 perform the factorization R(z) = z n D H 0(z)G 0(z), 3 assign the synthesis filters G 0(z) and G 1(z). The choice of R(z) and the particular factorization determine the properties of the filter bank. Choosing R(z) = H(z)H(z 1 ) requires the design of a single filter and leads to orthogonal or para-unitary filter banks. Choosing R(z) = H 0 (z)g 0 (z) requires the design of two filters and leads to the more general bi-orthogonal filter banks. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 28/38
29 Orthogonal FIR filter banks If H 0 (z) and H 1 (z) are FIR filters and m (z) = c 0 z n0, then the synthesis filters G 0 (z) and G 1 (z) will be FIR as well. For example, given an M th -order FIR filter H(z), one can define H 0 (z) = H(z) and H 1 (z) = z M H( z 1 ), which are called conjugate quadrature filters (CQFs). To ensure that the synthesis filters are also FIR we require that or, alternatively, H(z)H(z 1 ) + H( z)h( z 1 ) = 1, H(e jω ) 2 + H(e j(ω π) ) 2 = 1, since H(z) has real coefficients (implying H(e jω ) = H (e jω )). Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 29/38
30 Orthogonal FIR filter banks (cont.) In this case, m (z) = z M (suggesting c 0 = 1 and n 0 = M). If we also choose n D = M, the synthesis filters are given by G 0 (z) = 2H 1 ( z) = 2z M H( z 1 ), G 1 (z) = 2H 0 ( z) = 2H( z). We note that the product filter is given by in which case R(z) = z n D H 0 (z)g 0 (z) = 2H(z)H(z 1 ), R(z) + R( z) = 2H(z)H(z 1 ) + 2H( z)h( z 1 ) = 2. Thus, the proposed filter bank has the perfect reconstruction property. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 30/38
31 Orthogonal FIR filter banks (cont.) When the solution of the spectral factorization problem R(z) = = z n D H 0 (z)g 0 (z) = 2H(z)H(z 1 ) exists? It is straightforward to show that H(z) satisfying H(z)H(z 1 )+ +H( z)h( z 1 ) = 1 exists if R(e jω ) 0. Moreover, since R(z) is a half-band filter, we have h[k]h[n + k] = δ[n] k assuming that h[n] is normalized so that h 2 = 1. The above condition suggests that h[n] is orthogonal to its own even translates (with the exception of n = 0), thus leading to the name orthogonal filter banks. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 31/38
32 Design of a perfect reconstruction CQF bank 1 Design a lowpass zero-phase half-band FIR filter R 0 (z) of order 2M, where M is an odd integer. 2 If δ min = min ω R 0 (e jω ) < 0, replace R 0 (e jω ) by or, alternatively, replace r[n] by R + (e jω ) = R 0 (e jω ) + δ min, r + [n] = r 0 [n] + δ min δ[n]. 3 Scale R + (z) so that R + (e jπ/2 ) = 1/2, viz. R(z) = δ min R +(z). 4 Determine the minimum-phase filter H(z) by solving the spectral factorization problem R(z) = H(z)H(z 1 ). 5 Specify the remaining filters of the bank. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 32/38
33 Example First, we design an FIR prototype H(z) of order 2M = 14 using [r0,deltamin]=firpm(14,[ ],[ ],[1 1]). Note that, in this case, deltamin is negative. We can also obtain various mixed-phase filters. Yet, no linearphase filter is possible. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 33/38
34 Example (cont.) The figure below shows the magnitude responses of the original half-band filter R 0 (z) and the conjugate quadrature filters H 0 (z) and H 1 (z). Note that R 0 (z) and H 0 (z) have the same stopband edge frequency, but different stopband attenuations. Rule of thumb: to obtain H 0 (z) with A s db, one should start with R 0 (z) having stopband attenuation larger than (2A s + 6) db. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 34/38
35 Tree-structured filter banks Two-channel filter banks can be used to split a signal into multiple sub-bands. A simple way to design the required filters is by cascading twochannel filter banks. The above is the uniform-band tree structure. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 35/38
36 Tree-structured filter banks (cont.) Alternatively, we can use the octave-band tree structure. The octave-band tree structure provides a logarithmic decomposition of the spectrum, very useful in coding. This filter bank produces a multiresolution decomposition of x[n], which is closely related to the discrete wavelet transform. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 36/38
37 Example In many real-life applications, the sub-band signals x 1 [n], x 01 [n], x 001 [n],... tend to be sparse (i.e., dominated by small values). It means that the energy of x[n] is encoded in a relatively small number of relatively large values of its details (i.e., x 1 [n], x 01 [n], x 001 [n],...). Thus, retaining these large values while setting the rest to zero will not have any significant impact on x[n] after reconstruction. As the small-amplitude details are mainly associated with noises and inessential features, the procedure of thresholding has found numerous applications in such areas as 1 signal compression and coding, 2 signal enhancement (de-noising), 3 signal modelling and representation, 4 pattern recognition and machine learning. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 37/38
38 Example: Image de-noising Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 10 38/38
Multirate 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 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 informationBasic Multi-rate Operations: Decimation and Interpolation
1 Basic Multirate Operations 2 Interconnection of Building Blocks 1.1 Decimation and Interpolation 1.2 Digital Filter Banks Basic Multi-rate Operations: Decimation and Interpolation Building blocks for
More informationMULTIRATE DIGITAL SIGNAL PROCESSING
MULTIRATE DIGITAL SIGNAL PROCESSING Signal processing can be enhanced by changing sampling rate: Up-sampling before D/A conversion in order to relax requirements of analog antialiasing filter. Cf. audio
More informationDigital Speech Processing Lecture 10. Short-Time Fourier Analysis Methods - Filter Bank Design
Digital Speech Processing Lecture Short-Time Fourier Analysis Methods - Filter Bank Design Review of STFT j j ˆ m ˆ. X e x[ mw ] [ nˆ m] e nˆ function of nˆ looks like a time sequence function of ˆ looks
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 information! Introduction. ! Discrete Time Signals & Systems. ! Z-Transform. ! Inverse Z-Transform. ! Sampling of Continuous Time Signals
ESE 531: Digital Signal Processing Lec 25: April 24, 2018 Review Course Content! Introduction! Discrete Time Signals & Systems! Discrete Time Fourier Transform! Z-Transform! Inverse Z-Transform! Sampling
More informationMassachusetts Institute of Technology Department of Electrical Engineering and Computer Science. Fall Solutions for Problem Set 2
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science Issued: Tuesday, September 5. 6.: Discrete-Time Signal Processing Fall 5 Solutions for Problem Set Problem.
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 informationMultiresolution image processing
Multiresolution image processing Laplacian pyramids Some applications of Laplacian pyramids Discrete Wavelet Transform (DWT) Wavelet theory Wavelet image compression Bernd Girod: EE368 Digital Image Processing
More informationECSE 512 Digital Signal Processing I Fall 2010 FINAL EXAMINATION
FINAL EXAMINATION 9:00 am 12:00 pm, December 20, 2010 Duration: 180 minutes Examiner: Prof. M. Vu Assoc. Examiner: Prof. B. Champagne There are 6 questions for a total of 120 points. This is a closed book
More informationLecture 11: Two Channel Filter Bank
WAVELETS AND MULTIRATE DIGITAL SIGNAL PROCESSING Lecture 11: Two Channel Filter Bank Prof.V.M.Gadre, EE, IIT Bombay 1 Introduction In the previous lecture we studied Z domain analysis of two channel filter
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 informationModule 4 MULTI- RESOLUTION ANALYSIS. Version 2 ECE IIT, Kharagpur
Module MULTI- RESOLUTION ANALYSIS Version ECE IIT, Kharagpur Lesson Multi-resolution Analysis: Theory of Subband Coding Version ECE IIT, Kharagpur Instructional Objectives At the end of this lesson, the
More informationDiscrete Fourier Transform
Discrete Fourier Transform Virtually all practical signals have finite length (e.g., sensor data, audio records, digital images, stock values, etc). Rather than considering such signals to be zero-padded
More information! Downsampling/Upsampling. ! Practical Interpolation. ! Non-integer Resampling. ! Multi-Rate Processing. " Interchanging Operations
Lecture Outline ESE 531: Digital Signal Processing Lec 10: February 14th, 2017 Practical and Non-integer Sampling, Multirate Sampling! Downsampling/! Practical Interpolation! Non-integer Resampling! Multi-Rate
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 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 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 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 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 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 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 informationCourse and Wavelets and Filter Banks. Filter Banks (contd.): perfect reconstruction; halfband filters and possible factorizations.
Course 18.327 and 1.130 Wavelets and Filter Banks Filter Banks (contd.): perfect reconstruction; halfband filters and possible factorizations. Product Filter Example: Product filter of degree 6 P 0 (z)
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 informationESE 531: Digital Signal Processing
ESE 531: Digital Signal Processing Lec 9: February 13th, 2018 Downsampling/Upsampling and Practical Interpolation Lecture Outline! CT processing of DT signals! Downsampling! Upsampling 2 Continuous-Time
More informationThe basic structure of the L-channel QMF bank is shown below
-Channel QMF Bans The basic structure of the -channel QMF ban is shown below The expressions for the -transforms of various intermediate signals in the above structure are given by Copyright, S. K. Mitra
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 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 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 informationReview of Fundamentals of Digital Signal Processing
Chapter 2 Review of Fundamentals of Digital Signal Processing 2.1 (a) This system is not linear (the constant term makes it non linear) but is shift-invariant (b) This system is linear but not shift-invariant
More informationPerfect Reconstruction Two- Channel FIR Filter Banks
Perfect Reconstruction Two- Channel FIR Filter Banks A perfect reconstruction two-channel FIR filter bank with linear-phase FIR filters can be designed if the power-complementary requirement e jω + e jω
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 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 informationECE538 Final Exam Fall 2017 Digital Signal Processing I 14 December Cover Sheet
ECE58 Final Exam Fall 7 Digital Signal Processing I December 7 Cover Sheet Test Duration: hours. Open Book but Closed Notes. Three double-sided 8.5 x crib sheets allowed This test contains five problems.
More informationReview of Fundamentals of Digital Signal Processing
Solution Manual for Theory and Applications of Digital Speech Processing by Lawrence Rabiner and Ronald Schafer Click here to Purchase full Solution Manual at http://solutionmanuals.info Link download
More informationFilter Banks II. Prof. Dr.-Ing. G. Schuller. Fraunhofer IDMT & Ilmenau University of Technology Ilmenau, Germany
Filter Banks II Prof. Dr.-Ing. G. Schuller Fraunhofer IDMT & Ilmenau University of Technology Ilmenau, Germany Page Modulated Filter Banks Extending the DCT The DCT IV transform can be seen as modulated
More informationReview of Discrete-Time System
Review of Discrete-Time System Electrical & Computer Engineering University of Maryland, College Park Acknowledgment: ENEE630 slides were based on class notes developed by Profs. K.J. Ray Liu and Min Wu.
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 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 information2. Typical Discrete-Time Systems All-Pass Systems (5.5) 2.2. Minimum-Phase Systems (5.6) 2.3. Generalized Linear-Phase Systems (5.
. Typical Discrete-Time Systems.1. All-Pass Systems (5.5).. Minimum-Phase Systems (5.6).3. Generalized Linear-Phase Systems (5.7) .1. All-Pass Systems An all-pass system is defined as a system which has
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 information4.1 Introduction. 2πδ ω (4.2) Applications of Fourier Representations to Mixed Signal Classes = (4.1)
4.1 Introduction Two cases of mixed signals to be studied in this chapter: 1. Periodic and nonperiodic signals 2. Continuous- and discrete-time signals Other descriptions: Refer to pp. 341-342, textbook.
More informationECE 413 Digital Signal Processing Midterm Exam, Spring Instructions:
University of Waterloo Department of Electrical and Computer Engineering ECE 4 Digital Signal Processing Midterm Exam, Spring 00 June 0th, 00, 5:0-6:50 PM Instructor: Dr. Oleg Michailovich Student s name:
More informationEEL3135: Homework #4
EEL335: Homework #4 Problem : For each of the systems below, determine whether or not the system is () linear, () time-invariant, and (3) causal: (a) (b) (c) xn [ ] cos( 04πn) (d) xn [ ] xn [ ] xn [ 5]
More informationSubband Coding and Wavelets. National Chiao Tung University Chun-Jen Tsai 12/04/2014
Subband Coding and Wavelets National Chiao Tung Universit Chun-Jen Tsai /4/4 Concept of Subband Coding In transform coding, we use N (or N N) samples as the data transform unit Transform coefficients are
More informationEE482: 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 informationDiscrete-Time David Johns and Ken Martin University of Toronto
Discrete-Time David Johns and Ken Martin University of Toronto (johns@eecg.toronto.edu) (martin@eecg.toronto.edu) University of Toronto 1 of 40 Overview of Some Signal Spectra x c () t st () x s () t xn
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 informationMULTIRATE SYSTEMS. work load, lower filter order, lower coefficient sensitivity and noise,
MULIRAE SYSEMS ransfer signals between two systems that operate with different sample frequencies Implemented system functions more efficiently by using several sample rates (for example narrow-band filters).
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 informationESE 531: Digital Signal Processing
ESE 531: Digital Signal Processing Lec 8: February 7th, 2017 Sampling and Reconstruction Lecture Outline! Review " Ideal sampling " Frequency response of sampled signal " Reconstruction " Anti-aliasing
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 informationMultimedia Signals and Systems - Audio and Video. Signal, Image, Video Processing Review-Introduction, MP3 and MPEG2
Multimedia Signals and Systems - Audio and Video Signal, Image, Video Processing Review-Introduction, MP3 and MPEG2 Kunio Takaya Electrical and Computer Engineering University of Saskatchewan December
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 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 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 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 informationUNIVERSITY OF OSLO. Please make sure that your copy of the problem set is complete before you attempt to answer anything.
UNIVERSITY OF OSLO Faculty of mathematics and natural sciences Examination in INF3470/4470 Digital signal processing Day of examination: December 9th, 011 Examination hours: 14.30 18.30 This problem set
More informationIntroduction to Digital Signal Processing
Introduction to Digital Signal Processing 1.1 What is DSP? DSP is a technique of performing the mathematical operations on the signals in digital domain. As real time signals are analog in nature we need
More informationDirect Design of Orthogonal Filter Banks and Wavelets
Direct Design of Orthogonal Filter Banks and Wavelets W.-S. Lu T. Hinamoto Dept. of Electrical & Computer Engineering Graduate School of Engineering University of Victoria Hiroshima University Victoria,
More informationAnalog vs. discrete signals
Analog vs. discrete signals Continuous-time signals are also known as analog signals because their amplitude is analogous (i.e., proportional) to the physical quantity they represent. Discrete-time signals
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 informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science Discrete-Time Signal Processing Fall 2005
1 MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science 6.341 Discrete-Time Signal Processing Fall 2005 FINAL EXAM Friday, December 16, 2005 Walker (50-340) 1:30pm
More informationENEE630 ADSP RECITATION 5 w/ solution Ver
ENEE630 ADSP RECITATION 5 w/ solution Ver20209 Consider the structures shown in Fig RI, with input transforms and filter responses as indicated Sketch the quantities Y 0 (e jω ) and Y (e jω ) Figure RI:
More informationFilter Banks II. Prof. Dr.-Ing Gerald Schuller. Fraunhofer IDMT & Ilmenau Technical University Ilmenau, Germany
Filter Banks II Prof. Dr.-Ing Gerald Schuller Fraunhofer IDMT & Ilmenau Technical University Ilmenau, Germany Prof. Dr.-Ing. G. Schuller, shl@idmt.fraunhofer.de Page Modulated Filter Banks Extending the
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 information6.003: Signals and Systems. Sampling and Quantization
6.003: Signals and Systems Sampling and Quantization December 1, 2009 Last Time: Sampling and Reconstruction Uniform sampling (sampling interval T ): x[n] = x(nt ) t n Impulse reconstruction: x p (t) =
More informationEECS 123 Digital Signal Processing University of California, Berkeley: Fall 2007 Gastpar November 7, Exam 2
EECS 3 Digital Signal Processing University of California, Berkeley: Fall 7 Gastpar November 7, 7 Exam Last name First name SID You have hour and 45 minutes to complete this exam. he exam is closed-book
More informationECE-700 Review. Phil Schniter. January 5, x c (t)e jωt dt, x[n]z n, Denoting a transform pair by x[n] X(z), some useful properties are
ECE-7 Review Phil Schniter January 5, 7 ransforms Using x c (t) to denote a continuous-time signal at time t R, Laplace ransform: X c (s) x c (t)e st dt, s C Continuous-ime Fourier ransform (CF): ote that:
More informationNAME: 11 December 2013 Digital Signal Processing I Final Exam Fall Cover Sheet
NAME: December Digital Signal Processing I Final Exam Fall Cover Sheet Test Duration: minutes. Open Book but Closed Notes. Three 8.5 x crib sheets allowed Calculators NOT allowed. This test contains four
More informationMulti-rate Signal Processing 7. M-channel Maximally Decmiated Filter Banks
Multi-rate Signal Processing 7. M-channel Maximally Decmiated Filter Banks Electrical & Computer Engineering University of Maryland, College Park Acknowledgment: ENEE630 slides were based on class notes
More informationECE 301 Division 1 Final Exam Solutions, 12/12/2011, 3:20-5:20pm in PHYS 114.
ECE 301 Division 1 Final Exam Solutions, 12/12/2011, 3:20-5:20pm in PHYS 114. The exam for both sections of ECE 301 is conducted in the same room, but the problems are completely different. Your ID will
More informationFilter Banks with Variable System Delay. Georgia Institute of Technology. Atlanta, GA Abstract
A General Formulation for Modulated Perfect Reconstruction Filter Banks with Variable System Delay Gerald Schuller and Mark J T Smith Digital Signal Processing Laboratory School of Electrical Engineering
More informationELEN 4810 Midterm Exam
ELEN 4810 Midterm Exam Wednesday, October 26, 2016, 10:10-11:25 AM. One sheet of handwritten notes is allowed. No electronics of any kind are allowed. Please record your answers in the exam booklet. Raise
More informationMulti-rate Signal Processing 3. The Polyphase Representation
Multi-rate Signal Processing 3. The Polyphase Representation Electrical & Computer Engineering University of Maryland, College Park Acknowledgment: ENEE630 slides were based on class notes developed by
More informationTransform Analysis of Linear Time-Invariant Systems
Transform Analysis of Linear Time-Invariant Systems Discrete-Time Signal Processing Chia-Ping Chen Department of Computer Science and Engineering National Sun Yat-Sen University Kaohsiung, Taiwan ROC 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 informationDigital Signal Processing. Midterm 1 Solution
EE 123 University of California, Berkeley Anant Sahai February 15, 27 Digital Signal Processing Instructions Midterm 1 Solution Total time allowed for the exam is 8 minutes Some useful formulas: Discrete
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 informationDiscrete Time Signals and Systems Time-frequency Analysis. Gloria Menegaz
Discrete Time Signals and Systems Time-frequency Analysis Gloria Menegaz Time-frequency Analysis Fourier transform (1D and 2D) Reference textbook: Discrete time signal processing, A.W. Oppenheim and R.W.
More informationWavelets and Multiresolution Processing
Wavelets and Multiresolution Processing Wavelets Fourier transform has it basis functions in sinusoids Wavelets based on small waves of varying frequency and limited duration In addition to frequency,
More informationLecture 8: Signal Reconstruction, DT vs CT Processing. 8.1 Reconstruction of a Band-limited Signal from its Samples
EE518 Digital Signal Processing University of Washington Autumn 2001 Dept. of Electrical Engineering Lecture 8: Signal Reconstruction, D vs C Processing Oct 24, 2001 Prof: J. Bilmes
More informationLecture 16: Multiresolution Image Analysis
Lecture 16: Multiresolution Image Analysis Harvey Rhody Chester F. Carlson Center for Imaging Science Rochester Institute of Technology rhody@cis.rit.edu November 9, 2004 Abstract Multiresolution analysis
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 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 informationVoiced Speech. Unvoiced Speech
Digital Speech Processing Lecture 2 Homomorphic Speech Processing General Discrete-Time Model of Speech Production p [ n] = p[ n] h [ n] Voiced Speech L h [ n] = A g[ n] v[ n] r[ n] V V V p [ n ] = u [
More informationEC Signals and Systems
UNIT I CLASSIFICATION OF SIGNALS AND SYSTEMS Continuous time signals (CT signals), discrete time signals (DT signals) Step, Ramp, Pulse, Impulse, Exponential 1. Define Unit Impulse Signal [M/J 1], [M/J
More informationDesign and Application of Quincunx Filter Banks
Design and Application of Quincunx Filter Banks by Yi Chen B.Eng., Tsinghua University, China, A Thesis Submitted in Partial Fulfillment of the Requirements for the Degree of MASTER OF APPLIED SCIENCE
More informationPeriodic (Uniform) Sampling ELEC364 & ELEC442
M.A. Amer Concordia University Electrical and Computer Engineering Content and Figures are from: Periodic (Uniform) Sampling ELEC364 & ELEC442 Introduction to sampling Introduction to filter Ideal sampling:
More informationLecture 8 Finite Impulse Response Filters
Lecture 8 Finite Impulse Response Filters Outline 8. Finite Impulse Response Filters.......................... 8. oving Average Filter............................... 8.. Phase response...............................
More informationDiscrete-time Signals and Systems in
Discrete-time Signals and Systems in the Frequency Domain Chapter 3, Sections 3.1-39 3.9 Chapter 4, Sections 4.8-4.9 Dr. Iyad Jafar Outline Introduction The Continuous-Time FourierTransform (CTFT) The
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 information( ) John A. Quinn Lecture. ESE 531: Digital Signal Processing. Lecture Outline. Frequency Response of LTI System. Example: Zero on Real Axis
John A. Quinn Lecture ESE 531: Digital Signal Processing Lec 15: March 21, 2017 Review, Generalized Linear Phase Systems Penn ESE 531 Spring 2017 Khanna Lecture Outline!!! 2 Frequency Response of LTI System
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 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 informationQUESTION BANK SIGNALS AND SYSTEMS (4 th SEM ECE)
QUESTION BANK SIGNALS AND SYSTEMS (4 th SEM ECE) 1. For the signal shown in Fig. 1, find x(2t + 3). i. Fig. 1 2. What is the classification of the systems? 3. What are the Dirichlet s conditions of Fourier
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 informationDigital Signal Processing
Digital Signal Processing Multirate Signal Processing Dr. Manar Mohaisen Office: F28 Email: manar.subhi@kut.ac.kr School of IT Engineering Review of the Precedent ecture Introduced Properties of FIR Filters
More informationECE-314 Fall 2012 Review Questions for Midterm Examination II
ECE-314 Fall 2012 Review Questions for Midterm Examination II First, make sure you study all the problems and their solutions from homework sets 4-7. Then work on the following additional problems. Problem
More information