Implementing IIR/FIR Filters

More documents

Recommendations

Info

Difference Equation yn ( ) = 2{ α[ x( n) + 2x( n – 1) + xn ( – 2)] + γyn ( – 1) – βyn ( – 2) } X:(R0) 2α α γ – β α Data Structures Y:(R4) Y:(R5) Total Instruction Cycles 6 Icyc @ 20 MHz = 600ns x(n-1) x(n-2) y(n-1) y(n-2) DSP56001 Code ;Y1=x(n) (Input) ;X0=α MPY MAC MAC MAC MAC MOVE X0,Y1,A X0,Y0,A X0,Y0,A X0,Y0,A X0,Y0,A X:(R0)+,X0 X:(R0)+,X0 X:(R0)+,X0 X:(R0)+,X0 X:(R0)+,X0 A,X1 Y:(R4)+,Y0 Y:(R4),Y0 Y:(R5)+,Y0 Y:(R5),Y0 Y1,Y:(R4) A,Y:(R5) ;A=αx(n) ;A=A+2αx(n-1) ;A=A+αx(n-2) ;A=A+γy(n-1) ;A=A-βy(n-2) ;y(n-2)=2A (assumes scaling ;mode is set). ;X1 is Output. Figure 2-12 Figure 2-12 DSP56001 Code and Data Structures for Second-Order Direct-Form Implementation of a Lowpass <strong>IIR</strong> Filter MOTOROLA 2-15
Network Diagram Transfer Function Gain Phase x(n) Coefficients z -1 z -1 G( θ) Hz ( ) α 1– 2z 1 – z 2 – ( + ) 1/2 γz 1 – – βz 2 – = -------------------------------------------- + ( 1 – cosθ) ( 1 + cosθc) ( dsinθsinθc) 2 4( cosθ – cosθc) 2 [ + ] 1/2 = -------------------------------------------------------------------------------------------------- φθ ( ) α −2α α Σ – 1 π+ tan 2 θ θ ( cos – cos c) --------------------------------------dsinθsinθc tan 1 – 2 θ θ ⎧ ⎪ ⎪ = ⎨ ⎪ ( cos – cos c) ⎪ ---------------------------------------- ⎩ dsinθsinθc 1 β ⎛--⎞ 1 d 2 ⎝2⎠ ⁄ θ – sin c = --------------------------------- 1+ d ⁄ 2sinθc γ = ( 1/2 + β) cosθc α = ( 1/2 + β+ γ) ⁄ 4 for θ≤θ c for θ>θ c Figure 2-13 Direct-Form Implementation of Second Order Highpass <strong>IIR</strong> Filter and Analytical Formulas Relating Desired Response to Filter Coefficients 2-16 MOTOROLA 2 γ −β z -1 z -1 y(n)
Page 1 and 2: Motorola’s High-Performance DSP T
Page 3 and 4: © Motorola Inc. 1993 Motorola rese
Page 5 and 6: SECTION 6 Filter Design and Analysi
Page 7 and 8: Figure 2-3 Figure 2-4 Figure 2-5 Fi
Page 9 and 10: Figure 4-4 Figure 4-5 Figure 5-1 Fi
Page 11 and 12: Figure 7-9 Figure 7-10 Figure 7-11
Page 13 and 14: equivalent functions in the digital
Page 15 and 16: is taken to calculating the filter
Page 17 and 18: solving the differential equation i
Page 19 and 20: Gain Phase 1.5 1.0 0.5 0 -π/2 -π
Page 21 and 22: 1.3 Analog Highpass Filter The pass
Page 23 and 24: Gain Phase 1.5 1.0 0.5 π π/2 0 0
Page 25 and 26: V 0 R ----- = ---------------------
Page 27 and 28: V 0 R ----- = ---------------------
Page 29 and 30: 1-18 MOTOROLA
Page 31 and 32: In the digital domain, the continuo
Page 33 and 34: One mapping or transformation from
Page 35 and 36: Although this transformation was de
Page 37 and 38: property in particular as follows:
Page 39 and 40: Figure 2-10 are similar to the anal
Page 41 and 42: case but also in the code for the h
Page 43: Lowpass Gain Lowpass Phase 1.5 1.0
Page 47 and 48: MPY MAC MAC MAC MAC Difference Equa
Page 49 and 50: Network Diagram x(n) Transfer Funct
Page 51 and 52: BANDSTOP GAIN BANDSTOP PHASE 1.5 1.
Page 53 and 54: Network Diagram x(n) Transfer Funct
Page 55 and 56: BANDPASS GAIN BANDPASS PHASE 1.5 1.
Page 57 and 58: Pole Equation of H(z) Zp = rcosθp
Page 59 and 60: x(n) b 0 + + Figure 3-24 The Second
Page 61 and 62: The last step uses the definition f
Page 63 and 64: where: θ p is the angle (see Figur
Page 65 and 66: Since this result does not depend o
Page 67 and 68: MPY MAC MACR MAC MACR Difference Eq
Page 71 and 72: expression for the maximum gain at
Page 73 and 74: GAIN AT INTERNAL NODE, u(n) Transfe
Page 75 and 76: Gain Gain Gain 1.2 1.0 0.8 0.6 0.4
Page 79 and 80: Hs ( ) damping factors, dk, of each
Page 81 and 82: Using , Eqn. 5-7 becomes: Hs ( ) N/
Page 83 and 84: α k = β k = γ k = tan 2 ( θ c /
Page 85 and 86: pass Butterworth filter (three seco
Page 87 and 88: Network Diagram x(n) α1 z -1 z -1
Page 91 and 92: can be assembled by the DSP56001 as
Page 93 and 94: Phase (Radians) π π/2 0 -π/2 -π
Page 95 and 96:
Seconds 6.3 E-03 4.7 E-03 3.1 E-03
Page 97 and 98:
ZEROES OF TRANSFER FUNCTION Hd(z) R
Page 99 and 100:
1 COEFF ident 1,0 2 include 'head2.
Page 101 and 102:
77 P:0000 0C0040 jmp start ;jump to
Page 103 and 104:
6.2 Transpose Implementation (Direc
Page 105 and 106:
POLES OF TRANSFER FUNCTION Hd(z) Re
Page 107 and 108:
1 C0EFF ident 1,0 2 include 'head1.
Page 109 and 110:
75 76 P:0040 0003F8 ori #3,mr ;disa
Page 111 and 112:
6-22 MOTOROLA
Page 113 and 114:
7.19 Linear-Phase FIR Filter Struct
Page 115 and 116:
G( θ) hi ()e jiθ - N- 1 = i = 0 E
Page 117 and 118:
ut is twice the width of the origin
Page 119 and 120:
Therefore, a nonrecursive filter (F
Page 121 and 122:
sides of by e j2πkm/N and summing
Page 123 and 124:
hi () 1 j πr πkN 1 --- Ak ( )e N
Page 125 and 126:
given with AN ( ⁄ 2) = 0 and r= 1
Page 127 and 128:
Likewise, the y component (imaginar
Page 129 and 130:
h(n) 0.2 0.1 0 -0.1 0 16 31 Figure
Page 131 and 132:
The continuous frequency gain and p
Page 133 and 134:
Filter Gain, G(θ), (dB) Filter Gai
Page 135 and 136:
7-24 MOTOROLA Figure 7-40 FDAS Outp
Page 137 and 138:
If it is assumed that linear phase
Page 139 and 140:
Reference 5). Basically, a guess is
Page 141 and 142:
Figure 7-41 shows an example of an
Page 143 and 144:
R0 x(n) b 0 X:(R0) x(n) x(n-1) x(n-
Page 145 and 146:
7-34 MOTOROLA
Page 147:
16. McClellan, J. H. and T. W. Park
show all

Implementing IIR/FIR Filters

Create successful ePaper yourself

Delete template?

Save as template?