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

Similar documents
Types of Transfer Functions. Types of Transfer Functions. Ideal Filters. Ideal Filters. Ideal Filters

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

10. The Discrete-Time Fourier Transform (DTFT)

Stability Condition in Terms of the Pole Locations

5 Transform Analysis of LTI Systems

Slide 1. Slide 2. Slide 3 DIGITAL SIGNAL PROCESSING CLASSIFICATION OF SIGNALS

Exercise 1. Sketch the graph of the following function. (x 2

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

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

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

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

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

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

u 3 = u 3 (x 1, x 2, x 3 )

ANALYSIS IN THE FREQUENCY DOMAIN

The Matrix Exponential

Hydrogen Atom and One Electron Ions

Higher order derivatives

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

SECTION where P (cos θ, sin θ) and Q(cos θ, sin θ) are polynomials in cos θ and sin θ, provided Q is never equal to zero.

Search sequence databases 3 10/25/2016

The Matrix Exponential

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

Introduction to Arithmetic Geometry Fall 2013 Lecture #20 11/14/2013

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

BINOMIAL COEFFICIENTS INVOLVING INFINITE POWERS OF PRIMES. 1. Statement of results

BINOMIAL COEFFICIENTS INVOLVING INFINITE POWERS OF PRIMES

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

1 General boundary conditions in diffusion

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

A Propagating Wave Packet Group Velocity Dispersion

Bifurcation Theory. , a stationary point, depends on the value of α. At certain values

Discrete-Time Signal Processing

EEO 401 Digital Signal Processing Prof. Mark Fowler

Einstein Equations for Tetrad Fields

Communication Technologies

Introduction to Condensed Matter Physics

DISCRETE TIME FOURIER TRANSFORM (DTFT)

EE140 Introduction to Communication Systems Lecture 2

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

COMPUTER GENERATED HOLOGRAMS Optical Sciences 627 W.J. Dallas (Monday, April 04, 2005, 8:35 AM) PART I: CHAPTER TWO COMB MATH.

Problem Set 6 Solutions

Solution: APPM 1360 Final (150 pts) Spring (60 pts total) The following parts are not related, justify your answers:

Math 34A. Final Review

Addition of angular momentum

ECE 344 Microwave Fundamentals

Combinatorial Networks Week 1, March 11-12

Discrete Hilbert Transform. Numeric Algorithms

22/ Breakdown of the Born-Oppenheimer approximation. Selection rules for rotational-vibrational transitions. P, R branches.

Elements of Statistical Thermodynamics

Ch. 24 Molecular Reaction Dynamics 1. Collision Theory

Supplementary Materials

ECE 2210 / 00 Phasor Examples

EEO 401 Digital Signal Processing Prof. Mark Fowler

Quasi-Classical States of the Simple Harmonic Oscillator

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

COHORT MBA. Exponential function. MATH review (part2) by Lucian Mitroiu. The LOG and EXP functions. Properties: e e. lim.

1.2 Faraday s law A changing magnetic field induces an electric field. Their relation is given by:

First derivative analysis

Differential Equations

Collisions between electrons and ions

There is an arbitrary overall complex phase that could be added to A, but since this makes no difference we set it to zero and choose A real.

NEW APPLICATIONS OF THE ABEL-LIOUVILLE FORMULA

Thus, because if either [G : H] or [H : K] is infinite, then [G : K] is infinite, then [G : K] = [G : H][H : K] for all infinite cases.

cycle that does not cross any edges (including its own), then it has at least

10. Limits involving infinity

Addition of angular momentum

A. Limits and Horizontal Asymptotes ( ) f x f x. f x. x "±# ( ).

Title: Vibrational structure of electronic transition

5.80 Small-Molecule Spectroscopy and Dynamics

Calculus II (MAC )

CS 361 Meeting 12 10/3/18

LTI Discrete-Time Systems in Transform Domain Ideal Filters Zero Phase Transfer Functions Linear Phase Transfer Functions

1973 AP Calculus AB: Section I

Propositional Logic. Combinatorial Problem Solving (CPS) Albert Oliveras Enric Rodríguez-Carbonell. May 17, 2018

3 Finite Element Parametric Geometry

MA 262, Spring 2018, Final exam Version 01 (Green)

2. Finite Impulse Response Filters (FIR)

(Upside-Down o Direct Rotation) β - Numbers

Chapter 10. The singular integral Introducing S(n) and J(n)

Volterra Kernel Estimation for Nonlinear Communication Channels Using Deterministic Sequences

The Frequency Response of a Quarter-Wave Matching Network

That is, we start with a general matrix: And end with a simpler matrix:

Outline. Thanks to Ian Blockland and Randy Sobie for these slides Lifetimes of Decaying Particles Scattering Cross Sections Fermi s Golden Rule

Brief Notes on the Fermi-Dirac and Bose-Einstein Distributions, Bose-Einstein Condensates and Degenerate Fermi Gases Last Update: 28 th December 2008

4.2 Design of Sections for Flexure

EECE 301 Signals & Systems Prof. Mark Fowler

1 Isoparametric Concept

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

Equidistribution and Weyl s criterion

CHAPTER 1. Introductory Concepts Elements of Vector Analysis Newton s Laws Units The basis of Newtonian Mechanics D Alembert s Principle

2.3 Matrix Formulation

Alpha and beta decay equation practice

Objective Mathematics

4. Money cannot be neutral in the short-run the neutrality of money is exclusively a medium run phenomenon.

Contemporary, atomic, nuclear, and particle physics

Math 102. Rumbos Spring Solutions to Assignment #8. Solution: The matrix, A, corresponding to the system in (1) is

First order differential equation Linear equation; Method of integrating factors

6. The Interaction of Light and Matter

Digital Signal Processing, Fall 2006

Transcription:

Typs of Transfr Typs of Transfr x[n] X( LTI h[n] H( y[n] Y( y [ n] h[ k] x[ n k] k Y ( H ( X ( Th tim-domain classification of an LTI digital transfr function is basd on th lngth of its impuls rspons h[n]: - Finit impuls rspons (FIR transfr function - Infinit impuls rspons (IIR transfr function Typs of Transfr In th cas of digital transfr functions with frquncy-slctiv frquncy rsponss, thr ar two typs of classifications ( Classification basd on th shap of th magnitud function H ( j ( Classification basd on th th form of th phas function θ( Classification Basd on agnitud Charactristics On common classification is basd on an idal magnitud rspons A digital filtr dsignd to pass signal componnts of crtain frquncis without distortion should hav a magnitud rspons qual to on at ths frquncis, and should hav a magnitud rspons qual to zro at all othr frquncis 3 4 Th rang of frquncis whr th magnitud rspons taks th valu of on is calld th passband Th rang of frquncis whr th magnitud rspons taks th valu of zro is calld th stopband Frquncy rsponss of th four popular typs of idal digital filtrs with ral impuls rspons cofficints ar shown blow: π H LP c ( c Lowpass H BP ( π H HP ( π c c π Highpass H BS ( 5 6 π c c c c π π c c c c π Bandpass Bandstop

7 Lowpass filtr: Passband - Stopband - Highpass filtr: Passband - Stopband - c < π c π < c Bandpass filtr: Passband - c c Stopband - < c and c < π Bandstop filtr: Stopband - c < < c Passband - c and c π c 8 c Th frquncis c, c, and ar calld th cutoff frquncis An idal filtr has a magnitud rspons qual to on in th passband and zro in th stopband, and has a zro phas vrywhr 9 Earlir in th cours w drivd th invrs DTFT of th frquncy rspons H ( j LP of th idal lowpass filtr: sin n h n c LP [ ], < n < πn W hav also shown that th abov impuls rspons is not absolutly summabl, and hnc, th corrsponding transfr function is not BIBO stabl Also, h LP [n] is not causal and is of doubly infinit lngth Th rmaining thr idal filtrs ar also charactrizd by doubly infinit, noncausal impuls rsponss and ar not absolutly summabl Thus, th idal filtrs with th idal brick wall frquncy rsponss cannot b ralizd with finit dimnsional LTI filtr To dvlop stabl and ralizabl transfr functions, th idal frquncy rspons spcifications ar rlaxd by including a transition band btwn th passband and th stopband This prmits th magnitud rspons to dcay slowly from its maximum valu in th passband to th zro valu in th stopband orovr, th magnitud rspons is allowd to vary by a small amount both in th passband and th stopband Typical magnitud rspons spcifications of a lowpass filtr ar shown blow G( + δp _ δp δs Passband p s Transition band Stopband π

3 A causal stabl ral-cofficint transfr function H( is dfind as a boundd ral (BR transfr function if H ( for all valus of Lt x[n] and y[n] dnot, rspctivly, th input and output of a digital filtr charactrizd by a BR transfr function H( with X ( and Y ( dnoting thir DTFTs 4 Thn th condition H ( implis that Y ( Intgrating th abov from π to π, and applying Parsval s thorm w gt n y[ n] X ( n x[ n] Thus, for all finit-nrgy inputs, th output nrgy is lss than or qual to th input nrgy implying that a digital filtr charactrizd by a BR transfr function can b viwd as a passiv structur If H (, thn th output nrgy is qual to th input nrgy, and such a digital filtr is thrfor a losslss systm A causal stabl ral-cofficint transfr function H( with H ( is thus calld a losslss boundd ral (LBR transfr function Th BR and LBR transfr functions ar th kys to th ralization of digital filtrs with low cofficint snsitivity 5 6 Exampl Considr th causal stabl IIR transfr function H( K, < α < αz whr K and α ar ral constants Its squar-magnitud function is givn by K H( H( H( z z ( + α αcos Th maximum valu of H( is obtaind whn αcos in th dnominator is a maximum and th minimum valu is obtaind whn αcos is a minimum For α >, maximum valu of α cos is qual to α at, and minimum valu is α at π 7 8 3

Thus, for α >, th maximum valu of H( is qual to K /( α at and th minimum valu is qual to K /( + α at π On th othr hand, for α <, th maximum valu of αcos is qual to α at π and th minimum valu is qual to α at Hr, th maximum valu of H( is qual to K /( α at π and th minimum valu is qual to K / ( + α at Hnc, th maximum valu can b mad qual to by choosing K ± ( α, in which cas th minimum valu bcoms ( α /( + α 9 Hnc, K H( αz < α < is a BR function for K ± ( α Plots of th magnitud function for α ±.5 with valus of K chosn to mak H( a BR function ar shown on th nxt slid, agnitud.8.6.4 K ±.5, α.5...4.6.8 Lowpass Filtr agnitud.8.6.4 K ±.5, α -.5...4.6.8 Highpass Filtr 3 Dfinition An IIR transfr function A( with unity magnitud rspons for all frquncis, i.., A(, for all is calld an allpass transfr function An-th ordr causal ral-cofficint allpass transfr function is of th form... + d + d z d z z z + + + ± A (... + + d z + + d z + d z 4 If w dnot th dnominator polynomial of A ( as D (: + D ( + d z +... + d z + d z thn it follows that A ( can b writtn as: z D ( z A ( ± D ( Not from th abov that if z r jφ is a pol of a ral cofficint allpass transfr function, thn it has a zro at φ z j r 4

5 Th numrator of a ral-cofficint allpass transfr function is said to b th mirrorimag polynomial of th dnominator, and vic vrsa ~ W shall us th notation D ( to dnot th mirror-imag polynomial of a dgr- polynomial D (, i.., ~ D ( z D ( z 6 Th xprssion z D ( z A ( ± D ( implis that th pols and zros of a ralcofficint allpass function xhibit mirrorimag symmtry in th z-plan 3. +. 8z +. + z A3( 3 +. +. 8z. z Imaginary Part.5.5 -.5 - -.5-3 Ral Part 7 To show that A Thrfor A j A z D ( ( z ± D ( z ( w obsrv that z D ( z z D ( ( A ( z D ( D ( z Hnc A ( A ( A ( z z 8 Now, th pols of a causal stabl transfr function must li insid th unit circl in th z-plan Hnc, all zros of a causal stabl allpass transfr function must li outsid th unit circl in a mirror-imag symmtry with its pols situatd insid th unit circl Figur blow shows th principal valu of th phas of th 3rd-ordr allpass function 3. +. 8z +. + z A 3( 3 +. +. 8z. z Not th discontinuity by th amount of π in th phas θ( 4 Principal valu of phas If w unwrap th phas by rmoving th discontinuity, w arriv at th unwrappd phas function θ c ( indicatd blow Not: Th unwrappd phas function is a continuous function of Unwrappd phas - Phas, dgrs - Phas, dgrs -4-6 -8 9-4..4.6.8 3 -..4.6.8 5

3 Th unwrappd phas function of any arbitrary causal stabl allpass function is a continuous function of Proprtis (A causal stabl ral-cofficint allpass transfr function is a losslss boundd ral (LBR function or, quivalntly, a causal stabl allpass filtr is a losslss structur 3 (Th magnitud function of a stabl allpass function A( satisfis: <, for z > A(, for z >, for z < (3Lt τ( dnot th group dlay function of an allpass filtr A(, i.., d d c τ( [ θ ( ] 33 Th unwrappd phas function θ c ( of a stabl allpass function is a monotonically dcrasing function of so that τ( is vrywhr positiv in th rang < < π Th group dlay of an -th ordr stabl ral-cofficint allpass transfr function satisfis: π τ( d π 34 A Simpl Application A simpl but oftn usd application of an allpass filtr is as a dlay qualizr LtG( b th transfr function of a digital filtr dsignd to mt a prscribd magnitud rspons Th nonlinar phas rspons of G( can b corrctd by cascading it with an allpass filtr A( so that th ovrall cascad has a constant group dlay in th band of intrst G( j A( Sinc A(, w hav G( A( G( Ovrall group dlay is th givn by th sum of th group dlays of G( and A( Exampl Figur blow shows th group dlay of a 4 th ordr lliptic filtr with th following spcifications: p. 3π, δ p db, 35 db Group dlay, sampls 5 5 δ s Original Filtr..4.6.8 35 36 6

Figur blow shows th group dlay of th original lliptic filtr cascadd with an 8 th ordr allpass sction dsignd to qualiz th group dlay in th passband Group dlay, sampls 3 5 5 5 Group Dlay Equalizd Filtr..4.6.8 Classification Basd on Phas Charactristics A scond classification of a transfr function is with rspct to its phas charactristics In many applications, it is ncssary that th digital filtr dsignd dos not distort th phas of th input signal componnts with frquncis in th passband 37 38 Zro-Phas Transfr Zro-Phas Transfr On way to avoid any phas distortion is to mak th frquncy rspons of th filtr ral and nonngativ, i.., to dsign th filtr with a zro phas charactristic Howvr, it is not possibl to dsign a causal digital filtr with a zro phas For non-ral-tim procssing of ral-valud input signals of finit lngth, zro-phas filtring can b vry simply implmntd by rlaxing th causality rquirmnt On zro-phas filtring schm is sktchd blow x[n] H( v[n] u[n] H( w[n] 39 4 u[ n] v[ n], y[ n] w[ n] Zro-Phas Transfr Zro-Phas Transfr 4 It is asy to vrify th abov schm in th frquncy domain Lt X (, V ( j, U ( j, W ( j, and Y ( dnot th DTFTs of x[n], v[n], u[n], w[n], andy[n], rspctivly From th figur shown arlir and making us of th symmtry rlations w arriv at th rlations btwn various DTFTs as givn on th nxt slid 4 x[n] H( v[n] u[n] H( w[n] ( j u[ n] v[ n], y[ n] w[ n] ( j ( j W*( j V H ( X (, W H ( U ( ( j j U V*(, Y Combining th abov quations w gt Y ( W*( H *( U *( H ( j V ( H*( H( X( * H ( X ( 7

Zro-Phas Transfr 43 Th function filtfilt implmnts th abov zro-phas filtring schm In th cas of a causal transfr function with a nonzro phas rspons, th phas distortion can b avoidd by nsuring that th transfr function has a unity magnitud and a linar-phas charactristic in th frquncy band of intrst 44 Th most gnral typ of a filtr with a linar phas has a frquncy rspons givn by H ( D which has a linar phas from to π Not also H ( τ ( D 45 Th output y[n] of this filtr to an input x[ n] A j n is thn givn by y[ n] A D n A ( n D If x a (t and y a (t rprsnt th continuoustim signals whos sampld vrsions, sampld at t nt,arx[n] and y[n] givn abov, thn th dlay btwn x a (t and y a (t is prcisly th group dlay of amount D 46 If D is an intgr, thn y[n] is idntical to x[n], but dlayd by D sampls If D is not an intgr, y[n], bing dlayd by a fractional part, is not idntical to x[n] In th lattr cas, th wavform of th undrlying continuous-tim output is idntical to th wavform of th undrlying continuous-tim input and dlayd D units of tim If it is dsird to pass input signal componnts in a crtain frquncy rang undistortd in both magnitud and phas, thn th transfr function should xhibit a unity magnitud rspons and a linar-phas rspons in th band of intrst Figur blow shows th frquncy rspons if a lowpass filtr with a linar-phas charactristic in th passband H LP ( π c π c arg H LP ( 47 48 π c c π 8

49 Sinc th signal componnts in th stopband ar blockd, th phas rspons in th stopband can b of any shap Exampl-Dtrmin th impuls rspons of an idal lowpass filtr with a linar phas rspons: H ( LP no,, < < c π c 5 Applying th frquncy-shifting proprty of th DTFT to th impuls rspons of an idal zro-phas lowpass filtr w arriv at sin n n h n c( o LP [ ], < n < π( n n As bfor, th abov filtr is noncausal and of doubly infinit lngth, and hnc, unralizabl o By truncating th impuls rspons to a finit numbr of trms, a ralizabl FIR approximation to th idal lowpass filtr can b dvlopd Th truncatd approximation may or may not xhibit linar phas, dpnding on th valu of chosn n o n o If w choos N/ with N a positiv intgr, th truncatd and shiftd approximation ^ sin n N h n c( / LP[ ], n N π( n N / will b a lngth N+ causal linar-phas FIR filtr 5 5 53 Figur blow shows th filtr cofficints obtaind using th function sinc for two diffrnt valus of N Amplitud.6.4. N -. 4 6 8 Tim indx n Amplitud.6.4. N 3 -. 4 6 8 Tim indx n 54 Zro-Phas Rspons Bcaus of th symmtry of th impuls rspons cofficints as indicatd in th two figurs, th frquncy rspons of th truncatd approximation can b xprssd as: ^ H LP N ^ j n LP LP n ( h [ n] N/ H whr H LP(, calld th zro-phas rspons or amplitud rspons, is a ral function of ( 9

55 inimum-phas and aximum- Phas Transfr Considr th two st-ordr transfr functions: bz+ z z b, H ( < < + +, H (, a b a z+ a Both transfr functions hav a pol insid th unit circl at th sam location z a and ar stabl But th zro of H ( is insid th unit circl at z b, whras, th zro of H ( is at z situatd in a mirror-imag symmtry b 56 inimum-phas and aximum- Phas Transfr Figur blow shows th pol-zro plots of th two transfr functions H ( H ( inimum-phas and aximum- Phas Transfr inimum-phas and aximum- Phas Transfr 57 Howvr, both transfr functions hav an idntical magnitud function as H( H( z H( H( z Th corrsponding phas functions ar sin sin tan b+ cos a+ cos arg[ H( ] tan bsin sin tan + bcos a+ cos arg[ H( ] tan 58 Figur blow shows th unwrappd phas rsponss of th two transfr functions for a.8 and b. 5 Phas, dgrs - - -3 H ( H ( -4..4.6.8 59 inimum-phas and aximum- Phas Transfr From this figur it follows that H ( has an xcss phas lag with rspct to H ( Th xcss phas lag proprty of H ( with rspct to H ( can also b xplaind by obsrving that w can writ bz + z + b bz + H ( z + a z + a z + b H ( z A ( z 6 inimum-phas and aximum- Phas Transfr whr A( z ( bz + / ( z + b is a stabl allpass function Th phas functions of H ( and H ( ar thus rlatd through arg[ H( ] arg[ H( ] + arg[ A( ] As th unwrappd phas function of a stabl first-ordr allpass function is a ngativ function of, it follows from th abov that H ( has indd an xcss phas lag with rspct to H (

inimum-phas and aximum- Phas Transfr inimum-phas and aximum- Phas Transfr Gnralizing th abov rsult, lt H m ( b a causal stabl transfr function with all zros insid th unit circl and lt H( b anothr causal stabl transfr function j satisfying ( H Hm( Ths two transfr functions ar thn rlatd through H ( Hm( A( whr A( is a causal stabl allpass function Th unwrappd phas functions of H m ( and H( ar thus rlatd through arg[ H ( ] arg[ Hm( ] + arg[ A( ] H( has an xcss phas lag with rspct to H m ( A causal stabl transfr function with all zros insid th unit circl is calld a minimum-phas transfr function 6 6 63 inimum-phas and aximum- Phas Transfr A causal stabl transfr function with all zros outsid th unit circl is calld a maximum-phas transfr function A causal stabl transfr function with zros insid and outsid th unit circl is calld a mixd-phas transfr function inimum-phas and aximum- Phas Transfr Exampl Considr th mixd-phas transfr function ( +. 3z (. 4 z H ( (. z ( +. 5z W can rwrit H( as ( +. 3z (. H ( (. z ( +. 5z. 4 z. 64 inimum-phas function Allpass function