Digital Signal Processing, Fall 2006

Similar documents
2. Finite Impulse Response Filters (FIR)

The Transfer Function. The Transfer Function. The Transfer Function. The Transfer Function. The Transfer Function. The Transfer Function

Linear-Phase FIR Transfer Functions. Functions. Functions. Functions. Functions. Functions. Let

A Simple Method of Tuning PI Controllers for Interval Plant of Cold Rolling Mill

First order differential equation Linear equation; Method of integrating factors

Multiple Short Term Infusion Homework # 5 PHA 5127

Finite Element Analysis

DTFT Properties Using the differentiation property of the DTFT given in Table 3.2, we observe that the DTFT of nx[n] is given by

Discrete-Time Signal Processing

10/5/2012 S. THAI SUBHA CHAPTER-V

the output is Thus, the output lags in phase by θ( ωo) radians Rewriting the above equation we get

[ ] [ ] DFT: Discrete Fourier Transform ( ) ( ) ( ) ( ) Congruence (Integer modulo m) N-point signal

Problem Set #2 Due: Friday April 20, 2018 at 5 PM.

MA1506 Tutorial 2 Solutions. Question 1. (1a) 1 ) y x. e x. 1 exp (in general, Integrating factor is. ye dx. So ) (1b) e e. e c.

ENGR 7181 LECTURE NOTES WEEK 5 Dr. Amir G. Aghdam Concordia University

ANALYSIS IN THE FREQUENCY DOMAIN

EEO 401 Digital Signal Processing Prof. Mark Fowler

Lecture 2: Discrete-Time Signals & Systems. Reza Mohammadkhani, Digital Signal Processing, 2015 University of Kurdistan eng.uok.ac.

5 Transform Analysis of LTI Systems

Image Filtering: Noise Removal, Sharpening, Deblurring. Yao Wang Polytechnic University, Brooklyn, NY11201

2F1120 Spektrala transformer för Media Solutions to Steiglitz, Chapter 1

Y 0. Standing Wave Interference between the incident & reflected waves Standing wave. A string with one end fixed on a wall

Computing and Communications -- Network Coding

Types of Transfer Functions. Types of Transfer Functions. Types of Transfer Functions. Ideal Filters. Ideal Filters

Introduction to Medical Imaging. Lecture 4: Fourier Theory = = ( ) 2sin(2 ) Introduction

Thomas Whitham Sixth Form

Chapter 6. The Discrete Fourier Transform and The Fast Fourier Transform

Simple FIR Digital Filters. Simple FIR Digital Filters. Simple Digital Filters. Simple FIR Digital Filters. Simple FIR Digital Filters

Types of Transfer Functions. Types of Transfer Functions. Ideal Filters. Ideal Filters. Ideal Filters

3) Use the average steady-state equation to determine the dose. Note that only 100 mg tablets of aminophylline are available here.

Numerical methods, Mixed exercise 10

The Matrix Exponential

10. The Discrete-Time Fourier Transform (DTFT)

5. B To determine all the holes and asymptotes of the equation: y = bdc dced f gbd

The Matrix Exponential

2. Transfer function. Kanazawa University Microelectronics Research Lab. Akio Kitagawa

1 Isoparametric Concept

CS553 Lecture Register Allocation I 3

Addition of angular momentum

ECE 410 DIGITAL SIGNAL PROCESSING D. Munson University of Illinois Chapter 12

Addition of angular momentum

Numbering Systems Basic Building Blocks Scaling and Round-off Noise. Number Representation. Floating vs. Fixed point. DSP Design.

Sinusoidal Response Notes

2. Background Material

Digital Signal Processing

EECE 301 Signals & Systems Prof. Mark Fowler

The graph of y = x (or y = ) consists of two branches, As x 0, y + ; as x 0, y +. x = 0 is the

Junction Tree Algorithm 1. David Barber

ECE602 Exam 1 April 5, You must show ALL of your work for full credit.

SIGNIFICANCE OF SMITH CHART IN ANTENNA TECHNOLOGY

Steinberg s Conjecture is false

u x v x dx u x v x v x u x dx d u x v x u x v x dx u x v x dx Integration by Parts Formula

ECE 2210 / 00 Phasor Examples

Problem 22: Journey to the Center of the Earth

Case Study Vancomycin Answers Provided by Jeffrey Stark, Graduate Student

Math 61 : Discrete Structures Final Exam Instructor: Ciprian Manolescu. You have 180 minutes.

A Propagating Wave Packet Group Velocity Dispersion

a 1and x is any real number.

Thomas Whitham Sixth Form

Digital Signal Processing

Frequency Response & Digital Filters

Digital Signal Processing II Lecture 2: FIR & IIR Filter Design

Computer Engineering 4TL4: Digital Signal Processing (Fall 2003) Solutions to Final Exam

Lie Groups HW7. Wang Shuai. November 2015

READING ASSIGNMENTS. Signal Processing First. Problem Solving Skills LECTURE OBJECTIVES. x(t) = cos(αt 2 ) Fourier Series ANALYSIS.

Abstract Interpretation: concrete and abstract semantics

Modeling with first order equations (Sect. 2.3).

Digital Signal Processing Lecture 8 - Filter Design - IIR

nd the particular orthogonal trajectory from the family of orthogonal trajectories passing through point (0; 1).

Solutions to Homework 5

1 General boundary conditions in diffusion

EEO 401 Digital Signal Processing Prof. Mark Fowler

(2) If we multiplied a row of B by λ, then the value is also multiplied by λ(here lambda could be 0). namely

Responses of Digital Filters Chapter Intended Learning Outcomes:

Robust Digital Redesign of Continuous PID Controller for Power system Using Plant-Input-Mapping

Examples and applications on SSSP and MST

Design of IIR filters

ME 321 Kinematics and Dynamics of Machines S. Lambert Winter 2002

ECE Department Univ. of Maryland, College Park

u r du = ur+1 r + 1 du = ln u + C u sin u du = cos u + C cos u du = sin u + C sec u tan u du = sec u + C e u du = e u + C

The second condition says that a node α of the tree has exactly n children if the arity of its label is n.

Grade 12 (MCV4UE) AP Calculus Page 1 of 5 Derivative of a Function & Differentiability

ECE 407 Computer Aided Design for Electronic Systems. Instructor: Maria K. Michael. Overview. CAD tools for multi-level logic synthesis:

Sundials and Linear Algebra

Mathematics 1110H Calculus I: Limits, derivatives, and Integrals Trent University, Summer 2018 Solutions to the Actual Final Examination

Lecture 4: Laplace Transforms

Outline. Image processing includes. Edge detection. Advanced Multimedia Signal Processing #8:Image Processing 2 processing

Case Study 4 PHA 5127 Aminoglycosides Answers provided by Jeffrey Stark Graduate Student

(1) Then we could wave our hands over this and it would become:

CSE 373: AVL trees. Warmup: Warmup. Interlude: Exploring the balance invariant. AVL Trees: Invariants. AVL tree invariants review

AerE 344: Undergraduate Aerodynamics and Propulsion Laboratory. Lab Instructions

Introduction to the Fourier transform. Computer Vision & Digital Image Processing. The Fourier transform (continued) The Fourier transform (continued)

Time Delay Estimation by Bispectrum Interpolation

Math 34A. Final Review

as a derivative. 7. [3.3] On Earth, you can easily shoot a paper clip straight up into the air with a rubber band. In t sec

Integral Calculus What is integral calculus?

Fourier Transforms and the Wave Equation. Key Mathematics: More Fourier transform theory, especially as applied to solving the wave equation.

DSP-First, 2/e. LECTURE # CH2-3 Complex Exponentials & Complex Numbers TLH MODIFIED. Aug , JH McClellan & RW Schafer

Sensors and Actuators Sensor Physics

Mathematics. Complex Number rectangular form. Quadratic equation. Quadratic equation. Complex number Functions: sinusoids. Differentiation Integration

Transcription:

Digital Signal Prossing, Fall 006 Ltur 7: Filtr Dsign Zhng-ua an Dpartmnt of Eltroni Systms Aalborg Univrsity, Dnmar t@om.aau. Cours at a glan MM Disrt-tim signals an systms Systm MM Fourir-omain rprsntation Sampling an ronstrution Systm analysis MM5 Systm strutur MM6 MM4 Filtr sign -transform DF/FF MM7, MM8 MM3 MM9, MM0

Part I: Filtr sign Filtr sign IIR filtr sign Analog filtr sign IIR filtr sign by impuls invarian IIR filtr sign by bilinar transformation 3 Filtr sign pross Filtr, in broar sns, ovrs any systm. hr sign stps Problm Spifiations Prforman onstraints Approximations Systm funtion raliation Solution Magnitu rspons Phas rspons (frquny omain Complxity IIR or FIR Subtyp Strutur 4

Spifiations an xampl Spifiations for a isrt-tim lowpass filtr 0.0 ( ( 0.00, + 0.0, ω ω s 0 ω ω p δ 0.0 δ 0.00 5 Spifiations of frquny rspons ypial lowpass filtr spifiations in trms of tolrabl Passban istortion, as smallst as possibl Stopban attnuation, as gratst as possibl With of transition ban: as narrowst as possibl Improving on oftn worsns othrs a traoff Inrasing filtr orr improvs all 6 3

D filtr for C signals Disrt-tim filtr for th prossing of ontinuoustim signals Banlimit input signal igh nough sampling frquny hn, spifiations onvrsion is straightforwar jω (, Ω < π / ff ( jω 0, Ω > π / ω ( ff ( j, ω < π ω Ω Fig. 7. 7 Spifiations an xampl Spifiations for a ontinuous-tim lowpass filtr 0.0 ff ff ( jω + 0.0, 0 Ω π (000 ( jω 0.00, Ω π (3000 Fig 7.(a(b δ 0.0 δ 0.00 Ω Ω p s π (000 π (3000 8 4

Dsign a filtr Dsign goal: fin systm funtion to ma frquny rspons mt th spifiations (tolrans Infinit impuls rspons filtr Pols insir unit irl u to ausality an stability Rational funtion approximation Finit impuls rspons filtr Linar phas is oftn rquir Polynomial approximation 9 E.g. IIR filtr sign For rational systm funtion ( fin th systm offiints suh that th orrsponing frquny rspons ( ( provis a goo approximation to a sir rspons ( ( sir M 0 N b a ( Rational systm funtion Stabl ausal 0 5

FIR or IIR Eithr FIR or IIR is oftn pnnt on th phas rquirmnts Only FIR filtr an b at th sam tim stabl, ausal an GLP If ( is stabl an GLP, any non-trivial pol p insi th unit irl orrspons a pol /p outsi th unit irl, so that ( annot hav a ausal impuls rspons (as ROC is a ring inluing unit irl. Dsign prinipl If GLP is ssntial FIR If not IIR prfrabl (an mt spifiations with lowr omplxity FIR an IIR IIR Rational systm funtion Pols + ros Stabl/unstabl ar to ontrol phas Low orr (4-0 Dsign on th basis of analog filtr FIR Polynomial systm funtion Zros Stabl Easy to gt linar phas igh orr (0-000 Unrlat to analog filtr 6

Part II: IIR filtr sign Filtr sign IIR filtr sign Analog filtr sign IIR filtr sign by impuls invarian IIR filtr sign by bilinar transformation 3 Dsign IIR filtr bas on analog filtr h mapping is irt jω (, ff ( jω 0, ω ( ff ( j, Ω < π / Ω > π / ω < π Avan analog filtr sign thniqus Dsigning D filtr by transforming prototyp C filtr: ransform (map D spifiations to analog Dsign analog filtr Invrs-transform analog filtr to D 4 7

ransformation mtho ransform (map D spifiations to analog ω Ω Dsign analog filtr ( s or h ( t Invrs-transform to D ( or h[ n] h imaginary axis of th s-plan th unit irl of th -plan Pols in th lft half of th s-plan pols insi th unit irl in th - plan (stabl s σ + jω ( s r X ( ( jω X ( n n h( t st h( t x[ n] x[ n] jωt n t t n 5 Part III: Analog filtr sign Filtr sign IIR filtr sign Analog filtr sign IIR filtr sign by impuls invarian IIR filtr sign by bilinar transformation 6 8

Analog filtr sign Buttrworth Chbyshv I Chbyshv II Ellipial 7 Buttrworth lowpass filtrs h magnitu rspons Maximally flat in th passban Monotoni in both passban an stopban h squar magnitu rspons ( jω N + ( Ω / Ω FigB.. 8 9

Pols in s-plan 9 Chbyshv filtrs Chbyshv II filtrs: quirippl in passban, flat in th stopban Chbyshv II filtrs: quirippl in stopban, flat in th passban 0 0

Ellipti filtrs Equirippl both in stopban an in th passban Part IV: Dsign by impuls invarian Filtr sign IIR filtr sign Analog filtr sign IIR filtr sign by impuls invarian IIR filtr sign by bilinar transformation

Filtr sign by impuls invarian Impuls invarian: a mtho for obtaining a D systm whos ( is trmin by th ( jω of a C systm. h [ n] h ( n - 'sign' sampling intrval In D filtr sign, th spifiations ar provi in th isrt-tim, so has no rol. is inlu for isussion though. also has nothing to o with C/D an D/C onvrsion in Fig. 7. 3 Rlationship btw frquny rsponss Impuls rspons sampling: h [ n] h ( n Frquny rspons ω π ( ( j + j if th C filtr is banlimit ( jω 0, Ω π / thn ( ω ( j, ω π his is also th way to gt C filtr spifiations from ( by applying th rlation Ω ω / 4 x[ n] x ( n ω π X ( X ( j j

3 5 Aliasing in th impuls invarian sign ( ( j j j π ω ω + ( ] [ n h n h 6 Rlationship btw systm funtions h transform from C to D is asy to arry out as a transformation on th systm funtion Rational systm funtion, aftr partial fration xpansion N n s N n s n u A n u A n h n h ] [ ( ] [ ( ] [ < 0 0, 0, ( ( t t A t h s s A s N t s N N s A (

Impuls invarian with a Buttrworth filtr Spifiations 0.895 ( (, 0 ω 0.π 0.7783, 0.3π ω π Sin th sampling intrval anls in th impuls invarian prour, w hoos, so ω Ω Magnitu funtion for a C Buttrworth filtr 0.895 ( jω, 0 Ω 0.π ( jω 0.7783, 0.3π Ω π Du to th monotoni funtion of Buttrworth filtr ( j0.π 0.895 ( j0.3π 0.7783 7 Impuls invarian with a Buttrworth filtr ( (3 Squar magnitu funtion of a Buttrworth filtr ( Ω j N + ( Ω / Ω 0.π + ( Ω 0.3π + ( Ω N N N 5.8858 Ω ( ( 0.70474 0.895 0.7783 (4 (5 ( ( s N 6 Ω ( j0.π 0.895 ( j0.3π 0.7783 0.703 ( s + ( s / jω + ( s / j0.703 N 8 4

Impuls invarian with a Buttrworth filtr pols for th squar magnitu funtion h systm funtion has th thr pol pairs in th lft half of th s-plan 9 Impuls invarian with a Buttrworth filtr s ( s + 0.3640s + 0.4945( s 0.093 + 0.9945s + 0.4945( s ( 0.87 0.4466 ( (.97 + 0.6949.48 +.455 + (.069 + 0.3699.8557 0.6303 + ( 0.997 + 0.570 +.3585s + 0.4945 30 5

Impuls invarian with a Buttrworth filtr 3 Part V: Dsign by bilinar transformation Filtr sign IIR filtr sign Analog filtr sign IIR filtr sign by impuls invarian IIR filtr sign by bilinar transformation 3 6

Bilinar transformation By using impuls invarian, th rlation btwn C an D frquny is linar (xpt for aliasing, thus th shap of th frquny rspons is prsrv. But only propr for banlimit filtrs, problm for.g. highpass Bilinar transformation btwn s an s ( ( [ ( ] + + Invrs + ( ( / s / s 33 Bilinar transformation Givn s σ + jω if s jω + ( ( if σ < 0, if σ > 0, / s + σ / s σ < for any > for any / / Ω Ω + jω jω / / + jω jω so,, on th jω - axis i.. th / / for any s jω - axis maps onto th unit irl 34 7

Bilinar transformation frquny rlationship Consir frquny s ( + + s s jω jω jω jω ( + Ω tan( ω / ω artan( Ω [ / / / ( j sinω / j ] tan( ω / (osω / 35 Bilinar transformation h bilinar transformation maps th ntir jω -axis in th s-plan to on rvolution of th unit irl in th -plan. Compar with Ω ω / 36 8

Bilinar transformation of a Buttrworth filtr Spifiations 0.895 ( (, 0 ω 0.π 0.7783, 0.3π ω π Magnitu funtion for a C Buttrworth filtr 0.π 0.895 ( jω, 0 Ω tan( 0.3π ( jω 0.7783, tan( Ω Du to th monotoni funtion of Buttrworth filtr Choos ( j tan(0.π 0.895 ( j tan(0.5π 0.7783 37 Bilinar transformation of a Buttrworth filtr Squar magnitu funtion of a Buttrworth filtr ( Ω j N + ( Ω / Ω N 5.305 ( j tan(0.π 0.895 ( j tan(0.5π 0.7783 N 6 Ω 0.766 38 9

Bilinar transformation of a Buttrworth filtr 39 Summary Filtr sign IIR filtr sign Analog filtr sign IIR filtr sign by impuls invarian IIR filtr sign by bilinar transformation 40 0

Cours at a glan MM Disrt-tim signals an systms Systm MM Fourir-omain rprsntation Sampling an ronstrution Systm analysis MM5 Systm strutur MM6 MM4 Filtr sign -transform DF/FF MM7, MM8 MM3 MM9, MM0 4