Direct Design of Orthogonal Filter Banks and Wavelets
|
|
- Neil Wood
- 5 years ago
- Views:
Transcription
1 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, Canada Hiroshima, Japan 1
2 Outline Introduction Early and recent work Constrained linear updates and a convex QP formulation for least-squares design of conjugate quadrature (CQ) filters Constrained linear updates and an second-order cone programming formulation for minimax design of CQ filters Experimental results
3 Two-channel FIR filter bank 1. Introduction x(n) H (z) F (z) x(n) ^ H 1 (z) F 1 (z) 1 1 Xˆ ( z) = [ F( zh ) ( z) + F1( zh ) 1( z)] X( z) + [ F( zh ) ( z) + F1( zh ) 1( z)] X( z) Perfect reconstruction (PR) conditions F ( z) H ( z) + F( z) H ( z) = z l 1 1 F ( z) H ( z) + F( z) H ( z) = 1 1 3
4 A conjugate quadrature (CQ) filter bank assumes H ( z) = z H ( z ), F ( z) = z H ( z ), F( z) = z H ( z ) ( N 1) 1 ( N 1) 1 ( N 1) the nd PR condition is automatically satisfied (no aliasing) and the 1 st PR condition becomes H ( z) H ( z ) + H ( z) H ( z ) = (PS) 1 1 which is called the power symmetric (PS) condition because it implies 1 ( j( π θ) ) j( π + θ) ( ) H e + H e = for any θ 4
5 . Early and Recent Work Representative early and recent work include Smith and Barnwell (1984) Mintzer (1985) Vaidyanathan and Nguyen (1987) Rioul and Duhamel (1994) Lawton and Michelli (1997) Tuqan and Vaidyanathan (1998) Dumitrescu and Popeea () Tay (5, 6) 5
6 The most common design technique: A half-band filter P(z) is a zero-phase FIR filter satisfying P(z) + P( z) = Let P z = H z H z, then the PS condition 1 ( ) ( ) ( ) H ( z) H ( z ) + H ( z) H ( z ) = 1 1 becomes P(z) + P( z) = So P(z) is a half-band filer and it is nonnegative: ( ) P e = H ( e ) H ( e ) = H ( e ) (P) jω jω jω jω 6
7 Design steps: (a) Design a lowpass half-band FIR filter P(z) with j nonnegativity property Pe ( ω ) (b) Perform a spectral decomposition P z H z H z 1 ( ) = ( ) ( ) Vanishing moment (VM): the number of VMs equals to the number of zeros of H at ω = π : dh ( e ) l jω N 1 j l n n l h l n dω ω= π n= = ( ) ( 1) =, for l =, 1,, L 1 7
8 Problem Formulation Let 3. Least-Squares Design of CQ Filters H ( z) h z N 1 n = n with N even, and h = [h h 1 h N-1 ] T n= A direct approach: minimizing a least squares type objective function subject to the PS constraint: jω minimize H ( e ) dω h π ω a subject to: H ( z) H ( z ) + H ( z) H ( z ) = 1 1 8
9 The objective function is a positive definite quadratic form: π ω a ( jω ) T H e dω = h Qh with Q a symmetric positive definite Toeplitz matrix: Q 1 = toeplitz π ωa sinωa sin[( N 1) ωa] N 1 The constraint is the PS condition: H ( z) H ( z ) + H ( z) H ( z ) = (PS) 1 1 that is equivalent to N/ second-order equality constraints: 9
10 N 1 m n n+ m δ n= h h = ( m) for m=,1,,( N )/ with δ ( m) = 1 for m= and δ ( m) = for m. The design problem now becomes a polynomial optimization problem (POP): N 1 m n= minimize h n n+ m T 1/ hqh= Q h subject to: h h = δ ( m) for m ( N ) / 1
11 The POP can be modified to include VM requirement: minimize h N 1 m n= N 1 n= T 1/ hqh= Q h subject to: h h = δ ( m) for m ( N ) / n n+ m n l ( 1) nh = for l=,1,, L 1 n Features of these problems: All polynomials are of second-order. The objective function is convex Nonconvex problems because of the N/ 11
12 second-order equality constraints (the PS conditions). Examples of the PS constraints 1. N = 4 constraints: h + h + h + h = h h + h h = 1 3. N = 1 constraints: h + h + + h = 1 ( terms) 1 19 h h + h1 h3 + + h17 h19 h h4 + h1 h5 + + h15 h19 = = (18 terms) (16 terms) h h + h h + + h h = (14 terms) h h + h h + h h + h h = (4 terms) h h + h h = ( terms)
13 Constrained Linear Updates In the kth iteration of the algorithm we update filter coefficients from h (k) things: to h (k+1) = h (k) + d to achieve two to reduce the filter s stopband energy hqh T to better approximate constraints N 1 m n= h h = δ ( m), m N/ 1 n n+ m The stopband energy at h (k+1) is equal to 1/ ( k ) Q ( d + h ) The constraints at h (k+1) becomes h h + h d + d h + d d = δ ( m) ( k) ( k) ( k) ( k) n n+ m n n+ m n n+ m n n+ m n n n n 13
14 Imposing constraints on the smallness of increment vector d: di β for i = 1,,, N the nd -order constraints can be linearized: h h + h d + d h + d d ( k) ( k) ( k) ( k) n n+ m n n+ m n n+ m n n+ m n n n n ( k ) a known term, denoted by s ( m) = δ ( m) for m=,1,,( N )/ very small, drop ( k) ( k) ( k) ( k) hn hn+ m + hn dn+ m + dnhn+ m n n n linear updates 14
15 This leads to a set of (N )/ linear equations: n h d + d h = δ ( m) s ( m) u ( m) ( k) ( k) ( k) ( k) n n+ m n n+ m n which can be expressed as C d = u ( k) ( k) The smallness constraint on d is given by d β for 1 i n Ad b i The linear constraint on VMs is given by N 1 n= ( 1) n ( d + h ) = for l L 1 Dd = v n l ( k) ( k) n n 15
16 A Quadratic Programming (QP) Formulation Summarizing, the solution strategy is to iteratively update the filter coefficients from h (k) to h (k+1) = h (k) + d (k) with d (k) obtained by solving the QP problem 1/ ( k ) minimize Q ( d + h ) d subject to: Ad b ( k) ( k) C u d = ( k ) D v 16
17 Problem Formulation 4. Minimax Design of CQ Filters The formulation in this case is changed to N 1 m n= N 1 n= jω minimize maximize H ( e ) h subject to: h h = δ ( m) for m ( N ) / n n+ m a n l ( 1) nh = for l=,1,, L 1 n ω ω π 17
18 Constrained Linear Updates Like in the least squares design, the constrained linear update gives jω minimize maximize H ( e ) d ω ω π a subject to: Ad b ( k) ( k) C u d = ( k ) D v 18
19 Dealing with the objective function, we write N 1 jω = n jnω = T ω T n= H ( e ) h e h c( ) jh s( ω) T [ ] [ ] c( ω) = 1 cosω cos( N 1) ω, s( ω) = sinω sin( N 1) ω hence T jω T T c( ω) H( e ) = ( h c( ω) ) + ( h s( ω) ) = h T( ω) h T s( ω) T H (, ) ( ) ( ) ( ) e h + d = T ω h + d = T ω d + g jω ( k) ( k) ( k) ( k) ( k) ( k) 19
20 This converts the minimax problem into minimize η T ω d + g η ω ω π ( k) ( k) subject to: ( i) for { i} [ a, ] Ad b ( k) ( k) C u d = ( k ) D v which is an SOCP problem.
21 5. Experimental Results 5.1 Comparison with designs by Smith-Barnwell s method Filter H (z) Largest Eq. Error N = 8 N = 16 N = 3 H (z) of [] Refined H (z) < 1-15 H (z) of [] Refined H (z) < 1-15 H (z) of [] Refined H (z) <
22 1 N = 3, L =, ω a = Normalized frequency
23 5. LS and minimax designs with N = 96 and L =, 1,, 5 Least squares with N = 96, ω =.56π L Energy in Stopband Largest Equation Error < < < < < a 3
24 Minimax with N = 96, ω =.56π a L Instantaneous Energy in Stopband Largest Equation Error <
25 Minimax design with N = 96, L = 3, ω =.56π : a Normalized frequency 5
Direct Design of Orthogonal Filter Banks and Wavelets by Sequential Convex Quadratic Programming
Direct Design of Orthogonl Filter Bnks nd Wvelets by Sequentil Convex Qudrtic Progrmming W.-S. Lu Dept. of Electricl nd Computer Engineering University of Victori, Victori, Cnd September 8 1 Abstrct Two-chnnel
More informationTowards Global Design of Orthogonal Filter Banks and Wavelets
Towards Global Design of Orthogonal Filter Banks and Wavelets Jie Yan and Wu-Sheng Lu Department of Electrical and Computer Engineering University of Victoria Victoria, BC, Canada V8W 3P6 jyan@ece.uvic.ca,
More informationDesigns of Orthogonal Filter Banks and Orthogonal Cosine-Modulated Filter Banks
1 / 45 Designs of Orthogonal Filter Banks and Orthogonal Cosine-Modulated Filter Banks Jie Yan Department of Electrical and Computer Engineering University of Victoria April 16, 2010 2 / 45 OUTLINE 1 INTRODUCTION
More informationVariable Fractional Delay FIR Filters with Sparse Coefficients
Variable Fractional Delay FIR Filters with Sparse Coefficients W.-S. Lu T. Hinamoto Dept. of Electrical & Computer Engineering Graduate School of Engineering University of Victoria Hiroshima University
More informationEnhanced Steiglitz-McBride Procedure for. Minimax IIR Digital Filters
Enhanced Steiglitz-McBride Procedure for Minimax IIR Digital Filters Wu-Sheng Lu Takao Hinamoto University of Victoria Hiroshima University Victoria, Canada Higashi-Hiroshima, Japan May 30, 2018 1 Outline
More informationA Unified Approach to the Design of Interpolated and Frequency Response Masking FIR Filters
A Unified Approach to the Design of Interpolated and Frequency Response Masking FIR Filters Wu Sheng Lu akao Hinamoto University of Victoria Hiroshima University Victoria, Canada Higashi Hiroshima, Japan
More informationMinimax Design of Complex-Coefficient FIR Filters with Low Group Delay
Minimax Design of Complex-Coefficient FIR Filters with Low Group Delay Wu-Sheng Lu Takao Hinamoto Dept. of Elec. and Comp. Engineering Graduate School of Engineering University of Victoria Hiroshima University
More informationMultirate signal processing
Multirate signal processing Discrete-time systems with different sampling rates at various parts of the system are called multirate systems. The need for such systems arises in many applications, including
More informationSDP APPROXIMATION OF THE HALF DELAY AND THE DESIGN OF HILBERT PAIRS. Bogdan Dumitrescu
SDP APPROXIMATION OF THE HALF DELAY AND THE DESIGN OF HILBERT PAIRS Bogdan Dumitrescu Tampere International Center for Signal Processing Tampere University of Technology P.O.Box 553, 3311 Tampere, FINLAND
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 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 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 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 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 informationDesign of Orthonormal Wavelet Filter Banks Using the Remez Exchange Algorithm
Electronics and Communications in Japan, Part 3, Vol. 81, No. 6, 1998 Translated from Denshi Joho Tsushin Gakkai Ronbunshi, Vol. J80-A, No. 9, September 1997, pp. 1396 1402 Design of Orthonormal Wavelet
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 informationOptimal Design of Real and Complex Minimum Phase Digital FIR Filters
Optimal Design of Real and Complex Minimum Phase Digital FIR Filters Niranjan Damera-Venkata and Brian L. Evans Embedded Signal Processing Laboratory Dept. of Electrical and Computer Engineering The University
More informationReconstruction of Block-Sparse Signals by Using an l 2/p -Regularized Least-Squares Algorithm
Reconstruction of Block-Sparse Signals by Using an l 2/p -Regularized Least-Squares Algorithm Jeevan K. Pant, Wu-Sheng Lu, and Andreas Antoniou University of Victoria May 21, 2012 Compressive Sensing 1/23
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 informationDesign of Projection Matrix for Compressive Sensing by Nonsmooth Optimization
Design of Proection Matrix for Compressive Sensing by Nonsmooth Optimization W.-S. Lu T. Hinamoto Dept. of Electrical & Computer Engineering Graduate School of Engineering University of Victoria Hiroshima
More informationDesign of High-Performance Filter Banks for Image Coding
Design of High-Performance Filter Banks for Image Coding Di Xu and Michael D. Adams Dept. of Elec. and Comp. Engineering, University of Victoria PO Box 3055, STN CSC, Victoria, BC, V8W 3P6, Canada dixu@ece.uvic.ca
More informationTwo-Dimensional State-Space Digital Filters with Minimum Frequency-Weighted L 2 -Sensitivity
wo-dimensional State-Space Digital Filters with Minimum Frequency-Weighted L -Sensitivity under L -Scaling Constraints. Hinamoto,. Oumi, O. I. Omoifo W.-S. Lu Graduate School of Engineering Electrical
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 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 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 informationGröbner bases and wavelet design
Journal of Symbolic Computation 37 (24) 227 259 www.elsevier.com/locate/jsc Gröbner bases and wavelet design Jérôme Lebrun a,, Ivan Selesnick b a Laboratoire I3S, CNRS/UNSA, Sophia Antipolis, France b
More informationWavelet Bi-frames with Uniform Symmetry for Curve Multiresolution Processing
Wavelet Bi-frames with Uniform Symmetry for Curve Multiresolution Processing Qingtang Jiang Abstract This paper is about the construction of univariate wavelet bi-frames with each framelet being symmetric.
More informationTwo-Dimensional Orthogonal Filter Banks with Directional Vanishing Moments
Two-imensional Orthogonal Filter Banks with irectional Vanishing Moments Jianping Zhou and Minh N. o epartment of Electrical and Computer Engineering University of Illinois at Urbana-Champaign, Urbana,
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 informationMINIMUM PEAK IMPULSE FIR FILTER DESIGN
MINIMUM PEAK IMPULSE FIR FILTER DESIGN CHRISTINE S. LAW AND JON DATTORRO Abstract. Saturation pulses rf(t) are essential to many imaging applications. Criteria for desirable saturation profile are flat
More informationDigital Image Processing Lectures 15 & 16
Lectures 15 & 16, Professor Department of Electrical and Computer Engineering Colorado State University CWT and Multi-Resolution Signal Analysis Wavelet transform offers multi-resolution by allowing for
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 informationDESIGN OF ALIAS-FREE LINEAR PHASE QUADRATURE MIRROR FILTER BANKS USING EIGENVALUE-EIGENVECTOR APPROACH
DESIGN OF ALIAS-FREE LINEAR PHASE QUADRATURE MIRROR FILTER BANKS USING EIGENVALUE-EIGENVECTOR APPROACH ABSTRACT S. K. Agrawal 1# and O. P. Sahu 1& Electronics and Communication Engineering Department,
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 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 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 informationFast Wavelet/Framelet Transform for Signal/Image Processing.
Fast Wavelet/Framelet Transform for Signal/Image Processing. The following is based on book manuscript: B. Han, Framelets Wavelets: Algorithms, Analysis Applications. To introduce a discrete framelet transform,
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 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 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 informationFIR BAND-PASS DIGITAL DIFFERENTIATORS WITH FLAT PASSBAND AND EQUIRIPPLE STOPBAND CHARACTERISTICS. T. Yoshida, Y. Sugiura, N.
FIR BAND-PASS DIGITAL DIFFERENTIATORS WITH FLAT PASSBAND AND EQUIRIPPLE STOPBAND CHARACTERISTICS T. Yoshida, Y. Sugiura, N. Aikawa Tokyo University of Science Faculty of Industrial Science and Technology
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 informationResearch Article Design of Optimal Quincunx Filter Banks for Image Coding
Hindawi Publishing Corporation EURASIP Journal on Advances in Signal Processing Volume 27, Article ID 83858, 8 pages doi:.55/27/83858 Research Article Design of Optimal Quincunx Filter Banks for Image
More informationA Novel Fast Computing Method for Framelet Coefficients
American Journal of Applied Sciences 5 (11): 15-157, 008 ISSN 1546-939 008 Science Publications A Novel Fast Computing Method for Framelet Coefficients Hadeel N. Al-Taai Department of Electrical and Electronic
More informationMagnitude F y. F x. Magnitude
Design of Optimum Multi-Dimensional Energy Compaction Filters N. Damera-Venkata J. Tuqan B. L. Evans Imaging Technology Dept. Dept. of ECE Dept. of ECE Hewlett-Packard Labs Univ. of California Univ. of
More informationLecture 14: Windowing
Lecture 14: Windowing ECE 401: Signal and Image Analysis University of Illinois 3/29/2017 1 DTFT Review 2 Windowing 3 Practical Windows Outline 1 DTFT Review 2 Windowing 3 Practical Windows DTFT Review
More informationDesign of Nonuniform Filter Banks with Frequency Domain Criteria
Blekinge Institute of Technology Research Report No 24:3 Design of Nonuniform Filter Banks with Frequency Domain Criteria Jan Mark de Haan Sven Nordholm Ingvar Claesson School of Engineering Blekinge Institute
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 informationFrames. Hongkai Xiong 熊红凯 Department of Electronic Engineering Shanghai Jiao Tong University
Frames Hongkai Xiong 熊红凯 http://ivm.sjtu.edu.cn Department of Electronic Engineering Shanghai Jiao Tong University 2/39 Frames 1 2 3 Frames and Riesz Bases Translation-Invariant Dyadic Wavelet Transform
More information2D Wavelets for Different Sampling Grids and the Lifting Scheme
D Wavelets for Different Sampling Grids and the Lifting Scheme Miroslav Vrankić University of Zagreb, Croatia Presented by: Atanas Gotchev Lecture Outline 1D wavelets and FWT D separable wavelets D nonseparable
More informationNew Design of Orthogonal Filter Banks Using the Cayley Transform
New Design of Orthogonal Filter Banks Using the Cayley Transform Jianping Zhou, Minh N. Do and Jelena Kovačević Department of Electrical and Computer Engineering University of Illinois at Urbana-Champaign,
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 informationIterative reweighted l 1 design of sparse FIR filters
Iterative reweighted l 1 design of sparse FIR filters Cristian Rusu, Bogdan Dumitrescu Abstract Designing sparse 1D and 2D filters has been the object of research in recent years due mainly to the developments
More informationHilbert Transform Pairs of Wavelets. September 19, 2007
Hilbert Transform Pairs of Wavelets September 19, 2007 Outline The Challenge Hilbert Transform Pairs of Wavelets The Reason The Dual-Tree Complex Wavelet Transform The Dissapointment Background Characterization
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 informationQuadrature Prefilters for the Discrete Wavelet Transform. Bruce R. Johnson. James L. Kinsey. Abstract
Quadrature Prefilters for the Discrete Wavelet Transform Bruce R. Johnson James L. Kinsey Abstract Discrepancies between the Discrete Wavelet Transform and the coefficients of the Wavelet Series are known
More informationDesign of Biorthogonal FIR Linear Phase Filter Banks with Structurally Perfect Reconstruction
Electronics and Communications in Japan, Part 3, Vol. 82, No. 1, 1999 Translated from Denshi Joho Tsushin Gakkai Ronbunshi, Vol. J81-A, No. 1, January 1998, pp. 17 23 Design of Biorthogonal FIR Linear
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 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 informationOrganization of This Pile of Lecture Notes. Part V.F: Cosine-Modulated Filter Banks
Part V.F: Cosine-Modulated s This part shows how to efficiently implement a multichannel (M > ) analysis synthesis bank using a single prototype filter a proper cosine modulation scheme. It is mostly based
More informationESE 531: Digital Signal Processing
ESE 531: Digital Signal Processing Lec 22: April 10, 2018 Adaptive Filters Penn ESE 531 Spring 2018 Khanna Lecture Outline! Circular convolution as linear convolution with aliasing! Adaptive Filters Penn
More information3-D Directional Filter Banks and Surfacelets INVITED
-D Directional Filter Bans and Surfacelets INVITED Yue Lu and Minh N. Do Department of Electrical and Computer Engineering Coordinated Science Laboratory University of Illinois at Urbana-Champaign, Urbana
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 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 informationOptimization in Wireless Communication
Zhi-Quan (Tom) Luo Department of Electrical and Computer Engineering University of Minnesota 200 Union Street SE Minneapolis, MN 55455 2007 NSF Workshop Challenges Optimization problems from wireless applications
More informationx[n] = x a (nt ) x a (t)e jωt dt while the discrete time signal x[n] has the discrete-time Fourier transform x[n]e jωn
Sampling Let x a (t) be a continuous time signal. The signal is sampled by taking the signal value at intervals of time T to get The signal x(t) has a Fourier transform x[n] = x a (nt ) X a (Ω) = x a (t)e
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 11 Adaptive Filtering 14/03/04 http://www.ee.unlv.edu/~b1morris/ee482/
More informationLTI Discrete-Time Systems in Transform Domain Ideal Filters Zero Phase Transfer Functions Linear Phase Transfer Functions
LTI Discrete-Time Systems in Transform Domain Ideal Filters Zero Pase Transfer Functions Linear Pase Transfer Functions Tania Stataki 811b t.stataki@imperial.ac.uk Types of Transfer Functions Te time-domain
More informationWavelets and Filter Banks
Wavelets and Filter Banks Inheung Chon Department of Mathematics Seoul Woman s University Seoul 139-774, Korea Abstract We show that if an even length filter has the same length complementary filter in
More informationAvailable at ISSN: Vol. 2, Issue 2 (December 2007) pp (Previously Vol. 2, No.
Available at http://pvamu.edu.edu/pages/398.asp ISSN: 193-9466 Vol., Issue (December 007) pp. 136 143 (Previously Vol., No. ) Applications and Applied Mathematics (AAM): An International Journal A New
More informationDesign of Nearly Linear-Phase Recursive Digital Filters by Constrained Optimization
Design of Nearly Linear-Phase Recursive Digital Filters by Constrained Optimization by David Guindon B.Eng., University of Victoria, 2001 A Thesis Submitted in Partial Fulfillment of the Requirements for
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 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 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 informationEE364 Review Session 4
EE364 Review Session 4 EE364 Review Outline: convex optimization examples solving quasiconvex problems by bisection exercise 4.47 1 Convex optimization problems we have seen linear programming (LP) quadratic
More informationChapter 7 Wavelets and Multiresolution Processing. Subband coding Quadrature mirror filtering Pyramid image processing
Chapter 7 Wavelets and Multiresolution Processing Wavelet transform vs Fourier transform Basis functions are small waves called wavelet with different frequency and limited duration Multiresolution theory:
More informationLecture 13: Discrete Time Fourier Transform (DTFT)
Lecture 13: Discrete Time Fourier Transform (DTFT) ECE 401: Signal and Image Analysis University of Illinois 3/9/2017 1 Sampled Systems Review 2 DTFT and Convolution 3 Inverse DTFT 4 Ideal Lowpass Filter
More informationOptimum Design of Frequency-Response-Masking Filters Using Convex-Concave Procedure. Haiying Chen
Optimum Design of Frequency-Response-Masking Filters Using Convex-Concave Procedure by Haiying Chen B.Eng., University of Electrical Science and Technology of China, 2013 A Dissertation Submitted in partial
More informationLecture 15: Time and Frequency Joint Perspective
WAVELETS AND MULTIRATE DIGITAL SIGNAL PROCESSING Lecture 15: Time and Frequency Joint Perspective Prof.V.M.Gadre, EE, IIT Bombay Introduction In lecture 14, we studied steps required to design conjugate
More 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 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 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 informationMultidimensional digital signal processing
PSfrag replacements Two-dimensional discrete signals N 1 A 2-D discrete signal (also N called a sequence or array) is a function 2 defined over thex(n set 1 of, n 2 ordered ) pairs of integers: y(nx 1,
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 informationUniversity of Waterloo Department of Electrical and Computer Engineering ECE 413 Digital Signal Processing. Spring Home Assignment 2 Solutions
University of Waterloo Department of Electrical and Computer Engineering ECE 13 Digital Signal Processing Spring 017 Home Assignment Solutions Due on June 8, 017 Exercise 1 An LTI system is described by
More informationEE 4372 Tomography. Carlos E. Davila, Dept. of Electrical Engineering Southern Methodist University
EE 4372 Tomography Carlos E. Davila, Dept. of Electrical Engineering Southern Methodist University EE 4372, SMU Department of Electrical Engineering 86 Tomography: Background 1-D Fourier Transform: F(
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 informationFinite impulse response (FIR) filters have played a central role in digital
[ Timothy N. Davidson ] [ Efficient evaluation of inherent tradeoffs ] Finite impulse response (FIR) filters have played a central role in digital signal processing since its inception. As befits that
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 informationSymmetric Wavelet Tight Frames with Two Generators
Symmetric Wavelet Tight Frames with Two Generators Ivan W. Selesnick Electrical and Computer Engineering Polytechnic University 6 Metrotech Center, Brooklyn, NY 11201, USA tel: 718 260-3416, fax: 718 260-3906
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 information1 The Continuous Wavelet Transform The continuous wavelet transform (CWT) Discretisation of the CWT... 2
Contents 1 The Continuous Wavelet Transform 1 1.1 The continuous wavelet transform (CWT)............. 1 1. Discretisation of the CWT...................... Stationary wavelet transform or redundant wavelet
More informationCoefficients of Recursive Linear Time-Invariant First-Order Low-Pass and High-Pass Filters (v0.1)
Coefficients of Recursive Linear Time-Invariant First-Order Low-Pass and High-Pass Filters (v0. Cliff Sparks www.arpchord.com The following is a quick overview of recursive linear time-invariant first-order
More informationDIGITAL FILTERS Analysis, Design, and Applications by A. Antoniou ERRATA. Page 10, Table 1.1: The upper limit of the summation should be K.
DIGITAL FILTERS Analysis, Design, and Applications by A. Antoniou ERRATA Printing #1 Page vii, line 4 : Replace Geophysicists by Geoscientists. Page 1, Table 1.1: The upper limit of the summation should
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 informationYour solutions for time-domain waveforms should all be expressed as real-valued functions.
ECE-486 Test 2, Feb 23, 2017 2 Hours; Closed book; Allowed calculator models: (a) Casio fx-115 models (b) HP33s and HP 35s (c) TI-30X and TI-36X models. Calculators not included in this list are not permitted.
More informationAnalysis of Fractals, Image Compression and Entropy Encoding
Analysis of Fractals, Image Compression and Entropy Encoding Myung-Sin Song Southern Illinois University Edwardsville Jul 10, 2009 Joint work with Palle Jorgensen. Outline 1. Signal and Image processing,
More informationClosed-Form Design of Maximally Flat IIR Half-Band Filters
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: ANALOG AND DIGITAL SIGNAL PROCESSING, VOL. 49, NO. 6, JUNE 2002 409 Closed-Form Design of Maximally Flat IIR Half-B Filters Xi Zhang, Senior Member, IEEE,
More informationA Width-Recursive Depth-First Tree Search Approach for the Design of Discrete Coefficient Perfect Reconstruction Lattice Filter Bank
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: ANALOG AND DIGITAL SIGNAL PROCESSING, VOL. 50, NO. 6, JUNE 2003 257 A Width-Recursive Depth-First Tree Search Approach for the Design of Discrete Coefficient
More informationCCNY. BME I5100: Biomedical Signal Processing. Linear systems, Fourier and z-transform
BME I5100: Biomedical Signal Processing Linear systems, Fourier and z-transform Lucas C. Parra Biomedical Engineering Department CCNY 1 Schedule Week 1: Introduction Linear, stationary, normal - the stuff
More information